ValidRecMdcDedxAlg.cxx

Go to the documentation of this file.

#include "ValidRecMdcDedxAlg.h"
#include "GaudiKernel/IHistogramSvc.h"
#include "GaudiKernel/INTupleSvc.h"
#include "GaudiKernel/MsgStream.h"
#include "GaudiKernel/SmartDataPtr.h"
 
#include "Identifier/MdcID.h"
#include "MdcRecEvent/RecMdcDedx.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" << endmsg;
    return StatusCode::FAILURE;
  }
  else
  {
    log << MSG::DEBUG << "RecMdcDedxCol retrieved of size " << mdcDedx->size() << endmsg;
    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 << endmsg;
      log << MSG::INFO << "m_dedx=" << m_dedx << endmsg;
      log << MSG::INFO << " m_parttype=" << m_parttype << endmsg;
      StatusCode status1 = m_tuple1->write();
      if ( status1.isFailure() ) { log << MSG::ERROR << "Cannot fill Ntuple1" << endmsg; }
    }
  }
  return StatusCode::SUCCESS;
}
 
StatusCode ValidRecMdcDedxAlg::finalize() {
 
  MsgStream log( msgSvc(), name() );
  log << MSG::INFO << "Finalizing..." << endmsg;
  return StatusCode::SUCCESS;
}