EmcTimeRec Class Reference

#include <EmcTimeRec.h>

Inheritance diagram for EmcTimeRec:

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

Detailed Description

Definition at line 12 of file EmcTimeRec.h.

Constructor & Destructor Documentation

EmcTimeRec()

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

Definition at line 17 of file EmcTimeRec.cxx.

    : Algorithm( name, pSvcLocator ) {
  declareProperty( "Output", m_output = false );
}

Referenced by EmcTimeRec().

Member Function Documentation

execute()

StatusCode EmcTimeRec::execute

(

)

Definition at line 58 of file EmcTimeRec.cxx.

                               {
 
  MsgStream log( msgSvc(), name() );
  log << MSG::DEBUG << "in execute()" << endmsg;
 
  SmartDataPtr<RecEsTimeCol> estimeCol( eventSvc(), "/Event/Recon/RecEsTimeCol" );
 
  if ( !estimeCol || ( estimeCol->size() == 0 ) )
  {
    log << MSG::WARNING << "Could not find RecEsTimeCol!" << endmsg;
    return StatusCode::SUCCESS;
  }
 
  RecEsTimeCol::iterator iter_ESTime = estimeCol->begin();
  double                 t0          = ( *iter_ESTime )->getTest(); // get t0
 
  SmartDataPtr<RecEmcShowerCol> emcRecShowerCol( eventSvc(),
                                                 EventModel::Recon::RecEmcShowerCol );
  if ( !emcRecShowerCol )
  {
    log << MSG::WARNING << "Could not find RecEmcShowerCol!" << endmsg;
    return StatusCode::SUCCESS;
  }
 
  RecEmcShowerCol::iterator iShowerCol;
  for ( iShowerCol = emcRecShowerCol->begin(); iShowerCol != emcRecShowerCol->end();
        iShowerCol++ )
  {
 
    if ( ( *iShowerCol )->status() >= 100 ) continue; // estTime has been subtracted
 
    double tEmc = ( *iShowerCol )->time();
    ( *iShowerCol )->setTime( tEmc - (int)( t0 / 50 ) );
    ( *iShowerCol )->setStatus( ( *iShowerCol )->status() % 100 + 100 );
 
    if ( m_output )
    {
      m_nt_emcTime = tEmc;
      m_nt_estTime = (double)(int)( t0 / 50 );
      m_nt_corTime = ( *iShowerCol )->time();
      RecEmcID id( ( *iShowerCol )->getShowerId() );
      m_nt_npart  = EmcID::barrel_ec( id );
      m_nt_ntheta = EmcID::theta_module( id );
      m_nt_nphi   = EmcID::phi_module( id );
      m_nt_e5x5   = ( *iShowerCol )->e5x5();
      m_nt_status = ( *iShowerCol )->status();
      m_tuple->write();
    }
  }
 
  return StatusCode::SUCCESS;
}

finalize()

StatusCode EmcTimeRec::finalize

(

)

Definition at line 112 of file EmcTimeRec.cxx.

                                {
 
  MsgStream log( msgSvc(), name() );
  log << MSG::INFO << "in finalize()" << endmsg;
 
  return StatusCode::SUCCESS;
}

initialize()

StatusCode EmcTimeRec::initialize

(

)

Definition at line 23 of file EmcTimeRec.cxx.

                                  {
 
  MsgStream log( msgSvc(), name() );
  log << MSG::INFO << "in initialize()" << endmsg;
 
  if ( m_output )
  {
    NTuplePtr nt( ntupleSvc(), "FILE303/n1" );
    if ( nt ) m_tuple = nt;
    else
    {
      m_tuple = ntupleSvc()->book( "FILE303/n1", CLID_ColumnWiseTuple, "EmcTimeRec" );
      if ( m_tuple )
      {
        m_tuple->addItem( "emcTime", m_nt_emcTime );
        m_tuple->addItem( "estTime", m_nt_estTime );
        m_tuple->addItem( "corTime", m_nt_corTime );
        m_tuple->addItem( "npart", m_nt_npart );
        m_tuple->addItem( "ntheta", m_nt_ntheta );
        m_tuple->addItem( "nphi", m_nt_nphi );
        m_tuple->addItem( "e5x5", m_nt_e5x5 );
        m_tuple->addItem( "status", m_nt_status );
      }
      else
      { // did not manage to book the N tuple....
        log << MSG::ERROR << "Cannot book N-tuple:" << long( m_tuple ) << endmsg;
        return StatusCode::FAILURE;
      }
    }
  }
 
  return StatusCode::SUCCESS;
}

---

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

• EmcTimeRec.h
• EmcTimeRec.cxx