#include <MdcCalibFunSvc.h>

Inheritance diagram for MdcCalibFunSvc:

Public Member Functions
	MdcCalibFunSvc (const std::string &name, ISvcLocator *svcloc)

	~MdcCalibFunSvc ()

virtual StatusCode	initialize ()

virtual StatusCode	finalize ()

void	handle (const Incident &)

double	getVprop (int lay) const

double	getTprop (int lay, double z) const

double	driftTimeToDist (double drifttime, int layid, int cellid, int lr, double entrance=0.0) const

double	distToDriftTime (double dist, int layid, int cellid, int lr, double entrance=0.0) const

double	getSigma (int layid, int lr, double dist, double entrance=0.0, double tanlam=0.0, double z=0.0, double Q=1000.0) const

double	getSigmaLR (int layid, int lr, double dist, double entrance=0.0, double tanlam=0.0, double z=0.0, double Q=1000.0) const

double	getSigma1 (int layid, int lr, double dist, double entrance=0.0, double tanlam=0.0, double z=0.0, double Q=1000.0) const

double	getSigma2 (int layid, int lr, double dist, double entrance=0.0, double tanlam=0.0, double z=0.0, double Q=1000.0) const

double	getF (int layid, int lr, double dist, double entrance=0.0, double tanlam=0.0, double z=0.0, double Q=1000.0) const

double	getSigmaToT (int layid, int lr, double tdr, double entrance=0.0, double tanlam=0.0, double z=0.0, double Q=1000.0) const

double	getSigmaToTLR (int layid, int lr, double tdr, double entrance=0.0, double tanlam=0.0, double z=0.0, double Q=1000.0) const

void	setXtBegin ()

int	getNextXtpar (int &key, double &par)

void	getXtpar (int layid, int entr, int lr, double par[]) const

bool	getNewXtpar ()

TTree *	getNewXtparTree (int layid, int entr, int lr) const

double	getT0 (int layid, int cellid) const

double	getT0 (int wireid) const

double	getDelT0 (int wireid) const

double	getTimeWalk (int layid, double Q) const

double	getQtpar (int layid, int ord) const

double	getWireEff (int layid, int cellid) const

double	getWireEff (int wireid) const

void	setSdBegin ()

int	getNextSdpar (int &key, double &par)

void	getSdpar (int layid, int entr, int lr, double par[]) const

bool	getR2tpar ()

TTree *	getR2tTree (int layid) const

int	getXtEntrIndex (double entrance) const

int	getSdEntrIndex (double entrance) const

Public Attributes
int	m_run

Detailed Description

Definition at line 21 of file MdcCalibFunSvc.h.

Constructor & Destructor Documentation

◆ MdcCalibFunSvc()

MdcCalibFunSvc::MdcCalibFunSvc	(	const std::string &	name,
		ISvcLocator *	svcloc )

Definition at line 40 of file MdcCalibFunSvc.cxx.

                                                                       :
     base_class (name, svcloc), m_layInfSig(-1) {
 
     // declare properties
     declareProperty("CheckConst", m_checkConst = false);
     declareProperty("LayerInfSig", m_layInfSig);
     declareProperty("XtMode", m_xtMode = 1);
     declareProperty("NewXtFile", m_xtfile);
     declareProperty("ReadWireEffDb", m_readWireEffDb = true);
     declareProperty("WireEffFile", m_wireEffFile);
     declareProperty("LinearXT", m_linearXT = false);
     declareProperty("FixSigma", m_fixSigma = false);
     declareProperty("FixSigmaValue", m_fixSigmaValue = 130.0); // micron
     m_outputXtMode = true;
}

◆ ~MdcCalibFunSvc()

MdcCalibFunSvc::~MdcCalibFunSvc ( )

Definition at line 56 of file MdcCalibFunSvc.cxx.

56 {

57}

Member Function Documentation

◆ distToDriftTime()

double MdcCalibFunSvc::distToDriftTime	(	double	dist,
		int	layid,
		int	cellid,
		int	lr,
		double	entrance = 0.0 ) const

Definition at line 216 of file MdcCalibFunSvc.cxx.

                                                      {
     int i = 0;
     double time;
     int ord;
     double xtpar[8];
     double dxdtpar[5];
     double x;
     double dxdt;
     double deltax;
 
     int entr = getXtEntrIndex(entrance);
     getXtpar(layid, entr, lr, xtpar);
 
     double tm1 = xtpar[6];
     double tm2 = 2000.0;
     double dm1 = driftTimeToDist(tm1, layid, cellid, lr, entrance);
     double dm2 = driftTimeToDist(tm2, layid, cellid, lr, entrance);
 
     if(dist < 0){
      cout << "Warning in MdcCalibFunSvc: driftDist < 0" << endl;
      time = 0.0;
     } else if(dist < xtpar[0]){
      time = 0.0;
     } else if(dist < dm1){
      for(ord=0; ord<5; ord++){
           dxdtpar[ord] = (double)(ord+1) * xtpar[ord+1];
      }
      time = dist / 0.03;
      while(1){
           if( i > 50 ){
            cout << "Warning in MdcCalibFunSvc: "
             << "can not get the exact value in the dist-to-time conversion." 
             << endl;
            time = dist / 0.03;
            break;
           }
 
           x = 0.0;
           for(ord=0; ord<6; ord++)
            x += xtpar[ord] * pow(time, ord);
 
           deltax = x - dist;
           if( fabs(deltax) < 0.001 ){
            break;
           }
 
           dxdt = 0.0;
           for(ord=0; ord<5; ord++)
            dxdt += dxdtpar[ord] * pow(time, ord);
 
           time = time - (deltax / dxdt);
           i++;
      }
     } else if(dist < dm2){
      time = (dist - dm1) * (tm2 - tm1) / (dm2 - dm1) + tm1;
     } else{
      time = tm2;
     }
 
     if(m_linearXT) time = dist / 0.03;
     return time;
}

Referenced by EsTimeAlg::execute().

◆ driftTimeToDist()

double MdcCalibFunSvc::driftTimeToDist	(	double	drifttime,
		int	layid,
		int	cellid,
		int	lr,
		double	entrance = 0.0 ) const

Definition at line 151 of file MdcCalibFunSvc.cxx.

                                                      {
     double dist;
     if(0 == m_xtMode){
      dist = t2dPoly(drifttime, layid, cellid, lr, entrance);
     } else{
      if((0==lr) || (1==lr)) dist = t2dInter(drifttime, layid, cellid, lr, entrance);
      else{
           double dl = t2dInter(drifttime, layid, cellid, lr, entrance);
           double dr = t2dInter(drifttime, layid, cellid, lr, entrance);
           dist = (dl + dr) * 0.5;
      }
     }
//      cout << setw(15) << drifttime << setw(15) << dist << endl;
     if(m_linearXT) dist = 0.03 * drifttime;
     return dist;
}

Referenced by HoughHit::calDriftDist(), MdcxHit::d(), distToDriftTime(), MdcHit::driftDist(), TrkReco::execute(), and KalFitTrack::getDriftDist().

◆ finalize()

StatusCode MdcCalibFunSvc::finalize ( )

virtual

Definition at line 108 of file MdcCalibFunSvc.cxx.

                                   {
     MsgStream log(messageService(), name());
     log << MSG::INFO << "MdcCalibFunSvc::finalize()" << endreq;
 
     m_xtmap.clear();
     m_t0.clear();
     m_delt0.clear();
     m_qtpar0.clear();
     m_qtpar1.clear();
     m_sdmap.clear();
 
     return StatusCode::SUCCESS;
}

◆ getDelT0()

double MdcCalibFunSvc::getDelT0 ( int wireid ) const

inline

Definition at line 87 of file MdcCalibFunSvc.h.

87{ return m_delt0[wireid]; }

◆ getF()

double MdcCalibFunSvc::getF	(	int	layid,
		int	lr,
		double	dist,
		double	entrance = 0.0,
		double	tanlam = 0.0,
		double	z = 0.0,
		double	Q = 1000.0 ) const

Definition at line 354 of file MdcCalibFunSvc.cxx.

                                          {
 
     return 1.0;
}

◆ getNewXtpar()

bool MdcCalibFunSvc::getNewXtpar ( )

Definition at line 418 of file MdcCalibFunSvc.cxx.

                                 {
     MsgStream log(messageService(), name());
     log << MSG::INFO << "read calib const from TCDS" << endreq;
 
     for(int layid=0; layid<NLAYER; layid++){
      for(int entr=0; entr<NXTENTR; entr++){
           for(int lr=0; lr<2; lr++){
            double br_t,br_d;
            TTree* newXtTree = getNewXtparTree(layid,entr,lr);
            if(!newXtTree) return false;
            newXtTree -> SetBranchAddress("t", &br_t);
            newXtTree -> SetBranchAddress("d", &br_d);
            int nEntries = newXtTree -> GetEntries();
            if((nEntries<10) || (nEntries>=200)){
             log << MSG::ERROR << "wrong X-T constants: layer " << layid
                 << ", iEntr " << entr << ", lr " << lr << endreq;
             return false;
            }
            m_nNewXt[layid][entr][lr] = nEntries;
            for(int i=0; i<nEntries; i++){
             newXtTree->GetEntry(i);
             m_vt[layid][entr][lr][i] = br_t;
             m_vd[layid][entr][lr][i] = br_d;
            }//end loop entries
           }//end lr
      }//end entr
     }//end layid
 
     return true;
}

◆ getNewXtparTree()

TTree * MdcCalibFunSvc::getNewXtparTree	(	int	layid,
		int	entr,
		int	lr ) const

Definition at line 449 of file MdcCalibFunSvc.cxx.

                                                                       {
     MsgStream log(messageService(), name());
     string fullPath = "/Calib/MdcCal";
     SmartDataPtr<CalibData::MdcCalibData> calConst(m_pCalDataSvc, fullPath);
     if( ! calConst ){
      log << MSG::ERROR << "can not get MdcCalibConst via SmartPtr" << endreq;
      return NULL;
     }
 
     TTree* newXtTree = calConst->getNewXtpar(layid,entr,lr);
     return newXtTree;
}

Referenced by getNewXtpar().

◆ getNextSdpar()

int MdcCalibFunSvc::getNextSdpar	(	int &	key,
		double &	par )

Definition at line 558 of file MdcCalibFunSvc.cxx.

                                                     {
     if( m_sditer != m_sdmap.end() ){
      key = (*m_sditer).first;
      par = (*m_sditer).second;
      m_sditer++;
      return 1;
     }
     else return 0;
}

◆ getNextXtpar()

int MdcCalibFunSvc::getNextXtpar	(	int &	key,
		double &	par )

Definition at line 400 of file MdcCalibFunSvc.cxx.

                                                     {
     if( m_xtiter != m_xtmap.end() ){
      key = (*m_xtiter).first;
      par = (*m_xtiter).second;
      m_xtiter++;
      return 1;
     }
     else return 0;
}

◆ getQtpar()

double MdcCalibFunSvc::getQtpar	(	int	layid,
		int	ord ) const

Definition at line 538 of file MdcCalibFunSvc.cxx.

                                                        {
     if(0 == ord)       return m_qtpar0[layid];
     else if(1 == ord)  return m_qtpar1[layid];
     else {
      cout << "wrong order number" << endl;
      return 0.0;
     }
}

Referenced by getTimeWalk().

◆ getR2tpar()

bool MdcCalibFunSvc::getR2tpar ( )

Definition at line 462 of file MdcCalibFunSvc.cxx.

                               {
     for(int layid=0; layid<NLAYER; layid++){
      int br_iEntr,br_lr;
      double br_s,br_t;
      TTree* r2tTree = getR2tTree(layid);
      if(!r2tTree) return false;
      r2tTree -> SetBranchAddress("iEntr", &br_iEntr);
      r2tTree -> SetBranchAddress("lr", &br_lr);
      r2tTree -> SetBranchAddress("s", &br_s);
      r2tTree -> SetBranchAddress("t", &br_t);
      int nEntries = r2tTree -> GetEntries();
      for(int i=0; i<nEntries; i++){
           r2tTree->GetEntry(i);
           int bin = m_nR2t[layid][br_iEntr][br_lr];
           if(bin>=200){
            cout << "Error: number of sigma-time bins overflow" << endl;
            return false;
           }
           m_tR2t[layid][br_iEntr][br_lr][bin] = br_t;
           m_sR2t[layid][br_iEntr][br_lr][bin] = br_s;
           m_nR2t[layid][br_iEntr][br_lr]++;
      }
     }
     for(int layid=0; layid<NLAYER; layid++){
      for(int iEntr=0; iEntr<NXTENTR; iEntr++){
           for(int lr=0; lr<2; lr++){
            if((m_nR2t[layid][iEntr][lr]<10) || (m_nR2t[layid][iEntr][lr]>=200)) return false;
           }
      }
     }
     return true;
}

◆ getR2tTree()

TTree * MdcCalibFunSvc::getR2tTree ( int layid ) const

Definition at line 495 of file MdcCalibFunSvc.cxx.

                                                {
     MsgStream log(messageService(), name());
     string fullPath = "/Calib/MdcCal";
     SmartDataPtr<CalibData::MdcCalibData> calConst(m_pCalDataSvc, fullPath);
     if( ! calConst ){
      log << MSG::ERROR << "can not get MdcCalibConst via SmartPtr" << endreq;
      return NULL;
     }
 
     TTree* r2tTree = calConst->getR2tpar(layid);
     return r2tTree;
}

Referenced by getR2tpar().

◆ getSdEntrIndex()

int MdcCalibFunSvc::getSdEntrIndex ( double entrance ) const

Definition at line 599 of file MdcCalibFunSvc.cxx.

                                                        {
     int i;
     int index;
     int idmax = 5;
     double aglpi = 3.141592653;
     double aglmin = -1.570796327; // -90 degree
     double aglmax = 1.570796327; // 90 degree
     double delAngle = 0.523598776; // 30 degree
 
     MsgStream log(messageService(), name());
     if(entrance < aglmin){
      log << MSG::WARNING << "entrance angle < -pi/2" << endreq;
      while(1){
           entrance += aglpi;
           if(entrance >= aglmin) break;
      }
     } else if(entrance > aglmax){
      log << MSG::WARNING << "entrance angle > pi/2" << endreq;
      while(1){
           entrance -= aglpi;
           if(entrance <= aglmax) break;
      }
     }
 
     index = (int)((entrance-aglmin) / delAngle);
     if(index < 0) index = 0;
     else if(index > idmax) index = idmax;
 
     return index;
}

Referenced by getSigmaLR().

◆ getSdpar()

void MdcCalibFunSvc::getSdpar	(	int	layid,
		int	entr,
		int	lr,
		double	par[] ) const

Definition at line 547 of file MdcCalibFunSvc.cxx.

                                                                            {
     int parId;
     if( (entr < 0) || (entr >= 18) ){
      entr = 17;
     }
     for(int bin=0; bin<NSDBIN; bin++){
      parId = getSdparId(layid, entr, lr, bin);
      par[bin] = m_sdpar[parId];
     }
}

Referenced by getSigmaLR().

◆ getSigma()

double MdcCalibFunSvc::getSigma	(	int	layid,
		int	lr,
		double	dist,
		double	entrance = 0.0,
		double	tanlam = 0.0,
		double	z = 0.0,
		double	Q = 1000.0 ) const

Definition at line 280 of file MdcCalibFunSvc.cxx.

                                          {
     double sigma;
     if( (0 == lr) || (1 == lr) ){
      sigma = getSigmaLR(layid, lr, dist, entrance, tanlam, z, Q);
     } else{
      double sl = getSigmaLR(layid, 0, dist, entrance, tanlam, z, Q);
      double sr = getSigmaLR(layid, 1, dist, entrance, tanlam, z, Q);
      sigma = (sl + sr) * 0.5;
     }
 
     if(m_fixSigma) sigma = 0.001 * m_fixSigmaValue;
     if(layid == m_layInfSig) sigma = 9999.0;
     return sigma;
}

Referenced by TrkReco::execute(), KalFitTrack::getSigma(), getSigma1(), and MdcHit::sigma().

◆ getSigma1()

double MdcCalibFunSvc::getSigma1	(	int	layid,
		int	lr,
		double	dist,
		double	entrance = 0.0,
		double	tanlam = 0.0,
		double	z = 0.0,
		double	Q = 1000.0 ) const

Definition at line 340 of file MdcCalibFunSvc.cxx.

                                           {
     double sigma1 = getSigma(layid, lr, dist, entrance, tanlam, z, Q);
     return sigma1;
}

◆ getSigma2()

double MdcCalibFunSvc::getSigma2	(	int	layid,
		int	lr,
		double	dist,
		double	entrance = 0.0,
		double	tanlam = 0.0,
		double	z = 0.0,
		double	Q = 1000.0 ) const

Definition at line 347 of file MdcCalibFunSvc.cxx.

                                           {
 
     return 0.0;
}

◆ getSigmaLR()

double MdcCalibFunSvc::getSigmaLR	(	int	layid,
		int	lr,
		double	dist,
		double	entrance = 0.0,
		double	tanlam = 0.0,
		double	z = 0.0,
		double	Q = 1000.0 ) const

Definition at line 297 of file MdcCalibFunSvc.cxx.

                                            {
 
     double sigma = 9999.0;
     double par[NSDBIN];
 
     int entr = getSdEntrIndex(entrance);
     getSdpar(layid, entr, lr, par);
 
     int nmaxBin;
     double dw = 0.5;       // width of each distance bin
     double dmin = 0.25;    // mm
     if(layid < 8){
      nmaxBin = 20; // 11->20 2011-12-10
     } else{
      nmaxBin = 20; // 15->20 2011-12-10
     }
 
     double dref[2];
     double distAbs = fabs(dist);
     if(distAbs < dmin){
      sigma = par[0];
     } else{
      int bin = (int)((distAbs - dmin) / dw);
      if(bin >= nmaxBin){
           sigma = par[nmaxBin];
      } else if(bin < 0){
           sigma = par[0];
      } else{
           dref[0] = (double)bin * dw + 0.25;
           dref[1] = (double)(bin+1) * dw + 0.25;
           if((dref[1] - dref[0]) <= 0){
            sigma = 9999.0;
           } else{
            sigma = (par[bin+1] - par[bin]) * (distAbs - dref[0]) /
             (dref[1] - dref[0]) + par[bin];
           }
      }
     }
     return sigma;
}

Referenced by getSigma().

◆ getSigmaToT()

double MdcCalibFunSvc::getSigmaToT	(	int	layid,
		int	lr,
		double	tdr,
		double	entrance = 0.0,
		double	tanlam = 0.0,
		double	z = 0.0,
		double	Q = 1000.0 ) const

Definition at line 361 of file MdcCalibFunSvc.cxx.

                                                           {
     if(!m_fgR2t){
      cout << "ERROR: can not get sigma-time functions" << endl;
      return -999.0;
     } else if( (0 == lr) || (1 == lr) ){
      return getSigmaToTLR(layid, lr, tdr, entrance, tanlam, z, Q);
     } else{
      double sl = getSigmaToTLR(layid, 0, tdr, entrance, tanlam, z, Q);
      double sr = getSigmaToTLR(layid, 1, tdr, entrance, tanlam, z, Q);
      double sigma = (sl + sr) * 0.5;
      return sigma;
     }
}

◆ getSigmaToTLR()

double MdcCalibFunSvc::getSigmaToTLR	(	int	layid,
		int	lr,
		double	tdr,
		double	entrance = 0.0,
		double	tanlam = 0.0,
		double	z = 0.0,
		double	Q = 1000.0 ) const

Definition at line 376 of file MdcCalibFunSvc.cxx.

                                                             {
     double sigma;
     int iEntr = getXtEntrIndex(entrance);
     int nBin = m_nR2t[layid][iEntr][lr];
     if(tdr < m_tR2t[layid][iEntr][lr][0]){
      sigma = m_sR2t[layid][iEntr][lr][0];
     } else if(tdr < m_tR2t[layid][iEntr][lr][nBin-1]){
      for(int i=0; i<(nBin-1); i++){
           if((tdr>=m_tR2t[layid][iEntr][lr][i]) && (tdr<m_tR2t[layid][iEntr][lr][i+1])){
            double t1 = m_tR2t[layid][iEntr][lr][i];
            double t2 = m_tR2t[layid][iEntr][lr][i+1];
            double s1 = m_sR2t[layid][iEntr][lr][i];
            double s2 = m_sR2t[layid][iEntr][lr][i+1];
            sigma = (tdr-t1) * (s2-s1) / (t2-t1) + s1;
            break;
           }
      }
     } else{
      sigma = m_sR2t[layid][iEntr][lr][nBin-1];
     }
     return sigma;
}

Referenced by getSigmaToT().

◆ getT0() [1/2]

double MdcCalibFunSvc::getT0	(	int	layid,
		int	cellid ) const

Definition at line 509 of file MdcCalibFunSvc.cxx.

                                                        {
     int wireid = m_pMdcGeomSvc->Wire(layid, cellid)->Id();
     double t0 = getT0(wireid);
 
     return t0;
}

Referenced by HoughHit::driftTime(), HoughHit::driftTime(), EsTimeAlg::execute(), TrkReco::execute(), Hough2D::fit(), Hough3D::fit(), KalFitTrack::getDriftTime(), getT0(), MdcxHit::process(), and MdcHit::setCalibSvc().

◆ getT0() [2/2]

double MdcCalibFunSvc::getT0 ( int wireid ) const

inline

Definition at line 86 of file MdcCalibFunSvc.h.

86{ return m_t0[wireid]; }

◆ getTimeWalk()

double MdcCalibFunSvc::getTimeWalk	(	int	layid,
		double	Q ) const

Definition at line 516 of file MdcCalibFunSvc.cxx.

                                                            {
     double tw = 0.0;
     double qtpar[2];
     int ord;
 
     if(Q < 0.0001) Q = 0.0001;
 
     for(ord=0; ord<2; ord++){
      qtpar[ord] = getQtpar(layid, ord);
     }
 
     tw = qtpar[0] + qtpar[1] / sqrt( Q );
     if(m_run < 0) tw = 0.0;    // for MC
 
     return tw;
}

Referenced by HoughHit::driftTime(), HoughHit::driftTime(), EsTimeAlg::execute(), TrkReco::execute(), Hough2D::fit(), Hough3D::fit(), KalFitTrack::getDriftTime(), MdcxHit::process(), and MdcHit::setCalibSvc().

◆ getTprop()

double MdcCalibFunSvc::getTprop	(	int	lay,
		double	z ) const

Definition at line 145 of file MdcCalibFunSvc.cxx.

                                                      {
     double vp = getVprop(lay);
     double tp = fabs(z - m_zst[lay]) / vp;
     return tp;
}

Referenced by HoughHit::driftTime(), HoughHit::driftTime(), MdcHit::driftTime(), and MdcxHit::tcor().

◆ getVprop()

double MdcCalibFunSvc::getVprop ( int lay ) const

inline

Definition at line 215 of file MdcCalibFunSvc.h.

                                                   {
     if(lay<8) return 220.0;
     else return 240.0;
}

Referenced by getTprop().

◆ getWireEff() [1/2]

double MdcCalibFunSvc::getWireEff	(	int	layid,
		int	cellid ) const

Definition at line 533 of file MdcCalibFunSvc.cxx.

                                                             {
     int wireid = m_pMdcGeomSvc->Wire(layid, cellid)->Id();
     return m_wireEff[wireid];
}

◆ getWireEff() [2/2]

double MdcCalibFunSvc::getWireEff ( int wireid ) const

inline

Definition at line 93 of file MdcCalibFunSvc.h.

93{ return m_wireEff[wireid]; }

◆ getXtEntrIndex()

int MdcCalibFunSvc::getXtEntrIndex ( double entrance ) const

Definition at line 568 of file MdcCalibFunSvc.cxx.

                                                        {
     int i;
     int index;
     int idmax = 17;
     double aglpi = 3.141592653;
     double aglmin = -1.570796327; // -90 degree
     double aglmax = 1.570796327; // 90 degree
     double delAngle = 0.174532925; // 10 degree
 
     MsgStream log(messageService(), name());
     if(entrance < aglmin){
      log << MSG::WARNING << "entrance angle < -pi/2" << endreq;
      while(1){
           entrance += aglpi;
           if(entrance >= aglmin) break;
      }
     } else if(entrance > aglmax){
      log << MSG::WARNING << "entrance angle > pi/2" << endreq;
      while(1){
           entrance -= aglpi;
           if(entrance <= aglmax) break;
      }
     }
 
     index = (int)((entrance-aglmin) / delAngle);
     if(index < 0) index = 0;
     else if(index > idmax) index = idmax;
 
     return index;
}

Referenced by distToDriftTime(), and getSigmaToTLR().

◆ getXtpar()

void MdcCalibFunSvc::getXtpar	(	int	layid,
		int	entr,
		int	lr,
		double	par[] ) const

Definition at line 410 of file MdcCalibFunSvc.cxx.

                                                                            {
     int parId;
     for(int ord=0; ord<8; ord++){
      parId = getXtparId(layid, entr, lr, ord);
      par[ord] = m_xtpar[parId];
     }
}

Referenced by distToDriftTime().

◆ handle()

void MdcCalibFunSvc::handle ( const Incident & inc )

Definition at line 122 of file MdcCalibFunSvc.cxx.

                                              {
     MsgStream log( messageService(), name() );
     log << MSG::DEBUG << "handle: " << inc.type() << endreq;
 
     if ( inc.type() == "NewRun" ){
      log << MSG::DEBUG << "NewRun" << endreq;
 
      if( ! initCalibConst() ){
           log << MSG::ERROR 
           << "can not initilize Mdc Calib Constants" << endl
           << "  Please insert the following statement "
           << "in your \"jobOption.txt\" "
           << "before the include file of Mdc Reconstruction: "
           << endl << "        "
           << "#include \"$CALIBSVCROOT/share/job-CalibData.txt\""
           << endl
           << "  If still error, please contact with Wu Linghui "
           << "([email protected])."
           << endreq;
      }
     } 
}

◆ initialize()

StatusCode MdcCalibFunSvc::initialize ( )

virtual

Definition at line 68 of file MdcCalibFunSvc.cxx.

                                     {
     MsgStream log(messageService(), name());
     log << MSG::INFO << "MdcCalibFunSvc::initialize()" << endreq;
 
     StatusCode sc = Service::initialize();
     if( sc.isFailure() ) return sc;
 
     IIncidentSvc* incsvc;
     sc = service("IncidentSvc", incsvc);
     int priority = 100;
     if( sc.isSuccess() ){
      incsvc -> addListener(this, "NewRun", priority);
     }
 
     sc = service("CalibDataSvc", m_pCalDataSvc, true);
     if( sc == StatusCode::SUCCESS ){
      log << MSG::INFO << "Retrieve IDataProviderSvc" << endreq;
     }else{
      log << MSG::FATAL << "can not get IDataProviderSvc" << endreq;
     }
 
     sc = service("MdcGeomSvc", m_pMdcGeomSvc);
     if( sc != StatusCode::SUCCESS ){
      log << MSG::ERROR << "can not use MdcGeomSvc" << endreq;
      return StatusCode::FAILURE;
     }
 
     if(m_fixSigma) cout << "Fix MDC sigma to " << m_fixSigmaValue << " micron." << endl;
 
     m_updateNum = 0;
     for(int wir=0; wir<6796; wir++) m_wireEff[wir] = 1.0;
     for(int lay=0; lay<NLAYER; lay++){
      for(int iEntr=0; iEntr<NXTENTR; iEntr++){
           for(int lr=0; lr<2; lr++) m_nR2t[lay][iEntr][lr] = 0;
      }
     }
 
     return StatusCode::SUCCESS;
}

◆ setSdBegin()

void MdcCalibFunSvc::setSdBegin ( )

inline

Definition at line 95 of file MdcCalibFunSvc.h.

95{ m_sditer = m_sdmap.begin(); }

◆ setXtBegin()

void MdcCalibFunSvc::setXtBegin ( )

inline

Definition at line 79 of file MdcCalibFunSvc.h.

79{ m_xtiter = m_xtmap.begin(); }

Member Data Documentation

◆ m_run

int MdcCalibFunSvc::m_run

Definition at line 34 of file MdcCalibFunSvc.h.

Referenced by getTimeWalk().

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

Public Member Functions

Public Attributes

Detailed Description

Constructor & Destructor Documentation

◆ MdcCalibFunSvc()

◆ ~MdcCalibFunSvc()

Member Function Documentation

◆ distToDriftTime()

◆ driftTimeToDist()

◆ finalize()

◆ getDelT0()

◆ getF()

◆ getNewXtpar()

◆ getNewXtparTree()

◆ getNextSdpar()

◆ getNextXtpar()

◆ getQtpar()

◆ getR2tpar()

◆ getR2tTree()

◆ getSdEntrIndex()

◆ getSdpar()

◆ getSigma()

◆ getSigma1()

◆ getSigma2()

◆ getSigmaLR()

◆ getSigmaToT()

◆ getSigmaToTLR()

◆ getT0() [1/2]

◆ getT0() [2/2]

◆ getTimeWalk()

◆ getTprop()

◆ getVprop()

◆ getWireEff() [1/2]

◆ getWireEff() [2/2]

◆ getXtEntrIndex()

◆ getXtpar()

◆ handle()

◆ initialize()

◆ setSdBegin()

◆ setXtBegin()

Member Data Documentation

◆ m_run