Cluster.h

/**************************************************************************
 * BASF2 (Belle Analysis Framework 2)                                     *
 * Copyright(C) 2016 - Belle II Collaboration                             *
 *                                                                        *
 * Author: The Belle II Collaboration                                     *
 * Contributors: Torben Ferber                                            *
 *                                                                        *
 * This software is provided "as is" without any warranty.                *
 **************************************************************************/
 
#ifndef CLUSTER_H
#define CLUSTER_H
 
#include <framework/datastore/RelationsObject.h>
#include <framework/gearbox/Const.h>
 
namespace Belle2 {
  class Cluster : public RelationsObject {
 
  public:
 
    Cluster();
 
    void setLogLikelihood(Const::EDetector det,
                          const Const::Cluster& cluster,
                          float logl);
 
    bool isAvailable(Const::ClusterDetectorSet set) const {return  m_detectors.contains(set);}
 
    float getLogL(const Const::Cluster& cluster,
                  Const::ClusterDetectorSet set = Const::ClusterDetectorSet::set()) const;
 
    float getDeltaLogL(const Const::Cluster& c1,
                       const Const::Cluster& c2,
                       Const::ClusterDetectorSet set = Const::ClusterDetectorSet::set()) const
    {
      return getLogL(c1, set) - getLogL(c2, set);
    }
 
    double getProbability(const Const::Cluster& c1,
                          const Const::Cluster& c2,
                          Const::ClusterDetectorSet set = Const::ClusterDetectorSet::set()) const
    {
      return getProbability(c1, c2, 1.0, set);
    }
 
    double getProbability(const Const::Cluster& c1,
                          const Const::Cluster& c2,
                          double ratio,
                          Const::ClusterDetectorSet set = Const::ClusterDetectorSet::set()) const;
 
 
    double getProbability(const Const::Cluster& cluster,
                          const double* fractions = 0,
                          Const::ClusterDetectorSet set = Const::ClusterDetectorSet::set()) const;
 
    Const::Cluster getMostLikely(const double* fractions = 0,
                                 Const::ClusterDetectorSet set =
                                   Const::ClusterDetectorSet::set()) const;
    void printArray() const;
 
    std::string getInfoHTML() const;
 
  private:
 
    Const::DetectorSet m_detectors;   
    float m_logl[Const::ClusterDetectors::c_size][Const::Cluster::c_SetSize]; 
    void probability(double probabilities[],
                     const double* fractions,
                     Const::ClusterDetectorSet detSet) const;
 
 
    ClassDef(Cluster, 1); 
  };
 
} // end namespace Belle2
 
#endif