MdcCalibFunSvc Class Reference

#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 34 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;
}

Referenced by MdcCalibFunSvc().

~MdcCalibFunSvc()

MdcCalibFunSvc::~MdcCalibFunSvc

(

)

Definition at line 50 of file MdcCalibFunSvc.cxx.

{}

Member Function Documentation

distToDriftTime()

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

Definition at line 210 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;
}

driftTimeToDist()

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

Definition at line 142 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 distToDriftTime().

finalize()

StatusCode MdcCalibFunSvc::finalize ( )

virtual

Definition at line 100 of file MdcCalibFunSvc.cxx.

                                    {
  MsgStream log( msgSvc(), name() );
  log << MSG::INFO << "MdcCalibFunSvc::finalize()" << endmsg;
 
  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 80 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

Definition at line 342 of file MdcCalibFunSvc.cxx.

                                                        {
 
  return 1.0;
}

getNewXtpar()

bool MdcCalibFunSvc::getNewXtpar

(

)

Definition at line 413 of file MdcCalibFunSvc.cxx.

                                 {
  MsgStream log( msgSvc(), name() );
  log << MSG::INFO << "read calib const from TCDS" << endmsg;
 
  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 << endmsg;
          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( msgSvc(), name() );
  string                                fullPath = "/Calib/MdcCal";
  SmartDataPtr<CalibData::MdcCalibData> calConst( m_pCalDataSvc, fullPath );
  if ( !calConst )
  {
    log << MSG::ERROR << "can not get MdcCalibConst via SmartPtr" << endmsg;
    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 564 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 393 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 544 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 463 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 503 of file MdcCalibFunSvc.cxx.

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

Referenced by getR2tpar().

getSdEntrIndex()

int MdcCalibFunSvc::getSdEntrIndex

(

double

entrance

)

const

Definition at line 611 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( msgSvc(), name() );
  if ( entrance < aglmin )
  {
    log << MSG::WARNING << "entrance angle < -pi/2" << endmsg;
    while ( 1 )
    {
      entrance += aglpi;
      if ( entrance >= aglmin ) break;
    }
  }
  else if ( entrance > aglmax )
  {
    log << MSG::WARNING << "entrance angle > pi/2" << endmsg;
    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 554 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 269 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 getSigma1().

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 330 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 336 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 286 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 348 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 366 of file MdcCalibFunSvc.cxx.

                                                                                {
  double sigma=0;
  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 517 of file MdcCalibFunSvc.cxx.

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

Referenced by getT0().

getT0() [2/2]

double MdcCalibFunSvc::getT0 ( int wireid ) const

inline

Definition at line 79 of file MdcCalibFunSvc.h.

{ return m_t0[wireid]; }

getTimeWalk()

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

Definition at line 524 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;
}

getTprop()

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

Definition at line 136 of file MdcCalibFunSvc.cxx.

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

getVprop()

double MdcCalibFunSvc::getVprop ( int lay ) const

inline

Definition at line 207 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 539 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 86 of file MdcCalibFunSvc.h.

{ return m_wireEff[wireid]; }

getXtEntrIndex()

int MdcCalibFunSvc::getXtEntrIndex

(

double

entrance

)

const

Definition at line 575 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( msgSvc(), name() );
  if ( entrance < aglmin )
  {
    log << MSG::WARNING << "entrance angle < -pi/2" << endmsg;
    while ( 1 )
    {
      entrance += aglpi;
      if ( entrance >= aglmin ) break;
    }
  }
  else if ( entrance > aglmax )
  {
    log << MSG::WARNING << "entrance angle > pi/2" << endmsg;
    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 404 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 114 of file MdcCalibFunSvc.cxx.

                                                 {
  MsgStream log( msgSvc(), name() );
  log << MSG::DEBUG << "handle: " << inc.type() << endmsg;
 
  if ( inc.type() == "NewRun" )
  {
    log << MSG::DEBUG << "NewRun" << endmsg;
 
    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 "
          << "(wulh@mail.ihep.ac.cn)." << endmsg;
    }
  }
}

initialize()

StatusCode MdcCalibFunSvc::initialize ( )

virtual

Definition at line 61 of file MdcCalibFunSvc.cxx.

                                      {
  MsgStream log( msgSvc(), name() );
  log << MSG::INFO << "MdcCalibFunSvc::initialize()" << endmsg;
 
  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" << endmsg; }
  else { log << MSG::FATAL << "can not get IDataProviderSvc" << endmsg; }
 
  sc = service( "MdcGeomSvc", m_pMdcGeomSvc );
  if ( sc != StatusCode::SUCCESS )
  {
    log << MSG::ERROR << "can not use MdcGeomSvc" << endmsg;
    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 88 of file MdcCalibFunSvc.h.

{ m_sditer = m_sdmap.begin(); }

setXtBegin()

void MdcCalibFunSvc::setXtBegin ( )

inline

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