ManualGenerator Class Reference

#include <ManualGenerator.h>

Inheritance diagram for ManualGenerator:

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

Detailed Description

Definition at line 61 of file ManualGenerator.h.

Constructor & Destructor Documentation

ManualGenerator()

ManualGenerator::ManualGenerator	(	const string &	name,
		ISvcLocator *	pSvcLocator )

Definition at line 65 of file ManualGenerator.cxx.

                                                                            :Algorithm( name, pSvcLocator ){
 
 
    // get the name of the input file from the job options file
 
  declareProperty("InputFile", m_inputFileName = "");
 
}

Member Function Documentation

execute()

StatusCode ManualGenerator::execute

(

)

Definition at line 116 of file ManualGenerator.cxx.

                                   { 
 
 
    // read an event from the input file
 
  int nParticles;
  int idParticles[100];
  double pxParticle;
  double pyParticle;
  double pzParticle;
  double eParticle;
  CLHEP::HepLorentzVector pParticles[100];
 
  if (m_inputFile >> nParticles){
    for (int i = 0; i < nParticles; i++){
      m_inputFile >> idParticles[i];
      m_inputFile >> pxParticle;  pParticles[i].setPx(pxParticle);
      m_inputFile >> pyParticle;  pParticles[i].setPy(pyParticle);
      m_inputFile >> pzParticle;  pParticles[i].setPz(pzParticle);
      m_inputFile >> eParticle;   pParticles[i].setE(eParticle);
    }
  }
 
  else{
    cout << "ManualGenerator: RAN OUT OF INPUT EVENTS FROM FILE " 
         << m_inputFileName << endl;
    return StatusCode::SUCCESS;
  }
  
 
 
    // print out a message for this event to make sure all is well
 
  cout << "ManualGenerator: PROCESSING EVENT " << ++m_eventCounter << endl;
  cout << nParticles << endl;
  for (int i = 0; i < nParticles; i++){
    cout << idParticles[i] << "\t";
    cout << pParticles[i].px() << "\t";
    cout << pParticles[i].py() << "\t";
    cout << pParticles[i].pz() << "\t";
    cout << pParticles[i].e() << endl;
  }
 
 
    // retrieve the original MC event (HepMC::GenEvent) from the framework
 
  HepMC::GenEvent* genEvent = 0;
  SmartDataPtr<McGenEventCol> initialMcCol(eventSvc(), "/Event/Gen");
  if (initialMcCol == 0){
//    HepMC::GenEvent* genEvent = new GenEvent();  
//    genEvent->set_event_number(0);
 
    //create a Transient store 
//    McGenEventCol *mcColl = new McGenEventCol;
//    McGenEvent* mcEvent = new McGenEvent(genEvent);
//    mcColl->push_back(mcEvent);
//    StatusCode sc = eventSvc()->registerObject("/Event/Gen", mcColl);
//    if (sc != StatusCode::SUCCESS) return StatusCode::FAILURE;    
 
    cout << "ManualGenerator: COULD NOT FIND THE ORIGINAL MC EVENT COLLECTION" << endl;
    exit(0);
  }
//  else
//  {
    genEvent = (*(initialMcCol->begin()))->getGenEvt();
//  }
  //genEvent->print();
 
 
 
    // loop over MC particles and look for the particle we want to decay
 
  HepMC::GenParticle* decayParticle = NULL;
  for(HepMC::GenEvent::particle_const_iterator igenParticle = genEvent->particles_begin();
      igenParticle != genEvent->particles_end(); ++igenParticle ){
    decayParticle = *igenParticle;
    if (decayParticle->status() == 1 && abs(decayParticle->pdg_id()) != 22) break;
  }
  if (!decayParticle){
    cout << "ManualGenerator: COULD NOT FIND A PARTICLE TO DECAY" << endl;
    exit(0);
  }
 
 
    // create a new vertex from the decaying particle's production vertex
    //  and add it to the event
 
  HepMC::GenVertex* initialVertex = decayParticle->production_vertex();
  HepMC::GenVertex* newVertex = new HepMC::GenVertex(initialVertex->position());
  genEvent->add_vertex(newVertex);
  newVertex->add_particle_in(decayParticle);
 
 
 
    // boost the input four-vectors into the decaying particle's rest frame
 
  CLHEP::HepLorentzVector initialFourMomentum(decayParticle->momentum().px(),
                                              decayParticle->momentum().py(),
                                              decayParticle->momentum().pz(),
                                              decayParticle->momentum().e());
  CLHEP::HepLorentzVector finalFourMomentum(0,0,0,0);
  for (int i = 0; i < nParticles; i++){
  //  pParticles[i].boost(initialFourMomentum.boostVector());
    finalFourMomentum += pParticles[i];
  }
  if ((fabs(initialFourMomentum.px() - finalFourMomentum.px()) > 0.0001) || 
      (fabs(initialFourMomentum.py() - finalFourMomentum.py()) > 0.0001) ||
      (fabs(initialFourMomentum.pz() - finalFourMomentum.pz()) > 0.0001) ||
      (fabs(initialFourMomentum.e()  - finalFourMomentum.e())  > 0.01)){
    cout << "ManualGenerator: INITIAL AND FINAL STATES HAVE INCONSISTENT 4-MOMENTA" << endl;
    cout << "delta px = " << initialFourMomentum.px() - finalFourMomentum.px() << endl; 
    cout << "delta py = " << initialFourMomentum.py() - finalFourMomentum.py() << endl; 
    cout << "delta pz = " << initialFourMomentum.pz() - finalFourMomentum.pz() << endl; 
    cout << "delta e = " << initialFourMomentum.e() - finalFourMomentum.e() << endl; 
   // exit(0);
  }
 
 
    // add the input four-vectors to the new vertex
 
  for (int i = 0; i < nParticles; i++){
    newVertex->add_particle_out(new GenParticle(pParticles[i], idParticles[i], 1));
  }
  //genEvent->print();
 
 
 
  return StatusCode::SUCCESS; 
}

finalize()

StatusCode ManualGenerator::finalize

(

)

Definition at line 256 of file ManualGenerator.cxx.

{
  m_inputFile.close();
  return StatusCode::SUCCESS;
}

initialize()

StatusCode ManualGenerator::initialize

(

)

Definition at line 83 of file ManualGenerator.cxx.

                                      {
  //MsgStream log(messageService(), name());
  //log << MSG::INFO << "ManualGenerator initialize" << endreq;
 
 
    // open the input file
 
  m_inputFile.open(m_inputFileName.c_str());
  if (!m_inputFile){
    cout << "ManualGenerator: PROBLEMS OPENING FILE "
         << m_inputFileName << endl;
    exit(0);
  }
 
 
    // zero the event counter
 
  m_eventCounter = 0;
 
 
  return StatusCode::SUCCESS;
 
}

---

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

• ManualGenerator.h
• ManualGenerator.cxx