PreT0MdcCalib Class Reference

#include <PreT0MdcCalib.h>

Inheritance diagram for PreT0MdcCalib:

List of all members.

Public Member Functions
void	clear ()
void	clear ()
int	fillHist (MdcCalEvent *event)
int	fillHist (MdcCalEvent *event)
void	initialize (TObjArray *hlist, IMdcGeomSvc *mdcGeomSvc, IMdcCalibFunSvc *mdcFunSvc)
void	initialize (TObjArray *hlist, IMdcGeomSvc *mdcGeomSvc, IMdcCalibFunSvc *mdcFunSvc)
	PreT0MdcCalib ()
	PreT0MdcCalib ()
virtual void	read_file (std::vector< std::string > path)
virtual void	read_file (std::vector< std::string > path)
void	setParam (MdcCalParams &param)
void	setParam (MdcCalParams &param)
virtual void	settuple (std::string path)
virtual void	settuple (std::string path)
int	updateConst (MdcCalibConst *calconst)
int	updateConst (MdcCalibConst *calconst)
	~PreT0MdcCalib ()
	~PreT0MdcCalib ()
Static Private Member Functions
Double_t	funTmax (Double_t *x, Double_t *par)
Double_t	funTmax (Double_t *x, Double_t *par)
Double_t	funTmin (Double_t *x, Double_t *par)
Double_t	funTmin (Double_t *x, Double_t *par)
Double_t	xtfun (Double_t *x, Double_t *par)
Double_t	xtfun (Double_t *x, Double_t *par)
Private Attributes
TFolder *	m_fdTrec
TFolder *	m_fdTrec
TFolder *	m_fdTrecZ
TFolder *	m_fdTrecZ
TObjArray *	m_hlist
TObjArray *	m_hlist
TH1F *	m_hTrec [MdcCalNLayer][MdcCalLR]
TH1F *	m_hTrec [MdcCalNLayer][MdcCalLR]
TH1F *	m_hTrecCosm [MdcCalNLayer][2]
TH1F *	m_hTrecCosm [MdcCalNLayer][2]
TH1F *	m_hTrecZ [MdcCalNLayer][MdcCalLR][11]
TH1F *	m_hTrecZ [MdcCalNLayer][MdcCalLR][11]
IMdcCalibFunSvc *	m_mdcFunSvc
IMdcCalibFunSvc *	m_mdcFunSvc
IMdcGeomSvc *	m_mdcGeomSvc
IMdcGeomSvc *	m_mdcGeomSvc
int	m_nzbin
MdcCalParams	m_param
double	m_vp [MdcCalNLayer]
double	m_zst [MdcCalNLayer]
double	m_zwid [MdcCalNLayer]

---

Constructor & Destructor Documentation

PreT0MdcCalib::PreT0MdcCalib (   )

{ m_nzbin = 11; }

PreT0MdcCalib::~PreT0MdcCalib (   )

{ }

PreT0MdcCalib::PreT0MdcCalib (   )

PreT0MdcCalib::~PreT0MdcCalib (   )

---

Member Function Documentation

void PreT0MdcCalib::clear (   )  [virtual]

Implements MdcCalib.

void PreT0MdcCalib::clear (   )  [virtual]

Implements MdcCalib. { int lay; int lr; int bin; for(lay=0; lay<MdcCalNLayer; lay++){ for(lr=0; lr<MdcCalLR; lr++){ delete m_hTrec[lay][lr]; for(bin=0; bin<m_nzbin; bin++){ delete m_hTrecZ[lay][lr][bin]; } } } for(lay=0; lay<MdcCalNLayer; lay++){ for(bin=0; bin<2; bin++) delete m_hTrecCosm[lay][bin]; } delete m_fdTrec; delete m_fdTrecZ; MdcCalib::clear(); }

int PreT0MdcCalib::fillHist (  MdcCalEvent *  event  )  [virtual]

Implements MdcCalib.

int PreT0MdcCalib::fillHist (  MdcCalEvent *  event  )  [virtual]

Implements MdcCalib. { IMessageSvc* msgSvc; Gaudi::svcLocator() -> service("MessageSvc", msgSvc); MsgStream log(msgSvc, "PreT0MdcCalib"); log << MSG::DEBUG << "PreT0MdcCalib::fillHist()" << endreq; // MdcCalib::fillHist(event); int i; int k; int lay; int cel; int lr; int bin; double tdc; double traw; double trec; double tdr; double t0; double tp = 0.0; double tof = 0.0; double zhit; double zprop; MdcCalRecTrk* rectrk; MdcCalRecHit* rechit; IDataProviderSvc* eventSvc = NULL; Gaudi::svcLocator()->service("EventDataSvc", eventSvc); // get EsTimeCol double tes = event->getTes(); bool esCutFg = event->getEsCutFlag(); if( ! esCutFg ) return -1; int nhit; int ntrk = event -> getNTrk(); for(i=0; i<ntrk; i++){ rectrk = event->getRecTrk(i); nhit = rectrk -> getNHits(); for(k=0; k<nhit; k++){ rechit = rectrk -> getRecHit(k); lay = rechit -> getLayid(); cel = rechit -> getCellid(); lr = rechit -> getLR(); zhit = rechit -> getZhit(); tdc = rechit -> getTdc(); tdr = rechit -> getTdrift(); traw = tdc * MdcCalTdcCnv; // traw = tdc; // ns zprop = fabs(zhit - m_zst[lay]); bin = (int)(zprop / m_zwid[lay]); tp = zprop / m_vp[lay]; t0 = m_mdcFunSvc -> getT0(lay, cel); tof = rechit -> getTof(); const MdcGeoWire* pWire = m_mdcGeomSvc -> Wire(lay, cel); double y1 = pWire -> Backward().y(); double y2 = pWire -> Forward().y(); double ymid = (y1+y2)*0.5; if(5==m_param.particle){ // cosmic-ray if(ymid > 0) trec = traw - tp + tof - tes + m_param.timeShift; else trec = traw - tp - tof - tes + m_param.timeShift; } else{ trec = traw - tp - tof - tes + m_param.timeShift; } m_hTrec[lay][lr] -> Fill(trec); m_hTrec[lay][2] -> Fill(trec); // l-r in common if(ymid > 0) m_hTrecCosm[lay][0] -> Fill(trec); else m_hTrecCosm[lay][1] -> Fill(trec); if(bin < m_nzbin){ m_hTrecZ[lay][lr][bin] -> Fill(trec); m_hTrecZ[lay][2][bin] -> Fill(trec); } } } return 1; }

Double_t PreT0MdcCalib::funTmax (  Double_t *  x, Double_t *  par )  [static, private]

Double_t PreT0MdcCalib::funTmax (  Double_t *  x, Double_t *  par )  [static, private]

{ Double_t fitval; fitval = par[0] + par[1] / (1 + exp((x[0]-par[2])/par[3])); return fitval; }

Double_t PreT0MdcCalib::funTmin (  Double_t *  x, Double_t *  par )  [static, private]

Double_t PreT0MdcCalib::funTmin (  Double_t *  x, Double_t *  par )  [static, private]

{ Double_t fitval; fitval = par[0] + par[1]*exp( -par[2]*(x[0]-par[3]) ) / ( 1 + exp( -(x[0]-par[4])/par[5] )); return fitval; }

void PreT0MdcCalib::initialize (  TObjArray *  hlist, IMdcGeomSvc *  mdcGeomSvc, IMdcCalibFunSvc *  mdcFunSvc )  [virtual]

Implements MdcCalib.

void PreT0MdcCalib::initialize (  TObjArray *  hlist, IMdcGeomSvc *  mdcGeomSvc, IMdcCalibFunSvc *  mdcFunSvc )  [virtual]

Implements MdcCalib. { IMessageSvc* msgSvc; Gaudi::svcLocator() -> service("MessageSvc", msgSvc); MsgStream log(msgSvc, "PreT0MdcCalib"); log << MSG::INFO << "PreT0MdcCalib::initialize()" << endreq; m_hlist = hlist; m_mdcGeomSvc = mdcGeomSvc; m_mdcFunSvc = mdcFunSvc; // MdcCalib::initialize(m_hlist, m_mdcGeomSvc, m_mdcFunSvc); int lay; int lr; int bin; char hname[200]; m_fdTrec = new TFolder("Trec", "Trec"); m_hlist->Add(m_fdTrec); m_fdTrecZ = new TFolder("TrecZ", "TrecZ"); m_hlist->Add(m_fdTrecZ); for(lay=0; lay<MdcCalNLayer; lay++){ for(lr=0; lr<MdcCalLR; lr++){ if(0 == lr) sprintf(hname, "Trec%02d_L", lay); else if(1 == lr) sprintf(hname, "Trec%02d_R", lay); else sprintf(hname, "Trec%02d_C", lay); if(lay < 8) m_hTrec[lay][lr] = new TH1F(hname, "", 100, 0, 400); else m_hTrec[lay][lr] = new TH1F(hname, "", 125, 0, 500); m_fdTrec -> Add(m_hTrec[lay][lr]); } } for(lay=0; lay<MdcCalNLayer; lay++){ for(int iud=0; iud<2; iud++){ if(0 == iud) sprintf(hname, "TrecCosm%02d_Up", lay); else sprintf(hname, "TrecCosm%02d_Dw", lay); if(lay < 8) m_hTrecCosm[lay][iud] = new TH1F(hname, "", 100, 0, 400); else m_hTrecCosm[lay][iud] = new TH1F(hname, "", 125, 0, 500); m_fdTrec -> Add(m_hTrecCosm[lay][iud]); } } for(lay=0; lay<MdcCalNLayer; lay++){ for(lr=0; lr<MdcCalLR; lr++){ for(bin=0; bin<m_nzbin; bin++){ if(0 == lr) sprintf(hname, "Trec%02d_z%02d_L", lay, bin); else if(1 == lr) sprintf(hname, "Trec%02d_z%02d_R", lay, bin); else sprintf(hname, "Trec%02d_z%02d_C", lay, bin); if(lay < 8) m_hTrecZ[lay][lr][bin] = new TH1F(hname, "", 100, 0, 400); else m_hTrecZ[lay][lr][bin] = new TH1F(hname, "", 125, 0, 500); m_fdTrecZ -> Add(m_hTrecZ[lay][lr][bin]); } } } double zeast; double zwest; for(lay=0; lay<MdcCalNLayer; lay++){ zwest = m_mdcGeomSvc->Wire(lay, 0)->Forward().z(); zeast = m_mdcGeomSvc->Wire(lay, 0)->Backward().z(); m_zwid[lay] = (zeast - zwest) / (double)m_nzbin; if(lay < 8) m_vp[lay] = 220.0; // *10^9 mm/s else m_vp[lay] = 240.0; // *10^9 mm/s if( 0 == (lay % 2) ){ // west end m_zst[lay] = zwest; } else{ // east end m_zst[lay] = zeast; } } }

virtual void PreT0MdcCalib::read_file (  std::vector< std::string >  path  )  [virtual]

Reimplemented from MdcCalib.

void PreT0MdcCalib::read_file (  std::vector< std::string >  path  )  [virtual]

Reimplemented from MdcCalib. { MdcCalib::read_file(path); }

void PreT0MdcCalib::setParam (  MdcCalParams &  param  )  [virtual]

Implements MdcCalib.

void PreT0MdcCalib::setParam (  MdcCalParams &  param  )  [inline, virtual]

Implements MdcCalib. { MdcCalib::setParam(param); m_param = param; }

virtual void PreT0MdcCalib::settuple (  std::string  path  )  [virtual]

Reimplemented from MdcCalib.

void PreT0MdcCalib::settuple (  std::string  path  )  [virtual]

Reimplemented from MdcCalib. { MdcCalib::settuple(path); TFile *f1 = new TFile(path.c_str()); std::cout<<path<<std::endl; // // m_mdcGeomSvc = mdcGeomSvc; //m_mdcFunSvc = mdcFunSvc; // MdcCalib::initialize(m_hlist, m_mdcGeomSvc, m_mdcFunSvc); int lay; int lr; int bin; char hname[200]; m_fdTrec=(TFolder *)f1->Get("Trec;1"); m_fdTrecZ=(TFolder *)f1->Get("TrecZ;1"); for(lay=0; lay<MdcCalNLayer; lay++){ for(lr=0; lr<MdcCalLR; lr++){ if(0 == lr) sprintf(hname, "Trec%02d_L", lay); else if(1 == lr) sprintf(hname, "Trec%02d_R", lay); else sprintf(hname, "Trec%02d_C", lay); if(lay < 8) m_hTrec[lay][lr]=(TH1F *)m_fdTrec ->FindObject(hname); else m_hTrec[lay][lr]=(TH1F *)m_fdTrec ->FindObject(hname); } } for(lay=0; lay<MdcCalNLayer; lay++){ for(lr=0; lr<MdcCalLR; lr++){ for(bin=0; bin<m_nzbin; bin++){ if(0 == lr) sprintf(hname, "Trec%02d_z%02d_L", lay, bin); else if(1 == lr) sprintf(hname, "Trec%02d_z%02d_R", lay, bin); else sprintf(hname, "Trec%02d_z%02d_C", lay, bin); if(lay < 8) m_hTrecZ[lay][lr][bin]=(TH1F *)m_fdTrecZ ->FindObject(hname); else m_hTrecZ[lay][lr][bin]=(TH1F *)m_fdTrecZ ->FindObject(hname); } } } double zeast; double zwest; for(lay=0; lay<MdcCalNLayer; lay++){ zwest = m_mdcGeomSvc->Wire(lay, 0)->Forward().z(); zeast = m_mdcGeomSvc->Wire(lay, 0)->Backward().z(); m_zwid[lay] = (zeast - zwest) / (double)m_nzbin; if(lay < 8) m_vp[lay] = 220.0; // *10^9 mm/s else m_vp[lay] = 240.0; // *10^9 mm/s if( 0 == (lay % 2) ){ // west end m_zst[lay] = zwest; } else{ // east end m_zst[lay] = zeast; } } }

int PreT0MdcCalib::updateConst (  MdcCalibConst *  calconst  )  [virtual]

Implements MdcCalib.

int PreT0MdcCalib::updateConst (  MdcCalibConst *  calconst  )  [virtual]

Implements MdcCalib. { IMessageSvc* msgSvc; Gaudi::svcLocator() -> service("MessageSvc", msgSvc); MsgStream log(msgSvc, "PreT0MdcCalib"); log << MSG::DEBUG << "PreT0MdcCalib::updateConst()" << endreq; // MdcCalib::updateConst(calconst); // fit Tmin int lay; int wir; int lr; double t0Fit[MdcCalNLayer][MdcCalLR]; double t0Cal[MdcCalNLayer][MdcCalLR]; double tmax[MdcCalNLayer][MdcCalLR]; double initT0 = m_param.initT0; int fitTminFg[MdcCalNLayer][MdcCalLR]; int fitTmaxFg[MdcCalNLayer][MdcCalLR]; double chisq; double ndf; double chindfTmin[MdcCalNLayer][MdcCalLR]; double chindfTmax[MdcCalNLayer][MdcCalLR]; char funname[200]; // add for cosmic-ray TF1* ftminCosm[MdcCalNLayer][2]; double t0FitCosm[MdcCalNLayer][2]; TF1* ftmin[MdcCalNLayer][MdcCalLR]; for(lay=0; lay<MdcCalNLayer; lay++){ for(lr=0; lr<MdcCalLR; lr++){ fitTminFg[lay][lr] = 0; chindfTmin[lay][lr] = -1; sprintf(funname, "ftmin%02d_%d", lay, lr); ftmin[lay][lr] = new TF1(funname, funTmin, 0, 150, 6); if(1 == m_param.fgCalib[lay]){ ftmin[lay][lr] -> SetParameter(0, 0); ftmin[lay][lr] -> SetParameter(4, initT0); ftmin[lay][lr] -> SetParameter(5, 1); m_hTrec[lay][lr] -> Fit(funname, "Q", "", m_param.tminFitRange[lay][0], m_param.tminFitRange[lay][1]); gStyle -> SetOptFit(11); chisq = ftmin[lay][lr]->GetChisquare(); ndf = ftmin[lay][lr]->GetNDF(); chindfTmin[lay][lr] = chisq / ndf; if(chindfTmin[lay][lr] < m_param.tminFitChindf){ fitTminFg[lay][lr] = 1; t0Fit[lay][lr] = ftmin[lay][lr]->GetParameter(4); t0Fit[lay][lr] += m_param.t0Shift; t0Cal[lay][lr] = t0Fit[lay][lr] - m_param.timeShift; } } if(0 == fitTminFg[lay][lr]){ wir = m_mdcGeomSvc->Wire(lay, 0)->Id(); t0Cal[lay][lr] = calconst->getT0(wir); t0Fit[lay][lr] = t0Cal[lay][lr] + m_param.timeShift; } } for(int iud=0; iud<2; iud++){ sprintf(funname, "ftminCosm_%02d_%d", lay, iud); ftminCosm[lay][iud] = new TF1(funname, funTmin, 0, 150, 6); ftminCosm[lay][iud] -> SetParameter(0, 0); ftminCosm[lay][iud] -> SetParameter(4, initT0); ftminCosm[lay][iud] -> SetParameter(5, 1); m_hTrecCosm[lay][iud] -> Fit(funname, "Q", "", m_param.tminFitRange[lay][0], m_param.tminFitRange[lay][1]); gStyle -> SetOptFit(11); t0FitCosm[lay][iud] += m_param.t0Shift; t0FitCosm[lay][iud] = ftminCosm[lay][iud]->GetParameter(4); } } // fit Tmax TF1* ftmax[MdcCalNLayer][MdcCalLR]; for(lay=0; lay<MdcCalNLayer; lay++){ for(lr=0; lr<MdcCalLR; lr++){ fitTmaxFg[lay][lr] = 0; chindfTmax[lay][lr] = -1; sprintf(funname, "ftmax%02d_%d", lay, lr); ftmax[lay][lr] = new TF1(funname, funTmax, 250, 500, 4); if(1 == m_param.fgCalib[lay]){ ftmax[lay][lr] -> SetParameter(2, m_param.initTm[lay]); ftmax[lay][lr] -> SetParameter(3, 10); m_hTrec[lay][lr] -> Fit(funname, "Q+", "", m_param.tmaxFitRange[lay][0], m_param.tmaxFitRange[lay][1]); gStyle -> SetOptFit(11); chisq = ftmax[lay][lr]->GetChisquare(); ndf = ftmax[lay][lr]->GetNDF(); chindfTmax[lay][lr] = chisq / ndf; if(chindfTmax[lay][lr] < m_param.tmaxFitChindf){ fitTmaxFg[lay][lr] = 1; tmax[lay][lr] = ftmax[lay][lr]->GetParameter(2); } } if(0 == fitTmaxFg[lay][lr]){ tmax[lay][lr] = (calconst->getXtpar(lay, 0, lr, 6)) + t0Fit[lay][2]; } } } // output for check ofstream ft0("preT0.dat"); for(lay=0; lay<MdcCalNLayer; lay++){ ft0 << setw(5) << lay << setw(3) << fitTminFg[lay][2] << setw(15) << t0Cal[lay][2] << setw(15) << t0Fit[lay][2] << setw(15) << chindfTmin[lay][2] << endl; } ft0 << endl; for(lay=0; lay<MdcCalNLayer; lay++){ ft0 << setw(5) << lay << setw(3) << fitTmaxFg[lay][0] << setw(10) << tmax[lay][0] << setw(10) << chindfTmax[lay][0] << setw(3) << fitTmaxFg[lay][1] << setw(10) << tmax[lay][1] << setw(10) << chindfTmax[lay][1] << setw(3) << fitTmaxFg[lay][2] << setw(10) << tmax[lay][2] << setw(10) << chindfTmax[lay][2] << setw(10) << tmax[lay][0] - t0Fit[lay][2] << setw(10) << tmax[lay][1] - t0Fit[lay][2] << setw(10) << tmax[lay][2] - t0Fit[lay][2] << endl; } ft0.close(); cout << "preT0.dat was written." << endl; // output for cosmic T0 ofstream ft0cosm("cosmicT0.dat"); for(lay=0; lay<MdcCalNLayer; lay++){ ft0cosm << setw(5) << lay << setw(15) << t0Fit[lay][2] << setw(15) << t0FitCosm[lay][0] << setw(15) << t0FitCosm[lay][1] << endl; } ft0cosm.close(); // set T0 int i; int nwire = m_mdcGeomSvc -> getWireSize(); for(i=0; i<nwire; i++){ lay = m_mdcGeomSvc -> Wire(i) -> Layer(); if(1 == m_param.fgCalib[lay]){ calconst -> resetT0(i, t0Cal[lay][2]); calconst -> resetDelT0(i, 0.0); } } // set tm of X-T if(m_param.preT0SetTm){ int iEntr; double tm; for(lay=0; lay<MdcCalNLayer; lay++){ if(1 != m_param.fgCalib[lay]) continue; for(iEntr=0; iEntr<MdcCalNENTRXT; iEntr++){ for(lr=0; lr<MdcCalLR; lr++){ tm = tmax[lay][lr] - t0Fit[lay][2]; if( (tmax[lay][lr] > m_param.tmaxFitRange[lay][0]) && (tmax[lay][lr] < m_param.tmaxFitRange[lay][1]) ){ calconst -> resetXtpar(lay, iEntr, lr, 6, tm); } } } } } for(lay=0; lay<MdcCalNLayer; lay++){ for(lr=0; lr<MdcCalLR; lr++){ delete ftmin[lay][lr]; delete ftmax[lay][lr]; } } return 1; }

Double_t PreT0MdcCalib::xtfun (  Double_t *  x, Double_t *  par )  [static, private]

Double_t PreT0MdcCalib::xtfun (  Double_t *  x, Double_t *  par )  [static, private]

---

Member Data Documentation

TFolder* PreT0MdcCalib::m_fdTrec [private]

TFolder* PreT0MdcCalib::m_fdTrec [private]

TFolder* PreT0MdcCalib::m_fdTrecZ [private]

TFolder* PreT0MdcCalib::m_fdTrecZ [private]

TObjArray* PreT0MdcCalib::m_hlist [private]

Reimplemented from MdcCalib.

TObjArray* PreT0MdcCalib::m_hlist [private]

Reimplemented from MdcCalib.

TH1F* PreT0MdcCalib::m_hTrec[MdcCalNLayer][MdcCalLR] [private]

TH1F* PreT0MdcCalib::m_hTrec[MdcCalNLayer][MdcCalLR] [private]

TH1F* PreT0MdcCalib::m_hTrecCosm[MdcCalNLayer][2] [private]

TH1F* PreT0MdcCalib::m_hTrecCosm[MdcCalNLayer][2] [private]

TH1F* PreT0MdcCalib::m_hTrecZ[MdcCalNLayer][MdcCalLR][11] [private]

TH1F* PreT0MdcCalib::m_hTrecZ[MdcCalNLayer][MdcCalLR][11] [private]

IMdcCalibFunSvc* PreT0MdcCalib::m_mdcFunSvc [private]

Reimplemented from MdcCalib.

IMdcCalibFunSvc* PreT0MdcCalib::m_mdcFunSvc [private]

Reimplemented from MdcCalib.

IMdcGeomSvc* PreT0MdcCalib::m_mdcGeomSvc [private]

Reimplemented from MdcCalib.

IMdcGeomSvc* PreT0MdcCalib::m_mdcGeomSvc [private]

Reimplemented from MdcCalib.

int PreT0MdcCalib::m_nzbin [private]

MdcCalParams PreT0MdcCalib::m_param [private]

Reimplemented from MdcCalib.

double PreT0MdcCalib::m_vp [private]

double PreT0MdcCalib::m_zst [private]

double PreT0MdcCalib::m_zwid [private]

---

The documentation for this class was generated from the following files:

• /data0/mdestefa/bes3soft/Boss_6.5.5/dist/6.5.5/InstallArea/include/MdcCalibAlg/MdcCalibAlg/PreT0MdcCalib.h
• /data0/mdestefa/bes3soft/Boss_6.5.5/dist/6.5.5/Mdc/MdcCalibAlg/MdcCalibAlg-00-08-02/MdcCalibAlg/PreT0MdcCalib.h
• /data0/mdestefa/bes3soft/Boss_6.5.5/dist/6.5.5/Mdc/MdcCalibAlg/MdcCalibAlg-00-08-02/src/PreT0MdcCalib.cxx

---