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 41 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 );
}

Referenced by MdcNavigation().

Member Function Documentation

beginRun()

StatusCode MdcNavigation::beginRun

(

)

Definition at line 111 of file MdcNavigation.cxx.

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

execute()

StatusCode MdcNavigation::execute

(

)

Definition at line 122 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" << endmsg;
    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" << endmsg; }
    else
    {
      // For each McParticle ...
      t_mcTkNum                  = 0;
      McParticleCol::iterator it = mcParticles->begin();
      log << MSG::INFO << "mcParticles size = " << mcParticles->size()
          << endmsg; // 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" << endmsg;
    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 202 of file MdcNavigation.cxx.

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

initialize()

StatusCode MdcNavigation::initialize

(

)

Definition at line 62 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" << endmsg;
    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" << endmsg;
    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!" << endmsg;
    return StatusCode::FAILURE;
  }
 
  if ( m_hist )
  {
    sc = bookNTuple();
    if ( !sc.isSuccess() )
    {
      log << MSG::WARNING << " Could not book NTuple" << endmsg;
      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