ValidRecMdcDedxAlg Class Reference

#include <ValidRecMdcDedxAlg.h>

List of all members.

Public Member Functions
	ValidRecMdcDedxAlg (const std::string &name, ISvcLocator *pSvcLocator)
StatusCode	initialize ()
StatusCode	execute ()
StatusCode	finalize ()
Private Attributes
NTuple::Tuple *	m_tuple1
NTuple::Item< double >	m_dedx
NTuple::Item< double >	m_parttype
NTuple::Item< double >	m_chidedx
NTuple::Item< double >	m_chidedxE
NTuple::Item< double >	m_chidedxMu
NTuple::Item< double >	m_chidedxPi
NTuple::Item< double >	m_chidedxK
NTuple::Item< double >	m_chidedxP

---

Detailed Description

Definition at line 14 of file ValidRecMdcDedxAlg.h.

---

Constructor & Destructor Documentation

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

Definition at line 11 of file ValidRecMdcDedxAlg.cxx.

References m_tuple1.

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

---

Member Function Documentation

StatusCode ValidRecMdcDedxAlg::execute

(

)

Definition at line 49 of file ValidRecMdcDedxAlg.cxx.

References Bes_Common::DEBUG, electron, calibUtil::ERROR, genRecEmupikp::i, Bes_Common::INFO, kaon, m_chidedx, m_chidedxE, m_chidedxK, m_chidedxMu, m_chidedxP, m_chidedxPi, m_dedx, m_parttype, m_tuple1, msgSvc(), muon, pion, and proton.

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

StatusCode ValidRecMdcDedxAlg::finalize

(

)

Definition at line 86 of file ValidRecMdcDedxAlg.cxx.

References Bes_Common::INFO, and msgSvc().

                                        {

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

}

StatusCode ValidRecMdcDedxAlg::initialize

(

)

Definition at line 18 of file ValidRecMdcDedxAlg.cxx.

References calibUtil::ERROR, Bes_Common::INFO, m_chidedx, m_chidedxE, m_chidedxK, m_chidedxMu, m_chidedxP, m_chidedxPi, m_dedx, m_parttype, m_tuple1, msgSvc(), and ntupleSvc().

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


}

---

Member Data Documentation

NTuple::Item<double> ValidRecMdcDedxAlg::m_chidedx [private]

Definition at line 28 of file ValidRecMdcDedxAlg.h.

Referenced by execute(), and initialize().

NTuple::Item<double> ValidRecMdcDedxAlg::m_chidedxE [private]

Definition at line 29 of file ValidRecMdcDedxAlg.h.

Referenced by execute(), and initialize().

NTuple::Item<double> ValidRecMdcDedxAlg::m_chidedxK [private]

Definition at line 32 of file ValidRecMdcDedxAlg.h.

Referenced by execute(), and initialize().

NTuple::Item<double> ValidRecMdcDedxAlg::m_chidedxMu [private]

Definition at line 30 of file ValidRecMdcDedxAlg.h.

Referenced by execute(), and initialize().

NTuple::Item<double> ValidRecMdcDedxAlg::m_chidedxP [private]

Definition at line 33 of file ValidRecMdcDedxAlg.h.

Referenced by execute(), and initialize().

NTuple::Item<double> ValidRecMdcDedxAlg::m_chidedxPi [private]

Definition at line 31 of file ValidRecMdcDedxAlg.h.

Referenced by execute(), and initialize().

NTuple::Item<double> ValidRecMdcDedxAlg::m_dedx [private]

Definition at line 26 of file ValidRecMdcDedxAlg.h.

Referenced by execute(), and initialize().

NTuple::Item<double> ValidRecMdcDedxAlg::m_parttype [private]

Definition at line 27 of file ValidRecMdcDedxAlg.h.

Referenced by execute(), and initialize().

NTuple::Tuple* ValidRecMdcDedxAlg::m_tuple1 [private]

Definition at line 24 of file ValidRecMdcDedxAlg.h.

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

---