TruthDemo Class Reference

#include <TruthDemo.h>

Inheritance diagram for TruthDemo:

Public Member Functions
	TruthDemo (const std::string &name, ISvcLocator *pSvcLocator)
StatusCode	initialize ()
StatusCode	execute ()
StatusCode	finalize ()

Detailed Description

Definition at line 9 of file TruthDemo.h.

Constructor & Destructor Documentation

TruthDemo()

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

Definition at line 42 of file TruthDemo.cxx.

    : Algorithm( name, pSvcLocator ) {
  // Declare the algorithm's properties
  declareProperty( "HistogramFlag", m_produceHistogram = true );
  //  declareProperty("HistogramFlag", m_produceHistogram = false );
}

Referenced by TruthDemo().

Member Function Documentation

execute()

StatusCode TruthDemo::execute

(

)

Definition at line 67 of file TruthDemo.cxx.

                              {
  //   HepMC::IO_PDG_ParticleDataTable pdg_io("PDGTABLE");
  //   static HepMC::ParticleDataTable *pp = NULL;
  //   if (pp == NULL) pp = pdg_io.read_particle_data_table();
  MsgStream msglog( msgSvc(), name() );
  msglog << MSG::INFO << ">>> TruthDemo from execute" << endmsg;
  //
  // Px/Py Balance and total energy
  //
  float totenergy = 0.;
  float pxbalance = 0.;
  float pybalance = 0.;
  // Iterate over MC particles  We are using the IsGenStable predicate from
  IsGenStable                            ifs;
  std::vector<const HepMC::GenParticle*> particles;
  StatusCode                             stat = m_tesIO->getMC( particles, &ifs );
  for ( std::vector<const HepMC::GenParticle*>::iterator pitr = particles.begin();
        pitr != particles.end(); pitr++ )
  {
    pxbalance += ( *pitr )->momentum().x();
    pybalance += ( *pitr )->momentum().y();
    totenergy += ( *pitr )->momentum().e();
  }
  m_pxBalance->fill( pxbalance, 1. );
  m_pyBalance->fill( pybalance, 1. );
  m_totEnergy->fill( totenergy, 1. );
  // End of execution for each event
  return StatusCode::SUCCESS;
}

finalize()

StatusCode TruthDemo::finalize

(

)

Definition at line 97 of file TruthDemo.cxx.

                               {
  MsgStream msglog( msgSvc(), name() );
  msglog << MSG::INFO << ">>> TruthDemo from finalize" << endmsg;
  return StatusCode::SUCCESS;
}

initialize()

StatusCode TruthDemo::initialize

(

)

Definition at line 49 of file TruthDemo.cxx.

                                 {
  StatusCode result = StatusCode::SUCCESS;
  MsgStream  msglog( msgSvc(), name() );
  msglog << MSG::INFO << ">>> Truthdemo from Initialize" << endmsg;
  m_hgenerated = histoSvc()->book( "/stat/1Dhist/1", "Generated", 100, 0, 1200 );
  if ( 0 == m_hgenerated )
  {
    msglog << MSG::ERROR << " ERROR booking histogram" << endmsg;
    result = StatusCode::FAILURE;
  }
  m_pxBalance = histoSvc()->book( "/stat/1Dhist/25", "px balance", 50, -10., 10. );
  m_pyBalance = histoSvc()->book( "/stat/1Dhist/26", "py balance", 50, -10., 10. );
  m_totEnergy = histoSvc()->book( "/stat/1Dhist/27", "total energy", 50, 10000., 20000. );
  m_tesIO     = new GenAccessIO();
 
  //        return StatusCode::SUCCESS;
  return result;
}

---

The documentation for this class was generated from the following files:

• TruthDemo.h
• TruthDemo.cxx