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"
 
DECLARE_COMPONENT(ValidRecMdcDedxAlg)
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();
           m_parttype = (*it)->particleId();       
           //cout<<"m_parttype is = "<<m_parttype<<endl;
           if(m_parttype >5 || m_parttype<0) continue;
           m_chidedx = (*it)->chi(m_parttype-1);
               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;
 
}