Cluster Class Reference

Class to collect log likelihoods from Clusters from ECL and KLM aimed for output to mdst includes functions to return combined likelihood probability. More...

#include <Cluster.h>

Inheritance diagram for Cluster:

Public Member Functions
	Cluster ()
	Default constructor: log likelihoods and flags set to 0.
void	setLogLikelihood (Const::EDetector det, const Const::Cluster &cluster, float logl)
	Set log likelihood for a given detector and particle.
bool	isAvailable (Const::ClusterDetectorSet set) const
	Check whether cluster information from a given set of detectors is available.
float	getLogL (const Const::Cluster &cluster, Const::ClusterDetectorSet set=Const::ClusterDetectorSet::set()) const
	Return log likelihood for a given detector set and particle.
float	getDeltaLogL (const Const::Cluster &c1, const Const::Cluster &c2, Const::ClusterDetectorSet set=Const::ClusterDetectorSet::set()) const
	Return log likelihood difference for a given detector set and clusters.
double	getProbability (const Const::Cluster &c1, const Const::Cluster &c2, Const::ClusterDetectorSet set=Const::ClusterDetectorSet::set()) const
	Return combined likelihood probability for a cluster being c1 and not c2, assuming equal prior probabilities.
double	getProbability (const Const::Cluster &c1, const Const::Cluster &c2, double ratio, Const::ClusterDetectorSet set=Const::ClusterDetectorSet::set()) const
	Return combined likelihood probability for a cluster being c1 and not c2.
double	getProbability (const Const::Cluster &cluster, const double *fractions=0, Const::ClusterDetectorSet set=Const::ClusterDetectorSet::set()) const
	Return combined likelihood probability for a cluster according to clusterSet; if prior fractions not given equal prior probabilities assumed.
Const::Cluster	getMostLikely (const double *fractions=0, Const::ClusterDetectorSet set=Const::ClusterDetectorSet::set()) const
	Return most likely particle among clusterSet; if prior fractions not given equal prior probabilities assumed.
void	printArray () const
	Prints the content of a private array of log likelihoods.
std::string	getInfoHTML () const override
	Return HTML Info of cluster Likelihoods.
void	addRelationTo (const RelationsInterface< BASE > *object, float weight=1.0, const std::string &namedRelation="") const
	Add a relation from this object to another object (with caching).
void	addRelationTo (const TObject *object, float weight=1.0, const std::string &namedRelation="") const
	Add a relation from this object to another object (no caching, can be quite slow).
void	copyRelations (const RelationsInterface< BASE > *sourceObj)
	Copies all relations of sourceObj (pointing from or to sourceObj) to this object (including weights).
template<class TO >
RelationVector< TO >	getRelationsTo (const std::string &name="", const std::string &namedRelation="") const
	Get the relations that point from this object to another store array.
template<class FROM >
RelationVector< FROM >	getRelationsFrom (const std::string &name="", const std::string &namedRelation="") const
	Get the relations that point from another store array to this object.
template<class T >
RelationVector< T >	getRelationsWith (const std::string &name="", const std::string &namedRelation="") const
	Get the relations between this object and another store array.
template<class TO >
TO *	getRelatedTo (const std::string &name="", const std::string &namedRelation="") const
	Get the object to which this object has a relation.
template<class FROM >
FROM *	getRelatedFrom (const std::string &name="", const std::string &namedRelation="") const
	Get the object from which this object has a relation.
template<class T >
T *	getRelated (const std::string &name="", const std::string &namedRelation="") const
	Get the object to or from which this object has a relation.
template<class TO >
std::pair< TO *, float >	getRelatedToWithWeight (const std::string &name="", const std::string &namedRelation="") const
	Get first related object & weight of relation pointing to an array.
template<class FROM >
std::pair< FROM *, float >	getRelatedFromWithWeight (const std::string &name="", const std::string &namedRelation="") const
	Get first related object & weight of relation pointing from an array.
template<class T >
std::pair< T *, float >	getRelatedWithWeight (const std::string &name="", const std::string &namedRelation="") const
	Get first related object & weight of relation pointing from/to an array.
virtual std::string	getName () const
	Return a short name that describes this object, e.g.
std::string	getInfo () const
	Return a short summary of this object's contents in raw text format.
std::string	getArrayName () const
	Get name of array this object is stored in, or "" if not found.
int	getArrayIndex () const
	Returns this object's array index (in StoreArray), or -1 if not found.

Protected Member Functions
TClonesArray *	getArrayPointer () const
	Returns the pointer to the raw DataStore array holding this object (protected since these arrays are easy to misuse).

Private Member Functions
void	probability (double probabilities[], const double *fractions, Const::ClusterDetectorSet detSet) const
	Calculate likelihood probabilities.
	ClassDefOverride (Cluster, 1)
	Initial version.
	ClassDef (RelationsInterface, 0)
	defines interface for accessing relations of objects in StoreArray.

Private Attributes
Const::DetectorSet	m_detectors
	set of detectors with cluster information
float	m_logl [Const::ClusterDetectors::c_size][Const::Cluster::c_SetSize]
	log likelihoods
DataStore::StoreEntry *	m_cacheDataStoreEntry
	Cache of the data store entry to which this object belongs.
int	m_cacheArrayIndex
	Cache of the index in the TClonesArray to which this object belongs.

Detailed Description

Class to collect log likelihoods from Clusters from ECL and KLM aimed for output to mdst includes functions to return combined likelihood probability.

Definition at line 26 of file Cluster.h.

Constructor & Destructor Documentation

Cluster()

Cluster

(

)

Default constructor: log likelihoods and flags set to 0.

Definition at line 19 of file Cluster.cc.

{
  for (unsigned short i = 0; i < Const::ClusterDetectors::c_size; i++) {
    for (unsigned int k = 0; k < Const::Cluster::c_SetSize; k++) {
      m_logl[i][k] = 0.0;
    }
  }
}

Member Function Documentation

addRelationTo() [1/2]

void addRelationTo ( const RelationsInterface< BASE > *  object, float  weight = 1.0, const std::string &  namedRelation = ""  ) const

inlineinherited

Add a relation from this object to another object (with caching).

Parameters

object	The object to which the relation should point.
weight	The weight of the relation.
namedRelation	Additional name for the relation, or "" for the default naming

Definition at line 142 of file RelationsObject.h.

    {
      if (object)
        DataStore::Instance().addRelation(this, m_cacheDataStoreEntry, m_cacheArrayIndex,
                                          object, object->m_cacheDataStoreEntry, object->m_cacheArrayIndex, weight, namedRelation);
    }

addRelationTo() [2/2]

void addRelationTo ( const TObject *  object, float  weight = 1.0, const std::string &  namedRelation = ""  ) const

inlineinherited

Add a relation from this object to another object (no caching, can be quite slow).

Parameters

object	The object to which the relation should point.
weight	The weight of the relation.
namedRelation	Additional name for the relation, or "" for the default naming

Definition at line 155 of file RelationsObject.h.

    {
      StoreEntry* toEntry = nullptr;
      int toIndex = -1;
      DataStore::Instance().addRelation(this, m_cacheDataStoreEntry, m_cacheArrayIndex, object, toEntry, toIndex, weight, namedRelation);
    }

copyRelations()

void copyRelations ( const RelationsInterface< BASE > *  sourceObj )

inlineinherited

Copies all relations of sourceObj (pointing from or to sourceObj) to this object (including weights).

Useful if you want to make a complete copy of a StoreArray object to make modifications to it, but retain all information on linked objects.

Note: this only works if sourceObj inherits from the same base (e.g. RelationsObject), and only for related objects that also inherit from the same base.

Definition at line 170 of file RelationsObject.h.

    {
      if (!sourceObj)
        return;
      auto fromRels = sourceObj->getRelationsFrom<RelationsInterface<BASE>>("ALL");
      for (unsigned int iRel = 0; iRel < fromRels.size(); iRel++) {
        fromRels.object(iRel)->addRelationTo(this, fromRels.weight(iRel));
      }
 
      auto toRels = sourceObj->getRelationsTo<RelationsInterface<BASE>>("ALL");
      for (unsigned int iRel = 0; iRel < toRels.size(); iRel++) {
        this->addRelationTo(toRels.object(iRel), toRels.weight(iRel));
      }
    }

getArrayIndex()

int getArrayIndex ( ) const

inlineinherited

Returns this object's array index (in StoreArray), or -1 if not found.

Definition at line 385 of file RelationsObject.h.

    {
      DataStore::Instance().findStoreEntry(this, m_cacheDataStoreEntry, m_cacheArrayIndex);
      return m_cacheArrayIndex;
    }

getArrayName()

std::string getArrayName ( ) const

inlineinherited

Get name of array this object is stored in, or "" if not found.

Definition at line 377 of file RelationsObject.h.

    {
      DataStore::Instance().findStoreEntry(this, m_cacheDataStoreEntry, m_cacheArrayIndex);
      return m_cacheDataStoreEntry ? m_cacheDataStoreEntry->name : "";
    }

getArrayPointer()

TClonesArray * getArrayPointer ( ) const

inlineprotectedinherited

Returns the pointer to the raw DataStore array holding this object (protected since these arrays are easy to misuse).

Definition at line 418 of file RelationsObject.h.

    {
      DataStore::Instance().findStoreEntry(this, m_cacheDataStoreEntry, m_cacheArrayIndex);
      if (!m_cacheDataStoreEntry)
        return nullptr;
      return m_cacheDataStoreEntry->getPtrAsArray();
    }

getDeltaLogL()

float getDeltaLogL ( const Const::Cluster &  c1, const Const::Cluster &  c2, Const::ClusterDetectorSet  set = Const::ClusterDetectorSet::set()  ) const

inline

Return log likelihood difference for a given detector set and clusters.

Parameters

c1	cluster
c2	cluster
set	a set of cluster detectors to use

Returns

log likelihood difference logL_c1 - logL_c2

Definition at line 68 of file Cluster.h.

    {
      return getLogL(c1, set) - getLogL(c2, set);
    }

getInfo()

std::string getInfo ( ) const

inlineinherited

Return a short summary of this object's contents in raw text format.

Returns the contents of getInfoHTML() while translating line-breaks etc.

Note

: You don't need to implement this function (it's not virtual), getInfoHTML() is enough.

Definition at line 370 of file RelationsObject.h.

    {
      return _RelationsInterfaceImpl::htmlToPlainText(getInfoHTML());
    }

getInfoHTML()

std::string getInfoHTML ( ) const

overridevirtual

Return HTML Info of cluster Likelihoods.

Reimplemented from RelationsInterface< BASE >.

Definition at line 209 of file Cluster.cc.

{
  const string detectorName[Const::ClusterDetectors::c_size] = {"ECL", "KLM"};
 
  std::stringstream stream;
  stream << std::setprecision(4);
  stream << "<b>Likelihoods</b><br>";
  stream << "<table>";
  stream << "<tr><th>Cluster / Detector</th>";
  for (unsigned k = 0; k < Const::ClusterDetectors::c_size; k++)
    stream << "<th>" << detectorName[k] << "</th>";
  stream << "</tr>";
  for (unsigned i = 0; i < Const::Cluster::c_SetSize; i++) {
    stream << "<tr>";
    stream << "<td>" << Const::clusterSet.at(i).getPDGCode() << "</td>";
    for (unsigned k = 0; k < Const::ClusterDetectors::c_size; k++) {
      stream << "<td>" << m_logl[k][i] << "</td>";
    }
    stream << "</tr>";
  }
  stream << "</table>";
  stream << "<br>";
 
  return stream.str();
}

getLogL()

float getLogL	(	const Const::Cluster &	cluster,
		Const::ClusterDetectorSet	set = Const::ClusterDetectorSet::set()
	)		const

Return log likelihood for a given detector set and particle.

Parameters

cluster	cluster
set	a set of cluster detectors to use

Returns

log likelihood

Definition at line 49 of file Cluster.cc.

{
  float result = 0;
  for (Const::DetectorSet::Iterator it = Const::ClusterDetectorSet::set().begin();
       it != Const::ClusterDetectorSet::set().end(); ++it) {
    if (set.contains(it))
      result += m_logl[it.getIndex()][cluster.getIndex()];
  }
  return result;
}

getMostLikely()

Const::Cluster getMostLikely	(	const double *	fractions = 0,
		Const::ClusterDetectorSet	set = Const::ClusterDetectorSet::set()
	)		const

Return most likely particle among clusterSet; if prior fractions not given equal prior probabilities assumed.

Parameters

fractions	array of prior probabilities in the order defined in Const::Cluster
set	a set of PID detectors to use

Returns

particle type

Definition at line 102 of file Cluster.cc.

{
  const unsigned int n = Const::Cluster::c_SetSize;
  double prob[n];
  probability(prob, fractions, detSet);
 
  int k = 0;
  double maxProb = prob[k];
  for (unsigned i = 0; i < n; ++i) {
    if (prob[i] > maxProb) {maxProb = prob[i]; k = i;}
  }
  return Const::clusterSet.at(k);
 
}

getName()

virtual std::string getName ( ) const

inlinevirtualinherited

Return a short name that describes this object, e.g.

pi+ for an MCParticle.

Reimplemented in Particle, MCParticle, and SpacePoint.

Definition at line 344 of file RelationsObject.h.

{ return ""; }

getProbability() [1/3]

double getProbability ( const Const::Cluster &  c1, const Const::Cluster &  c2, Const::ClusterDetectorSet  set = Const::ClusterDetectorSet::set()  ) const

inline

Return combined likelihood probability for a cluster being c1 and not c2, assuming equal prior probabilities.

Parameters

c1	cluster
c2	cluster
set	a set of cluster detectors to use

Returns

likelihood probability

Definition at line 83 of file Cluster.h.

    {
      return getProbability(c1, c2, 1.0, set);
    }

getProbability() [2/3]

Double32_t getProbability	(	const Const::Cluster &	c1,
		const Const::Cluster &	c2,
		double	ratio,
		Const::ClusterDetectorSet	set = Const::ClusterDetectorSet::set()
	)		const

Return combined likelihood probability for a cluster being c1 and not c2.

Parameters

c1	cluster
c2	cluster
ratio	ratio of prior probabilities (c1/c2)
set	a set of cluster detectors to use

Returns

likelihood probability

Definition at line 62 of file Cluster.cc.

{
  if (ratio < 0) {
    B2ERROR("Cluster::probability argument 'ratio' is given with negative value");
    return 0;
  }
  if (ratio == 0) return 0;
 
  double dlogl = getLogL(c2, set) - getLogL(c1, set);
  double res;
  if (dlogl < 0) {
    double elogl = exp(dlogl);
    res = ratio / (ratio + elogl);
  } else {
    double elogl = exp(-dlogl) * ratio; // to prevent overflow for very large dlogl
    res = elogl / (1.0 + elogl);
  }
  //TODO: only necessary if one wants to use mcprod1405 MC sample. Remove when there's a good replacement.
  if (std::isfinite(res))
    return res;
  return 0.0;
}

getProbability() [3/3]

double getProbability	(	const Const::Cluster &	cluster,
		const double *	fractions = 0,
		Const::ClusterDetectorSet	set = Const::ClusterDetectorSet::set()
	)		const

Return combined likelihood probability for a cluster according to clusterSet; if prior fractions not given equal prior probabilities assumed.

Parameters

cluster	cluster
fractions	array of prior probabilities in the order defined in Const::Cluster
set	a set of cluster detectors to use

Returns

likelihood probability (a value btw. 0 and 1)

Definition at line 88 of file Cluster.cc.

{
  double prob[Const::Cluster::c_SetSize];
  probability(prob, fractions, detSet);
 
  int k = cluster.getIndex();
  if (k < 0) return 0;
 
  return prob[k];
 
}

getRelated()

T * getRelated ( const std::string &  name = "", const std::string &  namedRelation = ""  ) const

inlineinherited

Get the object to or from which this object has a relation.

Template Parameters

T	The class of objects to or from which the relation points.

Parameters

name	The name of the store array to or from which the relation points. If empty the default store array name for class T will be used. If the special name "ALL" is given all store arrays containing objects of type T are considered.
namedRelation	Additional name for the relation, or "" for the default naming

Returns

The first related object or a null pointer.

Definition at line 278 of file RelationsObject.h.

    {
      return static_cast<T*>(DataStore::Instance().getRelationWith(DataStore::c_BothSides, this, m_cacheDataStoreEntry, m_cacheArrayIndex,
                             T::Class(), name, namedRelation).object);
    }

getRelatedFrom()

FROM * getRelatedFrom ( const std::string &  name = "", const std::string &  namedRelation = ""  ) const

inlineinherited

Get the object from which this object has a relation.

Template Parameters

FROM	The class of objects from which the relation points.

Parameters

name	The name of the store array from which the relation points. If empty the default store array name for class FROM will be used. If the special name "ALL" is given all store arrays containing objects of type FROM are considered.
namedRelation	Additional name for the relation, or "" for the default naming

Returns

The first related object or a null pointer.

Definition at line 263 of file RelationsObject.h.

    {
      return static_cast<FROM*>(DataStore::Instance().getRelationWith(DataStore::c_FromSide, this, m_cacheDataStoreEntry,
                                m_cacheArrayIndex, FROM::Class(), name, namedRelation).object);
    }

getRelatedFromWithWeight()

std::pair< FROM *, float > getRelatedFromWithWeight ( const std::string &  name = "", const std::string &  namedRelation = ""  ) const

inlineinherited

Get first related object & weight of relation pointing from an array.

Template Parameters

FROM	The class of objects from which the relation points.

Parameters

name	The name of the store array from which the relation points. If empty the default store array name for class FROM will be used. If the special name "ALL" is given all store arrays containing objects of type FROM are considered.
namedRelation	Additional name for the relation, or "" for the default naming

Returns

Pair of first related object and the relation weight, or (NULL, 1.0) if none found.

Definition at line 314 of file RelationsObject.h.

    {
      RelationEntry entry = DataStore::Instance().getRelationWith(DataStore::c_FromSide, this, m_cacheDataStoreEntry, m_cacheArrayIndex,
                                                                  FROM::Class(), name, namedRelation);
      return std::make_pair(static_cast<FROM*>(entry.object), entry.weight);
    }

getRelatedTo()

TO * getRelatedTo ( const std::string &  name = "", const std::string &  namedRelation = ""  ) const

inlineinherited

Get the object to which this object has a relation.

Template Parameters

TO	The class of objects to which the relation points.

Parameters

name	The name of the store array to which the relation points. If empty the default store array name for class TO will be used. If the special name "ALL" is given all store arrays containing objects of type TO are considered.
namedRelation	Additional name for the relation, or "" for the default naming

Returns

The first related object or a null pointer.

Definition at line 248 of file RelationsObject.h.

    {
      return static_cast<TO*>(DataStore::Instance().getRelationWith(DataStore::c_ToSide, this, m_cacheDataStoreEntry, m_cacheArrayIndex,
                              TO::Class(), name, namedRelation).object);
    }

getRelatedToWithWeight()

std::pair< TO *, float > getRelatedToWithWeight ( const std::string &  name = "", const std::string &  namedRelation = ""  ) const

inlineinherited

Get first related object & weight of relation pointing to an array.

Template Parameters

TO	The class of objects to which the relation points.

Parameters

name	The name of the store array to which the relation points. If empty the default store array name for class TO will be used. If the special name "ALL" is given all store arrays containing objects of type TO are considered.
namedRelation	Additional name for the relation, or "" for the default naming

Returns

Pair of first related object and the relation weight, or (NULL, 1.0) if none found.

Definition at line 297 of file RelationsObject.h.

    {
      RelationEntry entry = DataStore::Instance().getRelationWith(DataStore::c_ToSide, this, m_cacheDataStoreEntry, m_cacheArrayIndex,
                                                                  TO::Class(), name, namedRelation);
      return std::make_pair(static_cast<TO*>(entry.object), entry.weight);
    }

getRelatedWithWeight()

std::pair< T *, float > getRelatedWithWeight ( const std::string &  name = "", const std::string &  namedRelation = ""  ) const

inlineinherited

Get first related object & weight of relation pointing from/to an array.

Template Parameters

T	The class of objects to or from which the relation points.

Parameters

name	The name of the store array to or from which the relation points. If empty the default store array name for class T will be used. If the special name "ALL" is given all store arrays containing objects of type T are considered.
namedRelation	Additional name for the relation, or "" for the default naming

Returns

Pair of first related object and the relation weight, or (NULL, 1.0) if none found.

Definition at line 331 of file RelationsObject.h.

    {
      RelationEntry entry = DataStore::Instance().getRelationWith(DataStore::c_BothSides, this, m_cacheDataStoreEntry, m_cacheArrayIndex,
                                                                  T::Class(), name, namedRelation);
      return std::make_pair(static_cast<T*>(entry.object), entry.weight);
    }

getRelationsFrom()

RelationVector< FROM > getRelationsFrom ( const std::string &  name = "", const std::string &  namedRelation = ""  ) const

inlineinherited

Get the relations that point from another store array to this object.

Template Parameters

FROM	The class of objects from which the relations point.

Parameters

name	The name of the store array from which the relations point. If empty the default store array name for class FROM will be used. If the special name "ALL" is given all store arrays containing objects of type FROM are considered.
namedRelation	Additional name for the relation, or "" for the default naming

Returns

A vector of relations.

Definition at line 212 of file RelationsObject.h.

    {
      return RelationVector<FROM>(DataStore::Instance().getRelationsWith(DataStore::c_FromSide, this, m_cacheDataStoreEntry,
                                  m_cacheArrayIndex, FROM::Class(), name, namedRelation));
    }

getRelationsTo()

RelationVector< TO > getRelationsTo ( const std::string &  name = "", const std::string &  namedRelation = ""  ) const

inlineinherited

Get the relations that point from this object to another store array.

Template Parameters

TO	The class of objects to which the relations point.

Parameters

name	The name of the store array to which the relations point. If empty the default store array name for class TO will be used. If the special name "ALL" is given all store arrays containing objects of type TO are considered.
namedRelation	Additional name for the relation, or "" for the default naming

Returns

A vector of relations.

Definition at line 197 of file RelationsObject.h.

    {
      return RelationVector<TO>(DataStore::Instance().getRelationsWith(DataStore::c_ToSide, this, m_cacheDataStoreEntry,
                                m_cacheArrayIndex, TO::Class(), name, namedRelation));
    }

getRelationsWith()

RelationVector< T > getRelationsWith ( const std::string &  name = "", const std::string &  namedRelation = ""  ) const

inlineinherited

Get the relations between this object and another store array.

Relations in both directions are returned.

Template Parameters

T	The class of objects to or from which the relations point.

Parameters

name	The name of the store array to or from which the relations point. If empty the default store array name for class T will be used. If the special name "ALL" is given all store arrays containing objects of type T are considered.
namedRelation	Additional name for the relation, or "" for the default naming

Returns

A vector of relations.

Definition at line 230 of file RelationsObject.h.

    {
      return RelationVector<T>(DataStore::Instance().getRelationsWith(DataStore::c_BothSides, this, m_cacheDataStoreEntry,
                               m_cacheArrayIndex, T::Class(), name, namedRelation));
    }

isAvailable()

bool isAvailable ( Const::ClusterDetectorSet  set ) const

inline

Check whether cluster information from a given set of detectors is available.

Parameters

set	a set of cluster detectors

Returns

true if the given set of detectors contributed to the cluster information

Definition at line 50 of file Cluster.h.

{return  m_detectors.contains(set);}

printArray()

void printArray

(

)

const

Prints the content of a private array of log likelihoods.

Definition at line 159 of file Cluster.cc.

{
 
  const string detectorName[Const::ClusterDetectors::c_size] = {"ECL", "KLM"};
  string hline("-------");
  for (unsigned i = 0; i < Const::Cluster::c_SetSize; i++)
    hline += "-----------";
  string Hline("=======");
  for (unsigned i = 0; i < Const::Cluster::c_SetSize; i++)
    Hline += "===========";
 
  cout << Hline << endl;
 
  cout << "PDGcode";
  for (unsigned i = 0; i < Const::Cluster::c_SetSize; i++)
    cout << setw(10) << Const::clusterSet.at(i).getPDGCode() << " ";
  cout << endl;
 
  cout << Hline << endl;
 
  float sum_logl[Const::Cluster::c_SetSize];
  for (unsigned i = 0; i < Const::Cluster::c_SetSize; i++) sum_logl[i] = 0;
 
  for (unsigned k = 0; k < Const::ClusterDetectors::c_size; k++) {
    cout << setw(7) << detectorName[k];
    for (unsigned i = 0; i < Const::Cluster::c_SetSize; i++) {
      cout << setw(10) << setprecision(4) << m_logl[k][i] << " ";
      sum_logl[i] += m_logl[k][i];
    }
    cout << endl;
  }
 
  cout << hline << endl;
 
  cout << setw(7) << "sum";
  for (unsigned i = 0; i < Const::Cluster::c_SetSize; i++)
    cout << setw(10) << setprecision(4) << sum_logl[i] << " ";
  cout << endl;
 
  if (isAvailable(Const::ECL) or isAvailable(Const::KLM)) {
    unsigned k = 0;
    for (unsigned i = 0; i < Const::Cluster::c_SetSize; i++)
      if (sum_logl[i] > sum_logl[k]) k = i;
    unsigned pos = 11 * (k + 1) + 3;
    Hline.replace(pos, 1, "^");
  }
  cout << Hline << endl;
 
}

probability()

void probability ( double  probabilities[], const double *  fractions, Const::ClusterDetectorSet  detSet  ) const

private

Calculate likelihood probabilities.

Parameters

fractions	array of prior fractions (not needed to be normalized)
probabilities	array of resulting probabilities
detSet	a set of cluster detectors to use

Definition at line 119 of file Cluster.cc.

{
  const unsigned int n = Const::Cluster::c_SetSize;
  double frac[n];
  if (!fractions) {
    for (unsigned int i = 0; i < n; ++i) frac[i] = 1.0; // normalization not needed
    fractions = frac;
  }
 
  double logL[n];
  double logLmax = 0;
  bool hasMax = false;
  for (unsigned i = 0; i < n; ++i) {
    logL[i] = 0;
    if (fractions[i] > 0) {
      logL[i] = getLogL(Const::clusterSet.at(i), detSet);
      if (!hasMax || logL[i] > logLmax) {
        logLmax = logL[i];
        hasMax = true;
      }
    }
  }
 
  double norm = 0;
  for (unsigned i = 0; i < n; ++i) {
    probabilities[i] = 0;
    if (fractions[i] > 0) probabilities[i] = exp(logL[i] - logLmax) * fractions[i];
    norm += probabilities[i];
  }
  if (norm == 0) return;
 
  for (unsigned i = 0; i < n; ++i) {
    probabilities[i] /= norm;
  }
 
}

setLogLikelihood()

void setLogLikelihood	(	Const::EDetector	det,
		const Const::Cluster &	cluster,
		float	logl
	)

Set log likelihood for a given detector and particle.

Parameters

det	detector enumerator
cluster	cluster
logl	log likelihood

Definition at line 29 of file Cluster.cc.

{
  int index = Const::ClusterDetectors::set().getIndex(det);
  if (index < 0) {
    B2ERROR("ClusterLikelihood::setLogLikelihood: detector is not a cluster device");
    return;
  }
  if (logl != logl or logl == INFINITY) {
    B2ERROR("ClusterLikelihood::setLogLikelihood: log-likelihood for detector " << det << " is " << logl <<
            " (i.e. +inf or NaN)! Ignoring this value. (" << Const::ECL << "=ECL, " << Const::KLM << "=KLM");
 
    return;
  }
  m_detectors += det;
  m_logl[index][part.getIndex()] = logl;
}

Member Data Documentation

m_cacheArrayIndex

int m_cacheArrayIndex

mutableprivateinherited

Cache of the index in the TClonesArray to which this object belongs.

Definition at line 432 of file RelationsObject.h.

m_cacheDataStoreEntry

DataStore::StoreEntry* m_cacheDataStoreEntry

mutableprivateinherited

Cache of the data store entry to which this object belongs.

Definition at line 429 of file RelationsObject.h.

m_detectors

Const::DetectorSet m_detectors

private

set of detectors with cluster information

Definition at line 138 of file Cluster.h.

m_logl

float m_logl[Const::ClusterDetectors::c_size][Const::Cluster::c_SetSize]

private

log likelihoods

Definition at line 139 of file Cluster.h.

---

The documentation for this class was generated from the following files:

• mdst/dataobjects/include/Cluster.h
• mdst/dataobjects/src/Cluster.cc