myPythia6ServiceWithCallback Class Reference

Inheritance diagram for myPythia6ServiceWithCallback:

Public Member Functions
	myPythia6ServiceWithCallback (const edm::ParameterSet &cfg)
void	setEvent (const HepMC::GenEvent *genEvt)
	~myPythia6ServiceWithCallback ()
Public Member Functions inherited from gen::Pythia6Service
void	closePYUPDA ()
void	closeSLHA ()
virtual void	enter ()
void	openPYUPDA (const char *, bool write_file)
void	openSLHA (const char *)
	Pythia6Service ()
	Pythia6Service (edm::ParameterSet const &)
CLHEP::HepRandomEngine *	randomEngine () const
void	setCSAParams ()
void	setGeneralParams ()
void	setPYUPDAParams (bool afterPyinit)
void	setRandomEngine (CLHEP::HepRandomEngine *v)
void	setSLHAFromHeader (const std::vector< std::string > &lines)
void	setSLHAParams ()
	~Pythia6Service ()
Public Member Functions inherited from gen::FortranInstance
void	call (void(&fn)())
template<typename T >
T	call (T(&fn)())
template<typename A >
void	call (void(&fn)(A), A a)
template<typename T , typename A >
T	call (T(&fn)(A), A a)
template<typename A1 , typename A2 >
void	call (void(&fn)(A1, A2), A1 a1, A2 a2)
template<typename T , typename A1 , typename A2 >
T	call (T(&fn)(A1, A2), A1 a1, A2 a2)
	FortranInstance ()
virtual void	leave ()
virtual	~FortranInstance ()

Private Member Functions
void	upEvnt ()
void	upInit ()
bool	upVeto ()

Private Attributes
double	beamEnergy_
const HepMC::GenEvent *	genEvt_
int	verbosity_

Additional Inherited Members
Static Public Member Functions inherited from gen::FortranInstance
template<typename T >
static T *	getInstance ()
Static Public Attributes inherited from gen::FortranInstance
static const std::string	kFortranInstance = "FortranInstance"

Detailed Description

Definition at line 32 of file GenMuonRadiationAlgorithm.cc.

Constructor & Destructor Documentation

myPythia6ServiceWithCallback::myPythia6ServiceWithCallback ( const edm::ParameterSet &  cfg )

inline

Definition at line 35 of file GenMuonRadiationAlgorithm.cc.

    : gen::Pythia6Service(cfg),
      beamEnergy_(cfg.getParameter<double>("beamEnergy")),      
      genEvt_(0),
      verbosity_(cfg.getParameter<int>("verbosity"))      
  {}

myPythia6ServiceWithCallback::~myPythia6ServiceWithCallback ( )

inline

Definition at line 41 of file GenMuonRadiationAlgorithm.cc.

{}

Member Function Documentation

void myPythia6ServiceWithCallback::setEvent ( const HepMC::GenEvent *  genEvt )

inline

Definition at line 43 of file GenMuonRadiationAlgorithm.cc.

References TtGenEvtProducer_cfi::genEvt, and genEvt_.

Referenced by GenMuonRadiationAlgorithm::compFSR().

  {  
    genEvt_ = genEvt;
  }

void myPythia6ServiceWithCallback::upEvnt ( )

inlineprivatevirtual

Reimplemented from gen::FortranInstance.

Definition at line 71 of file GenMuonRadiationAlgorithm.cc.

References lhef::HEPEUP::AQCDUP, lhef::HEPEUP::AQEDUP, assert(), gather_cfg::cout, lhef::CommonBlocks::fillHEPEUP(), genEvt_, lhef::HEPEUP::ICOLUP, lhef::HEPEUP::IDPRUP, lhef::HEPEUP::IDUP, reco::if(), lhef::HEPEUP::ISTUP, lhef::HEPEUP::MOTHUP, lhef::HEPEUP::NUP, SingleMuPt40Fwdv2_cfi_GEN_SIM::pdgId, lhef::HEPEUP::PUP, lhef::HEPEUP::resize(), lhef::HEPEUP::SCALUP, lhef::HEPEUP::SPINUP, verbosity_, lhef::HEPEUP::VTIMUP, lhef::HEPEUP::XPDWUP, and lhef::HEPEUP::XWGTUP.

  {     
    if ( verbosity_ ) std::cout << "<upEvnt>:" << std::endl;
    int numParticles = 0;
    for ( HepMC::GenEvent::particle_const_iterator genParticle = genEvt_->particles_begin();
      genParticle != genEvt_->particles_end(); ++genParticle ) {
      ++numParticles;
    }
    lhef::HEPEUP hepeup;
    hepeup.NUP = numParticles;
    if ( verbosity_ ) std::cout << " numParticles = " << numParticles << std::endl;
    hepeup.IDPRUP = 1;
    hepeup.XWGTUP = 1.;
    hepeup.XPDWUP.first = 0.;
    hepeup.XPDWUP.second = 0.;
    hepeup.SCALUP = -1.;
    hepeup.AQEDUP = -1.;
    hepeup.AQCDUP = -1.;
    hepeup.resize();
    std::map<int, int> barcodeToIdxMap;
    int genParticleIdx = 0;
    for ( HepMC::GenEvent::particle_const_iterator genParticle = genEvt_->particles_begin();
      genParticle != genEvt_->particles_end(); ++genParticle ) {
      barcodeToIdxMap[(*genParticle)->barcode()] = genParticleIdx + 1; // CV: conversion between C++ and Fortran array index conventions
      int pdgId = (*genParticle)->pdg_id();
      hepeup.IDUP[genParticleIdx] = pdgId;
      if ( pdgId == 2212 && !(*genParticle)->production_vertex() ) {
    hepeup.ISTUP[genParticleIdx] = -1;
      } else {
    hepeup.ISTUP[genParticleIdx] = (*genParticle)->status();
      }
      if ( (*genParticle)->production_vertex() ) {
    int firstMother = 0;
    bool firstMother_initialized = false;
    int lastMother = 0;
    bool lastMother_initialized = false;
    for ( HepMC::GenVertex::particles_in_const_iterator genMother = (*genParticle)->production_vertex()->particles_in_const_begin();
          genMother != (*genParticle)->production_vertex()->particles_in_const_end(); ++genMother ) {
      int genMotherBarcode = (*genMother)->barcode();
      assert(barcodeToIdxMap.find(genMotherBarcode) != barcodeToIdxMap.end());
      int genMotherIdx = barcodeToIdxMap[genMotherBarcode];
      if ( genMotherIdx < firstMother || !firstMother_initialized ) {
        firstMother = genMotherIdx;
        firstMother_initialized = true;
      }
      if ( genMotherIdx > lastMother || !lastMother_initialized ) {
        lastMother = genMotherIdx;
        lastMother_initialized = true;
      }
    }
    hepeup.MOTHUP[genParticleIdx].first = firstMother;
    hepeup.MOTHUP[genParticleIdx].second = ( lastMother > firstMother ) ? lastMother : 0;
      } else {
    hepeup.MOTHUP[genParticleIdx].first = 0;
    hepeup.MOTHUP[genParticleIdx].second = 0;
      }
      hepeup.ICOLUP[genParticleIdx].first = 0;
      hepeup.ICOLUP[genParticleIdx].second = 0;
      hepeup.PUP[genParticleIdx][0] = (*genParticle)->momentum().px();
      hepeup.PUP[genParticleIdx][1] = (*genParticle)->momentum().py();
      hepeup.PUP[genParticleIdx][2] = (*genParticle)->momentum().pz();
      hepeup.PUP[genParticleIdx][3] = (*genParticle)->momentum().e();
      hepeup.PUP[genParticleIdx][4] = (*genParticle)->momentum().m();
      if ( (*genParticle)->status() == 1 ) hepeup.VTIMUP[genParticleIdx] = 1.e+6;
      else hepeup.VTIMUP[genParticleIdx] = 0.;
      hepeup.SPINUP[genParticleIdx] = 9;
      if ( verbosity_ ) {
    std::cout << "adding genParticle #" << (genParticleIdx + 1)
          << " (pdgId = " << hepeup.IDUP[genParticleIdx] << ", status = " << hepeup.ISTUP[genParticleIdx] << "):" 
          << " En = " << hepeup.PUP[genParticleIdx][3] << "," 
          << " Px = " << hepeup.PUP[genParticleIdx][0] << "," 
          << " Py = " << hepeup.PUP[genParticleIdx][1] << "," 
          << " Pz = " << hepeup.PUP[genParticleIdx][2] << " (mass = " << hepeup.PUP[genParticleIdx][4] << ")" << std::endl;
    std::cout << "(mothers = " << hepeup.MOTHUP[genParticleIdx].first << ":" << hepeup.MOTHUP[genParticleIdx].second << ")" << std::endl;
      }
      ++genParticleIdx;
    }
    lhef::CommonBlocks::fillHEPEUP(&hepeup);
  }

void myPythia6ServiceWithCallback::upInit ( )

inlineprivatevirtual

Reimplemented from gen::FortranInstance.

Definition at line 49 of file GenMuonRadiationAlgorithm.cc.

References beamEnergy_, gather_cfg::cout, alignCSCRings::e, lhef::HEPRUP::EBMUP, lhef::CommonBlocks::fillHEPRUP(), lhef::HEPRUP::IDBMUP, lhef::HEPRUP::IDWTUP, lhef::HEPRUP::LPRUP, lhef::HEPRUP::NPRUP, lhef::HEPRUP::PDFGUP, lhef::HEPRUP::PDFSUP, verbosity_, lhef::HEPRUP::XERRUP, lhef::HEPRUP::XMAXUP, and lhef::HEPRUP::XSECUP.

  { 
    // initialize LEs Houches accord common block (see hep-ph/0109068)
    // needed to run PYTHIA in parton shower mode
    if ( verbosity_ ) std::cout << "<upInit:" << std::endl;
    lhef::HEPRUP heprup;
    heprup.IDBMUP.first  = 2212;
    heprup.IDBMUP.second = 2212;
    heprup.EBMUP.first   = beamEnergy_;
    heprup.EBMUP.second  = beamEnergy_;
    heprup.PDFGUP.first  = -9; // CV: LHE input event contains proton remnants
    heprup.PDFGUP.second = -9;
    heprup.PDFSUP.first  = 10042;
    heprup.PDFSUP.second = 10042;
    heprup.IDWTUP        = 3;
    heprup.NPRUP         = 1;
    heprup.XSECUP.push_back(1.e+3);
    heprup.XERRUP.push_back(0.);
    heprup.XMAXUP.push_back(1.);
    heprup.LPRUP.push_back(1);
    lhef::CommonBlocks::fillHEPRUP(&heprup);   
  }

bool myPythia6ServiceWithCallback::upVeto ( )

inlineprivatevirtual

Reimplemented from gen::FortranInstance.

Definition at line 150 of file GenMuonRadiationAlgorithm.cc.

References gather_cfg::cout, and verbosity_.

  { 
     if ( verbosity_ ) std::cout << "<upVeto>:" << std::endl;
    return false; // keep current event
  }

Member Data Documentation

double myPythia6ServiceWithCallback::beamEnergy_

private

Definition at line 156 of file GenMuonRadiationAlgorithm.cc.

Referenced by upInit().

const HepMC::GenEvent* myPythia6ServiceWithCallback::genEvt_

private

Definition at line 158 of file GenMuonRadiationAlgorithm.cc.

Referenced by setEvent(), and upEvnt().

int myPythia6ServiceWithCallback::verbosity_

private

Definition at line 160 of file GenMuonRadiationAlgorithm.cc.

Referenced by upEvnt(), upInit(), and upVeto().