PIDLikelihood.h

/**************************************************************************
 * basf2 (Belle II Analysis Software Framework)                           *
 * Author: The Belle II Collaboration                                     *
 *                                                                        *
 * See git log for contributors and copyright holders.                    *
 * This file is licensed under LGPL-3.0, see LICENSE.md.                  *
 **************************************************************************/
 
#pragma once
 
#include <framework/datastore/RelationsObject.h>
#include <framework/gearbox/Const.h>
 
#include <map>
#include <string>
 
namespace Belle2 {

 
  class PIDLikelihood : public RelationsObject {
 
  public:

    PIDLikelihood();

    void setLogLikelihood(Const::EDetector det,
                          const Const::ChargedStable& part,
                          float logl);

    void subtractMaximum();

    bool isAvailable(Const::PIDDetectorSet set = Const::PIDDetectorSet::set()) const;

    bool areAllAvailable(Const::PIDDetectorSet set = Const::PIDDetectorSet::set()) const
    {
      return  m_detectors.contains(set);
    }

    float getLogL(const Const::ChargedStable& part,
                  Const::PIDDetectorSet set = Const::PIDDetectorSet::set()) const;

    float getDeltaLogL(const Const::ChargedStable& p1,
                       const Const::ChargedStable& p2,
                       Const::PIDDetectorSet set = Const::PIDDetectorSet::set()) const
    {
      return getLogL(p1, set) - getLogL(p2, set);
    }

    double getProbability(const Const::ChargedStable& p1,
                          const Const::ChargedStable& p2,
                          Const::PIDDetectorSet set = Const::PIDDetectorSet::set()) const
    {
      return getProbability(p1, p2, 1.0, set);
    }

    double getProbability(const Const::ChargedStable& p1,
                          const Const::ChargedStable& p2,
                          double ratio,
                          Const::PIDDetectorSet set = Const::PIDDetectorSet::set()) const;
 

    double getProbability(const Const::ChargedStable& part,
                          const double* fractions = nullptr,
                          Const::PIDDetectorSet set = Const::PIDDetectorSet::set()) const;

    double getLogarithmicProbability(const Const::ChargedStable& part,
                                     const double* fractions = nullptr,
                                     Const::PIDDetectorSet set = Const::PIDDetectorSet::set()) const;

    Const::ChargedStable getMostLikely(const double* fractions = nullptr,
                                       Const::PIDDetectorSet set = Const::PIDDetectorSet::set()) const;
    void printArray() const;

    std::string getInfoHTML() const override;

    void addPreOfficialLikelihood(const std::string& preOfficialIdentifier,
                                  const double preOfficialLikelihood);

    double getPreOfficialLikelihood(const std::string& preOfficialIdentifier) const;
 
  private:

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

    std::map<std::string, double> m_preOfficialLikelihoods;
 
    float m_logl[Const::PIDDetectors::c_size][Const::ChargedStable::c_SetSize]; 
    Const::DetectorSet m_detectors;   
 
    ClassDefOverride(PIDLikelihood, 5); 
    // version 5: reorder class members to improve memory layout
 
  };

} // end namespace Belle2