TimeWalkCalibration.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/logging/Logger.h>
 
#include <TH1D.h>
#include <TH2D.h>
#include "string"
 
namespace Belle2 {
  namespace CDC {
    class TimeWalkCalibration {
    public:
      TimeWalkCalibration();
      virtual ~TimeWalkCalibration() {}
      virtual void setUseDB(bool useDB) {m_useDB = useDB; }
      virtual void setStoreHisto(bool storeHist) {m_storeHisto = storeHist;}
      virtual void setMinimumNDF(double minndf) {m_ndfmin = minndf;}
      virtual void setMinimumPval(double minPval) {m_Pvalmin = minPval;}
      virtual void InputFileNames(std::string inputname) {m_InputRootFileName.assign(inputname);}
      virtual void InputTWFileName(std::string tw_name) {m_InputTWFileName.assign(tw_name);}
      virtual void OutputTWFileName(std::string tw_out_name) {m_outputTWFileName.assign(tw_out_name);}
      virtual void OutputT0FileName(std::string t0_out_name) {m_outputT0FileName.assign(t0_out_name);}
      virtual void setMode(unsigned short mode)
      {
        m_twParamMode_new = mode;
        if (mode == 0)
          m_nTwParams_new = 1;
        else if (mode == 1)
          m_nTwParams_new = 2;
        else
          B2FATAL("Mode hasn't implemented yet");
      }
      virtual void execute()
      {
        calibrate();
      }
 
    protected:
      virtual bool calibrate();
      virtual void CreateHisto();
      virtual void Write();
      virtual void updateT0();
      virtual void readTW();
      virtual void fitToExponentialFunc(TH1D* h1);
 
    private:
      TH1D* m_h1[300];                               
      TH2D* m_h2[300];                               
      double m_xmin = 0.07;                          
      double m_ndfmin = 20;                          
      double m_Pvalmin = 0.;                         
      //      double m_oldTW[300][2] = {{0.}};       /**< Time Walk params of previous iteration*/
      //      double m_tw[300][3] = {{0.}};          /**< Time Walk params, result after fitting*
      double m_constTerm[300] = {0.};                
      std::vector<float> m_tw_old[300];              
      std::vector<float> m_tw_new[300];              
      bool m_storeHisto = true;                      
      bool m_useDB = false;                          
      unsigned short m_flag[300];                    
      std::string m_InputTWFileName = "tw.dat";      
      std::string m_InputT0FileName = "t0.dat";      
      std::string m_outputTWFileName = "tw_new.dat"; 
      std::string m_outputT0FileName = "t0.dat";     
      std::string m_InputRootFileName = "rootfile/output*"; 
      int m_firstExperiment = 0;                     
      int m_firstRun = 0;                            
      int m_lastExperiment = -1;                     
      int m_lastRun = -1;                            
      unsigned short m_twParamMode_old;              
      unsigned short m_twParamMode_new = 1;          
      unsigned short m_nTwParams_new = 2;            
      unsigned short m_nTwParams_old;                
    };
  }//namespace CDC
} // namespace Belle2