EsTimeAlg.h

Go to the documentation of this file.

#ifndef ESTIME_FLAG_
#define ESTIME_FLAG_
#include "GaudiKernel/Algorithm.h"
#include "GaudiKernel/DataSvc.h"
#include "GaudiKernel/NTuple.h"
#include "HepPDT/ParticleDataTable.hh"
#include <fstream>
#include <iostream>
#include <stdio.h>
#include <vector>
 
#include "CalibData/CalibBase1.h"
#include "CalibData/CalibModel.h"
#include "CalibData/EsTime/EsTimeCalibData.h"
#include "CalibData/Tof/TofCalibData.h"
#include "CalibDataSvc/ICalibRootSvc.h"
#include "EstTofCaliSvc/IEstTofCaliSvc.h"
#include "MdcCalibFunSvc/IMdcCalibFunSvc.h"
#include "MdcUtilitySvc/IMdcUtilitySvc.h"
#include "RawDataProviderSvc/IRawDataProviderSvc.h"
#include "TofQElecSvc/ITofQElecSvc.h"
 
#include "EstParameter.h"
 
typedef std::vector<double> Vdouble;
 
class EsTimeAlg : public Algorithm {
 
public:
  EsTimeAlg( const std::string& name, ISvcLocator* pSvcLocator );
 
  StatusCode initialize();
  StatusCode execute();
  StatusCode finalize();
 
public:
  int m_flag; // 1: use segment linear fitting method
  // 2: use fasttrk and mdc hit
 
  int m_nbunch; // 3: 3 bunches,0ns,8ns,16ns
  // 6: 6 bunches,0,4,8,12,16,20ns
  double m_bunchtime_MC;
  int    m_ntupleflag; // 0: donot draw ntuple
  // 1: add ntuple
 
  int m_optCosmic;    // optCosmic = 1: employ cosmic correction, 0: don't
  int m_cosmicScheme; // cosmicScheme to decide incoming and outgoing legs
  // 0: 2hits+TDCdiff  1: 2hits+Phi  2: All hits+Phi
 
  int    m_phase;          // Detector Version, 1 ETOF->96Scin, 2 ETOF->2MRPC+92Scin, 3->72MRPC
  int    m_userawtime_opt; // 1: use raw tdc time of TOF;  0: use TOF calibration TDC time
  double toffset_rawtime;  // offset for barrel using rawtime of TOF
  double toffset_rawtime_e; // offset for endcap using rawtime of TOF
  double offset_dt;
  double offset_dt_end;
  int    m_debug;
  int    m_beforrec;
  //  double m_t0shift;   //TOF event start time need abstract m_t0shift when do MDC
  //  calibration
  int    m_evtNo; // 1: cout event Number; 0: don't
  double m_ebeam; // beam energy;
  //  void doMDCMatching( int * ){};    // Simple hit matching
  // int   _doMatch;
  // int   _optCosmic, _cosmicScheme;
  // int   _SkipMC,_Skip;
  // double _ptCut, _drCut, _dzCut;
  // double _pathlCut, _ztofCutmin, _ztofCutmax;
  // double _timeOfs;
  // int   _MDC_Skip;
  // double _MDC_ptCut, _MDC_drCut, _MDC_dzCut;
  // int   _MDC_Tof, _MDC_Prop, _MDC_Walk;
  // int   _MDC_Debug;
  // int   _MDC_Xt;
  int _MDC_Inner;
 
  // private:
  //  static evtime * _evtprt;
  bool   m_useXT;
  bool   m_useT0cal;
  bool   m_useSw;
  bool   m_mdcopt;
  bool   m_TofOpt;
  double m_TofOpt_Cut;
  bool   m_useTimeFactor;
 
private:
  StatusCode storeTDS( double tEst, int tEstFlag, double quality ); // yzhang add
  double     Opt_new( const double arr[], const int size, const double sigma_cut );
  double     EST_Trimmer( double t0_original, double t0_offset, double raw_offset,
                          double offset_dt, double bunchTime );
  // double Opt_new(    vector<double>::iterator beg,vector<double>::iterator end);
  int ndriftt;
  int ntrkMC;
  int ntrk;
  int m_pass[5];
 
  IDataProviderSvc* m_pCalibDataSvc;
  ICalibRootSvc*    m_pRootSvc;
  IEstTofCaliSvc*   tofCaliSvc;
  ITofQElecSvc*     tofQElecSvc;
  IMdcCalibFunSvc*  m_mdcCalibFunSvc;
  IMdcUtilitySvc*   m_mdcUtilitySvc;
 
  HepPDT::ParticleDataTable* m_particleTable;
  IRawDataProviderSvc*       m_rawDataProviderSvc;
  NTuple::Tuple*             m_tuple2;
  // NTuple::Tuple*  m_tuple9;
 
  NTuple::Item<int> g_eventNo, g_runNo;
  // MCTruth
 
  NTuple::Item<int>     g_ntrkMC;
  NTuple::Array<double> g_theta0MC, g_phi0MC;
  NTuple::Array<double> g_pxMC, g_pyMC, g_pzMC, g_ptMC;
 
  // recon
  NTuple::Item<double>  g_nmatchbarrel, g_nmatchend;
  NTuple::Item<double>  g_nmatchbarrel_1, g_nmatchbarrel_2;
  NTuple::Item<int>     g_nmatch_tot;
  NTuple::Item<int>     g_ntrk, g_trigtiming;
  NTuple::Array<double> g_ttof;
  NTuple::Array<double> g_vel;
  NTuple::Array<double> g_abmom;
  NTuple::Array<int>    g_pid;
  NTuple::Array<double> g_t0for, g_t0back;
  NTuple::Item<double>  g_meant0, g_t0less, g_t0barrelTof;
  NTuple::Item<double>  g_ndriftt, g_nmatchmdc, g_EstimeMdc;
  NTuple::Item<double>  g_t0mean, g_T0, g_t0, g_mcTestime;
  NTuple::Item<double>  g_meantdc, g_Testime_fzisan, g_Testime;
  NTuple::Item<int>     g_ntofup, g_ntofdown, g_ntofup1, g_ntofdown1;
  NTuple::Item<double>  g_difftof_b, g_difftof_e;
  NTuple::Array<double> g_meantevup, g_meantevdown;
  NTuple::Item<int>     m_estFlag;
  NTuple::Item<double>  m_estTime;
 
  NTuple::Tuple*        m_tuple9;
  NTuple::Array<double> g_meantev;
  NTuple::Item<int>     g_nmatch;
 
  NTuple::Tuple*       m_tuple3;
  NTuple::Item<double> g_t0offset_b;
  NTuple::Item<double> g_t0offset_e;
  NTuple::Item<int>    g_bunchtime;
};
 
#endif