MdcNavigation Class Reference

#include <MdcNavigation.h>

Inheritance diagram for MdcNavigation:

Public Member Functions
	MdcNavigation (const std::string &name, ISvcLocator *pSvcLocator)
StatusCode	initialize ()
StatusCode	execute ()
StatusCode	finalize ()
StatusCode	beginRun ()

Detailed Description

Definition at line 23 of file MdcNavigation.h.

Constructor & Destructor Documentation

MdcNavigation()

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

Definition at line 44 of file MdcNavigation.cxx.

                                                                            : 
  Algorithm(name, pSvcLocator) {
    declareProperty("hist", m_hist = 0);
    declareProperty("nMcHit", m_nMcHit= 5);
    declareProperty("mc", m_mc = 1);
 
    declareProperty("maxMdcDigi",   m_maxMdcDigi= 0);
    declareProperty("keepBadTdc",   m_keepBadTdc= 0);
    declareProperty("dropHot",      m_dropHot= 0);
    declareProperty("keepUnmatch",  m_keepUnmatch= 0);
 
    declareProperty("poca", m_poca = false);
    declareProperty("doSag", m_doSag = false);
 
    declareProperty("d0Cut", m_d0Cut= 1.);
    declareProperty("z0Cut", m_z0Cut= 10.);
    declareProperty("debug", m_debug = 0);
 
  }

Member Function Documentation

beginRun()

StatusCode MdcNavigation::beginRun

(

)

Definition at line 114 of file MdcNavigation.cxx.

                                  {
  MsgStream log(msgSvc(), name());
  log << MSG::INFO << "in beginRun()" << endreq;
 
  m_gm = MdcDetector::instance(m_doSag);
  if(NULL == m_gm) return StatusCode::FAILURE;
 
  return StatusCode::SUCCESS;
}

execute()

StatusCode MdcNavigation::execute

(

)

Definition at line 126 of file MdcNavigation.cxx.

                                  {
  setFilterPassed(false);
  MsgStream log(msgSvc(), name());
  StatusCode sc = StatusCode::SUCCESS;
 
  // Get EventNavigator from the TDS
  SmartDataPtr<EventNavigator> navigator (eventSvc(),"/Event/Navigator");
  if( ! navigator ) {
    log << MSG::WARNING<< " Unable to retrieve EventNavigator" << endreq;
    m_rawData = true;
  }
  SmartDataPtr<RecMdcTrackCol> recMdcTrackCol(eventSvc(), "/Event/Recon/RecMdcTrackCol");
  SmartDataPtr<RecMdcHitCol> recMdcHitCol(eventSvc(), "/Event/Recon/RecMdcHitCol");
 
  // get eventNo, t0 and MdcDigi
  if(m_hist){
    sc = fillInit();
    if ( sc!=StatusCode::SUCCESS ) {
      return StatusCode::FAILURE;
    }
  }
  if (m_mc){
    // Get McParticleCol 
    SmartDataPtr<Event::McParticleCol> mcParticles(eventSvc(),"/Event/MC/McParticleCol");
    SmartDataPtr<Event::MdcMcHitCol> mcHit(eventSvc(),"/Event/MC/MdcMcHitCol");
    if( ! mcParticles ) {
      log << MSG::WARNING<< " Unable to retrieve McParticleCol" << endreq;
    }else{
      // For each McParticle ...       
      t_mcTkNum = 0;
      McParticleCol::iterator it= mcParticles->begin();
      log <<MSG::INFO << "mcParticles size = "<<mcParticles->size() << endreq;//yzhang debug
      for(; it!= mcParticles->end(); it++ ) {     
    //int tkId = (*it)->trackIndex();
    t_mcTkNum++;
      } 
    }
  }
  t_mcTkNum = 0;
  t_recTkNum = 0;
  //for each rec track
 
  if (!recMdcTrackCol){
    log << MSG::WARNING<< " Unable to retrieve recMdcTrackCol" << endreq;
    return StatusCode::SUCCESS;
  }
  t_recTkNum = recMdcTrackCol->size();
 
  //=============loop over tracks==============
  RecMdcTrackCol::iterator it = recMdcTrackCol->begin(); 
  for(;it != recMdcTrackCol->end(); it++ ) {   
    if ( m_mc && navigator ) {
      McParticleVector particles = navigator->getMcParticles(*it);
      t_mcTkNum = particles.size();
      RecMdcTrackVector tracks = navigator->getMdcTracks(particles[0]);
      // for FIRST parent particle
      if ( sc!=StatusCode::SUCCESS ) return StatusCode::FAILURE;
    }
    sc = fillRecTrack(*it, t_mcTkNum, t_recTkNum);
    t_nTk++;
    if ( sc!=StatusCode::SUCCESS ) return StatusCode::FAILURE;
  }//end of loop over tracks
 
  //=============loop over hits==============
  fillRecHits(*recMdcHitCol);
 
 
  if(m_hist){fillEvent();}
 
  return StatusCode::SUCCESS;
}

finalize()

StatusCode MdcNavigation::finalize

(

)

Definition at line 200 of file MdcNavigation.cxx.

                                   {
  MsgStream log(msgSvc(), name());
  log << MSG::INFO << "in finalize()" << endreq;
  std::cout<< "nTk ==  "<<t_nTk << std::endl;//yzhang debug
  delete keepedParticles; //FIXME
  return StatusCode::SUCCESS;
}

initialize()

StatusCode MdcNavigation::initialize

(

)

Definition at line 67 of file MdcNavigation.cxx.

                                    {
  MsgStream log(msgSvc(), name());
  StatusCode sc = StatusCode::SUCCESS;
 
  t_nTk = 0;
  //Initailize magnetic filed 
  sc = service ("MagneticFieldSvc",m_pIMF);
  if(sc!=StatusCode::SUCCESS) {
    log << MSG::ERROR << "Unable to open Magnetic field service"<<endreq;
    return StatusCode::FAILURE;
  }
 
 
  // Get the Particle Properties Service
  IPartPropSvc* p_PartPropSvc;
  static const bool CREATEIFNOTTHERE(true);
  sc = service("PartPropSvc", p_PartPropSvc, CREATEIFNOTTHERE);
  if (!sc.isSuccess() || 0 == p_PartPropSvc) {
    log << MSG::ERROR << " Could not initialize PartPropSvc" << endreq;
    return sc;
  }
 
  m_particleTable = p_PartPropSvc->PDT();
 
  IRawDataProviderSvc* irawDataProviderSvc;
  sc = service ("RawDataProviderSvc", irawDataProviderSvc);
  m_rawDataProviderSvc = dynamic_cast<RawDataProviderSvc*> (irawDataProviderSvc);
  if ( sc.isFailure() ){
    log << MSG::FATAL << "Could not load RawDataProviderSvc!" << endreq;
    return StatusCode::FAILURE;
  }
 
 
  if (m_hist){ 
    sc = bookNTuple();
    if (!sc.isSuccess()) { 
      log << MSG::WARNING << " Could not book NTuple" << endreq; 
      m_hist = 0;
    }
  }
 
  keepedParticles = new int(211);
 
 
  return StatusCode::SUCCESS;
}

---

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

• MdcNavigation.h
• MdcNavigation.cxx