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 43 of file TruthDemo.cxx.

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

Member Function Documentation

execute()

StatusCode TruthDemo::execute

(

)

Definition at line 68 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(messageService(), name());
  msglog << MSG::INFO << ">>> TruthDemo from execute" << endreq;
    //
  // 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(messageService(), name());
        msglog << MSG::INFO << ">>> TruthDemo from finalize" << endreq;
        return StatusCode::SUCCESS;
}

initialize()

StatusCode TruthDemo::initialize

(

)

Definition at line 51 of file TruthDemo.cxx.

                                {
  StatusCode result = StatusCode::SUCCESS;
  MsgStream msglog(messageService(), name());
  msglog << MSG::INFO << ">>> Truthdemo from Initialize" << endreq;
  m_hgenerated = histoSvc()->book("/stat/1Dhist/1","Generated",100,0,1200);
  if (0 == m_hgenerated) {
    msglog << MSG::ERROR << " ERROR booking histogram" << endreq;
    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