#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 &)

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

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_

int	angularspecselector_

double	bfixed_

double	bmax_

double	bmin_

int	cflag_

double	comenergy

double	cosphi0_

bool	docollisionalenloss_
	DEFAULT = true. More...

bool	doradiativeenloss_

bool	embedding_

HepMC::GenEvent *	evt

double	fracsoftmult_
	DEFAULT = true. More...

double	hadfreeztemp_

std::string	hymode_

unsigned int	maxEventsToPrint_

double	maxlongy_

double	maxtrany_

int	nhard_

int	nmultiplicity_

unsigned int	nquarkflavor_

int	nsoft_

int	nsub_

double	phi0_

edm::ParameterSet	pset_

Pythia6Service *	pythia6Service_

unsigned int	pythiaPylistVerbosity_

double	qgpt0_

double	qgptau0_

bool	rotate_

unsigned int	shadowingswitch_

double	signn_

double	sinphi0_

edm::InputTag	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< 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 35 of file HydjetHadronizer.h.

Constructor & Destructor Documentation

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

Definition at line 64 of file HydjetHadronizer.cc.

References embedding_, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), LogDebug, maxEventsToPrint_, pythiaPylistVerbosity_, 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<bool>("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")),
       pythia6Service_(new Pythia6Service(pset)) {
   // Default constructor
 
   // 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_)
     src_ = pset.getParameter<edm::InputTag>("backgroundLabel");
 }

HydjetHadronizer::~HydjetHadronizer ( )

override

Definition at line 115 of file HydjetHadronizer.cc.

References pythia6Service_.

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

Member Function Documentation

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

private

Definition at line 125 of file HydjetHadronizer.cc.

References 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;
 }

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

private

Definition at line 151 of file HydjetHadronizer.cc.

References cosphi0_, GenParticle::GenParticle, hyjets, gen::p, and sinphi0_.

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;
 }

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

private

Definition at line 173 of file HydjetHadronizer.cc.

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

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;
 }

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

private

Definition at line 343 of file HydjetHadronizer.cc.

References HYINIT.

Referenced by initializeForInternalPartons().

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

const char * HydjetHadronizer::classname ( ) const

Definition at line 484 of file HydjetHadronizer.cc.

484 { return "gen::HydjetHadronizer"; }

bool HydjetHadronizer::decay ( )

Definition at line 476 of file HydjetHadronizer.cc.

476 { return true; }

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

inline

Definition at line 48 of file HydjetHadronizer.h.

References data-class-funcs::classname, myMessageLogger_cff::statistics, and findQualityFiles::v.

48 { return true; }

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

Definition at line 453 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;
 }

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

overrideprivatevirtual

Reimplemented from gen::BaseHadronizer.

Definition at line 122 of file HydjetHadronizer.cc.

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

122 { pythia6Service_->setRandomEngine(v); }

gen::v

double v[5][pyjets_maxn]

Definition: Cascade2Hadronizer.cc:74

gen::HydjetHadronizer::pythia6Service_

Pythia6Service * pythia6Service_

Definition: HydjetHadronizer.h:124

gen::Pythia6Service::setRandomEngine

void setRandomEngine(CLHEP::HepRandomEngine *v)

Definition: Pythia6Service.h:46

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

inlineoverrideprivatevirtual

Reimplemented from gen::BaseHadronizer.

Definition at line 56 of file HydjetHadronizer.h.

56 { return theSharedResources; }

gen::HydjetHadronizer::theSharedResources

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

Definition: HydjetHadronizer.h:58

void HydjetHadronizer::finalizeEvent ( )

Definition at line 480 of file HydjetHadronizer.cc.

480 {}

bool HydjetHadronizer::generatePartonsAndHadronize ( )

Definition at line 189 of file HydjetHadronizer.cc.

References add_heavy_ion_rec(), bfixed_, cflag_, funct::cos(), cosphi0_, MillePedeFileConverter_cfg::e, embedding_, gen::BaseHadronizer::event(), edm::errors::EventCorruption, evt, get_particles(), edm::Event::getByLabel(), gen::BaseHadronizer::getEDMEvent(), edm::HepMCProduct::GetEvent(), HYEVNT, hyflow, hyfpar, hyjpar, input, nhard_, nsoft_, nsub_, phi0_, pypars, pyqpar, pythia6Service_, rotate_, rotateEvtPlane(), funct::sin(), sinphi0_, and src_.

                                                    {
   Pythia6Service::InstanceWrapper guard(pythia6Service_);
 
   // generate single event
   if (embedding_) {
     cflag_ = 0;
     const edm::Event& e = getEDMEvent();
     Handle<HepMCProduct> input;
     e.getByLabel(src_, input);
     const HepMC::GenEvent* inev = input->GetEvent();
     const HepMC::HeavyIon* hi = inev->heavy_ion();
     if (hi) {
       bfixed_ = hi->impact_parameter();
       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;
 
   // generate a HYDJET event
   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();
       nsoft_ = hyfpar.nhyd;
       nsub_ = hyjpar.njet;
       nhard_ = hyfpar.npyt;
       ++ntry;
     }
   }
 
   if (hyjpar.nhsel < 3)
     nsub_++;
 
   // event information
   HepMC::GenEvent* evt = new HepMC::GenEvent();
 
   if (nhard_ > 0 || nsoft_ > 0)
     get_particles(evt);
 
   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);
 
   event().reset(evt);
   return true;
 }

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

private

Definition at line 260 of file HydjetHadronizer.cc.

References build_hyjet(), build_hyjet_vertex(), GenParticle::GenParticle, hyjets, hyjpar, mps_fire::i, LogDebug, nhard_, nsoft_, nsub_, and ecalTrigSettings_cff::particles.

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("SubEvent") << "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_;
 
   vector<HepMC::GenVertex*> sub_vertices(nsub_);
 
   int ihy = 0;
   for (int isub = 0; isub < nsub_; isub++) {
     LogDebug("SubEvent") << "Sub Event ID : " << isub;
 
     int sub_up = (isub + 1) * 50000;  // Upper limit in mother index, determining the range of Sub-Event
     vector<HepMC::GenParticle*> particles;
     vector<int> mother_ids;
     vector<HepMC::GenVertex*> prods;
 
     sub_vertices[isub] = new HepMC::GenVertex(HepMC::FourVector(0, 0, 0, 0), isub);
     evt->add_vertex(sub_vertices[isub]);
     if (!evt->signal_process_vertex())
       evt->set_signal_process_vertex(sub_vertices[isub]);
 
     while (ihy < nhard_ + nsoft_ && (hyjets.khj[2][ihy] < sub_up || ihy > nhard_)) {
       particles.push_back(build_hyjet(ihy, ihy + 1));
       prods.push_back(build_hyjet_vertex(ihy, isub));
       mother_ids.push_back(hyjets.khj[2][ihy]);
       LogDebug("DecayChain") << "Mother index : " << hyjets.khj[2][ihy];
 
       ihy++;
     }
 
     //Produce Vertices and add them to the GenEvent. Remember that GenParticles are adopted by
     //GenVertex and GenVertex is adopted by GenEvent.
 
     LogDebug("Hydjet") << "Number of particles in vector " << particles.size();
 
     for (unsigned int i = 0; i < particles.size(); i++) {
       HepMC::GenParticle* part = particles[i];
 
       //The Fortran code is modified to preserve mother id info, by seperating the beginning
       //mother indices of successive subevents by 5000
       int mid = mother_ids[i] - isub * 50000 - 1;
       LogDebug("DecayChain") << "Particle " << i;
       LogDebug("DecayChain") << "Mother's ID " << mid;
       LogDebug("DecayChain") << "Particle's PDG ID " << part->pdg_id();
 
       if (mid <= 0) {
         sub_vertices[isub]->add_particle_out(part);
         continue;
       }
 
       if (mid > 0) {
         HepMC::GenParticle* mother = particles[mid];
         LogDebug("DecayChain") << "Mother's PDG ID " << mother->pdg_id();
 
         HepMC::GenVertex* prod_vertex = mother->end_vertex();
         if (!prod_vertex) {
           prod_vertex = prods[i];
           prod_vertex->add_particle_in(mother);
           evt->add_vertex(prod_vertex);
           prods[i] = nullptr;  // mark to protect deletion
         }
         prod_vertex->add_particle_out(part);
       }
     }
     // cleanup vertices not assigned to evt
     for (unsigned int i = 0; i < prods.size(); i++) {
       if (prods[i])
         delete prods[i];
     }
   }
   return true;
 }

bool HydjetHadronizer::hadronize ( )

Definition at line 474 of file HydjetHadronizer.cc.

474 { return false; }

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

private

Definition at line 352 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;
 }

bool gen::HydjetHadronizer::initializeForExternalPartons ( )

bool HydjetHadronizer::initializeForInternalPartons ( )

Definition at line 433 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;
 }

double HydjetHadronizer::nuclear_radius ( ) const

inlineprivate

Definition at line 128 of file HydjetHadronizer.h.

References funct::pow().

Referenced by add_heavy_ion_rec(), and initializeForInternalPartons().

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

bool HydjetHadronizer::readSettings ( int )

Definition at line 424 of file HydjetHadronizer.cc.

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

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

bool HydjetHadronizer::residualDecay ( )

Definition at line 478 of file HydjetHadronizer.cc.

478 { return true; }

void HydjetHadronizer::rotateEvtPlane ( )

private

Definition at line 466 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_);
 }

void HydjetHadronizer::statistics ( )

Definition at line 482 of file HydjetHadronizer.cc.

482 {}

Member Data Documentation

double gen::HydjetHadronizer::abeamtarget_

private

Definition at line 71 of file HydjetHadronizer.h.

Referenced by initializeForInternalPartons().

int gen::HydjetHadronizer::angularspecselector_

private

Definition at line 72 of file HydjetHadronizer.h.

Referenced by hydjet_init().

double gen::HydjetHadronizer::bfixed_

private

Definition at line 73 of file HydjetHadronizer.h.

Referenced by generatePartonsAndHadronize(), and initializeForInternalPartons().

double gen::HydjetHadronizer::bmax_

private

Definition at line 74 of file HydjetHadronizer.h.

Referenced by initializeForInternalPartons().

double gen::HydjetHadronizer::bmin_

private

Definition at line 76 of file HydjetHadronizer.h.

Referenced by initializeForInternalPartons().

int gen::HydjetHadronizer::cflag_

private

Definition at line 78 of file HydjetHadronizer.h.

Referenced by generatePartonsAndHadronize(), and initializeForInternalPartons().

double gen::HydjetHadronizer::comenergy

private

Definition at line 82 of file HydjetHadronizer.h.

Referenced by initializeForInternalPartons().

double gen::HydjetHadronizer::cosphi0_

private

Definition at line 116 of file HydjetHadronizer.h.

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

bool gen::HydjetHadronizer::docollisionalenloss_

private

DEFAULT = true.

Definition at line 84 of file HydjetHadronizer.h.

Referenced by hydjet_init().

bool gen::HydjetHadronizer::doradiativeenloss_

private

Definition at line 83 of file HydjetHadronizer.h.

Referenced by hydjet_init().

bool gen::HydjetHadronizer::embedding_

private

Definition at line 81 of file HydjetHadronizer.h.

Referenced by generatePartonsAndHadronize(), and HydjetHadronizer().

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

private

Definition at line 69 of file HydjetHadronizer.h.

Referenced by generatePartonsAndHadronize().

double gen::HydjetHadronizer::fracsoftmult_

private

DEFAULT = true.

Definition at line 85 of file HydjetHadronizer.h.

Referenced by hydjet_init().

double gen::HydjetHadronizer::hadfreeztemp_

private

Definition at line 91 of file HydjetHadronizer.h.

Referenced by hydjet_init().

std::string gen::HydjetHadronizer::hymode_

private

Definition at line 93 of file HydjetHadronizer.h.

Referenced by hydjet_init().

unsigned int gen::HydjetHadronizer::maxEventsToPrint_

private

Definition at line 94 of file HydjetHadronizer.h.

Referenced by HydjetHadronizer().

double gen::HydjetHadronizer::maxlongy_

private

Definition at line 95 of file HydjetHadronizer.h.

Referenced by hydjet_init().

double gen::HydjetHadronizer::maxtrany_

private

Definition at line 98 of file HydjetHadronizer.h.

Referenced by hydjet_init().

int gen::HydjetHadronizer::nhard_

private

Definition at line 102 of file HydjetHadronizer.h.

Referenced by generatePartonsAndHadronize(), and get_particles().

int gen::HydjetHadronizer::nmultiplicity_

private

Definition at line 103 of file HydjetHadronizer.h.

Referenced by initializeForInternalPartons().

unsigned int gen::HydjetHadronizer::nquarkflavor_

private

Definition at line 106 of file HydjetHadronizer.h.

Referenced by hydjet_init().

int gen::HydjetHadronizer::nsoft_

private

Definition at line 105 of file HydjetHadronizer.h.

Referenced by generatePartonsAndHadronize(), and get_particles().

int gen::HydjetHadronizer::nsub_

private

Definition at line 101 of file HydjetHadronizer.h.

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

double gen::HydjetHadronizer::phi0_

private

Definition at line 114 of file HydjetHadronizer.h.

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

edm::ParameterSet gen::HydjetHadronizer::pset_

private

Definition at line 70 of file HydjetHadronizer.h.

Referenced by initializeForInternalPartons().

Pythia6Service* gen::HydjetHadronizer::pythia6Service_

private

Definition at line 124 of file HydjetHadronizer.h.

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

unsigned int gen::HydjetHadronizer::pythiaPylistVerbosity_

private

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

Definition at line 108 of file HydjetHadronizer.h.

Referenced by HydjetHadronizer().

double gen::HydjetHadronizer::qgpt0_

private

Definition at line 109 of file HydjetHadronizer.h.

Referenced by hydjet_init().

double gen::HydjetHadronizer::qgptau0_

private

Definition at line 111 of file HydjetHadronizer.h.

Referenced by hydjet_init().

bool gen::HydjetHadronizer::rotate_

private

Definition at line 117 of file HydjetHadronizer.h.

Referenced by generatePartonsAndHadronize().

unsigned int gen::HydjetHadronizer::shadowingswitch_

private

Definition at line 119 of file HydjetHadronizer.h.

Referenced by hydjet_init().

double gen::HydjetHadronizer::signn_

private

Definition at line 121 of file HydjetHadronizer.h.

Referenced by hydjet_init().

double gen::HydjetHadronizer::sinphi0_

private

Definition at line 115 of file HydjetHadronizer.h.

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

edm::InputTag gen::HydjetHadronizer::src_

private

Definition at line 125 of file HydjetHadronizer.h.

Referenced by generatePartonsAndHadronize(), and HydjetHadronizer().

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

staticprivate

Initial value:

= {edm::SharedResourceNames::kPythia6,

gen::FortranInstance::kFortranInstance}

Definition at line 58 of file HydjetHadronizer.h.

Public Member Functions

Private Member Functions

Private Attributes

Static Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation