#include <ReggeGribovPartonMC.h>

Inheritance diagram for gen::ReggeGribovPartonMCHadronizer:

Public Member Functions
const char *	classname () const

bool	decay ()

bool	declareSpecialSettings (const std::vector< std::string > &)

bool	declareStableParticles (const std::vector< int > &)

void	finalizeEvent ()

bool	generatePartonsAndHadronize ()

bool	hadronize ()

bool	initializeForExternalPartons ()

bool	initializeForInternalPartons ()

bool	initializeTablePaths ()

bool	readSettings (int)

	ReggeGribovPartonMCHadronizer (const edm::ParameterSet &)

	ReggeGribovPartonMCHadronizer (const edm::ParameterSet &)

bool	residualDecay ()

void	statistics ()

	~ReggeGribovPartonMCHadronizer () override

	~ReggeGribovPartonMCHadronizer () 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 Member Functions
void	doSetRandomEngine (CLHEP::HepRandomEngine *v) override

Private Attributes
int	m_BeamID

double	m_BeamMomentum

double	m_bMax

double	m_bMin

int	m_HEModel

double	m_ImpactParameter

bool	m_IsInitialized

int	m_NEvent

int	m_NParticles

edm::FileInPath	m_ParamFileName

double	m_PartEnergy [99990]

int	m_PartID [99990]

double	m_PartMass [99990]

double	m_PartPx [99990]

double	m_PartPy [99990]

double	m_PartPz [99990]

int	m_PartStatus [99990]

bool	m_SkipNuclFrag

int	m_TargetID

double	m_TargetMomentum

edm::ParameterSet	pset_

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::BaseHadronizer
std::string	lheFile_

int	randomIndex_

Detailed Description

Definition at line 20 of file ReggeGribovPartonMC.h.

Constructor & Destructor Documentation

◆ ReggeGribovPartonMCHadronizer() [1/2]

ReggeGribovPartonMCHadronizer::ReggeGribovPartonMCHadronizer ( const edm::ParameterSet & pset )

Definition at line 58 of file ReggeGribovPartonMCHadronizer.cc.

     : BaseHadronizer(pset),
       pset_(pset),
       m_BeamMomentum(pset.getParameter<double>("beammomentum")),
       m_TargetMomentum(pset.getParameter<double>("targetmomentum")),
       m_BeamID(pset.getParameter<int>("beamid")),
       m_TargetID(pset.getParameter<int>("targetid")),
       m_HEModel(pset.getParameter<int>("model")),
       m_bMin(pset.getParameter<double>("bmin")),
       m_bMax(pset.getParameter<double>("bmax")),
       m_ParamFileName(pset.getUntrackedParameter<string>("paramFileName")),
       m_SkipNuclFrag(pset.getParameter<bool>("skipNuclFrag")),
       m_NEvent(0),
       m_NParticles(0),
       m_ImpactParameter(0.),
       m_IsInitialized(false) {
   int nevet = 1;      //needed for CS
   int noTables = 0;   //don't calculate tables
   int LHEoutput = 0;  //no lhe
   int dummySeed = 123;
   char dummyName[] = "dummy";
   crmc_set_f_(nevet,
               dummySeed,
               m_BeamMomentum,
               m_TargetMomentum,
               m_BeamID,
               m_TargetID,
               m_HEModel,
               noTables,
               LHEoutput,
               dummyName,
               m_ParamFileName.fullPath().c_str());
  
   //additionally initialise tables
   initializeTablePaths();
  
   //change impact paramter
   nucl2_.bminim = float(m_bMin);
   nucl2_.bmaxim = float(m_bMax);
 }

References crmc_set_f_(), dqmMemoryStats::float, edm::FileInPath::fullPath(), initializeTablePaths(), m_BeamID, m_BeamMomentum, m_bMax, m_bMin, m_HEModel, m_ParamFileName, m_TargetID, m_TargetMomentum, and nucl2_.

◆ ~ReggeGribovPartonMCHadronizer() [1/2]

ReggeGribovPartonMCHadronizer::~ReggeGribovPartonMCHadronizer ( )

override

Definition at line 100 of file ReggeGribovPartonMCHadronizer.cc.

                                                               {
   // destructor
 }

◆ ReggeGribovPartonMCHadronizer() [2/2]

gen::ReggeGribovPartonMCHadronizer::ReggeGribovPartonMCHadronizer ( const edm::ParameterSet & )

◆ ~ReggeGribovPartonMCHadronizer() [2/2]

gen::ReggeGribovPartonMCHadronizer::~ReggeGribovPartonMCHadronizer ( )

override

Member Function Documentation

◆ classname()

const char * ReggeGribovPartonMCHadronizer::classname ( ) const

Definition at line 210 of file ReggeGribovPartonMCHadronizer.cc.

210 { return "gen::ReggeGribovPartonMCHadronizer"; }

◆ decay()

bool ReggeGribovPartonMCHadronizer::decay ( )

Definition at line 202 of file ReggeGribovPartonMCHadronizer.cc.

202 { return true; }

◆ declareSpecialSettings()

bool gen::ReggeGribovPartonMCHadronizer::declareSpecialSettings ( const std::vector< std::string > & )

inline

Definition at line 193 of file ReggeGribovPartonMCHadronizer.h.

193 { return true; }

◆ declareStableParticles()

bool ReggeGribovPartonMCHadronizer::declareStableParticles ( const std::vector< int > & )

Definition at line 212 of file ReggeGribovPartonMCHadronizer.cc.

212 { return true; }

◆ doSetRandomEngine()

void ReggeGribovPartonMCHadronizer::doSetRandomEngine ( CLHEP::HepRandomEngine * v )

overrideprivatevirtual

Reimplemented from gen::BaseHadronizer.

Definition at line 105 of file ReggeGribovPartonMCHadronizer.cc.

105 { reggeGribovRandomEngine = v; }

References reggeGribovRandomEngine, and gen::v.

◆ finalizeEvent()

void ReggeGribovPartonMCHadronizer::finalizeEvent ( )

Definition at line 206 of file ReggeGribovPartonMCHadronizer.cc.

206 { return; }

◆ generatePartonsAndHadronize()

bool ReggeGribovPartonMCHadronizer::generatePartonsAndHadronize ( )

Definition at line 108 of file ReggeGribovPartonMCHadronizer.cc.

                                                                 {
   int iout = 0, ievent = 0;
   crmc_f_(iout,
           ievent,
           m_NParticles,
           m_ImpactParameter,
           m_PartID[0],
           m_PartPx[0],
           m_PartPy[0],
           m_PartPz[0],
           m_PartEnergy[0],
           m_PartMass[0],
           m_PartStatus[0]);
   LogDebug("ReggeGribovPartonMCInterface") << "event generated" << endl;
  
   const bool isHeavyIon = (m_TargetID + m_BeamID > 2);
  
   if (isHeavyIon)
     conv.set_trust_beam_particles(false);
  
   conv.set_skip_nuclear_fragments(m_SkipNuclFrag);
  
   HepMC::GenEvent* evt = conv.read_next_event();
  
   evt->set_event_number(m_NEvent++);
   int sig_id = -1;
   switch (int(c2evt_.typevt))  // if negative typevt mini plasma was created by event (except -4)
   {
     case 0:
       break;  //unknown for qgsjetII
     case 1:
       sig_id = 101;
       break;
     case -1:
       sig_id = 101;
       break;
     case 2:
       sig_id = 105;
       break;
     case -2:
       sig_id = 105;
       break;
     case 3:
       sig_id = 102;
       break;
     case -3:
       sig_id = 102;
       break;
     case 4:
       sig_id = 103;
       break;
     case -4:
       sig_id = 104;
       break;
     default:
       LogDebug("ReggeGribovPartonMCInterface") << "Signal ID not recognised for setting HEPEVT" << endl;
   }
   evt->set_signal_process_id(sig_id);  //an integer ID uniquely specifying the signal process (i.e. MSUB in Pythia)
  
 #ifdef HEPMC_HAS_CROSS_SECTION
   // set cross section information for this event
   HepMC::GenCrossSection theCrossSection;
   theCrossSection.set_cross_section(double(hadr5_.sigineaa) * 1e9);  //required in pB
   evt->set_cross_section(theCrossSection);
 #endif
  
   if (isHeavyIon)  //other than pp
   {
     HepMC::HeavyIon ion(cevt_.kohevt,                 // Number of hard scatterings
                         cevt_.npjevt,                 // Number of projectile participants
                         cevt_.ntgevt,                 // Number of target participants
                         cevt_.kolevt,                 // Number of NN (nucleon-nucleon) collisions
                         cevt_.npnevt + cevt_.ntnevt,  // Number of spectator neutrons
                         cevt_.nppevt + cevt_.ntpevt,  // Number of spectator protons
                         -1,                           // Number of N-Nwounded collisions
                         -1,                           // Number of Nwounded-N collisons
                         -1,                           // Number of Nwounded-Nwounded collisions
                         cevt_.bimevt,                 // Impact Parameter(fm) of collision
                         cevt_.phievt,                 // Azimuthal angle of event plane
                         c2evt_.fglevt,                // eccentricity of participating nucleons
                         hadr5_.sigine * 1e9);         // nucleon-nucleon inelastic (in pB)
     evt->set_heavy_ion(ion);
   }
  
   event().reset(evt);
   //evt->print();
   //EPOS::EPOS_Wrapper::print_hepcom();
  
   return true;
 }

References c2evt_, cevt_, conv, crmc_f_(), gen::BaseHadronizer::event(), hadr5_, hcal_dqm_sourceclient-live_cfg::isHeavyIon, LogDebug, m_BeamID, m_ImpactParameter, m_NEvent, m_NParticles, m_PartEnergy, m_PartID, m_PartMass, m_PartPx, m_PartPy, m_PartPz, m_PartStatus, m_SkipNuclFrag, m_TargetID, EPOS::IO_EPOS::set_skip_nuclear_fragments(), and EPOS::IO_EPOS::set_trust_beam_particles().

◆ hadronize()

bool ReggeGribovPartonMCHadronizer::hadronize ( )

Definition at line 200 of file ReggeGribovPartonMCHadronizer.cc.

200 { return false; }

◆ initializeForExternalPartons()

bool gen::ReggeGribovPartonMCHadronizer::initializeForExternalPartons ( )

inline

Definition at line 190 of file ReggeGribovPartonMCHadronizer.h.

190 { return true; }

◆ initializeForInternalPartons()

bool ReggeGribovPartonMCHadronizer::initializeForInternalPartons ( )

Definition at line 214 of file ReggeGribovPartonMCHadronizer.cc.

                                                                  {
   if (!m_IsInitialized) {
     //use set parameters to init models
     crmc_init_f_();
     m_IsInitialized = true;
   }
   return true;
 }

References crmc_init_f_(), and m_IsInitialized.

◆ initializeTablePaths()

bool ReggeGribovPartonMCHadronizer::initializeTablePaths ( )

Definition at line 223 of file ReggeGribovPartonMCHadronizer.cc.

                                                          {
   //epos
   string path_fnii(FileInPath("GeneratorInterface/ReggeGribovPartonMCInterface/data/epos.initl").fullPath());
   string path_fnie(FileInPath("GeneratorInterface/ReggeGribovPartonMCInterface/data/epos.iniev").fullPath());
   string path_fnrj(FileInPath("GeneratorInterface/ReggeGribovPartonMCInterface/data/epos.inirj").fullPath());
   string path_fncs(FileInPath("GeneratorInterface/ReggeGribovPartonMCInterface/data/epos.inics").fullPath());
  
   if (path_fnii.length() >= 500)
     LogError("ReggeGribovPartonMCInterface") << "table path too long" << endl;
   else
     nfname_.nfnii = path_fnii.length();
   if (path_fnie.length() >= 500)
     LogError("ReggeGribovPartonMCInterface") << "table path too long" << endl;
   else
     nfname_.nfnie = path_fnie.length();
   if (path_fnrj.length() >= 500)
     LogError("ReggeGribovPartonMCInterface") << "table path too long" << endl;
   else
     nfname_.nfnrj = path_fnrj.length();
   if (path_fncs.length() >= 500)
     LogError("ReggeGribovPartonMCInterface") << "table path too long" << endl;
   else
     nfname_.nfncs = path_fncs.length();
  
   strcpy(fname_.fnii, path_fnii.c_str());
   strcpy(fname_.fnie, path_fnie.c_str());
   strcpy(fname_.fnrj, path_fnrj.c_str());
   strcpy(fname_.fncs, path_fncs.c_str());
  
   //qgsjet
   string path_fndat(FileInPath("GeneratorInterface/ReggeGribovPartonMCInterface/data/qgsjet.dat").fullPath());
   string path_fnncs(FileInPath("GeneratorInterface/ReggeGribovPartonMCInterface/data/qgsjet.ncs").fullPath());
  
   if (path_fndat.length() >= 500)
     LogError("ReggeGribovPartonMCInterface") << "table path too long" << endl;
   else
     qgsnfname_.nfndat = path_fndat.length();
   if (path_fnncs.length() >= 500)
     LogError("ReggeGribovPartonMCInterface") << "table path too long" << endl;
   else
     qgsnfname_.nfnncs = path_fnncs.length();
  
   strcpy(qgsfname_.fndat, path_fndat.c_str());
   strcpy(qgsfname_.fnncs, path_fnncs.c_str());
  
   qgsfname_.ifdat = 1;  //option to overwrite the normal path
   qgsfname_.ifncs = 2;
  
   //qgsjetII
   string path_fniidat(FileInPath("GeneratorInterface/ReggeGribovPartonMCInterface/data/qgsdat-II-04.lzma").fullPath());
   string path_fniincs(FileInPath("GeneratorInterface/ReggeGribovPartonMCInterface/data/sectnu-II-04").fullPath());
  
   if (path_fniidat.length() >= 500)
     LogError("ReggeGribovPartonMCInterface") << "table path too long" << endl;
   else
     qgsiinfname_.nfniidat = path_fniidat.length();
   if (path_fniincs.length() >= 500)
     LogError("ReggeGribovPartonMCInterface") << "table path too long" << endl;
   else
     qgsiinfname_.nfniincs = path_fniincs.length();
  
   strcpy(qgsiifname_.fniidat, path_fniidat.c_str());
   strcpy(qgsiifname_.fniincs, path_fniincs.c_str());
  
   qgsiifname_.ifiidat = 1;  //option to overwrite the normal path
   qgsiifname_.ifiincs = 2;
  
   return true;
 }