/home/bes3soft/bes3soft/Boss/7.0.2/dist/7.0.2/Validation/ValidRecMdcDedx/ValidRecMdcDedx-00-00-07-p01/src/ValidRecMdcDedxAlg.cxx

Go to the documentation of this file.

#include "ValidRecMdcDedx/ValidRecMdcDedxAlg.h"
#include "GaudiKernel/MsgStream.h"
#include "GaudiKernel/SmartDataPtr.h"
#include "GaudiKernel/IHistogramSvc.h"
#include "GaudiKernel/INTupleSvc.h"

#include "MdcRecEvent/RecMdcDedx.h"
#include "Identifier/MdcID.h"


ValidRecMdcDedxAlg::ValidRecMdcDedxAlg(const std::string& name, ISvcLocator* pSvcLocator) :
   Algorithm(name, pSvcLocator) {
    
     m_tuple1 = 0; 
       
   }

StatusCode ValidRecMdcDedxAlg::initialize(){
   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;


}


StatusCode ValidRecMdcDedxAlg::execute() {
   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() {

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

}

---