CrudeT0CalibrationAlgorithm.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 <TH1D.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;    
    };
  }

}