#include <TruthDemo.h>

Inheritance diagram for TruthDemo:

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

StatusCode	initialize ()

StatusCode	execute ()

StatusCode	finalize ()

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

StatusCode	initialize ()

StatusCode	execute ()

StatusCode	finalize ()

Detailed Description

Definition at line 9 of file Generator/GenAnalysisTools/TruthExamples/TruthExamples-00-00-07/TruthExamples/TruthDemo.h.

Constructor & Destructor Documentation

◆ TruthDemo() [1/2]

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

◆ TruthDemo() [2/2]

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

Member Function Documentation

◆ execute() [1/2]

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

◆ execute() [2/2]

StatusCode TruthDemo::execute ( )

◆ finalize() [1/2]

StatusCode TruthDemo::finalize ( )

Definition at line 96 of file TruthDemo.cxx.

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

◆ finalize() [2/2]

StatusCode TruthDemo::finalize ( )

◆ initialize() [1/2]

StatusCode TruthDemo::initialize ( )

Definition at line 50 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;
}

◆ initialize() [2/2]