ARICHReconstructionPar.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 <TObject.h>
#include <vector>
#include <TF1.h>
#include <TF2.h>
 
namespace Belle2 {
 
  class ARICHReconstructionPar : public TObject {
 
  public:
 
    ARICHReconstructionPar() {};
 
    void initializeDefault();
 
    void setBackgroundPDF(TF1* bkgPDF)
    {
      m_bkgPDF = bkgPDF;
    }
 
    void setBackgroundPhiPDF(TF2* bkgPhiPDF)
    {
      m_bkgPhiPDF = bkgPhiPDF;
    }
 
 
    void setThcResolution(TF1* thcRes)
    {
      m_thcResolution = thcRes;
    }
 
    void setAerogelFOM(const std::vector<float>& aerogelFOM)
    {
      m_aerogelFOM =  aerogelFOM;
    }
 
    void setFlatBkgPerPad(const float flatBkgPerPad)
    {
      m_flatBkgPerPad = flatBkgPerPad;
    }
 
    void setParameters(const std::vector<float>& pars)
    {
      m_pars =  pars;
    }
 
    const std::vector<float>& getParameters() const
    {
      return m_pars;
    }
 
    const TF1* getBackgroundPDF() const
    {
      return m_bkgPDF;
    }
 
    const TF2* getBackgroundPhiPDF() const
    {
      return m_bkgPhiPDF;
    }
 
    double getThcResolution(double momentum) const
    {
      return  m_thcResolution->Eval(momentum);
    }
 
    double getAerogelFOM(unsigned iLayer) const
    {
      return m_aerogelFOM.at(iLayer);
    }
 
    float getFlatBkgPerPad() const
    {
      return m_flatBkgPerPad;
    }
 
    double getBackgroundPerPad(double th_cer, const std::vector<double>& pars) const;
 
    double getPhiCorrectedBackgroundPerPad(double fi_cer_trk, double th_cer, const std::vector<double>& pars) const;
 
 
    double getExpectedBackgroundHits(const std::vector<double>& pars, double minThc = 0.1, double maxThc = 0.5) const;
 
    double getNPadsInRing(double maxThc, double minThc = 0.0, double trackTh = 0.45) const;
 
    void print() const;
 
  private:
 
    mutable TF2* m_bkgPhiPDF = NULL; 
    mutable TF1* m_bkgPDF = NULL; 
    TF1* m_thcResolution = NULL; 
    std::vector<float> m_pars; 
    std::vector<float> m_aerogelFOM; 
    float m_flatBkgPerPad = 0.0; 
    ClassDef(ARICHReconstructionPar, 2);  
  };
 
} // namespace