#include <Hydjet2Hadronizer.h>

Inheritance diagram for gen::Hydjet2Hadronizer:

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 ()

	Hydjet2Hadronizer (const edm::ParameterSet &, edm::ConsumesCollector &&)

bool	initializeForExternalPartons ()

bool	initializeForInternalPartons ()

bool	readSettings (int)

bool	residualDecay ()

void	statistics ()

	~Hydjet2Hadronizer () 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< HepMC3::GenEvent >	getGenEvent3 ()

std::unique_ptr< GenEventInfoProduct >	getGenEventInfo ()

std::unique_ptr< GenEventInfoProduct3 >	getGenEventInfo3 ()

virtual std::unique_ptr< GenLumiInfoHeader >	getGenLumiInfoHeader () const

GenRunInfoProduct &	getGenRunInfo ()

std::unique_ptr< lhef::LHEEvent >	getLHEEvent ()

const std::shared_ptr< lhef::LHERunInfo > &	getLHERunInfo () const

unsigned int	getVHepMC ()

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	resetEvent3 (std::unique_ptr< HepMC3::GenEvent > event3)

void	resetEventInfo (std::unique_ptr< GenEventInfoProduct > eventInfo)

void	resetEventInfo3 (std::unique_ptr< GenEventInfoProduct3 > 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	add_heavy_ion_rec (HepMC::GenEvent *evt)

HepMC::GenParticle *	build_hyjet2 (int index, int barcode)

HepMC::GenVertex *	build_hyjet2_vertex (int i, int id)

int	convertStatus (int)

int	convertStatusForComponents (int, int, int)

void	doSetRandomEngine (CLHEP::HepRandomEngine *v) override

std::vector< std::string > const &	doSharedResources () const override

bool	get_particles (HepMC::GenEvent *evt)

double	nuclear_radius () const

void	rotateEvtPlane ()

Private Attributes
double	cosphi0_

bool	ev = false

HepMC::GenEvent *	evt

InitialParamsHydjet_t	fParams

HepMC::FourVector *	fVertex_

Hydjet2 *	hj2

unsigned int	maxEventsToPrint_

int	nhard_

int	nsoft_

int	nsub_

double	phi0_

edm::ParameterSet	pset

Pythia6Service *	pythia6Service_

unsigned int	pythiaPylistVerbosity_

bool	rotate_

bool	separateHydjetComponents_

double	Sigin

double	Sigjet

std::vector< double >	signalVtx_

double	sinphi0_

edm::EDGetTokenT< CrossingFrame< edm::HepMCProduct > >	src_

Static Private Attributes
static const std::vector< std::string >	theSharedResources = {edm::SharedResourceNames::kPythia6}

Additional Inherited Members
Protected Member Functions inherited from gen::BaseHadronizer
std::unique_ptr< HepMC::GenEvent > &	event ()

std::unique_ptr< HepMC3::GenEvent > &	event3 ()

std::unique_ptr< GenEventInfoProduct > &	eventInfo ()

std::unique_ptr< GenEventInfoProduct3 > &	eventInfo3 ()

lhef::LHEEvent *	lheEvent ()

lhef::LHERunInfo *	lheRunInfo ()

GenRunInfoProduct &	runInfo ()

Protected Attributes inherited from gen::BaseHadronizer
unsigned int	ivhepmc = 2

std::string	lheFile_

int	randomIndex_

Detailed Description

Definition at line 44 of file Hydjet2Hadronizer.h.

Constructor & Destructor Documentation

◆ Hydjet2Hadronizer()

Hydjet2Hadronizer::Hydjet2Hadronizer	(	const edm::ParameterSet &	pset,
		edm::ConsumesCollector &&	iC
	)

Definition at line 82 of file Hydjet2Hadronizer.cc.

References edm::ParameterSet::exists(), DeadROC_duringRun::f1, DeadROC_duringRun::f2, fParams, fVertex_, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), HLT_2022v15_cff::InputTag, LogDebug, maxEventsToPrint_, pset, pythiaPylistVerbosity_, separateHydjetComponents_, signalVtx_, and src_.

     : BaseHadronizer(pset),
       rotate_(pset.getParameter<bool>("rotateEventPlane")),
       evt(nullptr),
       nsub_(0),
       nhard_(0),
       nsoft_(0),
       phi0_(0.),
       sinphi0_(0.),
       cosphi0_(1.),
       fVertex_(nullptr),
       pythia6Service_(new Pythia6Service(pset))
 
 {
   fParams.doPrintInfo = false;
   fParams.allowEmptyEvent = false;
   fParams.fNevnt = 0;                                      //not used in CMSSW
   fParams.femb = pset.getParameter<int>("embeddingMode");  //
   fParams.fSqrtS = pset.getParameter<double>("fSqrtS");    // C.m.s. energy per nucleon pair
   fParams.fAw = pset.getParameter<double>("fAw");          // Atomic weigth of nuclei, fAw
   fParams.fIfb = pset.getParameter<int>(
       "fIfb");  // Flag of type of centrality generation, fBfix (=0 is fixed by fBfix, >0 distributed [fBfmin, fBmax])
   fParams.fBmin = pset.getParameter<double>("fBmin");  // Minimum impact parameter in units of nuclear radius, fBmin
   fParams.fBmax = pset.getParameter<double>("fBmax");  // Maximum impact parameter in units of nuclear radius, fBmax
   fParams.fBfix = pset.getParameter<double>("fBfix");  // Fixed impact parameter in units of nuclear radius, fBfix
   fParams.fT = pset.getParameter<double>("fT");        // Temperature at chemical freeze-out, fT [GeV]
   fParams.fMuB = pset.getParameter<double>("fMuB");    // Chemical baryon potential per unit charge, fMuB [GeV]
   fParams.fMuS = pset.getParameter<double>("fMuS");    // Chemical strangeness potential per unit charge, fMuS [GeV]
   fParams.fMuC = pset.getParameter<double>(
       "fMuC");  // Chemical charm potential per unit charge, fMuC [GeV] (used if charm production is turned on)
   fParams.fMuI3 = pset.getParameter<double>("fMuI3");  // Chemical isospin potential per unit charge, fMuI3 [GeV]
   fParams.fThFO = pset.getParameter<double>("fThFO");  // Temperature at thermal freeze-out, fThFO [GeV]
   fParams.fMu_th_pip =
       pset.getParameter<double>("fMu_th_pip");  // Chemical potential of pi+ at thermal freeze-out, fMu_th_pip [GeV]
   fParams.fTau = pset.getParameter<double>(
       "fTau");  // Proper time proper at thermal freeze-out for central collisions, fTau [fm/c]
   fParams.fSigmaTau = pset.getParameter<double>(
       "fSigmaTau");  // Duration of emission at thermal freeze-out for central collisions, fSigmaTau [fm/c]
   fParams.fR = pset.getParameter<double>(
       "fR");  // Maximal transverse radius at thermal freeze-out for central collisions, fR [fm]
   fParams.fYlmax =
       pset.getParameter<double>("fYlmax");  // Maximal longitudinal flow rapidity at thermal freeze-out, fYlmax
   fParams.fUmax = pset.getParameter<double>(
       "fUmax");  // Maximal transverse flow rapidity at thermal freeze-out for central collisions, fUmax
   fParams.frhou2 = pset.getParameter<double>("fRhou2");  //parameter to swich ON/OFF = 0) rhou2
   fParams.frhou3 = pset.getParameter<double>("fRhou3");  //parameter to swich ON/OFF(0) rhou3
   fParams.frhou4 = pset.getParameter<double>("fRhou4");  //parameter to swich ON/OFF(0) rhou4
   fParams.fDelta =
       pset.getParameter<double>("fDelta");  // Momentum azimuthal anizotropy parameter at thermal freeze-out, fDelta
   fParams.fEpsilon =
       pset.getParameter<double>("fEpsilon");  // Spatial azimuthal anisotropy parameter at thermal freeze-out, fEpsilon
   fParams.fv2 = pset.getParameter<double>("fKeps2");  //parameter to swich ON/OFF(0) epsilon2 fluctuations
   fParams.fv3 = pset.getParameter<double>("fKeps3");  //parameter to swich ON/OFF(0) epsilon3 fluctuations
   fParams.fIfDeltaEpsilon = pset.getParameter<double>(
       "fIfDeltaEpsilon");  // Flag to specify fDelta and fEpsilon values, fIfDeltaEpsilon (=0 user's ones, >=1 calculated)
   fParams.fDecay =
       pset.getParameter<int>("fDecay");  // Flag to switch on/off hadron decays, fDecay (=0 decays off, >=1 decays on)
   fParams.fWeakDecay = pset.getParameter<double>(
       "fWeakDecay");  // Low decay width threshold fWeakDecay[GeV]: width<fWeakDecay decay off, width>=fDecayWidth decay on; can be used to switch off weak decays
   fParams.fEtaType = pset.getParameter<double>(
       "fEtaType");  // Flag to choose longitudinal flow rapidity distribution, fEtaType (=0 uniform, >0 Gaussian with the dispersion Ylmax)
   fParams.fTMuType = pset.getParameter<double>(
       "fTMuType");  // Flag to use calculated T_ch, mu_B and mu_S as a function of fSqrtS, fTMuType (=0 user's ones, >0 calculated)
   fParams.fCorrS = pset.getParameter<double>(
       "fCorrS");  // Strangeness supression factor gamma_s with fCorrS value (0<fCorrS <=1, if fCorrS <= 0 then it is calculated)
   fParams.fCharmProd = pset.getParameter<int>(
       "fCharmProd");  // Flag to include thermal charm production, fCharmProd (=0 no charm production, >=1 charm production)
   fParams.fCorrC = pset.getParameter<double>(
       "fCorrC");  // Charmness enhancement factor gamma_c with fCorrC value (fCorrC >0, if fCorrC<0 then it is calculated)
   fParams.fNhsel = pset.getParameter<int>(
       "fNhsel");  //Flag to include jet (J)/jet quenching (JQ) and hydro (H) state production, fNhsel (0 H on & J off, 1 H/J on & JQ off, 2 H/J/HQ on, 3 J on & H/JQ off, 4 H off & J/JQ on)
   fParams.fPyhist = pset.getParameter<int>(
       "fPyhist");  // Flag to suppress the output of particle history from PYTHIA, fPyhist (=1 only final state particles; =0 full particle history from PYTHIA)
   fParams.fIshad = pset.getParameter<int>(
       "fIshad");  // Flag to switch on/off nuclear shadowing, fIshad (0 shadowing off, 1 shadowing on)
   fParams.fPtmin =
       pset.getParameter<double>("fPtmin");  // Minimal pt of parton-parton scattering in PYTHIA event, fPtmin [GeV/c]
   fParams.fT0 = pset.getParameter<double>(
       "fT0");  // Initial QGP temperature for central Pb+Pb collisions in mid-rapidity, fT0 [GeV]
   fParams.fTau0 = pset.getParameter<double>("fTau0");  // Proper QGP formation time in fm/c, fTau0 (0.01<fTau0<10)
   fParams.fNf = pset.getParameter<int>("fNf");         // Number of active quark flavours in QGP, fNf (0, 1, 2 or 3)
   fParams.fIenglu = pset.getParameter<int>(
       "fIenglu");  // Flag to fix type of partonic energy loss, fIenglu (0 radiative and collisional loss, 1 radiative loss only, 2 collisional loss only)
   fParams.fIanglu = pset.getParameter<int>(
       "fIanglu");  // Flag to fix type of angular distribution of in-medium emitted gluons, fIanglu (0 small-angular, 1 wide-angular, 2 collinear).
 
   edm::FileInPath f1("externals/hydjet2/particles.data");
   strcpy(fParams.partDat, (f1.fullPath()).c_str());
 
   edm::FileInPath f2("externals/hydjet2/tabledecay.txt");
   strcpy(fParams.tabDecay, (f2.fullPath()).c_str());
 
   fParams.fPythiaTune = false;
 
   if (pset.exists("signalVtx"))
     signalVtx_ = pset.getUntrackedParameter<std::vector<double>>("signalVtx");
 
   if (signalVtx_.size() == 4) {
     if (!fVertex_)
       fVertex_ = new HepMC::FourVector();
     LogDebug("EventSignalVertex") << "Setting event signal vertex "
                                   << " x = " << signalVtx_.at(0) << " y = " << signalVtx_.at(1)
                                   << "  z= " << signalVtx_.at(2) << " t = " << signalVtx_.at(3) << endl;
     fVertex_->set(signalVtx_.at(0), signalVtx_.at(1), signalVtx_.at(2), signalVtx_.at(3));
   }
 
   // PYLIST Verbosity Level
   // Valid PYLIST arguments are: 1, 2, 3, 5, 7, 11, 12, 13
   pythiaPylistVerbosity_ = pset.getUntrackedParameter<int>("pythiaPylistVerbosity", 0);
   LogDebug("PYLISTverbosity") << "Pythia PYLIST verbosity level = " << pythiaPylistVerbosity_;
   //Max number of events printed on verbosity level
   maxEventsToPrint_ = pset.getUntrackedParameter<int>("maxEventsToPrint", 0);
   LogDebug("Events2Print") << "Number of events to be printed = " << maxEventsToPrint_;
   if (fParams.femb == 1) {
     fParams.fIfb = 0;
     src_ = iC.consumes<CrossingFrame<edm::HepMCProduct>>(
         pset.getUntrackedParameter<edm::InputTag>("backgroundLabel", edm::InputTag("mix", "generatorSmeared")));
   }
 
   separateHydjetComponents_ = pset.getUntrackedParameter<bool>("separateHydjetComponents", false);
 }

◆ ~Hydjet2Hadronizer()

Hydjet2Hadronizer::~Hydjet2Hadronizer ( )

override

Definition at line 204 of file Hydjet2Hadronizer.cc.

References pythia6Service_.

                                       {
   call_pystat(1);
   delete pythia6Service_;
 }

Member Function Documentation

◆ add_heavy_ion_rec()

void Hydjet2Hadronizer::add_heavy_ion_rec ( HepMC::GenEvent * evt )

private

Definition at line 516 of file Hydjet2Hadronizer.cc.

References evt, hj2, nsub_, nuclear_radius(), phi0_, and Sigin.

Referenced by generatePartonsAndHadronize().

                                                             {
   // heavy ion record in the final CMSSW Event
   int nproj = static_cast<int>((hj2->GetNpart()) / 2);
   int ntarg = static_cast<int>((hj2->GetNpart()) - nproj);
 
   HepMC::HeavyIon *hi = new HepMC::HeavyIon(nsub_,                              // Ncoll_hard/N of SubEvents
                                             nproj,                              // Npart_proj
                                             ntarg,                              // Npart_targ
                                             hj2->GetNbcol(),                    // Ncoll
                                             0,                                  // spectator_neutrons
                                             0,                                  // spectator_protons
                                             0,                                  // N_Nwounded_collisions
                                             0,                                  // Nwounded_N_collisions
                                             0,                                  // Nwounded_Nwounded_collisions
                                             hj2->GetBgen() * nuclear_radius(),  // impact_parameter in [fm]
                                             phi0_,                              // event_plane_angle
                                             hj2->GetPsiv3(),                    // eccentricity <<<---- psi for v3!!!
                                             Sigin                               // sigma_inel_NN
   );
 
   evt->set_heavy_ion(*hi);
   delete hi;
 }

◆ build_hyjet2()

HepMC::GenParticle * Hydjet2Hadronizer::build_hyjet2	(	int	index,
		int	barcode
	)

private

Definition at line 475 of file Hydjet2Hadronizer.cc.

References convertStatusForComponents(), cosphi0_, GenParticle::GenParticle, hj2, gen::p, multPhiCorr_741_25nsDY_cfi::px, multPhiCorr_741_25nsDY_cfi::py, and sinphi0_.

Referenced by get_particles().

                                                                         {
   // Build particle object corresponding to index in hyjets (soft+hard)
   double px0 = (hj2->GetPx()).at(index);
   double py0 = (hj2->GetPy()).at(index);
 
   double px = px0 * cosphi0_ - py0 * sinphi0_;
   double py = py0 * cosphi0_ + px0 * sinphi0_;
 
   HepMC::GenParticle *p = new HepMC::GenParticle(
       HepMC::FourVector(px,                        // px
                         py,                        // py
                         (hj2->GetPz()).at(index),  // pz
                         (hj2->GetE()).at(index)),  // E
       (hj2->GetPdg()).at(index),                   // id
       convertStatusForComponents(
           (hj2->GetFinal()).at(index), (hj2->GetType()).at(index), (hj2->GetPythiaStatus()).at(index))  // status
   );
 
   p->suggest_barcode(barcode);
   return p;
 }

◆ build_hyjet2_vertex()

HepMC::GenVertex * Hydjet2Hadronizer::build_hyjet2_vertex	(	int	i,
		int	id
	)

private

Definition at line 498 of file Hydjet2Hadronizer.cc.

References cosphi0_, MillePedeFileConverter_cfg::e, hj2, mps_fire::i, sinphi0_, submitPVValidationJobs::t, bphysicsOniaDQM_cfi::vertex, and x.

Referenced by get_particles().

                                                                     {
   // build verteces for the hyjets stored events
   double x0 = (hj2->GetX()).at(i);
   double y0 = (hj2->GetY()).at(i);
 
   // convert to mm (as in PYTHIA6)
   const double fm_to_mm = 1e-12;
   double x = fm_to_mm * (x0 * cosphi0_ - y0 * sinphi0_);
   double y = fm_to_mm * (y0 * cosphi0_ + x0 * sinphi0_);
   double z = fm_to_mm * (hj2->GetZ()).at(i);
   double t = fm_to_mm * (hj2->GetT()).at(i);
 
   HepMC::GenVertex *vertex = new HepMC::GenVertex(HepMC::FourVector(x, y, z, t), id);
 
   return vertex;
 }

◆ classname()

const char * Hydjet2Hadronizer::classname ( ) const

Definition at line 385 of file Hydjet2Hadronizer.cc.

385 { return "gen::Hydjet2Hadronizer"; }

◆ convertStatus()

int Hydjet2Hadronizer::convertStatus ( int st )

private

Definition at line 66 of file Hydjet2Hadronizer.cc.

Referenced by get_particles().

                                            {
   if (st <= 0)
     return 0;
   if (st <= 10)
     return 1;
   if (st <= 20)
     return 2;
   if (st <= 30)
     return 3;
   else
     return -1;
 }

◆ convertStatusForComponents()

int Hydjet2Hadronizer::convertStatusForComponents	(	int	sta,
		int	typ,
		int	pySta
	)

private

Definition at line 43 of file Hydjet2Hadronizer.cc.

References Exception.

Referenced by build_hyjet2(), and get_particles().

                                                                              {
   int st = -1;
   if (typ == 0)  //soft
     st = 2 - sta;
   else if (typ == 1)
     st = convertStatus(pySta);
 
   if (st == -1)
     throw cms::Exception("ConvertStatus") << "Wrong status code!" << endl;
 
   if (separateHydjetComponents_) {
     if (st == 1 && typ == 0)
       return 6;
     if (st == 1 && typ == 1)
       return 7;
     if (st == 2 && typ == 0)
       return 16;
     if (st == 2 && typ == 1)
       return 17;
   }
   return st;
 }

◆ decay()

bool Hydjet2Hadronizer::decay ( )

Definition at line 381 of file Hydjet2Hadronizer.cc.

381 { return true; }

◆ declareSpecialSettings()

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

inline

Definition at line 50 of file Hydjet2Hadronizer.h.

50 { return true; }

◆ declareStableParticles()

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

Definition at line 368 of file Hydjet2Hadronizer.cc.

References gen::call_pygive(), gather_cfg::cout, mps_fire::i, and gen::pycomp_().

                                                                          {
   std::vector<int> pdg = _pdg;
   for (size_t i = 0; i < pdg.size(); i++) {
     int pyCode = pycomp_(pdg[i]);
     std::ostringstream pyCard;
     pyCard << "MDCY(" << pyCode << ",1)=0";
     std::cout << pyCard.str() << std::endl;
     call_pygive(pyCard.str());
   }
   return true;
 }

◆ doSetRandomEngine()

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

overrideprivatevirtual

Reimplemented from gen::BaseHadronizer.

Definition at line 210 of file Hydjet2Hadronizer.cc.

References hjRandomEngine, pythia6Service_, gen::Pythia6Service::setRandomEngine(), and gen::v.

                                                                  {
   pythia6Service_->setRandomEngine(v);
   hjRandomEngine = v;
 }

◆ doSharedResources()

std::vector<std::string> const& gen::Hydjet2Hadronizer::doSharedResources ( ) const

inlineoverrideprivatevirtual

Reimplemented from gen::BaseHadronizer.

Definition at line 71 of file Hydjet2Hadronizer.h.

References theSharedResources.

71 { return theSharedResources; }

gen::Hydjet2Hadronizer::theSharedResources

static const std::vector< std::string > theSharedResources

Definition: Hydjet2Hadronizer.h:72

◆ finalizeEvent()

void Hydjet2Hadronizer::finalizeEvent ( )

Definition at line 383 of file Hydjet2Hadronizer.cc.

383 {}

◆ generatePartonsAndHadronize()

bool Hydjet2Hadronizer::generatePartonsAndHadronize ( )

Definition at line 244 of file Hydjet2Hadronizer.cc.

References add_heavy_ion_rec(), funct::cos(), cosphi0_, MillePedeFileConverter_cfg::e, ev, gen::BaseHadronizer::event(), edm::errors::EventCorruption, evt, Exception, fParams, fVertex_, get_particles(), gen::BaseHadronizer::getEDMEvent(), edm::HepMCProduct::GetEvent(), hj2, edm::HepMCProduct::isVtxGenApplied(), LogDebug, GeneratorMix_cff::mix, eostools::move(), nhard_, nsoft_, nsub_, nuclear_radius(), phi0_, edm::Handle< T >::product(), pypars, pythia6Service_, rotate_, rotateEvtPlane(), Sigin, Sigjet, funct::sin(), sinphi0_, and src_.

                                                     {
   Pythia6Service::InstanceWrapper guard(pythia6Service_);
 
   // generate single event
   if (fParams.femb == 1) {
     const edm::Event &e = getEDMEvent();
     HepMC::GenVertex *genvtx = nullptr;
     const HepMC::GenEvent *inev = nullptr;
     Handle<CrossingFrame<HepMCProduct>> cf;
     e.getByToken(src_, cf);
     MixCollection<HepMCProduct> mix(cf.product());
     if (mix.size() < 1) {
       throw cms::Exception("MatchVtx") << "Mixing has " << mix.size() << " sub-events, should have been at least 1"
                                        << endl;
     }
     const HepMCProduct &bkg = mix.getObject(0);
     if (!(bkg.isVtxGenApplied())) {
       throw cms::Exception("MatchVtx") << "Input background does not have smeared vertex!" << endl;
     } else {
       inev = bkg.GetEvent();
     }
 
     genvtx = inev->signal_process_vertex();
 
     if (!genvtx)
       throw cms::Exception("MatchVtx") << "Input background does not have signal process vertex!" << endl;
 
     double aX, aY, aZ, aT;
 
     aX = genvtx->position().x();
     aY = genvtx->position().y();
     aZ = genvtx->position().z();
     aT = genvtx->position().t();
 
     if (!fVertex_) {
       fVertex_ = new HepMC::FourVector();
     }
     LogInfo("MatchVtx") << " setting vertex "
                         << " aX " << aX << " aY " << aY << " aZ " << aZ << " aT " << aT << endl;
     fVertex_->set(aX, aY, aZ, aT);
 
     const HepMC::HeavyIon *hi = inev->heavy_ion();
 
     if (hi) {
       fParams.fBfix = (hi->impact_parameter()) / nuclear_radius();
       phi0_ = hi->event_plane_angle();
       sinphi0_ = sin(phi0_);
       cosphi0_ = cos(phi0_);
     } else {
       LogWarning("EventEmbedding") << "Background event does not have heavy ion record!";
     }
 
   } else if (rotate_)
     rotateEvtPlane();
 
   nsoft_ = 0;
   nhard_ = 0;
 
   // generate one HYDJET event
   int ntry = 0;
 
   while (nsoft_ == 0 && nhard_ == 0) {
     if (ntry > 100) {
       LogError("Hydjet2EmptyEvent") << "##### HYDJET2: No Particles generated, Number of tries =" << ntry;
       // Throw an exception. Use the EventCorruption exception since it maps onto SkipEvent
       // which is what we want to do here.
       std::ostringstream sstr;
       sstr << "Hydjet2HadronizerProducer: No particles generated after " << ntry << " tries.\n";
       edm::Exception except(edm::errors::EventCorruption, sstr.str());
       throw except;
     } else {
       hj2->GenerateEvent(fParams.fBfix);
 
       if (hj2->IsEmpty()) {
         continue;
       }
 
       nsoft_ = hj2->GetNhyd();
       nsub_ = hj2->GetNjet();
       nhard_ = hj2->GetNpyt();
 
       //100 trys
       ++ntry;
     }
   }
 
   if (ev == 0) {
     Sigin = hj2->GetSigin();
     Sigjet = hj2->GetSigjet();
   }
   ev = true;
 
   if (fParams.fNhsel < 3)
     nsub_++;
 
   // event information
   std::unique_ptr<HepMC::GenEvent> evt = std::make_unique<HepMC::GenEvent>();
   std::unique_ptr<edm::HepMCProduct> HepMCEvt = std::make_unique<edm::HepMCProduct>();
 
   if (nhard_ > 0 || nsoft_ > 0)
     get_particles(evt.get());
 
   evt->set_signal_process_id(pypars.msti[0]);  // type of the process
   evt->set_event_scale(pypars.pari[16]);       // Q^2
   add_heavy_ion_rec(evt.get());
 
   if (fVertex_) {
     // generate new vertex & apply the shift
     // Copy the HepMC::GenEvent
     HepMCEvt = std::make_unique<edm::HepMCProduct>(evt.get());
     HepMCEvt->applyVtxGen(fVertex_);
     evt = std::make_unique<HepMC::GenEvent>(*HepMCEvt->GetEvent());
   }
 
   HepMC::HEPEVT_Wrapper::check_hepevt_consistency();
   LogDebug("HEPEVT_info") << "Ev numb: " << HepMC::HEPEVT_Wrapper::event_number()
                           << " Entries number: " << HepMC::HEPEVT_Wrapper::number_entries() << " Max. entries "
                           << HepMC::HEPEVT_Wrapper::max_number_entries() << std::endl;
 
   event() = std::move(evt);
   return kTRUE;
 }

◆ get_particles()

bool Hydjet2Hadronizer::get_particles ( HepMC::GenEvent * evt )

private

Definition at line 396 of file Hydjet2Hadronizer.cc.

References funct::abs(), build_hyjet2(), build_hyjet2_vertex(), convertStatus(), convertStatusForComponents(), evt, GenParticle::GenParticle, hj2, LogDebug, nhard_, nsoft_, and nsub_.

Referenced by generatePartonsAndHadronize().

                                                         {
   LogDebug("Hydjet2") << " Number of sub events " << nsub_;
   LogDebug("Hydjet2") << " Number of hard events " << hj2->GetNjet();
   LogDebug("Hydjet2") << " Number of hard particles " << nhard_;
   LogDebug("Hydjet2") << " Number of soft particles " << nsoft_;
   LogDebug("Hydjet2") << " nhard_ + nsoft_ = " << nhard_ + nsoft_ << " Ntot = " << hj2->GetNtot() << endl;
 
   int ihy = 0;
   int isub_l = -1;
   int stab = 0;
 
   vector<HepMC::GenParticle *> particle(hj2->GetNtot());
   HepMC::GenVertex *sub_vertices = nullptr;
 
   while (ihy < hj2->GetNtot()) {
     if ((hj2->GetiJet().at(ihy)) != isub_l) {
       sub_vertices = new HepMC::GenVertex(HepMC::FourVector(0, 0, 0, 0), hj2->GetiJet().at(ihy));
       evt->add_vertex(sub_vertices);
       if (!evt->signal_process_vertex())
         evt->set_signal_process_vertex(sub_vertices);
       isub_l = hj2->GetiJet().at(ihy);
     }
 
     if ((convertStatusForComponents(
             (hj2->GetFinal()).at(ihy), (hj2->GetType()).at(ihy), (hj2->GetPythiaStatus().at(ihy)))) == 1)
       stab++;
     LogDebug("Hydjet2_array") << ihy << " MULTin ev.:" << hj2->GetNtot() << " SubEv.#" << hj2->GetiJet().at(ihy)
                               << " Part #" << ihy + 1 << ", PDG: " << hj2->GetPdg().at(ihy) << " (st. "
                               << convertStatus(hj2->GetPythiaStatus().at(ihy))
                               << ") mother=" << hj2->GetMotherIndex().at(ihy) + 1 << ", childs ("
                               << hj2->GetFirstDaughterIndex().at(ihy) + 1 << "-"
                               << hj2->GetLastDaughterIndex().at(ihy) + 1 << "), vtx (" << hj2->GetX().at(ihy) << ","
                               << hj2->GetY().at(ihy) << "," << hj2->GetZ().at(ihy) << ") " << std::endl;
 
     if ((hj2->GetMotherIndex().at(ihy)) <= 0) {
       particle.at(ihy) = build_hyjet2(ihy, ihy + 1);
       if (!sub_vertices)
         LogError("Hydjet2_array") << "##### HYDJET2: Vertex not initialized!";
       else
         sub_vertices->add_particle_out(particle.at(ihy));
       LogDebug("Hydjet2_array") << " ---> " << ihy + 1 << std::endl;
     } else {
       particle.at(ihy) = build_hyjet2(ihy, ihy + 1);
       int mid = hj2->GetMotherIndex().at(ihy);
 
       while (((mid + 1) < ihy) && (std::abs(hj2->GetPdg().at(mid)) < 100) &&
              ((hj2->GetFirstDaughterIndex().at(mid + 1)) <= ihy)) {
         mid++;
         LogDebug("Hydjet2_array") << "======== MID changed to " << mid
                                   << " ======== PDG(mid) = " << hj2->GetPdg().at(mid) << std::endl;
       }
 
       if (std::abs(hj2->GetPdg().at(mid)) < 100) {
         mid = hj2->GetMotherIndex().at(ihy);
         LogDebug("Hydjet2_array") << "======== MID changed BACK to " << mid
                                   << " ======== PDG(mid) = " << hj2->GetPdg().at(mid) << std::endl;
       }
 
       HepMC::GenParticle *mother = particle.at(mid);
       HepMC::GenVertex *prod_vertex = mother->end_vertex();
 
       if (!prod_vertex) {
         prod_vertex = build_hyjet2_vertex(ihy, (hj2->GetiJet().at(ihy)));
         prod_vertex->add_particle_in(mother);
         LogDebug("Hydjet2_array") << " <--- " << mid + 1 << std::endl;
         evt->add_vertex(prod_vertex);
       }
       prod_vertex->add_particle_out(particle.at(ihy));
       LogDebug("Hydjet2_array") << " ---" << mid + 1 << "---> " << ihy + 1 << std::endl;
     }
     ihy++;
   }
 
   LogDebug("Hydjet2_array") << " MULTin ev.:" << hj2->GetNtot() << ", last index: " << ihy - 1
                             << ", stable particles: " << stab << std::endl;
   return kTRUE;
 }

◆ hadronize()

bool Hydjet2Hadronizer::hadronize ( )

Definition at line 380 of file Hydjet2Hadronizer.cc.

380 { return false; }

◆ initializeForExternalPartons()

bool gen::Hydjet2Hadronizer::initializeForExternalPartons ( )

◆ initializeForInternalPartons()

bool Hydjet2Hadronizer::initializeForInternalPartons ( )

Definition at line 228 of file Hydjet2Hadronizer.cc.

References fParams, hj2, nuclear_radius(), and pythia6Service_.

                                                      {
   Pythia6Service::InstanceWrapper guard(pythia6Service_);
 
   // the input impact parameter (bxx_) is in [fm]; transform in [fm/RA] for hydjet usage
   const double ra = nuclear_radius();
   LogInfo("Hydjet2Hadronizer|RAScaling") << "Nuclear radius(RA) =  " << ra;
   fParams.fBmin /= ra;
   fParams.fBmax /= ra;
   fParams.fBfix /= ra;
 
   hj2 = new Hydjet2(fParams);
 
   return kTRUE;
 }

◆ nuclear_radius()

double gen::Hydjet2Hadronizer::nuclear_radius ( ) const

inlineprivate

Definition at line 107 of file Hydjet2Hadronizer.h.

References fParams, and funct::pow().

Referenced by add_heavy_ion_rec(), generatePartonsAndHadronize(), and initializeForInternalPartons().

                                                  {
     // Return the nuclear radius derived from the
     // beam/target atomic mass number.
 
     return 1.15 * pow((double)fParams.fAw, 1. / 3.);
   }

◆ readSettings()

bool Hydjet2Hadronizer::readSettings ( int )

Definition at line 216 of file Hydjet2Hadronizer.cc.

References fParams, hjRandomEngine, pythia6Service_, and gen::Pythia6Service::setGeneralParams().

                                         {
   Pythia6Service::InstanceWrapper guard(pythia6Service_);
   pythia6Service_->setGeneralParams();
 
   fParams.fSeed = hjRandomEngine->CLHEP::HepRandomEngine::getSeed();
   LogInfo("Hydjet2Hadronizer|GenSeed") << "Seed for random number generation: "
                                        << hjRandomEngine->CLHEP::HepRandomEngine::getSeed();
 
   return kTRUE;
 }

◆ residualDecay()

bool Hydjet2Hadronizer::residualDecay ( )

Definition at line 382 of file Hydjet2Hadronizer.cc.

382 { return true; }

◆ rotateEvtPlane()

void Hydjet2Hadronizer::rotateEvtPlane ( )

private

Definition at line 388 of file Hydjet2Hadronizer.cc.

References funct::cos(), cosphi0_, phi0_, pi, gen::pyr_(), funct::sin(), and sinphi0_.

Referenced by generatePartonsAndHadronize().

                                        {
   const double pi = 3.14159265358979;
   phi0_ = 2. * pi * gen::pyr_(nullptr) - pi;
   sinphi0_ = sin(phi0_);
   cosphi0_ = cos(phi0_);
 }

◆ statistics()

void Hydjet2Hadronizer::statistics ( )

Definition at line 384 of file Hydjet2Hadronizer.cc.

384 {}

Member Data Documentation

◆ cosphi0_

double gen::Hydjet2Hadronizer::cosphi0_

private

Definition at line 91 of file Hydjet2Hadronizer.h.

Referenced by build_hyjet2(), build_hyjet2_vertex(), generatePartonsAndHadronize(), and rotateEvtPlane().

◆ ev

bool gen::Hydjet2Hadronizer::ev = false

private

Definition at line 82 of file Hydjet2Hadronizer.h.

Referenced by generatePartonsAndHadronize().

◆ evt

HepMC::GenEvent* gen::Hydjet2Hadronizer::evt

private

Definition at line 85 of file Hydjet2Hadronizer.h.

Referenced by add_heavy_ion_rec(), generatePartonsAndHadronize(), and get_particles().

◆ fParams

InitialParamsHydjet_t gen::Hydjet2Hadronizer::fParams

private

Definition at line 79 of file Hydjet2Hadronizer.h.

Referenced by generatePartonsAndHadronize(), Hydjet2Hadronizer(), initializeForInternalPartons(), nuclear_radius(), and readSettings().

◆ fVertex_

HepMC::FourVector* gen::Hydjet2Hadronizer::fVertex_

private

Definition at line 99 of file Hydjet2Hadronizer.h.

Referenced by generatePartonsAndHadronize(), and Hydjet2Hadronizer().

◆ hj2

Hydjet2* gen::Hydjet2Hadronizer::hj2

private

Definition at line 80 of file Hydjet2Hadronizer.h.

Referenced by add_heavy_ion_rec(), build_hyjet2(), build_hyjet2_vertex(), generatePartonsAndHadronize(), get_particles(), and initializeForInternalPartons().

◆ maxEventsToPrint_

unsigned int gen::Hydjet2Hadronizer::maxEventsToPrint_

private

Definition at line 94 of file Hydjet2Hadronizer.h.

Referenced by Hydjet2Hadronizer().

◆ nhard_

int gen::Hydjet2Hadronizer::nhard_

private

Definition at line 87 of file Hydjet2Hadronizer.h.

Referenced by generatePartonsAndHadronize(), and get_particles().

◆ nsoft_

int gen::Hydjet2Hadronizer::nsoft_

private

Definition at line 88 of file Hydjet2Hadronizer.h.

Referenced by generatePartonsAndHadronize(), and get_particles().

◆ nsub_

int gen::Hydjet2Hadronizer::nsub_

private

Definition at line 86 of file Hydjet2Hadronizer.h.

Referenced by add_heavy_ion_rec(), generatePartonsAndHadronize(), and get_particles().

◆ phi0_

double gen::Hydjet2Hadronizer::phi0_

private

Definition at line 89 of file Hydjet2Hadronizer.h.

Referenced by add_heavy_ion_rec(), generatePartonsAndHadronize(), and rotateEvtPlane().

◆ pset

edm::ParameterSet gen::Hydjet2Hadronizer::pset

private

Definition at line 96 of file Hydjet2Hadronizer.h.

Referenced by Hydjet2Hadronizer().

◆ pythia6Service_

Pythia6Service* gen::Hydjet2Hadronizer::pythia6Service_

private

Definition at line 103 of file Hydjet2Hadronizer.h.

Referenced by doSetRandomEngine(), generatePartonsAndHadronize(), initializeForInternalPartons(), readSettings(), and ~Hydjet2Hadronizer().

◆ pythiaPylistVerbosity_

unsigned int gen::Hydjet2Hadronizer::pythiaPylistVerbosity_

private

Definition at line 93 of file Hydjet2Hadronizer.h.

Referenced by Hydjet2Hadronizer().

◆ rotate_

bool gen::Hydjet2Hadronizer::rotate_

private

Definition at line 84 of file Hydjet2Hadronizer.h.

Referenced by generatePartonsAndHadronize().

◆ separateHydjetComponents_

bool gen::Hydjet2Hadronizer::separateHydjetComponents_

private

Definition at line 83 of file Hydjet2Hadronizer.h.

Referenced by Hydjet2Hadronizer().

◆ Sigin

double gen::Hydjet2Hadronizer::Sigin

private

Definition at line 97 of file Hydjet2Hadronizer.h.

Referenced by add_heavy_ion_rec(), and generatePartonsAndHadronize().

◆ Sigjet

double gen::Hydjet2Hadronizer::Sigjet

private

Definition at line 97 of file Hydjet2Hadronizer.h.

Referenced by generatePartonsAndHadronize().

◆ signalVtx_

std::vector<double> gen::Hydjet2Hadronizer::signalVtx_

private

Definition at line 101 of file Hydjet2Hadronizer.h.

Referenced by Hydjet2Hadronizer().

◆ sinphi0_

double gen::Hydjet2Hadronizer::sinphi0_

private

Definition at line 90 of file Hydjet2Hadronizer.h.

Referenced by build_hyjet2(), build_hyjet2_vertex(), generatePartonsAndHadronize(), and rotateEvtPlane().

◆ src_

edm::EDGetTokenT<CrossingFrame<edm::HepMCProduct> > gen::Hydjet2Hadronizer::src_

private

Definition at line 104 of file Hydjet2Hadronizer.h.

Referenced by generatePartonsAndHadronize(), and Hydjet2Hadronizer().

◆ theSharedResources

const std::vector< std::string > Hydjet2Hadronizer::theSharedResources = {edm::SharedResourceNames::kPythia6}

staticprivate

Definition at line 72 of file Hydjet2Hadronizer.h.

Referenced by doSharedResources().

Public Member Functions

Private Member Functions

Private Attributes

Static Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ Hydjet2Hadronizer()

◆ ~Hydjet2Hadronizer()

Member Function Documentation

◆ add_heavy_ion_rec()

◆ build_hyjet2()

◆ build_hyjet2_vertex()

◆ classname()

◆ convertStatus()

◆ convertStatusForComponents()

◆ decay()

◆ declareSpecialSettings()

◆ declareStableParticles()

◆ doSetRandomEngine()

◆ doSharedResources()

◆ finalizeEvent()

◆ generatePartonsAndHadronize()

◆ get_particles()

◆ hadronize()

◆ initializeForExternalPartons()

◆ initializeForInternalPartons()

◆ nuclear_radius()

◆ readSettings()

◆ residualDecay()

◆ rotateEvtPlane()

◆ statistics()

Member Data Documentation

◆ cosphi0_

◆ ev

◆ evt

◆ fParams

◆ fVertex_

◆ hj2

◆ maxEventsToPrint_

◆ nhard_

◆ nsoft_

◆ nsub_

◆ phi0_

◆ pset

◆ pythia6Service_

◆ pythiaPylistVerbosity_

◆ rotate_

◆ separateHydjetComponents_

◆ Sigin

◆ Sigjet

◆ signalVtx_

◆ sinphi0_

◆ src_

◆ theSharedResources