#include <HydjetHadronizer.h>

Inheritance diagram for gen::HydjetHadronizer:

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

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

bool	initializeForExternalPartons ()

bool	initializeForInternalPartons ()

bool	readSettings (int)

bool	residualDecay ()

void	statistics ()

	~HydjetHadronizer () 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_hyjet (int index, int barcode)

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

bool	call_hyinit (double energy, double a, int ifb, double bmin, double bmax, double bfix, int nh)

int	convertStatus (int st)

void	doSetRandomEngine (CLHEP::HepRandomEngine *v) override

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

bool	get_particles (HepMC::GenEvent *evt)

bool	hydjet_init (const edm::ParameterSet &pset)

double	nuclear_radius () const

void	rotateEvtPlane ()

Private Attributes
double	abeamtarget_
	beam/target atomic mass number More...

int	angularspecselector_
	angular emitted gluon spectrum selection More...

double	bfixed_
	fixed impact param (fm); valid only if cflag_=0 More...

double	bmax_

double	bmin_

int	cflag_

double	comenergy
	collision energy More...

double	cosphi0_

bool	docollisionalenloss_
	DEFAULT = true. More...

bool	doradiativeenloss_
	DEFAULT = true. More...

bool	embedding_
	Switch for embedding mode. More...

HepMC::GenEvent *	evt

double	fracsoftmult_

HepMC::FourVector *	fVertex_
	Event signal vertex. More...

double	hadfreeztemp_

std::string	hymode_
	Hydjet running mode. More...

unsigned int	maxEventsToPrint_
	Events to print if verbosity. More...

double	maxlongy_

double	maxtrany_

int	nhard_
	multiplicity of PYTHIA(+PYQUEN)-induced particles in event More...

int	nmultiplicity_

unsigned int	nquarkflavor_

int	nsoft_
	multiplicity of HYDRO-induced particles in event More...

int	nsub_
	number of sub-events More...

double	phi0_
	Event plane angle. More...

edm::ParameterSet	pset_

Pythia6Service *	pythia6Service_

unsigned int	pythiaPylistVerbosity_
	pythia verbosity; def=1 More...

double	qgpt0_

double	qgptau0_

bool	rotate_
	Switch to rotate event plane. More...

unsigned int	shadowingswitch_

std::vector< double >	signalVtx_
	Pset double vector to set event signal vertex. More...

double	signn_

double	sinphi0_

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

Static Private Attributes
static const std::vector< std::string >	theSharedResources

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 40 of file HydjetHadronizer.h.

Constructor & Destructor Documentation

◆ HydjetHadronizer()

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

Definition at line 55 of file HydjetHadronizer.cc.

References cflag_, embedding_, fVertex_, HYJVER, ProducerED_cfi::InputTag, LogDebug, maxEventsToPrint_, muonDTDigis_cfi::pset, pythiaPylistVerbosity_, signalVtx_, and src_.

     : BaseHadronizer(pset),
       evt(nullptr),
       pset_(pset),
       abeamtarget_(pset.getParameter<double>("aBeamTarget")),
       angularspecselector_(pset.getParameter<int>("angularSpectrumSelector")),
       bfixed_(pset.getParameter<double>("bFixed")),
       bmax_(pset.getParameter<double>("bMax")),
       bmin_(pset.getParameter<double>("bMin")),
       cflag_(pset.getParameter<int>("cFlag")),
       embedding_(pset.getParameter<int>("embeddingMode")),
       comenergy(pset.getParameter<double>("comEnergy")),
       doradiativeenloss_(pset.getParameter<bool>("doRadiativeEnLoss")),
       docollisionalenloss_(pset.getParameter<bool>("doCollisionalEnLoss")),
       fracsoftmult_(pset.getParameter<double>("fracSoftMultiplicity")),
       hadfreeztemp_(pset.getParameter<double>("hadronFreezoutTemperature")),
       hymode_(pset.getParameter<string>("hydjetMode")),
       maxEventsToPrint_(pset.getUntrackedParameter<int>("maxEventsToPrint", 1)),
       maxlongy_(pset.getParameter<double>("maxLongitudinalRapidity")),
       maxtrany_(pset.getParameter<double>("maxTransverseRapidity")),
       nsub_(0),
       nhard_(0),
       nmultiplicity_(pset.getParameter<int>("nMultiplicity")),
       nsoft_(0),
       nquarkflavor_(pset.getParameter<int>("qgpNumQuarkFlavor")),
       pythiaPylistVerbosity_(pset.getUntrackedParameter<int>("pythiaPylistVerbosity", 0)),
       qgpt0_(pset.getParameter<double>("qgpInitialTemperature")),
       qgptau0_(pset.getParameter<double>("qgpProperTimeFormation")),
       phi0_(0.),
       sinphi0_(0.),
       cosphi0_(1.),
       rotate_(pset.getParameter<bool>("rotateEventPlane")),
       shadowingswitch_(pset.getParameter<int>("shadowingSwitch")),
       signn_(pset.getParameter<double>("sigmaInelNN")),
       fVertex_(nullptr),
       pythia6Service_(new Pythia6Service(pset)) {
   // Default constructor
 
   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 (embedding_) {
     cflag_ = 0;
     src_ = iC.consumes<CrossingFrame<edm::HepMCProduct> >(
         pset.getUntrackedParameter<edm::InputTag>("backgroundLabel", edm::InputTag("mix", "generatorSmeared")));
   }
 
   int cm = 1, va, vb, vc;
   HYJVER(cm, va, vb, vc);
   HepMC::HEPEVT_Wrapper::set_max_number_entries(4000);
 }

◆ ~HydjetHadronizer()

HydjetHadronizer::~HydjetHadronizer ( )

override

Definition at line 126 of file HydjetHadronizer.cc.

References pythia6Service_.

                                     {
   // destructor
   call_pystat(1);
   delete pythia6Service_;
 }

Member Function Documentation

◆ add_heavy_ion_rec()

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

private

Definition at line 136 of file HydjetHadronizer.cc.

References evt, hyfpar, hyjpar, npart, nsub_, nuclear_radius(), and phi0_.

Referenced by generatePartonsAndHadronize().

                                                            {
   // heavy ion record in the final CMSSW Event
   double npart = hyfpar.npart;
   int nproj = static_cast<int>(npart / 2);
   int ntarg = static_cast<int>(npart - nproj);
 
   HepMC::HeavyIon* hi = new HepMC::HeavyIon(nsub_,                           // Ncoll_hard/N of SubEvents
                                             nproj,                           // Npart_proj
                                             ntarg,                           // Npart_targ
                                             static_cast<int>(hyfpar.nbcol),  // Ncoll
                                             0,                               // spectator_neutrons
                                             0,                               // spectator_protons
                                             0,                               // N_Nwounded_collisions
                                             0,                               // Nwounded_N_collisions
                                             0,                               // Nwounded_Nwounded_collisions
                                             hyfpar.bgen * nuclear_radius(),  // impact_parameter in [fm]
                                             phi0_,                           // event_plane_angle
                                             0,  //hypsi3.psi3,                                   // eccentricity
                                             hyjpar.sigin  // sigma_inel_NN
   );
 
   evt->set_heavy_ion(*hi);
   delete hi;
 }

◆ build_hyjet()

HepMC::GenParticle * HydjetHadronizer::build_hyjet	(	int	index,
		int	barcode
	)

private

Definition at line 162 of file HydjetHadronizer.cc.

References convertStatus(), cosphi0_, GenParticle::GenParticle, hyjets, gen::p, sinphi0_, and x.

Referenced by get_particles().

                                                                       {
   // Build particle object corresponding to index in hyjets (soft+hard)
   double x0 = hyjets.phj[0][index];
   double y0 = hyjets.phj[1][index];
 
   double x = x0 * cosphi0_ - y0 * sinphi0_;
   double y = y0 * cosphi0_ + x0 * sinphi0_;
 
   HepMC::GenParticle* p = new HepMC::GenParticle(HepMC::FourVector(x,                      // px
                                                                    y,                      // py
                                                                    hyjets.phj[2][index],   // pz
                                                                    hyjets.phj[3][index]),  // E
                                                  hyjets.khj[1][index],                     // id
                                                  convertStatus(hyjets.khj[0][index]        // status
                                                                ));
 
   p->suggest_barcode(barcode);
   return p;
 }

◆ build_hyjet_vertex()

HepMC::GenVertex * HydjetHadronizer::build_hyjet_vertex	(	int	i,
		int	id
	)

private

Definition at line 183 of file HydjetHadronizer.cc.

References cosphi0_, hyjets, mps_fire::i, sinphi0_, submitPVValidationJobs::t, bphysicsOniaDQM_cfi::vertex, and x.

Referenced by get_particles().

                                                                   {
   // build verteces for the hyjets stored events
   double x0 = hyjets.vhj[0][i];
   double y0 = hyjets.vhj[1][i];
   double x = x0 * cosphi0_ - y0 * sinphi0_;
   double y = y0 * cosphi0_ + x0 * sinphi0_;
   double z = hyjets.vhj[2][i];
   double t = hyjets.vhj[4][i];
 
   HepMC::GenVertex* vertex = new HepMC::GenVertex(HepMC::FourVector(x, y, z, t), id);
   return vertex;
 }

◆ call_hyinit()

bool HydjetHadronizer::call_hyinit	(	double	energy,
		double	a,
		int	ifb,
		double	bmin,
		double	bmax,
		double	bfix,
		int	nh
	)

private

Definition at line 424 of file HydjetHadronizer.cc.

References a, ReggeGribovPartonMC_AdvancedParameters_cfi::bmax, ReggeGribovPartonMC_AdvancedParameters_cfi::bmin, hcalRecHitTable_cff::energy, HYINIT, and cms::cuda::nh.

Referenced by initializeForInternalPartons().

                                                                                                                   {
   // initialize hydjet
 
   pydatr.mrpy[2] = 1;
   HYINIT(energy, a, ifb, bmin, bmax, bfix, nh);
   return true;
 }

◆ classname()

const char * HydjetHadronizer::classname ( ) const

Definition at line 565 of file HydjetHadronizer.cc.

565 { return "gen::HydjetHadronizer"; }

◆ convertStatus()

int HydjetHadronizer::convertStatus ( int st )

private

Definition at line 38 of file HydjetHadronizer.cc.

Referenced by build_hyjet(), and get_particles().

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

◆ decay()

bool HydjetHadronizer::decay ( )

Definition at line 557 of file HydjetHadronizer.cc.

557 { return true; }

◆ declareSpecialSettings()

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

inline

Definition at line 53 of file HydjetHadronizer.h.

53 { return true; }

◆ declareStableParticles()

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

Definition at line 534 of file HydjetHadronizer.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 HydjetHadronizer::doSetRandomEngine ( CLHEP::HepRandomEngine * v )

overrideprivatevirtual

Reimplemented from gen::BaseHadronizer.

Definition at line 133 of file HydjetHadronizer.cc.

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

133 { pythia6Service_->setRandomEngine(v); }

gen::v

double v[5][pyjets_maxn]

Definition: Pythia6Hadronizer.cc:89

gen::HydjetHadronizer::pythia6Service_

Pythia6Service * pythia6Service_

Definition: HydjetHadronizer.h:132

gen::Pythia6Service::setRandomEngine

void setRandomEngine(CLHEP::HepRandomEngine *v)

Definition: Pythia6Service.h:46

◆ doSharedResources()

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

inlineoverrideprivatevirtual

Reimplemented from gen::BaseHadronizer.

Definition at line 61 of file HydjetHadronizer.h.

References theSharedResources.

61 { return theSharedResources; }

gen::HydjetHadronizer::theSharedResources

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

Definition: HydjetHadronizer.h:63

◆ finalizeEvent()

void HydjetHadronizer::finalizeEvent ( )

Definition at line 561 of file HydjetHadronizer.cc.

561 {}

◆ generatePartonsAndHadronize()

bool HydjetHadronizer::generatePartonsAndHadronize ( )

Definition at line 198 of file HydjetHadronizer.cc.

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

                                                    {
   Pythia6Service::InstanceWrapper guard(pythia6Service_);
 
   // generate single event
   if (embedding_) {
     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) {
       bfixed_ = (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;
 
   edm::LogInfo("HYDJETmode") << "##### HYDJET  nhsel = " << hyjpar.nhsel;
   edm::LogInfo("HYDJETfpart") << "##### HYDJET fpart = " << hyflow.fpart;
   edm::LogInfo("HYDJETtf") << "##### HYDJET hadron freez-out temp, Tf = " << hyflow.Tf;
   edm::LogInfo("HYDJETinTemp") << "##### HYDJET: QGP init temperature, T0 =" << pyqpar.T0u;
   edm::LogInfo("HYDJETinTau") << "##### HYDJET: QGP formation time,tau0 =" << pyqpar.tau0u;
 
   int ntry = 0;
   while (nsoft_ == 0 && nhard_ == 0) {
     if (ntry > 100) {
       edm::LogError("HydjetEmptyEvent") << "##### HYDJET: 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 << "HydjetHadronizerProducer: No particles generated after " << ntry << " tries.\n";
       edm::Exception except(edm::errors::EventCorruption, sstr.str());
       throw except;
     } else {
       HYEVNT(bfixed_);
       nsoft_ = hyfpar.nhyd;
       nsub_ = hyjpar.njet;
       nhard_ = hyfpar.npyt;
       ++ntry;
     }
   }
 
   if (hyjpar.nhsel < 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 true;
 }

◆ get_particles()

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

private

Definition at line 315 of file HydjetHadronizer.cc.

References funct::abs(), build_hyjet(), build_hyjet_vertex(), ALPAKA_ACCELERATOR_NAMESPACE::brokenline::constexpr(), convertStatus(), evt, Exception, GenParticle::GenParticle, hyjets, hyjpar, ALPAKA_ACCELERATOR_NAMESPACE::caPixelDoublets::if(), LogDebug, nhard_, nsoft_, and nsub_.

Referenced by generatePartonsAndHadronize().

                                                        {
   // Hard particles. The first nhard_ lines from hyjets array.
   // Pythia/Pyquen sub-events (sub-collisions) for a given event
   // Return T/F if success/failure
   // Create particles from lujet entries, assign them into vertices and
   // put the vertices in the GenEvent, for each SubEvent
   // The SubEvent information is kept by storing indeces of main vertices
   // of subevents as a vector in GenHIEvent.
 
   LogDebug("Hydjet") << " Number of sub events " << nsub_;
   LogDebug("Hydjet") << " Number of hard events " << hyjpar.njet;
   LogDebug("Hydjet") << " Number of hard particles " << nhard_;
   LogDebug("Hydjet") << " Number of soft particles " << nsoft_;
   LogDebug("Hydjet") << " nhard_ + nsoft_ = " << nhard_ + nsoft_ << " hyjets.nhj = " << hyjets.nhj << endl;
 
   int ihy = 0;
   int isub = -1;
   int isub_l = -1;
   int stab = 0;
 
   vector<HepMC::GenParticle*> primary_particle(hyjets.nhj);
   vector<HepMC::GenParticle*> particle(hyjets.nhj);
 
   HepMC::GenVertex* sub_vertices = nullptr;  // just initialization
 
   // contain the last index in for each subevent
   vector<int> index(nsub_);
 
   while (ihy < hyjets.nhj) {
     constexpr int kMaxMultiplicity = 200000;
     isub = std::floor((hyjets.khj[2][ihy] / kMaxMultiplicity));
     int hjoffset = isub * kMaxMultiplicity;
 
     if (isub != isub_l) {
       sub_vertices = new HepMC::GenVertex(HepMC::FourVector(0, 0, 0, 0), isub);
       evt->add_vertex(sub_vertices);
       if (!evt->signal_process_vertex())
         evt->set_signal_process_vertex(sub_vertices);
       if (isub >= static_cast<int>(index.size())) {
         index.resize(isub + 1);
       }
       index[isub] = ihy - 1;
       isub_l = isub;
     }
 
     if (convertStatus(hyjets.khj[0][ihy]) == 1)
       stab++;
     LogDebug("Hydjet_array") << ihy << " MULTin ev.:" << hyjets.nhj << " SubEv.#" << isub << " Part #" << ihy + 1
                              << ", PDG: " << hyjets.khj[1][ihy] << " (st. " << convertStatus(hyjets.khj[0][ihy])
                              << ") mother=" << hyjets.khj[2][ihy] - (isub * kMaxMultiplicity) + index[isub] + 1 << " ("
                              << hyjets.khj[2][ihy] << "), childs ("
                              << hyjets.khj[3][ihy] - (isub * kMaxMultiplicity) + index[isub] + 1 << "-"
                              << hyjets.khj[4][ihy] - (isub * kMaxMultiplicity) + index[isub] + 1 << "), vtx ("
                              << hyjets.vhj[0][ihy] << "," << hyjets.vhj[1][ihy] << "," << hyjets.vhj[2][ihy] << ") "
                              << std::endl;
 
     if (hyjets.khj[2][ihy] == 0) {
       primary_particle[ihy] = build_hyjet(ihy, ihy + 1);
       if (!sub_vertices)
         LogError("Hydjet_array") << "##### HYDJET2: Vertex not initialized!";
       else
         sub_vertices->add_particle_out(primary_particle[ihy]);
       LogDebug("Hydjet_array") << " ---> " << ihy + 1 << std::endl;
     } else {
       particle[ihy] = build_hyjet(ihy, ihy + 1);
       int mid = hyjets.khj[2][ihy] - hjoffset + index[isub];
       if (mid > ihy) {
         throw cms::Exception("BadHydjetParent") << "Parent ID " << mid << " is greater than child id " << ihy
                                                 << " this is not allowed.\n When this happened in the past it was "
                                                    "caused the sub event multiplicity being > "
                                                 << kMaxMultiplicity << " which caused an offset overflow.";
       }
       int mid_t = mid;
       while (((mid + 1) < ihy) && (std::abs(hyjets.khj[1][mid]) < 100) &&
              (hyjets.khj[3][mid + 1] - hjoffset + index[isub] <= ihy))
         mid++;
       if (std::abs(hyjets.khj[1][mid]) < 100)
         mid = mid_t;
 
       HepMC::GenParticle* mother = primary_particle.at(mid);
       HepMC::GenVertex* prods = build_hyjet_vertex(ihy, isub);
 
       if (!mother) {
         mother = particle[mid];
         primary_particle[mid] = mother;
       }
 
       HepMC::GenVertex* prod_vertex = mother->end_vertex();
       if (!prod_vertex) {
         prod_vertex = prods;
         prod_vertex->add_particle_in(mother);
         LogDebug("Hydjet_array") << " <--- " << mid + 1 << std::endl;
         evt->add_vertex(prod_vertex);
         prods = nullptr;
       }
 
       prod_vertex->add_particle_out(particle[ihy]);
       LogDebug("Hydjet_array") << " ---" << mid + 1 << "---> " << ihy + 1 << std::endl;
       delete prods;
     }
     ihy++;
   }
   LogDebug("Hydjet_array") << " MULTin ev.:" << hyjets.nhj << ", last index: " << ihy - 1
                            << ", Sub events: " << isub + 1 << ", stable particles: " << stab << std::endl;
 
   return true;
 }

◆ hadronize()

bool HydjetHadronizer::hadronize ( )

Definition at line 555 of file HydjetHadronizer.cc.

555 { return false; }

◆ hydjet_init()

bool HydjetHadronizer::hydjet_init ( const edm::ParameterSet & pset )

private

Definition at line 433 of file HydjetHadronizer.cc.

References angularspecselector_, docollisionalenloss_, doradiativeenloss_, fracsoftmult_, hadfreeztemp_, hyflow, hyjpar, hymode_, maxlongy_, maxtrany_, nquarkflavor_, pyqpar, qgpt0_, qgptau0_, shadowingswitch_, and signn_.

Referenced by initializeForInternalPartons().

                                                            {
   // set hydjet options
 
   // hydjet running mode mode
   // kHydroOnly --- nhsel=0 jet production off (pure HYDRO event), nhsel=0
   // kHydroJets --- nhsle=1 jet production on, jet quenching off (HYDRO+njet*PYTHIA events)
   // kHydroQJet --- nhsel=2 jet production & jet quenching on (HYDRO+njet*PYQUEN events)
   // kJetsOnly  --- nhsel=3 jet production on, jet quenching off, HYDRO off (njet*PYTHIA events)
   // kQJetsOnly --- nhsel=4 jet production & jet quenching on, HYDRO off (njet*PYQUEN events)
 
   if (hymode_ == "kHydroOnly")
     hyjpar.nhsel = 0;
   else if (hymode_ == "kHydroJets")
     hyjpar.nhsel = 1;
   else if (hymode_ == "kHydroQJets")
     hyjpar.nhsel = 2;
   else if (hymode_ == "kJetsOnly")
     hyjpar.nhsel = 3;
   else if (hymode_ == "kQJetsOnly")
     hyjpar.nhsel = 4;
   else
     hyjpar.nhsel = 2;
 
   // fraction of soft hydro induced multiplicity
   hyflow.fpart = fracsoftmult_;
 
   // hadron freez-out temperature
   hyflow.Tf = hadfreeztemp_;
 
   // maximum longitudinal collective rapidity
   hyflow.ylfl = maxlongy_;
 
   // maximum transverse collective rapidity
   hyflow.ytfl = maxtrany_;
 
   // shadowing on=1, off=0
   hyjpar.ishad = shadowingswitch_;
 
   // set inelastic nucleon-nucleon cross section
   hyjpar.sigin = signn_;
 
   // angular emitted gluon spectrum selection
   pyqpar.ianglu = angularspecselector_;
 
   // number of active quark flavors in qgp
   pyqpar.nfu = nquarkflavor_;
 
   // initial temperature of QGP
   pyqpar.T0u = qgpt0_;
 
   // proper time of QGP formation
   pyqpar.tau0u = qgptau0_;
 
   // type of medium induced partonic energy loss
   if (doradiativeenloss_ && docollisionalenloss_) {
     edm::LogInfo("HydjetEnLoss") << "##### Radiative AND Collisional partonic energy loss ON ####";
     pyqpar.ienglu = 0;
   } else if (doradiativeenloss_) {
     edm::LogInfo("HydjetenLoss") << "##### Only RADIATIVE partonic energy loss ON ####";
     pyqpar.ienglu = 1;
   } else if (docollisionalenloss_) {
     edm::LogInfo("HydjetEnLoss") << "##### Only COLLISIONAL partonic energy loss ON ####";
     pyqpar.ienglu = 2;
   } else {
     edm::LogInfo("HydjetEnLoss") << "##### Radiative AND Collisional partonic energy loss ON ####";
     pyqpar.ienglu = 0;
   }
   return true;
 }

◆ initializeForExternalPartons()

bool gen::HydjetHadronizer::initializeForExternalPartons ( )

◆ initializeForInternalPartons()

bool HydjetHadronizer::initializeForInternalPartons ( )

Definition at line 514 of file HydjetHadronizer.cc.

References abeamtarget_, bfixed_, bmax_, bmin_, call_hyinit(), cflag_, comenergy, hydjet_init(), nmultiplicity_, nuclear_radius(), pset_, and pythia6Service_.

                                                     {
   Pythia6Service::InstanceWrapper guard(pythia6Service_);
   // pythia6Service_->setGeneralParams();
 
   // the input impact parameter (bxx_) is in [fm]; transform in [fm/RA] for hydjet usage
   const float ra = nuclear_radius();
   LogInfo("RAScaling") << "Nuclear radius(RA) =  " << ra;
   bmin_ /= ra;
   bmax_ /= ra;
   bfixed_ /= ra;
 
   // hydjet running options
   hydjet_init(pset_);
   // initialize hydjet
   LogInfo("HYDJETinAction") << "##### Calling HYINIT(" << comenergy << "," << abeamtarget_ << "," << cflag_ << ","
                             << bmin_ << "," << bmax_ << "," << bfixed_ << "," << nmultiplicity_ << ") ####";
   call_hyinit(comenergy, abeamtarget_, cflag_, bmin_, bmax_, bfixed_, nmultiplicity_);
   return true;
 }

◆ nuclear_radius()

double HydjetHadronizer::nuclear_radius ( ) const

inlineprivate

Return the nuclear radius derived from the beam/target atomic mass number.

Definition at line 140 of file HydjetHadronizer.h.

References abeamtarget_, and funct::pow().

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

140 { return 1.15 * pow((double)abeamtarget_, 1. / 3.); }

gen::HydjetHadronizer::abeamtarget_

double abeamtarget_

beam/target atomic mass number

Definition: HydjetHadronizer.h:77

funct::pow

Power< A, B >::type pow(const A &a, const B &b)

Definition: Power.h:29

◆ readSettings()

bool HydjetHadronizer::readSettings ( int )

Definition at line 505 of file HydjetHadronizer.cc.

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

                                        {
   Pythia6Service::InstanceWrapper guard(pythia6Service_);
   pythia6Service_->setGeneralParams();
 
   return true;
 }

◆ residualDecay()

bool HydjetHadronizer::residualDecay ( )

Definition at line 559 of file HydjetHadronizer.cc.

559 { return true; }

◆ rotateEvtPlane()

void HydjetHadronizer::rotateEvtPlane ( )

private

Definition at line 547 of file HydjetHadronizer.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 HydjetHadronizer::statistics ( )

Definition at line 563 of file HydjetHadronizer.cc.

563 {}

Member Data Documentation

◆ abeamtarget_

double gen::HydjetHadronizer::abeamtarget_

private

beam/target atomic mass number

Definition at line 77 of file HydjetHadronizer.h.

Referenced by initializeForInternalPartons(), and nuclear_radius().

◆ angularspecselector_

int gen::HydjetHadronizer::angularspecselector_

private

angular emitted gluon spectrum selection

Definition at line 78 of file HydjetHadronizer.h.

Referenced by hydjet_init().

◆ bfixed_

double gen::HydjetHadronizer::bfixed_

private

fixed impact param (fm); valid only if cflag_=0

Definition at line 79 of file HydjetHadronizer.h.

Referenced by generatePartonsAndHadronize(), and initializeForInternalPartons().

◆ bmax_

double gen::HydjetHadronizer::bmax_

private

max impact param; units of nucl radius

Definition at line 80 of file HydjetHadronizer.h.

Referenced by initializeForInternalPartons().

◆ bmin_

double gen::HydjetHadronizer::bmin_

private

min impact param; units of nucl radius

Definition at line 82 of file HydjetHadronizer.h.

Referenced by initializeForInternalPartons().

◆ cflag_

int gen::HydjetHadronizer::cflag_

private

centrality flag = 0 fixed impact param, <> 0 between bmin and bmax

Definition at line 84 of file HydjetHadronizer.h.

Referenced by HydjetHadronizer(), and initializeForInternalPartons().

◆ comenergy

double gen::HydjetHadronizer::comenergy

private

collision energy

Definition at line 88 of file HydjetHadronizer.h.

Referenced by initializeForInternalPartons().

◆ cosphi0_

double gen::HydjetHadronizer::cosphi0_

private

Definition at line 122 of file HydjetHadronizer.h.

Referenced by build_hyjet(), build_hyjet_vertex(), generatePartonsAndHadronize(), and rotateEvtPlane().

◆ docollisionalenloss_

bool gen::HydjetHadronizer::docollisionalenloss_

private

DEFAULT = true.

Definition at line 90 of file HydjetHadronizer.h.

Referenced by hydjet_init().

◆ doradiativeenloss_

bool gen::HydjetHadronizer::doradiativeenloss_

private

DEFAULT = true.

Definition at line 89 of file HydjetHadronizer.h.

Referenced by hydjet_init().

◆ embedding_

bool gen::HydjetHadronizer::embedding_

private

Switch for embedding mode.

Definition at line 87 of file HydjetHadronizer.h.

Referenced by generatePartonsAndHadronize(), and HydjetHadronizer().

◆ evt

HepMC::GenEvent* gen::HydjetHadronizer::evt

private

Definition at line 75 of file HydjetHadronizer.h.

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

◆ fracsoftmult_

double gen::HydjetHadronizer::fracsoftmult_

private

fraction of soft hydro induced hadronic multiplicity proportional to no of nucleon participants (1-fracsoftmult_)— fraction of soft multiplicity proportional to the numebr of nucleon-nucleon binary collisions DEFAULT=1., allowed range [0.01,1]

Definition at line 91 of file HydjetHadronizer.h.

Referenced by hydjet_init().

◆ fVertex_

HepMC::FourVector* gen::HydjetHadronizer::fVertex_

private

Event signal vertex.

Definition at line 129 of file HydjetHadronizer.h.

Referenced by generatePartonsAndHadronize(), and HydjetHadronizer().

◆ hadfreeztemp_

double gen::HydjetHadronizer::hadfreeztemp_

private

hadron freez-out temperature DEFAULT=0.14MeV, allowed ranges [0.08,0.2]MeV

Definition at line 97 of file HydjetHadronizer.h.

Referenced by hydjet_init().

◆ hymode_

std::string gen::HydjetHadronizer::hymode_

private

Hydjet running mode.

Definition at line 99 of file HydjetHadronizer.h.

Referenced by hydjet_init().

◆ maxEventsToPrint_

unsigned int gen::HydjetHadronizer::maxEventsToPrint_

private

Events to print if verbosity.

Definition at line 100 of file HydjetHadronizer.h.

Referenced by HydjetHadronizer().

◆ maxlongy_

double gen::HydjetHadronizer::maxlongy_

private

max longitudinal collective rapidity: controls width of eta-spectra DEFAULT=4, allowed range [0.01,7.0]

Definition at line 101 of file HydjetHadronizer.h.

Referenced by hydjet_init().

◆ maxtrany_

double gen::HydjetHadronizer::maxtrany_

private

max transverse collective rapidity: controls slope of low-pt spectra DEFAULT=1.5, allowed range [0.01,3.0]

Definition at line 104 of file HydjetHadronizer.h.

Referenced by hydjet_init().

◆ nhard_

int gen::HydjetHadronizer::nhard_

private

multiplicity of PYTHIA(+PYQUEN)-induced particles in event

Definition at line 108 of file HydjetHadronizer.h.

Referenced by generatePartonsAndHadronize(), and get_particles().

◆ nmultiplicity_

int gen::HydjetHadronizer::nmultiplicity_

private

mean soft multiplicity in central PbPb automatically calculated for other centralitie and beams

Definition at line 109 of file HydjetHadronizer.h.

Referenced by initializeForInternalPartons().

◆ nquarkflavor_

unsigned int gen::HydjetHadronizer::nquarkflavor_

private

number of active quark flavors in qgp DEFAULT=0; allowed values: 0,1,2,3.

Definition at line 112 of file HydjetHadronizer.h.

Referenced by hydjet_init().

◆ nsoft_

int gen::HydjetHadronizer::nsoft_

private

multiplicity of HYDRO-induced particles in event

Definition at line 111 of file HydjetHadronizer.h.

Referenced by generatePartonsAndHadronize(), and get_particles().

◆ nsub_

int gen::HydjetHadronizer::nsub_

private

number of sub-events

Definition at line 107 of file HydjetHadronizer.h.

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

◆ phi0_

double gen::HydjetHadronizer::phi0_

private

Event plane angle.

Definition at line 120 of file HydjetHadronizer.h.

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

◆ pset_

edm::ParameterSet gen::HydjetHadronizer::pset_

private

Definition at line 76 of file HydjetHadronizer.h.

Referenced by initializeForInternalPartons().

◆ pythia6Service_

Pythia6Service* gen::HydjetHadronizer::pythia6Service_

private

Definition at line 132 of file HydjetHadronizer.h.

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

◆ pythiaPylistVerbosity_

unsigned int gen::HydjetHadronizer::pythiaPylistVerbosity_

private

pythia verbosity; def=1

Definition at line 114 of file HydjetHadronizer.h.

Referenced by HydjetHadronizer().

◆ qgpt0_

double gen::HydjetHadronizer::qgpt0_

private

initial temperature of QGP DEFAULT = 1GeV; allowed range [0.2,2.0]GeV;

Definition at line 115 of file HydjetHadronizer.h.

Referenced by hydjet_init().

◆ qgptau0_

double gen::HydjetHadronizer::qgptau0_

private

proper time of QGP formation DEFAULT = 0.1 fm/c; allowed range [0.01,10.0]fm/

Definition at line 117 of file HydjetHadronizer.h.

Referenced by hydjet_init().

◆ rotate_

bool gen::HydjetHadronizer::rotate_

private

Switch to rotate event plane.

Definition at line 123 of file HydjetHadronizer.h.

Referenced by generatePartonsAndHadronize().

◆ shadowingswitch_

unsigned int gen::HydjetHadronizer::shadowingswitch_

private

shadowing switcher 1-ON, 0-OFF

Definition at line 125 of file HydjetHadronizer.h.

Referenced by hydjet_init().

◆ signalVtx_

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

private

Pset double vector to set event signal vertex.

Definition at line 130 of file HydjetHadronizer.h.

Referenced by HydjetHadronizer().

◆ signn_

double gen::HydjetHadronizer::signn_

private

inelastic nucleon nucleon cross section [mb] DEFAULT= 58 mb

Definition at line 127 of file HydjetHadronizer.h.

Referenced by hydjet_init().

◆ sinphi0_

double gen::HydjetHadronizer::sinphi0_

private

Definition at line 121 of file HydjetHadronizer.h.

Referenced by build_hyjet(), build_hyjet_vertex(), generatePartonsAndHadronize(), and rotateEvtPlane().

◆ src_

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

private

Definition at line 133 of file HydjetHadronizer.h.

Referenced by generatePartonsAndHadronize(), and HydjetHadronizer().

◆ theSharedResources

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

staticprivate

Initial value:

= {edm::SharedResourceNames::kPythia6,

gen::FortranInstance::kFortranInstance}

Definition at line 63 of file HydjetHadronizer.h.

Referenced by doSharedResources().

Public Member Functions

Private Member Functions

Private Attributes

Static Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ HydjetHadronizer()

◆ ~HydjetHadronizer()

Member Function Documentation

◆ add_heavy_ion_rec()

◆ build_hyjet()

◆ build_hyjet_vertex()

◆ call_hyinit()

◆ classname()

◆ convertStatus()

◆ decay()

◆ declareSpecialSettings()

◆ declareStableParticles()

◆ doSetRandomEngine()

◆ doSharedResources()

◆ finalizeEvent()

◆ generatePartonsAndHadronize()

◆ get_particles()

◆ hadronize()

◆ hydjet_init()

◆ initializeForExternalPartons()

◆ initializeForInternalPartons()

◆ nuclear_radius()

◆ readSettings()

◆ residualDecay()

◆ rotateEvtPlane()

◆ statistics()

Member Data Documentation

◆ abeamtarget_

◆ angularspecselector_

◆ bfixed_

◆ bmax_

◆ bmin_

◆ cflag_

◆ comenergy

◆ cosphi0_

◆ docollisionalenloss_

◆ doradiativeenloss_

◆ embedding_

◆ evt

◆ fracsoftmult_

◆ fVertex_

◆ hadfreeztemp_

◆ hymode_

◆ maxEventsToPrint_

◆ maxlongy_

◆ maxtrany_

◆ nhard_

◆ nmultiplicity_

◆ nquarkflavor_

◆ nsoft_

◆ nsub_

◆ phi0_

◆ pset_

◆ pythia6Service_

◆ pythiaPylistVerbosity_

◆ qgpt0_

◆ qgptau0_

◆ rotate_

◆ shadowingswitch_

◆ signalVtx_

◆ signn_

◆ sinphi0_

◆ src_

◆ theSharedResources