gen::Py8MassGun Class Reference

Inheritance diagram for gen::Py8MassGun:

Public Member Functions
const char *	classname () const override
bool	generatePartonsAndHadronize () override
	Py8MassGun (edm::ParameterSet const &)
	~Py8MassGun () override
Public Member Functions inherited from gen::Py8GunBase
void	evtGenDecay ()
void	finalizeEvent () override
bool	initializeForInternalPartons () override
	Py8GunBase (edm::ParameterSet const &ps)
virtual bool	residualDecay ()
void	setRandomEngine (CLHEP::HepRandomEngine *v)
std::vector< std::string > const &	sharedResources () const
void	statistics () override
	~Py8GunBase () override
Public Member Functions inherited from gen::Py8InterfaceBase
bool	decay ()
bool	declareSpecialSettings (const std::vector< std::string > &)
bool	declareStableParticles (const std::vector< int > &)
void	makeTmpSLHA (const std::string &)
void	p8SetRandomEngine (CLHEP::HepRandomEngine *v)
	Py8InterfaceBase (edm::ParameterSet const &ps)
P8RndmEngine &	randomEngine ()
bool	readSettings (int)
	~Py8InterfaceBase () override
Public Member Functions inherited from gen::BaseHadronizer
	BaseHadronizer (edm::ParameterSet const &ps)
void	cleanLHE ()
void	generateLHE (edm::LuminosityBlock const &lumi, CLHEP::HepRandomEngine *rengine, unsigned int ncpu)
edm::Event &	getEDMEvent () const
std::unique_ptr< HepMC::GenEvent >	getGenEvent ()
std::unique_ptr< GenEventInfoProduct >	getGenEventInfo ()
virtual std::unique_ptr< GenLumiInfoHeader >	getGenLumiInfoHeader () const
GenRunInfoProduct &	getGenRunInfo ()
std::unique_ptr< lhef::LHEEvent >	getLHEEvent ()
const std::shared_ptr< lhef::LHERunInfo > &	getLHERunInfo () const
const std::string &	gridpackPath () const
int	randomIndex () const
const std::string &	randomInitConfigDescription () const
void	randomizeIndex (edm::LuminosityBlock const &lumi, CLHEP::HepRandomEngine *rengine)
void	resetEvent (std::unique_ptr< HepMC::GenEvent > event)
void	resetEventInfo (std::unique_ptr< GenEventInfoProduct > eventInfo)
virtual bool	select (HepMC::GenEvent *) const
void	setEDMEvent (edm::Event &event)
void	setLHEEvent (std::unique_ptr< lhef::LHEEvent > event)
void	setLHERunInfo (std::unique_ptr< lhef::LHERunInfo > runInfo)
void	setRandomEngine (CLHEP::HepRandomEngine *v)
std::vector< std::string > const &	sharedResources () const
virtual	~BaseHadronizer () noexcept(false)

Private Attributes
double	fMaxEta
double	fMaxM
double	fMaxP
double	fMaxPt
double	fMinEta
double	fMinM
double	fMinP
double	fMinPt
int	fMomMode

Additional Inherited Members
Protected Member Functions inherited from gen::BaseHadronizer
std::unique_ptr< HepMC::GenEvent > &	event ()
std::unique_ptr< GenEventInfoProduct > &	eventInfo ()
lhef::LHEEvent *	lheEvent ()
lhef::LHERunInfo *	lheRunInfo ()
GenRunInfoProduct &	runInfo ()
Protected Attributes inherited from gen::Py8GunBase
double	fMaxPhi
double	fMinPhi
std::vector< int >	fPartIDs
Protected Attributes inherited from gen::Py8InterfaceBase
HepMC::IO_AsciiParticles *	ascii_io
std::shared_ptr< Pythia8::EvtGenDecays >	evtgenDecays
std::string	evtgenDecFile
std::string	evtgenPdlFile
std::vector< std::string >	evtgenUserFiles
std::unique_ptr< Pythia8::Pythia >	fDecayer
std::unique_ptr< Pythia8::Pythia >	fMasterGen
edm::ParameterSet	fParameters
unsigned int	maxEventsToPrint
bool	pythiaHepMCVerbosity
bool	pythiaHepMCVerbosityParticles
unsigned int	pythiaPylistVerbosity
std::string	slhafile_
HepMC::Pythia8ToHepMC	toHepMC
bool	useEvtGen
Protected Attributes inherited from gen::BaseHadronizer
std::string	lheFile_
int	randomIndex_

Detailed Description

Definition at line 11 of file Py8MassGun.cc.

Constructor & Destructor Documentation

Py8MassGun()

gen::Py8MassGun::Py8MassGun

(

edm::ParameterSet const &

ps

)

Definition at line 34 of file Py8MassGun.cc.

References fMaxEta, fMaxM, fMaxP, fMaxPt, fMinEta, fMinM, fMinP, fMinPt, fMomMode, and edm::ParameterSet::getParameter().

                                                  : Py8GunBase(ps) {
    // ParameterSet defpset ;
    edm::ParameterSet pgun_params = ps.getParameter<edm::ParameterSet>("PGunParameters");  // , defpset ) ;
    fMinEta = pgun_params.getParameter<double>("MinEta");                                  // ,-2.2);
    fMaxEta = pgun_params.getParameter<double>("MaxEta");                                  // , 2.2);
    fMinP = pgun_params.getParameter<double>("MinP");                                      // ,  0.);
    fMaxP = pgun_params.getParameter<double>("MaxP");                                      // ,  0.);
    fMinPt = pgun_params.getParameter<double>("MinPt");                                    // ,  0.);
    fMaxPt = pgun_params.getParameter<double>("MaxPt");                                    // ,  0.);
    fMinM = pgun_params.getParameter<double>("MinM");                                      // ,  0.);
    fMaxM = pgun_params.getParameter<double>("MaxM");                                      // ,  0.);
    fMomMode = pgun_params.getParameter<int>("MomMode");                                   // ,  1);
  }

~Py8MassGun()

gen::Py8MassGun::~Py8MassGun ( )

inlineoverride

Definition at line 14 of file Py8MassGun.cc.

{}

Member Function Documentation

classname()

const char * gen::Py8MassGun::classname ( ) const

overridevirtual

Implements gen::Py8InterfaceBase.

Definition at line 120 of file Py8MassGun.cc.

{ return "Py8MassGun"; }

generatePartonsAndHadronize()

bool gen::Py8MassGun::generatePartonsAndHadronize ( )

overridevirtual

Implements gen::Py8InterfaceBase.

Definition at line 48 of file Py8MassGun.cc.

References mps_setup::append, funct::cos(), StorageManager_cfg::e1, PVValHelper::eta, gen::BaseHadronizer::event(), JetChargeProducer_cfi::exp, gen::P8RndmEngine::flat(), gen::Py8InterfaceBase::fMasterGen, fMaxEta, fMaxM, fMaxP, gen::Py8GunBase::fMaxPhi, fMaxPt, fMinEta, fMinM, fMinP, gen::Py8GunBase::fMinPhi, fMinPt, fMomMode, gen::Py8GunBase::fPartIDs, callgraph::m2, M_PI, createTree::pp, PVValHelper::pT, gen::Py8InterfaceBase::randomEngine(), funct::sin(), mathSSE::sqrt(), and gen::Py8InterfaceBase::toHepMC.

                                               {
    fMasterGen->event.reset();
    size_t pSize = fPartIDs.size();
    if (pSize > 2)
      return false;

    int NTotParticles = fPartIDs.size();

    // energy below is dummy, it is not used
    (fMasterGen->event).append(990, -11, 0, 0, 2, 1 + NTotParticles, 0, 0, 0., 0., 0., 15000., 15000.);

    // Pick a flat mass range
    double phi, eta, the, ee, pp;
    double m0 = (fMaxM - fMinM) * randomEngine().flat() + fMinM;
    // Global eta
    eta = (fMaxEta - fMinEta) * randomEngine().flat() + fMinEta;

    if (pSize == 2) {
      // Masses.
      double m1 = fMasterGen->particleData.m0(fPartIDs[0]);
      double m2 = fMasterGen->particleData.m0(fPartIDs[1]);

      // Energies and absolute momentum in the rest frame.
      if (m1 + m2 > m0)
        return false;
      double e1 = 0.5 * (m0 * m0 + m1 * m1 - m2 * m2) / m0;
      double e2 = 0.5 * (m0 * m0 + m2 * m2 - m1 * m1) / m0;
      double pAbs = 0.5 * sqrt((m0 - m1 - m2) * (m0 + m1 + m2) * (m0 + m1 - m2) * (m0 - m1 + m2)) / m0;
      // Isotropic angles in rest frame give three-momentum.
      double cosTheta = 2. * randomEngine().flat() - 1.;
      double sinTheta = sqrt(1. - cosTheta * cosTheta);
      phi = 2. * M_PI * randomEngine().flat();

      double pX = pAbs * sinTheta * cos(phi);
      double pY = pAbs * sinTheta * sin(phi);
      double pZ = pAbs * cosTheta;

      (fMasterGen->event).append(fPartIDs[0], 1, 1, 0, 0, 0, 0, 0, pX, pY, pZ, e1, m1);
      (fMasterGen->event).append(fPartIDs[1], 1, 1, 0, 0, 0, 0, 0, -pX, -pY, -pZ, e2, m2);
    } else {
      (fMasterGen->event).append(fPartIDs[0], 1, 1, 0, 0, 0, 0, 0, 0.0, 0.0, 0.0, m0, m0);
    }

    //now the boost (from input params)
    if (fMomMode == 0) {
      pp = (fMaxP - fMinP) * randomEngine().flat() + fMinP;
    } else {
      double pT = (fMaxPt - fMinPt) * randomEngine().flat() + fMinPt;
      pp = pT * cosh(eta);
    }
    ee = sqrt(m0 * m0 + pp * pp);

    //the boost direction (from input params)
    //
    phi = (fMaxPhi - fMinPhi) * randomEngine().flat() + fMinPhi;
    the = 2. * atan(exp(-eta));

    double betaX = pp / ee * std::sin(the) * std::cos(phi);
    double betaY = pp / ee * std::sin(the) * std::sin(phi);
    double betaZ = pp / ee * std::cos(the);

    // boost all particles
    //
    (fMasterGen->event).bst(betaX, betaY, betaZ);

    if (!fMasterGen->next())
      return false;

    event() = std::make_unique<HepMC::GenEvent>();
    return toHepMC.fill_next_event(fMasterGen->event, event().get());
  }

Member Data Documentation

fMaxEta

double gen::Py8MassGun::fMaxEta

private

Definition at line 22 of file Py8MassGun.cc.

Referenced by generatePartonsAndHadronize(), and Py8MassGun().

fMaxM

double gen::Py8MassGun::fMaxM

private

Definition at line 28 of file Py8MassGun.cc.

Referenced by generatePartonsAndHadronize(), and Py8MassGun().

fMaxP

double gen::Py8MassGun::fMaxP

private

Definition at line 24 of file Py8MassGun.cc.

Referenced by generatePartonsAndHadronize(), and Py8MassGun().

fMaxPt

double gen::Py8MassGun::fMaxPt

private

Definition at line 26 of file Py8MassGun.cc.

Referenced by generatePartonsAndHadronize(), and Py8MassGun().

fMinEta

double gen::Py8MassGun::fMinEta

private

Definition at line 21 of file Py8MassGun.cc.

Referenced by generatePartonsAndHadronize(), and Py8MassGun().

fMinM

double gen::Py8MassGun::fMinM

private

Definition at line 27 of file Py8MassGun.cc.

Referenced by generatePartonsAndHadronize(), and Py8MassGun().

fMinP

double gen::Py8MassGun::fMinP

private

Definition at line 23 of file Py8MassGun.cc.

Referenced by generatePartonsAndHadronize(), and Py8MassGun().

fMinPt

double gen::Py8MassGun::fMinPt

private

Definition at line 25 of file Py8MassGun.cc.

Referenced by generatePartonsAndHadronize(), and Py8MassGun().

fMomMode

int gen::Py8MassGun::fMomMode

private

Definition at line 29 of file Py8MassGun.cc.

Referenced by generatePartonsAndHadronize(), and Py8MassGun().