---

EvtGenProducer.h

Go to the documentation of this file.

// -*- C++ -*-
//
// Package:    EvtGenInterface
// Class:      EvtGenProducer
// 
//
// Original Author:  Nello Nappi
//         Created:  Fri May 11 15:33:30 CEST 2007
// $Id: EvtGenProducer.h,v 1.4 2008/03/12 16:28:35 covarell Exp $
//
//

// system include files
#include <memory>

// user include files
#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/Framework/interface/EDProducer.h"

#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/MakerMacros.h"

#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "CLHEP/Random/RandFlat.h"

#include "EvtGen/EvtGen.hh"     
#include "EvtGenBase/EvtId.hh"
#include "EvtGenBase/EvtPDL.hh"
#include "EvtGenBase/EvtDecayTable.hh"
#include "EvtGenBase/EvtSpinType.hh"
#include "EvtGenBase/EvtVector4R.hh"
#include "EvtGenBase/EvtParticle.hh"
#include "EvtGenBase/EvtScalarParticle.hh"
#include "EvtGenBase/EvtStringParticle.hh"
#include "EvtGenBase/EvtDiracParticle.hh"
#include "EvtGenBase/EvtVectorParticle.hh"
#include "EvtGenBase/EvtRaritaSchwingerParticle.hh"
#include "EvtGenBase/EvtTensorParticle.hh"
#include "EvtGenBase/EvtHighSpinParticle.hh"
#include "EvtGenBase/EvtStdHep.hh"
#include "EvtGenBase/EvtSecondary.hh"
#include "EvtGenModels/EvtPythia.hh"

namespace CLHEP {
  class HepRandomEngine;
  class RandFlat;
}

namespace HepMC {
  class GenParticle;
  class GenEvent;
}

class EvtGenProducer : public edm::EDProducer
{
public:

  explicit EvtGenProducer(edm::ParameterSet const & p);
  virtual ~EvtGenProducer();

private:                                                
  virtual void beginJob(const edm::EventSetup & c);                
  virtual void endJob();                                           
  virtual void produce(edm::Event & e, const edm::EventSetup & c); 
  void decay(HepMC::GenParticle* partHep, EvtId idEvt, HepMC::GenEvent* theEvent, bool del_daug);
  void go_through_daughters(EvtParticle* part);
  void update_candlist( int theIndex, HepMC::GenParticle *thePart );
  
  // from Pythia
  void call_pygive(const std::string& iParm );

private:
  CLHEP::RandFlat* m_flat;   
  EvtGen *m_EvtGen;
  std::vector<EvtId> forced_Evt;     // EvtId's of particels with forced decay
  std::vector<int> forced_Hep;       // HepId's of particels with forced decay
  int nforced;                       // number of particels with forced decay
  int ntotal, npartial, nevent;      // generic counters

  int nPythia;
  std::vector<std::string> pythia_params;  // Pythia stuff

  int nlist; 
  HepMC::GenParticle *listp[10]; 
  int index[10];                     // list of candidates to be forced 

};

---