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

Referenced by ValidRecMdcDedxAlg().

Member Function Documentation

execute()

StatusCode ValidRecMdcDedxAlg::execute

(

)

Definition at line 48 of file ValidRecMdcDedxAlg.cxx.

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

finalize()

StatusCode ValidRecMdcDedxAlg::finalize

(

)

Definition at line 83 of file ValidRecMdcDedxAlg.cxx.

                                        {
 
  MsgStream log( msgSvc(), name() );
  log << MSG::INFO << "Finalizing..." << endmsg;
  return StatusCode::SUCCESS;
}

initialize()

StatusCode ValidRecMdcDedxAlg::initialize

(

)

Definition at line 17 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