MrpcRec Class Reference

#include <MrpcRec.h>

Inheritance diagram for MrpcRec:

Public Member Functions
	MrpcRec (const std::string &name, ISvcLocator *pSvcLocator)
StatusCode	initialize ()
StatusCode	beginRun ()
StatusCode	execute ()
StatusCode	finalize ()
void	clearTofTrackVec (std::vector< MrpcTrack * > *&tofTrackVec)

Detailed Description

Definition at line 43 of file MrpcRec.h.

Constructor & Destructor Documentation

MrpcRec()

MrpcRec::MrpcRec	(	const std::string &	name,
		ISvcLocator *	pSvcLocator
	)

Definition at line 48 of file MrpcRec.cxx.

                                                                :
  Algorithm(name, pSvcLocator)
{   
  declareProperty( "AcceleratorStatus", m_acceleratorStatus );
  declareProperty( "MagneticField",     m_magneticField     );
  declareProperty( "FirstIteration",    m_firstIteration    );
  declareProperty( "PrintOutInfo",      m_printOutInfo      );
  declareProperty( "neighborhood",      m_neighborhood=6    );//For the time being: 6 is the standardvalue -> It yield the highest reconstruction efficiency! 
}

Member Function Documentation

beginRun()

StatusCode MrpcRec::beginRun

(

)

Definition at line 98 of file MrpcRec.cxx.

                            {
  MsgStream log(msgSvc(), name());
  log << MSG::INFO <<"MrpcRec begin_run!"<< endreq;
  return StatusCode::SUCCESS;
}

clearTofTrackVec()

void MrpcRec::clearTofTrackVec

(

std::vector< MrpcTrack * > *&

tofTrackVec

)

Definition at line 250 of file MrpcRec.cxx.

                                                                   {
  if( tofTrackVec ) {
    std::vector<MrpcTrack*>::iterator iter = tofTrackVec->begin();
    for( ; iter < tofTrackVec->end(); iter++ ) {
      delete (*iter);
    }
    tofTrackVec->clear();
    delete tofTrackVec;
  }
  return;
}

Referenced by execute().

execute()

StatusCode MrpcRec::execute

(

)

Definition at line 106 of file MrpcRec.cxx.

                            {
 
  MsgStream log(msgSvc(), name());
  log << MSG::INFO << "MrpcRec in execute()!" << endreq;
 
  SmartDataPtr<Event::EventHeader> eventHeader(eventSvc(),"/Event/EventHeader");
  if( !eventHeader ) {
    log << MSG::FATAL << "MrpcRec could not find Event Header!" << endreq;
    return StatusCode::FAILURE;
  }
  int run = eventHeader->runNumber();
  int event = eventHeader->eventNumber();
  if( ( event % 1000 == 0 ) && m_printOutInfo ) {
    std::cout << "run:" << run << "  event: " << event << std::endl;
  }
  log << MSG::INFO << "run= " << run << "  event= " << event << endreq;
 
  MrpcRecTDS tds;
  tds.RegisterNullRecTofTrackCol();
 
  // Magnetic Field: Colliding data and Cosmic Ray
  if( m_acceleratorStatus == "Colliding" ) {
    // Retrieve Event Start Time
    SmartDataPtr<RecEsTimeCol> estimeCol(eventSvc(),"/Event/Recon/RecEsTimeCol");
    if( !estimeCol || ( estimeCol->size() == 0 ) ) { 
      log << MSG::WARNING << "MrpcRec Could not find RecEsTimeCol!   Run = " << run << " Event = " << event << endreq;
      return StatusCode::SUCCESS;
    } 
    RecEsTimeCol::iterator iter_ESTime=estimeCol->begin();
    double t0  = (*iter_ESTime)->getTest();
    int t0Stat = (*iter_ESTime)->getStat();
    log << MSG::INFO << "t0= " << t0 << "  t0Stat= " << t0Stat << endreq;
 
 
 
    // Kalman Filter Track
    SmartDataPtr<RecMdcKalTrackCol> mdcKalTrackCol(eventSvc(),"/Event/Recon/RecMdcKalTrackCol");
    if( !mdcKalTrackCol ) {
      log << MSG::WARNING << "No MdcKalTrackCol in TDS!   Run = " << run << "  Event = " << event << endreq;
      return StatusCode::SUCCESS;
    }
 
 
 
    // Tof Get Extrapolated Track
    SmartDataPtr<RecExtTrackCol> extTrackCol(eventSvc(),"/Event/Recon/RecExtTrackCol");
    if( !extTrackCol ) {
      log << MSG::WARNING << "No ExtTrackCol in TDS!   Run = " << run << "  Event = " << event << endreq;
      return StatusCode::SUCCESS;
    }
    else {
      if( m_printOutInfo ) { m_printOut->setExtTrackNum( extTrackCol->size() ); }
    }
 
 
 
    TofTrackVec* tofTrackVec = new TofTrackVec;
    RecExtTrackCol::iterator iter_track = extTrackCol->begin();
    for( ; iter_track < extTrackCol->end(); iter_track++ ) {
      RecMdcKalTrackCol::iterator iter_kal = mdcKalTrackCol->begin();
      for( ; iter_kal != mdcKalTrackCol->end(); iter_kal++ ) {
    if( (*iter_kal)->trackId() == (*iter_track)->trackId() ) break;
      }
      int kal[5] = {-1};
      if( iter_kal != mdcKalTrackCol->end() ) {
    for( unsigned int i=0; i<5; i++ ) {
      kal[i] = (*iter_kal)->getStat( 1, i );
    }
      }
      MrpcTrack* tof = new MrpcTrack;
      tof->setExtTrack( (*iter_track), kal, t0, t0Stat );
 
      if( tofTrackVec->size()>0 ) {
    std::vector<MrpcTrack*>::iterator iterExt = tofTrackVec->begin();
    for( ; iterExt < tofTrackVec->end(); iterExt++ ) {
      if( (*iterExt)->isNoHit() ) continue;
      tof->getMultiHit( *iterExt );
    }
      }
 
      tofTrackVec->push_back( tof );
    }
 
    if( m_printOutInfo ) {
      m_printOut->setTrack1Col( tofTrackVec );
    }
 
    // Retrieve Tof Digi Data
    TofDataMap tofDataMap = tofDigiSvc->tofDataMapTof( t0 );
    if( tofDataMap.empty() ) {
      log << MSG::WARNING << "No Tof Data Map in RawDataProviderSvc!   Run=" << run << " Event=" << event << endreq;
    }
 
 
    std::vector<MrpcTrack*>::iterator iter = tofTrackVec->begin();
    for( ; iter < tofTrackVec->end(); iter++ ) {
      if( (*iter)->isNoHit() ) continue;
      (*iter)->setTofData( tofDataMap,m_neighborhood );
      if( m_printOutInfo ) { m_printOut->setTrack2( (*iter) ); }
      if( (*iter)->isNoHit() ) continue;
      (*iter)->match(tofTrackVec );
      if( m_printOutInfo ) { m_printOut->setTrack3( (*iter) ); }
    }
 
    iter = tofTrackVec->begin();
    for( ; iter < tofTrackVec->end(); iter++ ) 
      {
    (*iter)->setCalibration();
      }
 
    tds.RegisterTDS( eventHeader->runNumber(), eventHeader->eventNumber(), tofTrackVec );
 
    clearTofTrackVec( tofTrackVec );
 
// Check RecTofTrackCol Registered
    SmartDataPtr<RecTofTrackCol> tofTrackCol(eventSvc(),"/Event/Recon/RecTofTrackCol");  
    if (!tofTrackCol) { 
      log << MSG::FATAL << "MrpcRec could not find RecTofTrackCol!" << endreq;
      return StatusCode::FAILURE;
    }
 
  }
  else {
    log << MSG::FATAL << "In MrpcRec: AcceleratorStatus is NOT correct! m_acceleratorStatus = " << m_acceleratorStatus << endreq;
    return StatusCode::FAILURE;
  }
 
  return StatusCode::SUCCESS;
 
}

finalize()

StatusCode MrpcRec::finalize

(

)

Definition at line 239 of file MrpcRec.cxx.

                             {
  if( m_printOutInfo ) {
    m_printOut->final();
    delete m_printOut;
  }
  return StatusCode::SUCCESS;
}

initialize()

StatusCode MrpcRec::initialize

(

)

Definition at line 60 of file MrpcRec.cxx.

                              {
 
  MsgStream log(msgSvc(), name());
  log << MSG::INFO << "MrpcRec in initialize()" << endreq;   
 
  //Get TOF Geomertry Service
  StatusCode scg = service("TofGeomSvc", tofGeomSvc);
  if (scg==  StatusCode::SUCCESS) {
    log << MSG::INFO << "MrpcRec Get Geometry Service Sucessfully!" << endreq;
  }else {
    log << MSG::ERROR << "MrpcRec Get Geometry Service Failed !" << endreq;
    return StatusCode::FAILURE;
  }
 
  //Get TOF Calibtration Service
  StatusCode scc = service("TofCaliSvc", tofCaliSvc);
  if (scc ==  StatusCode::SUCCESS) { 
    log << MSG::INFO << "MrpcRec Get Calibration Service Sucessfully!" << endreq;
  } else {
    log << MSG::ERROR << "MrpcRec Get Calibration Service Failed !" << endreq;
    return StatusCode::FAILURE;
  }
 
  //Get TOF Raw Data Provider Service
  StatusCode scd = service("RawDataProviderSvc", tofDigiSvc);
  if (scd ==  StatusCode::SUCCESS) { 
    log << MSG::INFO << "MrpcRec Get Tof Raw Data Service Sucessfully!" << endreq;
  } else {
    log << MSG::ERROR << "MrpcRec Get Tof Raw Data Service Failed !" << endreq;
    return StatusCode::FAILURE;
  }
 
 
  if( m_printOutInfo ) { m_printOut = new MrpcCount; }
 
  return StatusCode::SUCCESS;
}

---

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

• MrpcRec.h
• MrpcRec.cxx