TofEnergyRec Class Reference

#include <TofEnergyRec.h>

Inheritance diagram for TofEnergyRec:

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

Detailed Description

Definition at line 15 of file TofEnergyRec.h.

Constructor & Destructor Documentation

TofEnergyRec()

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

Definition at line 42 of file TofEnergyRec.cxx.

    : Algorithm( name, pSvcLocator ) {
  declareProperty( "Output", m_output = false );
  declareProperty( "EventNb", m_eventNb = 0 );
  declareProperty( "CalibConst", m_calibConst = 0.0085 );
  declareProperty( "IsData", m_isData = true );
}

Referenced by TofEnergyRec().

Member Function Documentation

execute()

StatusCode TofEnergyRec::execute

(

)

Definition at line 95 of file TofEnergyRec.cxx.

                                 {
 
  MsgStream log( msgSvc(), name() );
  log << MSG::INFO << "TofEnergyRec in execute()!!" << endmsg;
 
  SmartDataPtr<Event::EventHeader> eventHeader( eventSvc(), "/Event/EventHeader" );
  if ( !eventHeader )
  {
    log << MSG::FATAL << "TofEnergyRec could not find Event Header!" << endmsg;
    return StatusCode::FAILURE;
  }
  int run   = eventHeader->runNumber();
  int event = eventHeader->eventNumber();
  if ( m_eventNb != 0 && m_event % m_eventNb == 0 )
    std::cout << "--------event: " << m_event << ", run= " << run << "  event= " << event
              << std::endl;
  log << MSG::INFO << "run= " << run << "  event= " << event << endmsg;
  m_event++;
 
  IRawDataProviderSvc* tofDigiSvc;
  StatusCode           sc = service( "RawDataProviderSvc", tofDigiSvc );
  if ( sc != StatusCode::SUCCESS )
  {
    log << MSG::FATAL << "TofEnergyRec could not find RawDataProviderSvc!" << endmsg;
    return StatusCode::SUCCESS;
  }
 
  // check whether the Recon has been already registered
  DataObject* aReconEvent;
  eventSvc()->findObject( "/Event/Recon", aReconEvent );
  if ( aReconEvent == NULL )
  {
    log << MSG::INFO << "Register ReconEvent" << endmsg;
    // then register Recon
    aReconEvent = new ReconEvent();
    sc          = eventSvc()->registerObject( "/Event/Recon", aReconEvent );
    if ( sc != StatusCode::SUCCESS )
    {
      log << MSG::FATAL << "Could not register ReconEvent" << endmsg;
      return StatusCode::FAILURE;
    }
  }
 
  // Register RecTofTrackCol
  SmartDataPtr<RecTofTrackCol> tofTrackCol( eventSvc(), EventModel::Recon::RecTofTrackCol );
  if ( !tofTrackCol )
  {
    log << MSG::INFO << "Register RecTofTrackCol!" << endmsg;
    tofTrackCol = new RecTofTrackCol;
    sc          = eventSvc()->registerObject( EventModel::Recon::RecTofTrackCol, tofTrackCol );
    if ( sc != StatusCode::SUCCESS )
    { log << MSG::FATAL << "Could not register RecTofTrackCol!" << endmsg; }
  }
 
  // Retrieve RecEsTimeCol
  SmartDataPtr<RecEsTimeCol> estimeCol( eventSvc(), "/Event/Recon/RecEsTimeCol" );
  if ( !estimeCol || ( estimeCol->size() == 0 ) )
  {
    log << MSG::WARNING << "TofRec Could not find RecEsTimeCol!   Run = " << run
        << " Event = " << event << endmsg;
    return StatusCode::SUCCESS;
  }
  RecEsTimeCol::iterator iter_ESTime = estimeCol->begin();
 
  double t0  = 0;
  t0         = ( *iter_ESTime )->getTest();
  int t0Stat = ( *iter_ESTime )->getStat();
 
  std::vector<TofData*> tofDataVec;
  tofDataVec = tofDigiSvc->tofDataVectorTof();
 
  m_tofShower->findShower( tofDataVec, tofTrackCol, t0 );
 
  return StatusCode::SUCCESS;
}

finalize()

StatusCode TofEnergyRec::finalize

(

)

Definition at line 173 of file TofEnergyRec.cxx.

                                  {
  cout << "Total event:" << m_event << endl;
  delete m_tofShower;
  return StatusCode::SUCCESS;
}

initialize()

StatusCode TofEnergyRec::initialize

(

)

Definition at line 52 of file TofEnergyRec.cxx.

                                    {
 
  MsgStream log( msgSvc(), name() );
  log << MSG::INFO << "TofEnergyRec in initialize()" << endmsg;
 
  m_event     = 0;
  m_tofShower = new TofShower;
  m_tofShower->setCalibConst( m_calibConst );
  m_tofShower->setIsData( m_isData );
  // m_tofShower->readCalibPar();
 
  if ( m_output )
  {
    cout << "--------Book Tof Energy Ntuple!" << endl;
    NTuplePtr nt( ntupleSvc(), "FILE209/sci" );
    NTuplePtr nt1( ntupleSvc(), "FILE209/shower" );
    NTuplePtr nt2( ntupleSvc(), "FILE209/seed" );
    if ( nt && nt1 && nt2 )
    {
      m_tuple  = nt;
      m_tuple1 = nt1;
      m_tuple2 = nt2;
    }
    else
    {
      m_tuple  = ntupleSvc()->book( "FILE209/sci", CLID_ColumnWiseTuple, "TofEnergyRec" );
      m_tuple1 = ntupleSvc()->book( "FILE209/shower", CLID_ColumnWiseTuple, "TofEnergyRec" );
      m_tuple2 = ntupleSvc()->book( "FILE209/seed", CLID_ColumnWiseTuple, "TofEnergyRec" );
      if ( m_tuple && m_tuple1 && m_tuple2 )
      { m_tofShower->BookNtuple( m_tuple, m_tuple1, m_tuple2 ); }
      else
      {
        log << MSG::ERROR << " ----- TofShower cannot book N-tuple:" << long( m_tuple )
            << endmsg;
      }
    }
  }
 
  return StatusCode::SUCCESS;
}

---

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

• TofEnergyRec.h
• TofEnergyRec.cxx