TofCorrPID.h

Go to the documentation of this file.

#ifndef ParticleID_TofCorrPID_H
#define ParticleID_TofCorrPID_H
//
// TofCorrPID package: particle Identification with TOF
//
// The offset for hadrons are corrected.
//
// Author: SUN Shengsen, LV Meng     date: 01/16/2013, created
//
 
 
#include "ParticleID/ParticleIDBase.h"
 
class TofCorrPID : public ParticleIDBase {
 
 public:
  
  static TofCorrPID * instance();
  ~TofCorrPID(){;} 
  
  void init();
  void calculate();
  bool IsPidInfoValid() const {return (m_ndof > 0); }
  double chi(int n) const {return m_chi[n];}
  double prob(int n) const {return m_prob[n];}
  double offset(int n) const {return m_offset[n];}
  double sigma(int n) const {return m_sigma[n];}
  int ndof() const {return m_ndof;}
 
 protected:
 
  int particleIDCalculation();
  int neuronPIDCalculation() { return -1;}
  int LikelihoodCalculation() {return -1;}
 
  void   inputParameter( int run );
  double offsetTof( unsigned int ispecies, bool barrel, unsigned int ipmt, double betaGamma, int charge, double zrhit, double dt );
  double offsetTof( unsigned int ispecies, int tofid, double zrhit, double dt );
  double offsetTof( unsigned int ispecies, int tofid1, int tofid2, double zrhit1, double zrhit2, double dt );
  double sigmaTof( unsigned int ispecies, int charge, bool barrel, unsigned int ipmt, int tofid[2], double zrhit[2], double betaGamma );
  double sigmaTof( unsigned int ispecies, int charge, bool barrel, unsigned int ipmt, double zrhit, int ibgbin );
 
  double qCurveFunc( unsigned int layer, double betaGamma );
  double bSigma( unsigned int end, int tofid, double zrhit );
  double bSigma( int tofid[2], double zrhit[2] );
  double eSigma( int tofid, double zrhit );
 
  bool correlationCheck();
 
 private:
  double m_chi[5];
  double m_offset[5];
  double m_prob[5];
  double m_sigma[5];
  double m_chimin;
  double m_pdfmin;
  int    m_ndof;
 
 private:
  double  m_deltaT[7][5]; 
 
 private:
  int     m_runBegin, m_runEnd;
  double  m_p_weight[176][3][5], m_p_common, m_ec_sigma[96][3];
  double  m_q0_bg[3][5], m_par_ab[5][2][5], m_par_pbar_ab[5][2][5];
  double  m_par_sigma[4][8][9], m_par_sig_mom[3][7][3];
  double  m_p_offset[2][4][10][20], m_p_sigma[2][7][10][20];
 
 private:
  TofCorrPID();
  static TofCorrPID *m_pointer;
};
 
#endif