T0CalibrationAlgorithm.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 <calibration/CalibrationAlgorithm.h>
#include <TH1F.h>
#include <cdc/dbobjects/CDCGeometry.h>
#include <framework/database/DBObjPtr.h>
#include "string"
namespace Belle2 {
  namespace CDC {
    class T0CalibrationAlgorithm: public CalibrationAlgorithm {
    public:
      T0CalibrationAlgorithm();
      ~T0CalibrationAlgorithm() {}
      void storeHisto(bool storeHist = false) {m_storeHisto = storeHist;}
      void setMinimumNDF(double minndf) {m_ndfmin = minndf;}
      void setMinimumPval(double minPval) {m_Pvalmin = minPval;}
      void setMaxRMSDt(double maxRMSDt) {m_maxRMSDt = maxRMSDt;}
      void setMaxMeanDt(double maxMeanDt) {m_maxMeanDt = maxMeanDt;}
      void setOffsetMeanDt(double offsetMeanDt) {m_offsetMeanDt = offsetMeanDt;}
      void setMaxBadChannel(double max_bad_channel) {m_maxBadChannel = max_bad_channel;}
 
      void setCommonT0(double commonT0) {m_commonT0 = commonT0;}
 
      void enableTextOutput(bool output = true) {m_textOutput = output;}
 
      void setOutputFileName(std::string outputname) {m_outputT0FileName.assign(outputname);}
 
      void setHistFileName(const std::string& name) {m_histName = "histT0_" + name + ".root";}
 
    protected:
      EResult calibrate() override;
      void createHisto();
      int write();
      double getMeanT0(TH1F* h1);
    private:
      TH1F* m_hTotal;       
      TH1F* m_h1[56][385];    
      TH1F* m_hT0b[300];      
      double m_xmin = 0.07;   
      double m_ndfmin = 5;    
      double m_Pvalmin = 0.;  
      /*Condition to stop iterate minDt <m_maxDt and rmsDt<m_maxRMS*/
      double m_maxMeanDt = 0.2;   
      double m_offsetMeanDt = -0.4;   
      double m_maxRMSDt = 0.22;   
      double m_maxBadChannel = 50;   
      double dt[56][385] = {{0.}};     
      double err_dt[56][385] = {{0.}}; 
      double dtb[300] = {0.};        
      double err_dtb[300] =  {0.};    
      double m_commonT0 = 4825.;  
      bool m_storeHisto = false; 
      bool  m_textOutput = false; 
      std::string m_outputT0FileName = "t0_new.dat"; 
      std::string m_histName = "histT0.root"; 
      DBObjPtr<CDCGeometry> m_cdcGeo; 
    };
  }// name space CDC
} // namespace Belle2