CrudeT0CalibrationAlgorithm.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 <TH1F.h>
#include "vector"
namespace Belle2 {
  namespace CDC {
    class CrudeT0CalibrationAlgorithm: public CalibrationAlgorithm {
    public:
      CrudeT0CalibrationAlgorithm();
      virtual ~CrudeT0CalibrationAlgorithm() {}
 
      void setFitWindow(std::vector<unsigned short> window)
      {
        if (window.size() != 2) {
          B2FATAL("Please specify [TDC_min, TDC_max]");
        }
        m_tdcMin = window.at(0);
        m_tdcMax = window.at(1);
      }
      void setInitialValue(float t)
      {
        m_initT0 = t;
      }
 
      void setZOffset(float z)
      {
        m_zOffset = z;
      }
 
      void setCosmics(bool b)
      {
        m_cosmic = b;
      }
 
      void saveHisto();
    protected:
      EResult calibrate() override;
      virtual void createHisto(StoreObjPtr<EventMetaData>& evtPtr);
      virtual void write(StoreObjPtr<EventMetaData>& evtPtr);
    private:
 
      double m_t0b[300];    
      double m_t0[56][400]; 
      bool m_flag[56][400] = {{false}}; 
      TH1D* m_hTDC[56][400];    
      TH1D* m_hTDCBoard[300];   
      TH1D* m_hT0All;           
      unsigned short m_tdcMin = 4500; 
      unsigned short m_tdcMax = 5000; 
      float m_initT0 = 4887.;   
      unsigned short m_minEntries = 100; 
      float m_zOffset = 0.0;        
      bool m_cosmic = true;    
    };
  }
}