ValidRecMdcDedxAlg Class Reference

#include <ValidRecMdcDedxAlg.h>

Inheritance diagram for ValidRecMdcDedxAlg:

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

Detailed Description

Definition at line 14 of file ValidRecMdcDedxAlg.h.

Constructor & Destructor Documentation

ValidRecMdcDedxAlg()

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

Definition at line 11 of file ValidRecMdcDedxAlg.cxx.

                                                                                      :
   Algorithm(name, pSvcLocator) {
    
     m_tuple1 = 0; 
       
   }

Member Function Documentation

execute()

StatusCode ValidRecMdcDedxAlg::execute

(

)

Definition at line 49 of file ValidRecMdcDedxAlg.cxx.

                                       {
   MsgStream log(msgSvc(), name());
   SmartDataPtr<RecMdcDedxCol> mdcDedx(eventSvc(),"/Event/Recon/RecMdcDedxCol");
   if( ! mdcDedx )
      {
         log << MSG::ERROR << "Unable to retrieve RecMdcDedxCol" << endreq;
         return  StatusCode::FAILURE;
      } else {
         log << MSG::DEBUG << "RecMdcDedxCol retrieved of size "<< mdcDedx->size() << endreq;
         int i = 0;         
     for(RecMdcDedxCol::iterator it=mdcDedx->begin(); it!=mdcDedx->end(); it++,i++)
       {  m_dedx = (*it)->probPH();
          pid_dedx m_pid_dedx = (*it)->particleType();       
          m_parttype = (int)m_pid_dedx;
           //cout<<"m_parttype is = "<<m_parttype<<endl;
          if((m_pid_dedx != electron) && (m_pid_dedx != muon) && (m_pid_dedx != pion)&& (m_pid_dedx != kaon)&& (m_pid_dedx != proton) ) continue;
           m_chidedx = (*it)->chi(m_parttype);
               m_chidedxE = (*it)->chiE();
           m_chidedxMu = (*it)->chiMu();
           m_chidedxPi = (*it)->chiPi();
           m_chidedxK = (*it)->chiK();
           m_chidedxP = (*it)->chiP();
         
               log << MSG::INFO << "particle chi="<< m_chidedx <<endreq;
           log << MSG::INFO << "m_dedx="<<m_dedx<<endreq;
           log << MSG::INFO << " m_parttype="<< m_parttype<<endreq;
           StatusCode status1 = m_tuple1->write();
           if ( status1.isFailure() )  {
         log << MSG::ERROR << "Cannot fill Ntuple1" << endreq;
           }
               
        }
      }
     return StatusCode::SUCCESS;
}

finalize()

StatusCode ValidRecMdcDedxAlg::finalize

(

)

Definition at line 86 of file ValidRecMdcDedxAlg.cxx.

                                        {
 
   MsgStream log(msgSvc(), name());
   log << MSG::INFO << "Finalizing..." << endreq;
 
}

initialize()

StatusCode ValidRecMdcDedxAlg::initialize

(

)

Definition at line 18 of file ValidRecMdcDedxAlg.cxx.

                                         {
   StatusCode status;
   MsgStream log( msgSvc(), name() );
 
   //  Book N-tuple 1
   NTuplePtr nt1(ntupleSvc(), "DedxTuples/Mdc_dedx");
   if ( nt1 ) m_tuple1 = nt1;
   else {
     m_tuple1 = ntupleSvc()->book ("DedxTuples/Mdc_dedx", CLID_RowWiseTuple, "Dedx parameter");
     if ( m_tuple1 )    {
       status = m_tuple1->addItem("dedx",m_dedx);
       status = m_tuple1->addItem("parttype",m_parttype);
       status = m_tuple1->addItem("chidedx",m_chidedx);
       status = m_tuple1->addItem("chidedxE",m_chidedxE);
       status = m_tuple1->addItem("chidedxMu",m_chidedxMu);
       status = m_tuple1->addItem("chidedxPi",m_chidedxPi);
       status = m_tuple1->addItem("chidedxK",m_chidedxK);
       status = m_tuple1->addItem("chidedxP",m_chidedxP);
     }
     else    {   // did not manage to book the N tuple....
       log << MSG::ERROR << "    Cannot book N-tuple:" << long(m_tuple1) << endmsg;
       return StatusCode::FAILURE;
     }
   }
   log << MSG::INFO << "Finished booking NTuples" << endmsg;
  return StatusCode::SUCCESS;
 
 
}

---

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

• ValidRecMdcDedxAlg.h
• ValidRecMdcDedxAlg.cxx