TimeWalkCalibrationAlgorithm.h

/**************************************************************************
 * BASF2 (Belle Analysis Framework 2)                                     *
 * Copyright(C) 2010 - Belle II Collaboration                             *
 *                                                                        *
 * Author: The Belle II Collaboration                                     *
 * Contributors:  Makoto Uchida                                           *
 *                                                                        *
 * This software is provided "as is" without any warranty.                *
 **************************************************************************/
 
#pragma once
#include <calibration/CalibrationAlgorithm.h>
#include <cdc/dbobjects/CDCGeometry.h>
#include <framework/database/DBObjPtr.h>
#include <TH1D.h>
#include <TH2D.h>
#include "string"
 
namespace Belle2 {
  namespace CDC {
    class TimeWalkCalibrationAlgorithm  :  public CalibrationAlgorithm {
    public:
      TimeWalkCalibrationAlgorithm();
 
      ~TimeWalkCalibrationAlgorithm() {}
 
      void setDebug(bool debug) {m_debug = debug; }
 
      void setStoreHisto(bool storeHist) {m_storeHisto = storeHist;}
 
      void setMinimumNDF(double ndf) {m_minNdf = ndf;}
 
      void setMinimumPval(double pval) {m_minPval = pval;}
 
      void enableTextOutput(bool output = true) {m_textOutput = output;}
 
      void outputTWFileName(std::string tw_out_name) {m_outputTWFileName.assign(tw_out_name);}
      void outputT0FileName(std::string t0_out_name) {m_outputT0FileName.assign(t0_out_name);}
 
      void setHistFileName(const std::string& name) {m_histName = "histTW_" + name + ".root";}
 
      void setMode(unsigned short mode)
      {
        m_twParamMode = mode;
        if (mode == 0)
          m_nTwParams = 1;
        else if (mode == 1)
          m_nTwParams = 2;
        else
          B2FATAL("Mode hasn't implemented yet");
      }
 
    protected:
      EResult calibrate() override;
 
      void createHisto();
 
      void write();
 
      void updateT0();
 
      void fitToExponentialFunc(TH1D* h1);
 
      void prepare();
 
      EResult checkConvergence();
 
      void storeHist();
 
    private:
 
      TH1D* m_h1[300];                               
      TH2D* m_h2[300];                               
      double m_xmin = 0.07;                          
      double m_minNdf = 20;                          
      double m_minPval = 0.;                         
      double m_constTerm[300] = {0.};                
      unsigned short m_flag[300] = {0};                    
      std::vector<float> m_tw_old[300];              
      std::vector<float> m_tw_new[300];              
      bool m_storeHisto = true;                      
      std::string m_inputTWFileName = "tw.dat";      
      std::string m_inputT0FileName = "t0.dat";      
      std::string m_outputTWFileName = "tw_new.dat"; 
      std::string m_outputT0FileName = "t0_new.dat"; 
      std::string m_histName = "histTW.root";        
      bool m_debug = false;                          
      bool m_textOutput = false;                     
      unsigned short m_twParamMode = 1;              
      unsigned short m_nTwParams = 2;                
      DBObjPtr<CDCGeometry> m_cdcGeo;                
    };
  }
}