EmcTimeRec Class Reference

#include <EmcTimeRec.h>

List of all members.

Public Member Functions
	EmcTimeRec (const std::string &name, ISvcLocator *pSvcLocator)
StatusCode	initialize ()
StatusCode	execute ()
StatusCode	finalize ()
Private Attributes
bool	m_output
NTuple::Tuple *	m_tuple
NTuple::Item< float >	m_nt_emcTime
NTuple::Item< float >	m_nt_estTime
NTuple::Item< float >	m_nt_corTime
NTuple::Item< long >	m_nt_npart
NTuple::Item< long >	m_nt_ntheta
NTuple::Item< long >	m_nt_nphi
NTuple::Item< float >	m_nt_e5x5
NTuple::Item< long >	m_nt_status

---

Detailed Description

Definition at line 12 of file EmcTimeRec.h.

---

Constructor & Destructor Documentation

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

Definition at line 16 of file EmcTimeRec.cxx.

References m_output.

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

---

Member Function Documentation

StatusCode EmcTimeRec::execute

(

)

Definition at line 54 of file EmcTimeRec.cxx.

References EmcID::barrel_ec(), Bes_Common::DEBUG, m_nt_corTime, m_nt_e5x5, m_nt_emcTime, m_nt_estTime, m_nt_npart, m_nt_nphi, m_nt_ntheta, m_nt_status, m_output, m_tuple, msgSvc(), EmcID::phi_module(), EventModel::Recon::RecEmcShowerCol, EmcID::theta_module(), and Bes_Common::WARNING.

                               {

  MsgStream log(msgSvc(), name());
  log << MSG::DEBUG << "in execute()" << endreq;   

  SmartDataPtr<RecEsTimeCol> estimeCol(eventSvc(),"/Event/Recon/RecEsTimeCol");

  if( !estimeCol || ( estimeCol->size() == 0 ) ) {
    log << MSG::WARNING << "Could not find RecEsTimeCol!" << endreq;
    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!" << endreq;
    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;
}

StatusCode EmcTimeRec::finalize

(

)

Definition at line 104 of file EmcTimeRec.cxx.

References Bes_Common::INFO, and msgSvc().

                                {

  MsgStream log(msgSvc(), name());
  log << MSG::INFO << "in finalize()" << endreq;

  return StatusCode::SUCCESS;
}

StatusCode EmcTimeRec::initialize

(

)

Definition at line 23 of file EmcTimeRec.cxx.

References calibUtil::ERROR, Bes_Common::INFO, m_nt_corTime, m_nt_e5x5, m_nt_emcTime, m_nt_estTime, m_nt_npart, m_nt_nphi, m_nt_ntheta, m_nt_status, m_output, m_tuple, msgSvc(), and ntupleSvc().

                                  {

  MsgStream log(msgSvc(), name());
  log << MSG::INFO << "in initialize()" << endreq;

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

---

Member Data Documentation

NTuple::Item<float> EmcTimeRec::m_nt_corTime [private]

Definition at line 25 of file EmcTimeRec.h.

Referenced by execute(), and initialize().

NTuple::Item<float> EmcTimeRec::m_nt_e5x5 [private]

Definition at line 29 of file EmcTimeRec.h.

Referenced by execute(), and initialize().

NTuple::Item<float> EmcTimeRec::m_nt_emcTime [private]

Definition at line 23 of file EmcTimeRec.h.

Referenced by execute(), and initialize().

NTuple::Item<float> EmcTimeRec::m_nt_estTime [private]

Definition at line 24 of file EmcTimeRec.h.

Referenced by execute(), and initialize().

NTuple::Item<long> EmcTimeRec::m_nt_npart [private]

Definition at line 26 of file EmcTimeRec.h.

Referenced by execute(), and initialize().

NTuple::Item<long> EmcTimeRec::m_nt_nphi [private]

Definition at line 28 of file EmcTimeRec.h.

Referenced by execute(), and initialize().

NTuple::Item<long> EmcTimeRec::m_nt_ntheta [private]

Definition at line 27 of file EmcTimeRec.h.

Referenced by execute(), and initialize().

NTuple::Item<long> EmcTimeRec::m_nt_status [private]

Definition at line 30 of file EmcTimeRec.h.

Referenced by execute(), and initialize().

bool EmcTimeRec::m_output [private]

Definition at line 21 of file EmcTimeRec.h.

Referenced by EmcTimeRec(), execute(), and initialize().

NTuple::Tuple* EmcTimeRec::m_tuple [private]

Definition at line 22 of file EmcTimeRec.h.

Referenced by execute(), and initialize().

---