Mdc/MdcCalibAlg/share/distcalib/src/include/CalibBase.h

Go to the documentation of this file.

#ifndef CALIBBASE_H
#define CALIBBASE_H
 
#include <cstring>
#include <fstream>
#include <iomanip>
#include <iostream>
 
#include "MdcCalibConst.h"
#include "MdcCosGeom.h"
#include "fun.h"
 
#include "TFile.h"
#include "TFolder.h"
#include "TGraph.h"
#include "TH1F.h"
#include "TObjArray.h"
 
using namespace std;
 
class CalibBase {
 
public:
  CalibBase();
  virtual ~CalibBase();
  virtual void init( TObjArray* hlist, MdcCosGeom* pGeom )                                = 0;
  virtual void mergeHist( TFile* fhist )                                                  = 0;
  virtual void calib( MdcCalibConst* calconst, TObjArray* newXtList, TObjArray* r2tList ) = 0;
 
private:
  void renameHist();
 
  static const int NPhiBin   = 20;
  static const int NThetaBin = 9;
 
  TFolder* m_fdcom;
  TH1F*    m_effNtrk;
  TH1F*    m_effNtrkRecHit;
 
  TH1F* m_hresAllInc;
  TH1F* m_hresAllExc;
  TH1F* m_hresAllAve;
 
  TH1F* m_hresInnInc;
  TH1F* m_hresInnExc;
 
  TH1F* m_hresStpInc;
  TH1F* m_hresStpExc;
 
  TH1F* m_hresOutInc;
  TH1F* m_hresOutExc;
 
  TFolder* m_fdResQ;
  TH1F*    m_hresAveAllQ[14];
  TH1F*    m_hresAveOutQ[14];
  TH1F*    m_hresAveLayQ[43][14];
 
  TH1F* m_hratio;
 
  /*      TH1F* m_hTes[10]; */
  TH1F* m_hbbTrkFlg;
  TH1F* m_hTesAll;
  TH1F* m_hTesGood;
  TH1F* m_hTesAllFlag;
  TH1F* m_hTesRec;
  TH1F* m_hTesCalFlag;
  TH1F* m_hTesCalUse;
 
  TH1F* m_hnRawHit;
 
  TH1F* m_hpt;
  TH1F* m_hpMax;
  TH1F* m_hpMaxCms;
  TH1F* m_hptPos;
  TH1F* m_hptNeg;
 
  TH1F* m_hp;
  TH1F* m_hp_cms;
  TH1F* m_hpPos;
  TH1F* m_hpNeg;
  TH1F* m_hpPoscms;
  TH1F* m_hpNegcms;
 
  // results after four momentum cut
  TH1F* m_hp_cut;
 
  TH1F* m_hnTrk;
  TH1F* m_hnTrkCal;
  TH1F* m_hnhitsRec;
  TH1F* m_hnhitsRecInn;
  TH1F* m_hnhitsRecStp;
  TH1F* m_hnhitsRecOut;
  TH1F* m_hnhitsCal;
  TH1F* m_hnhitsCalInn;
  TH1F* m_hnhitsCalStp;
  TH1F* m_hnhitsCalOut;
  TH1F* m_wirehitmap;
  TH1F* m_layerhitmap;
  TH1F* m_hchisq;
 
  // histogram for distribution of noise
  TH1F* m_hnoisephi;
  TH1F* m_hnoiselay;
  TH1F* m_hnoisenhits;
 
  TH1F* m_hdr;
  TH1F* m_hphi0;
  TH1F* m_hkap;
  TH1F* m_hdz;
  TH1F* m_htanl;
  TH1F* m_hcosthe;
  TH1F* m_hcostheNeg;
  TH1F* m_hcosthePos;
 
  TH1F* m_hx0;
  TH1F* m_hy0;
  TH1F* m_hdelZ0;
 
  TGraph* m_grX0Y0;
  int     m_nGrPoint;
 
  TH1F* m_hitEffAll;
  TH1F* m_hitEffRaw;
  TH1F* m_hitEffRec;
 
  /* histograms drift time */
  TFolder* m_fdTime;
  TH1F*    m_htraw[NLAYER];
  TH1F*    m_htdr[NLAYER];
  TH1F*    m_htdrlr[NLAYER][2];
 
  /* histograms for spatial resolution */
  TFolder* m_fdres;
  TH1F*    m_hresInc[NLAYER];
  TH1F*    m_hreslrInc[NLAYER][2];
  TH1F*    m_hresExc[NLAYER];
  TH1F*    m_hreslrExc[NLAYER][2];
 
  TFolder* m_fdmomPhi;
  // in experimental reference frame
  TH1F* m_ppPhi[NPhiBin];
  TH1F* m_pnPhi[NPhiBin];
  TH1F* m_ppThe[NThetaBin];
  TH1F* m_pnThe[NThetaBin];
  TH1F* m_ppThePhi[NThetaBin][NPhiBin];
  TH1F* m_pnThePhi[NThetaBin][NPhiBin];
 
  // in CMS
  TH1F* m_ppPhiCms[NPhiBin];
  TH1F* m_pnPhiCms[NPhiBin];
  TH1F* m_ppTheCms[NThetaBin];
  TH1F* m_pnTheCms[NThetaBin];
  TH1F* m_ppThePhiCms[NThetaBin][NPhiBin];
  TH1F* m_pnThePhiCms[NThetaBin][NPhiBin];
 
  TFolder* m_fdr2d;
  TH1F*    m_hr2dInc[43][2][2][NSDBIN];
  TH1F*    m_hr2dExc[43][2][2][NSDBIN];
 
  TFolder* m_fdr2t;
  TH1F*    m_hr2t[43][18][2][45];
};
 
#endif /* CALIBBASE_H */