CgemCalibFunSvc Class Reference

#include <CgemCalibFunSvc.h>

Inheritance diagram for CgemCalibFunSvc:

Public Member Functions
	CgemCalibFunSvc (const std::string &name, ISvcLocator *svcloc)
	~CgemCalibFunSvc ()
virtual StatusCode	queryInterface (const InterfaceID &riid, void **ppvUnknown)
virtual StatusCode	initialize ()
virtual StatusCode	finalize ()
void	handle (const Incident &)
double	getSigma (int layer, int xvFlag, int readoutMode, double angle, double Q, double T) const
double	getTimeRising (int layer, int xvFlag, int sheet, int stripID, double Q=100., double z=0.) const
double	getTimeFalling (int layer, int xvFlag, int sheet, int stripID, double Q=100., double z=0.) const
double	getTimeWalk (int layer, int xvFlag, int sheet, int stripID, double Q) const
double	getTimeWalk (double Q, double threshold) const
Public Member Functions inherited from ICgemCalibFunSvc

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

Detailed Description

Definition at line 23 of file CgemCalibFunSvc.h.

Constructor & Destructor Documentation

CgemCalibFunSvc()

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

Definition at line 37 of file CgemCalibFunSvc.cxx.

                                                                         :Service (name, svcloc){
     declareProperty("PrintInfo", m_printInfo=false);
     declareProperty("TimeFitFile", m_timeFitFile);
     declareProperty("TimeFitFilePath", m_timeFitFilePath);
     declareProperty("TimeWalkCalibFile", m_timeWalkFileName);
     declareProperty("LUTFile", m_fileLUT);
}

~CgemCalibFunSvc()

CgemCalibFunSvc::~CgemCalibFunSvc

(

)

Definition at line 45 of file CgemCalibFunSvc.cxx.

                                 {
}

Member Function Documentation

finalize()

StatusCode CgemCalibFunSvc::finalize ( )

virtual

Definition at line 118 of file CgemCalibFunSvc.cxx.

                                    {
     MsgStream log(messageService(), name());
     log << MSG::INFO << "CgemCalibFunSvc::finalize()" << endreq;
 
     return StatusCode::SUCCESS;
}

getSigma()

double CgemCalibFunSvc::getSigma ( int layer, int xvFlag, int readoutMode, double angle, double Q, double T ) const

virtual

Implements ICgemCalibFunSvc.

Definition at line 271 of file CgemCalibFunSvc.cxx.

                                                                                                              {
     double sigma = 0.13;   // unit: mm
     return sigma;
}

Referenced by CgemHitOnTrack::getFitStuff(), and CgemHitOnTrack::getFitStuff().

getTimeFalling()

double CgemCalibFunSvc::getTimeFalling ( int layer, int xvFlag, int sheet, int stripID, double Q = 100., double z = 0. ) const

virtual

Implements ICgemCalibFunSvc.

Definition at line 280 of file CgemCalibFunSvc.cxx.

                                                                                                             {
     return m_tFalling[layer][xvFlag][sheet];
}

getTimeRising()

double CgemCalibFunSvc::getTimeRising ( int layer, int xvFlag, int sheet, int stripID, double Q = 100., double z = 0. ) const

virtual

Implements ICgemCalibFunSvc.

Definition at line 276 of file CgemCalibFunSvc.cxx.

                                                                                                            {
     return m_tRising[layer][xvFlag][sheet];
}

getTimeWalk() [1/2]

double CgemCalibFunSvc::getTimeWalk ( double Q, double threshold ) const

virtual

Implements ICgemCalibFunSvc.

Definition at line 300 of file CgemCalibFunSvc.cxx.

                                                                   {
     MsgStream log(messageService(), name());
 
     double tw = 0.0;
     int thresholdId = -1;
     double minval = 1000.;
     for(int i=0; i<m_valThreshold.size(); i++){
      double val = fabs(threshold-m_valThreshold[i]);
      if(val<minval){
           minval = val;
           thresholdId = i;
      }
     }
 
     int binMin = m_thIdFrom[thresholdId];
     int binMax = m_thIdEnd[thresholdId];
     double qMin = m_charge[binMin];
     double qMax = m_charge[binMax];
 
     if(minval>1.0){
      log << MSG::WARNING << "can not find time walk calibration constants for threshold="
          << threshold << " , use the parameters with threshold=" << m_valThreshold[thresholdId] << endreq;
      if(m_printInfo){
           cout << "threshold_binId = " << thresholdId << "   qMin = " << qMin
            << "    qMax = " << qMax << endl;
      }
     }
 
     if(Q < 0.0){
      log << MSG::WARNING << "charge<0, calculate time walk with charge=0" << endreq;
      tw = (0.0-m_charge[binMin])*(m_tw[binMin]-m_tw[binMin+1])/(m_charge[binMin]-m_charge[binMin+1]) + m_tw[binMin];
     } else if(Q < qMin){
      log << MSG::WARNING << "charge < minimumCharge in time walk calibration function " << qMin << endreq;
      tw = (Q-m_charge[binMin])*(m_tw[binMin]-m_tw[binMin+1])/(m_charge[binMin]-m_charge[binMin+1]) + m_tw[binMin];
     } else if(Q > qMax){
      log << MSG::WARNING << "charge > maximumCharge in time walk calibration function " << qMax << endreq;
      if(m_printInfo) cout << "binMin = " << binMin << "   binMax = " << binMax << endl;
      tw = (Q-m_charge[binMax])*(m_tw[binMax]-m_tw[binMax-1])/(m_charge[binMax]-m_charge[binMax-1]) + m_tw[binMax];
     } else{
      for(int i=binMin; i<binMax; i++){
           if((Q>=m_charge[i]) && (Q<m_charge[i+1])){
            tw = (m_tw[i+1]-m_tw[i])*(Q-m_charge[i])/(m_charge[i+1]-m_charge[i]) + m_tw[i];
           }
      }
     }
 
     return tw;
}

getTimeWalk() [2/2]

double CgemCalibFunSvc::getTimeWalk ( int layer, int xvFlag, int sheet, int stripID, double Q ) const

virtual

Implements ICgemCalibFunSvc.

Definition at line 284 of file CgemCalibFunSvc.cxx.

                                                                                                {
     double tw = 0.0;
     if(!m_readLUT){
      cout << "ERROR: can not calculate time walk with LUT, return tw=0." << endl;
      return tw;
     }
     double threshold = m_cgemLUTReader->Get_thr_T_fC(layer, sheet, xvFlag, stripID);
     if(m_printInfo){
      cout << setw(8)<<layer <<setw(8)<<sheet<<setw(8)<<xvFlag<<setw(8)<<stripID
           <<setw(15)<<threshold<<endl;
     }
 
     tw = getTimeWalk(Q, threshold);
     return tw;
}

Referenced by getTimeWalk().

handle()

void CgemCalibFunSvc::handle

(

const Incident &

inc

)

Definition at line 125 of file CgemCalibFunSvc.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 Cgem Calib Constants" << endl
               << "  Please insert the following statement "
               << "in your \"jobOption.txt\" "
               << "before the include file of Cgem Reconstruction: "
               << endl << "        "
               << "#include \"$CALIBSVCROOT/share/job-CalibData.txt\""
               << endreq;
          }
     } 
}

initialize()

StatusCode CgemCalibFunSvc::initialize ( )

virtual

Definition at line 57 of file CgemCalibFunSvc.cxx.

                                      {
     MsgStream log(messageService(), name());
     log << MSG::INFO << "CgemCalibFunSvc::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);
     // }
 
     m_getRunNo = false;
     static IJobOptionsSvc* jobSvc = 0;
     if ( jobSvc == 0 ) {
      sc = service("JobOptionsSvc", jobSvc);
      if ( sc.isFailure() ) {
           std::cout << "Can't get the JobOptionsSvc @ DistBoss::GetPropertyValue()" << std::endl;
           return sc;
      }
     }
     const vector<const Property*>* properties = jobSvc->getProperties("ReadCosmicRayData");
     if ( properties == NULL ) {
      log << MSG::WARNING << "can not get runNo from ReadCosmicRayData" << endreq;
     } else{
      for ( unsigned int i = 0; i < properties->size(); ++i ) {
           if ( properties->at(i)->name() == "runNo" ) {
            cout << "name of property: " << properties->at(i)->name() << endl;
            string strRnd = properties->at(i)->toString();
            sscanf(strRnd.c_str(), "%d", &m_run);
            // cout << "runNo: " << m_run << endl;
            break;
           }
      }
      m_getRunNo = true;
     }
     if( !initCalibConst() ){
          log << MSG::ERROR << "can not initialize CGEM calibration constants" << endreq;
      return StatusCode::FAILURE;
     }
 
     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;
     }
 
     m_readLUT = true;
     m_cgemLUTReader = new CgemLUTReader(m_fileLUT.c_str());
     cout << "LUT file: " << m_fileLUT << endl;
     if(!m_cgemLUTReader->ReadLUT()){
      m_readLUT = false;
      log << MSG::FATAL << "can not read LUT " << m_fileLUT << endl;
     }
 
     return StatusCode::SUCCESS;
}

queryInterface()

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

virtual

Definition at line 48 of file CgemCalibFunSvc.cxx.

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

---

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

• CgemCalibFunSvc.h
• CgemCalibFunSvc.cxx