TOP1Dpdf.h

/**************************************************************************
 * BASF2 (Belle Analysis Framework 2)                                     *
 * Copyright(C) 2018 - Belle II Collaboration                             *
 *                                                                        *
 * Author: The Belle II Collaboration                                     *
 * Contributors: Marko Staric                                             *
 *                                                                        *
 * This software is provided "as is" without any warranty.                *
 **************************************************************************/
 
#pragma once
 
#include <vector>
#include <string>
#include <TH1F.h>
#include <top/reconstruction/TOPreco.h>
#include <framework/datastore/StoreArray.h>
#include <top/dataobjects/TOPDigit.h>
 
 
namespace Belle2 {
  namespace TOP {
 
    class TOP1Dpdf {
 
    public:
 
      TOP1Dpdf(TOPreco& reco, const StoreArray<TOPDigit>& digits,
               int moduleID, double binSize);
 
      double getBinSize() const {return m_binSize;}
 
      int getModuleID() const {return m_moduleID;}
 
      TH1F getHistogram(std::string name, std::string title) const;
 
      const std::vector<double>& getPhotonTimes() const {return m_times;}
 
      int getNumOfPhotons() const {return m_times.size();}
 
      double getExpectedSignal() const {return m_expectedSignal;}
 
      double getExpectedBG() const {return m_expectedBG;}
 
      double getTminPDF() const {return m_tminPDF;}
 
      double getTmaxPDF() const {return m_tmaxPDF;}
 
      double getTminFot() const {return m_tminFot;}
 
      double getTmaxFot() const {return m_tmaxFot;}
 
      double getMinT0() const {return m_minT0;}
 
      double getMaxT0() const {return m_maxT0;}
 
      int getNumBinsT0() const {return m_numBinsT0;}
 
      double getLogL(double timeShift) const;
 
 
    private:
 
      int m_moduleID = 0;   
      int m_numBins = 0;    
      double m_minTime = 0; 
      double m_maxTime = 0; 
      std::vector<double> m_times;  
      double m_tminPDF = 0; 
      double m_tmaxPDF = 0; 
      double m_tminFot = 0; 
      double m_tmaxFot = 0; 
      double m_minT0 = 0;   
      double m_maxT0 = 0;   
      int m_numBinsT0 = 0;  
      double m_expectedBG = 0; 
      double m_expectedSignal = 0; 
      double m_bkg = 0;     
      double m_binSize = 0;        
      std::vector<double> m_logF;  
      double m_logBkg = 0;         
    };
 
  } // namespace TOP
} // namespace Belle2