MdcCalibFunSvc Class Reference

#include <MdcCalibFunSvc.h>

Inheritance diagram for MdcCalibFunSvc:

Public Member Functions
	MdcCalibFunSvc (const std::string &name, ISvcLocator *svcloc)
	~MdcCalibFunSvc ()
virtual StatusCode	queryInterface (const InterfaceID &riid, void **ppvUnknown)
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 Member Functions inherited from IMdcCalibFunSvc

Public Attributes
int	m_run

Additional Inherited Members
Static Public Member Functions inherited from IMdcCalibFunSvc
static const InterfaceID &	interfaceID ()

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 38 of file MdcCalibFunSvc.cxx.

                                                                       :
     Service (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 54 of file MdcCalibFunSvc.cxx.

                               {
}

Member Function Documentation

distToDriftTime()

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

virtual

Implements IMdcCalibFunSvc.

Definition at line 214 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

virtual

Implements IMdcCalibFunSvc.

Definition at line 149 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 DotsHelixFitter::calculateNewHelix(), DotsHelixFitter::calculateRoughDD(), MdcxHit::d(), distToDriftTime(), MdcHit::driftDist(), TrkReco::execute(), KalFitTrack::getDriftDist(), and DotsHelixFitter::updateDcDigiInfo().

finalize()

StatusCode MdcCalibFunSvc::finalize ( )

virtual

Definition at line 106 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

inlinevirtual

Implements IMdcCalibFunSvc.

Definition at line 87 of file MdcCalibFunSvc.h.

{ 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

virtual

Implements IMdcCalibFunSvc.

Definition at line 352 of file MdcCalibFunSvc.cxx.

                                          {
 
     return 1.0;
}

getNewXtpar()

bool MdcCalibFunSvc::getNewXtpar ( )

virtual

Implements IMdcCalibFunSvc.

Definition at line 416 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

virtual

Implements IMdcCalibFunSvc.

Definition at line 447 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 )

virtual

Implements IMdcCalibFunSvc.

Definition at line 556 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 )

virtual

Implements IMdcCalibFunSvc.

Definition at line 398 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

virtual

Implements IMdcCalibFunSvc.

Definition at line 536 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 ( )

virtual

Implements IMdcCalibFunSvc.

Definition at line 460 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

virtual

Implements IMdcCalibFunSvc.

Definition at line 493 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

virtual

Implements IMdcCalibFunSvc.

Definition at line 597 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

virtual

Implements IMdcCalibFunSvc.

Definition at line 545 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

virtual

Implements IMdcCalibFunSvc.

Definition at line 278 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 DotsHelixFitter::calculateNewHelix(), TrkReco::execute(), KalFitTrack::getSigma(), getSigma1(), MdcHit::sigma(), and DotsHelixFitter::updateDcDigiInfo().

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

virtual

Implements IMdcCalibFunSvc.

Definition at line 338 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

virtual

Implements IMdcCalibFunSvc.

Definition at line 345 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

virtual

Implements IMdcCalibFunSvc.

Definition at line 295 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

virtual

Implements IMdcCalibFunSvc.

Definition at line 359 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

virtual

Implements IMdcCalibFunSvc.

Definition at line 374 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

virtual

Implements IMdcCalibFunSvc.

Definition at line 507 of file MdcCalibFunSvc.cxx.

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

Referenced by DotsHelixFitter::calculateNewHelix(), DotsHelixFitter::calculateRoughDD(), HoughHit::driftTime(), EsTimeAlg::execute(), TrkReco::execute(), KalFitTrack::getDriftTime(), getT0(), MdcxHit::process(), MdcHit::setCalibSvc(), and DotsHelixFitter::updateDcDigiInfo().

getT0() [2/2]

double MdcCalibFunSvc::getT0 ( int wireid ) const

inlinevirtual

Implements IMdcCalibFunSvc.

Definition at line 86 of file MdcCalibFunSvc.h.

{ return m_t0[wireid]; }

getTimeWalk()

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

virtual

Implements IMdcCalibFunSvc.

Definition at line 514 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 DotsHelixFitter::calculateNewHelix(), DotsHelixFitter::calculateRoughDD(), HoughHit::driftTime(), EsTimeAlg::execute(), TrkReco::execute(), KalFitTrack::getDriftTime(), MdcxHit::process(), MdcHit::setCalibSvc(), and DotsHelixFitter::updateDcDigiInfo().

getTprop()

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

virtual

Implements IMdcCalibFunSvc.

Definition at line 143 of file MdcCalibFunSvc.cxx.

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

Referenced by DotsHelixFitter::calculateNewHelix(), DotsHelixFitter::calculateRoughDD(), HoughHit::driftTime(), MdcHit::driftTime(), MdcxHit::tcor(), and DotsHelixFitter::updateDcDigiInfo().

getVprop()

double MdcCalibFunSvc::getVprop ( int lay ) const

inlinevirtual

Implements IMdcCalibFunSvc.

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

virtual

Implements IMdcCalibFunSvc.

Definition at line 531 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

inlinevirtual

Implements IMdcCalibFunSvc.

Definition at line 93 of file MdcCalibFunSvc.h.

{ return m_wireEff[wireid]; }

getXtEntrIndex()

int MdcCalibFunSvc::getXtEntrIndex ( double entrance ) const

virtual

Implements IMdcCalibFunSvc.

Definition at line 566 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

virtual

Implements IMdcCalibFunSvc.

Definition at line 408 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 120 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 66 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;
}

queryInterface()

StatusCode MdcCalibFunSvc::queryInterface ( const InterfaceID & riid, void ** ppvUnknown )

virtual

Definition at line 57 of file MdcCalibFunSvc.cxx.

                                                                                     {
     if( IID_IMdcCalibFunSvc.versionMatch(riid) ){
      *ppvInterface = static_cast<IMdcCalibFunSvc*> (this);
     } else{
      return Service::queryInterface(riid, ppvInterface);
     }
     return StatusCode::SUCCESS;
}

setSdBegin()

void MdcCalibFunSvc::setSdBegin ( )

inlinevirtual

Implements IMdcCalibFunSvc.

Definition at line 95 of file MdcCalibFunSvc.h.

{ m_sditer = m_sdmap.begin(); }

setXtBegin()

void MdcCalibFunSvc::setXtBegin ( )

inlinevirtual

Implements IMdcCalibFunSvc.

Definition at line 79 of file MdcCalibFunSvc.h.

{ 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:

• MdcCalibFunSvc.h
• MdcCalibFunSvc.cxx