#include <HijingHadronizer.h>

Inheritance diagram for gen::HijingHadronizer:

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

	HijingHadronizer (const edm::ParameterSet &)

bool	initializeForExternalPartons ()

bool	initializeForInternalPartons ()

bool	readSettings (int)

bool	residualDecay ()

void	statistics ()

	~HijingHadronizer () override

Public Member Functions inherited from gen::BaseHadronizer
	BaseHadronizer (edm::ParameterSet const &ps)

void	cleanLHE ()

void	generateLHE (edm::LuminosityBlock const &lumi, CLHEP::HepRandomEngine *rengine, unsigned int ncpu)

edm::Event &	getEDMEvent () const

std::unique_ptr< HepMC::GenEvent >	getGenEvent ()

std::unique_ptr< GenEventInfoProduct >	getGenEventInfo ()

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

GenRunInfoProduct &	getGenRunInfo ()

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

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

const std::string &	gridpackPath () const

int	randomIndex () const

const std::string &	randomInitConfigDescription () const

void	randomizeIndex (edm::LuminosityBlock const &lumi, CLHEP::HepRandomEngine *rengine)

void	resetEvent (std::unique_ptr< HepMC::GenEvent > event)

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

virtual bool	select (HepMC::GenEvent *) const

void	setEDMEvent (edm::Event &event)

void	setLHEEvent (std::unique_ptr< lhef::LHEEvent > event)

void	setLHERunInfo (std::unique_ptr< lhef::LHERunInfo > runInfo)

void	setRandomEngine (CLHEP::HepRandomEngine *v)

std::vector< std::string > const &	sharedResources () const

virtual	~BaseHadronizer () noexcept(false)

Private Member Functions
void	add_heavy_ion_rec (HepMC::GenEvent *evt)

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

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

bool	call_hijset (double efrm, std::string frame, std::string proj, std::string targ, int iap, int izp, int iat, int izt)

void	doSetRandomEngine (CLHEP::HepRandomEngine *v) override

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

bool	get_particles (HepMC::GenEvent *evt)

void	rotateEvtPlane ()

Private Attributes
double	bmax_

double	bmin_

double	cosphi0_

double	efrm_

HepMC::GenEvent *	evt

std::string	frame_

int	iap_

int	iat_

int	izp_

int	izt_

double	phi0_

std::string	proj_

edm::ParameterSet	pset_

bool	rotate_

double	sinphi0_

std::string	targ_

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< 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 27 of file HijingHadronizer.h.

Constructor & Destructor Documentation

◆ HijingHadronizer()

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

Definition at line 54 of file HijingHadronizer.cc.

     : BaseHadronizer(pset),
       evt(nullptr),
       pset_(pset),
       bmax_(pset.getParameter<double>("bMax")),
       bmin_(pset.getParameter<double>("bMin")),
       efrm_(pset.getParameter<double>("comEnergy")),
       frame_(pset.getParameter<string>("frame")),
       proj_(pset.getParameter<string>("proj")),
       targ_(pset.getParameter<string>("targ")),
       iap_(pset.getParameter<int>("iap")),
       izp_(pset.getParameter<int>("izp")),
       iat_(pset.getParameter<int>("iat")),
       izt_(pset.getParameter<int>("izt")),
       phi0_(0.),
       sinphi0_(0.),
       cosphi0_(1.),
       rotate_(pset.getParameter<bool>("rotateEventPlane")) {
   // Default constructor
 }

◆ ~HijingHadronizer()

HijingHadronizer::~HijingHadronizer ( )

override

Definition at line 76 of file HijingHadronizer.cc.

                                     {
   // destructor
 }

Member Function Documentation

◆ add_heavy_ion_rec()

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

private

Definition at line 84 of file HijingHadronizer.cc.

References evt, himain1, hiparnt, and phi0_.

Referenced by generatePartonsAndHadronize().

                                                            {
   // heavy ion record in the final CMSSW Event
   HepMC::HeavyIon* hi = new HepMC::HeavyIon(himain1.jatt,  // Ncoll_hard/N of SubEvents
                                             himain1.np,    // Npart_proj
                                             himain1.nt,    // Npart_targ
                                             himain1.n0 + himain1.n01 + himain1.n10 + himain1.n11,  // Ncoll
                                             0,                                                     // spectator_neutrons
                                             0,                                                     // spectator_protons
                                             himain1.n01,  // N_Nwounded_collisions
                                             himain1.n10,  // Nwounded_N_collisions
                                             himain1.n11,  // Nwounded_Nwounded_collisions
                                             //gsfs Changed from 19 to 18 (Fortran counts from 1 , not 0)
                                             hiparnt.hint1[18],  // impact_parameter in [fm]
                                             phi0_,              // event_plane_angle
                                             0,                  // eccentricity
                                             //gsfs Changed from 12 to 11 (Fortran counts from 1 , not 0)
                                             hiparnt.hint1[11]  // sigma_inel_NN
   );
   evt->set_heavy_ion(*hi);
   delete hi;
 }

◆ build_hijing()

HepMC::GenParticle * HijingHadronizer::build_hijing	(	int	index,
		int	barcode
	)

private

Definition at line 107 of file HijingHadronizer.cc.

References cosphi0_, GenParticle::GenParticle, himain2, gen::p, sinphi0_, and x.

Referenced by get_particles().

                                                                        {
   // Build particle object corresponding to index in hijing
 
   double x0 = himain2.patt[0][index];
   double y0 = himain2.patt[1][index];
 
   double x = x0 * cosphi0_ - y0 * sinphi0_;
   double y = y0 * cosphi0_ + x0 * sinphi0_;
 
   // Hijing gives V0's status=4, they need to have status=1 to be decayed in geant
   // also change status=11 to status=2
   if (himain2.katt[3][index] <= 10 && himain2.katt[3][index] > 0)
     himain2.katt[3][index] = 1;
   if (himain2.katt[3][index] <= 20 && himain2.katt[3][index] > 10)
     himain2.katt[3][index] = 2;
 
   HepMC::GenParticle* p = new HepMC::GenParticle(HepMC::FourVector(x,                        // px
                                                                    y,                        // py
                                                                    himain2.patt[2][index],   // pz
                                                                    himain2.patt[3][index]),  // E
                                                  himain2.katt[0][index],                     // id
                                                  himain2.katt[3][index]                      // status
   );
   p->suggest_barcode(barcode);
 
   return p;
 }

◆ build_hijing_vertex()

HepMC::GenVertex * HijingHadronizer::build_hijing_vertex	(	int	i,
		int	id
	)

private

Definition at line 136 of file HijingHadronizer.cc.

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

Referenced by get_particles().

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

◆ call_hijset()

bool HijingHadronizer::call_hijset	(	double	efrm,
		std::string	frame,
		std::string	proj,
		std::string	targ,
		int	iap,
		int	izp,
		int	iat,
		int	izt
	)

private

Definition at line 242 of file HijingHadronizer.cc.

References amptDefault_cfi::frame, HIJSET, amptDefault_cfi::iap, amptDefault_cfi::iat, amptDefault_cfi::izp, amptDefault_cfi::izt, amptDefault_cfi::proj, and amptDefault_cfi::targ.

Referenced by initializeForInternalPartons().

                                                                                                       {
   float ef = efrm;
   // initialize hydjet
   frame.resize(4, ' ');
   proj.resize(4, ' ');
   targ.resize(4, ' ');
   HIJSET(ef, frame.data(), proj.data(), targ.data(), iap, izp, iat, izt, frame.length(), proj.length(), targ.length());
   return true;
 }

◆ classname()

const char * HijingHadronizer::classname ( ) const

Definition at line 288 of file HijingHadronizer.cc.

288 { return "gen::HijingHadronizer"; }

◆ decay()

bool HijingHadronizer::decay ( )

Definition at line 280 of file HijingHadronizer.cc.

280 { return true; }

◆ declareSpecialSettings()

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

inline

Definition at line 40 of file HijingHadronizer.h.

40 { return true; }

◆ declareStableParticles()

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

Definition at line 268 of file HijingHadronizer.cc.

268 { return true; }

◆ doSetRandomEngine()

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

overrideprivatevirtual

Reimplemented from gen::BaseHadronizer.

Definition at line 81 of file HijingHadronizer.cc.

References hijRandomEngine, and gen::v.

81 { hijRandomEngine = v; }

gen::v

double v[5][pyjets_maxn]

Definition: Pythia6Hadronizer.cc:89

hijRandomEngine

static CLHEP::HepRandomEngine * hijRandomEngine

Definition: HijingHadronizer.cc:30

◆ doSharedResources()

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

inlineoverrideprivatevirtual

Reimplemented from gen::BaseHadronizer.

Definition at line 48 of file HijingHadronizer.h.

References theSharedResources.

48 { return theSharedResources; }

gen::HijingHadronizer::theSharedResources

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

Definition: HijingHadronizer.h:82

◆ finalizeEvent()

void HijingHadronizer::finalizeEvent ( )

Definition at line 284 of file HijingHadronizer.cc.

284 { return; }

◆ generatePartonsAndHadronize()

bool HijingHadronizer::generatePartonsAndHadronize ( )

Definition at line 149 of file HijingHadronizer.cc.

References add_heavy_ion_rec(), bmax_, bmin_, gen::BaseHadronizer::event(), evt, frame_, get_particles(), HIJING, rotate_, and rotateEvtPlane().

                                                    {
   // generate single event
   if (rotate_)
     rotateEvtPlane();
 
   // generate a HIJING event
 
   float f_bmin = bmin_;
   float f_bmax = bmax_;
   HIJING(frame_.c_str(), f_bmin, f_bmax, frame_.length());
 
   // event information
   HepMC::GenEvent* evt = new HepMC::GenEvent();
   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;
 }

◆ get_particles()

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

private

Definition at line 174 of file HijingHadronizer.cc.

References build_hijing(), build_hijing_vertex(), evt, GenParticle::GenParticle, himain1, himain2, mps_fire::i, LogDebug, ecalTrigSettings_cff::particles, and mathSSE::sqrt().

Referenced by generatePartonsAndHadronize().

                                                        {
   HepMC::GenVertex* vertice;
 
   vector<HepMC::GenParticle*> particles;
   vector<int> mother_ids;
   vector<HepMC::GenVertex*> prods;
 
   vertice = new HepMC::GenVertex(HepMC::FourVector(0, 0, 0, 0), 0);
   evt->add_vertex(vertice);
   if (!evt->signal_process_vertex())
     evt->set_signal_process_vertex(vertice);
 
   const unsigned int knumpart = himain1.natt;
 
   for (unsigned int ipart = 0; ipart < knumpart; ipart++) {
     int mid = himain2.katt[2][ipart] - 1;  // careful of fortan to c++ array index
 
     particles.push_back(build_hijing(ipart, ipart + 1));
     prods.push_back(build_hijing_vertex(ipart, 0));
     mother_ids.push_back(mid);
     LogDebug("DecayChain") << "Mother index : " << mid;
   }
 
   LogDebug("Hijing") << "Number of particles in vector " << particles.size();
 
   for (unsigned int ipart = 0; ipart < particles.size(); ipart++) {
     HepMC::GenParticle* part = particles[ipart];
 
     int mid = mother_ids[ipart];
     LogDebug("DecayChain") << "Particle " << ipart;
     LogDebug("DecayChain") << "Mother's ID " << mid;
     LogDebug("DecayChain") << "Particle's PDG ID " << part->pdg_id();
 
     // remove zero pT particles from list, protection for fastJet against pt=0 jets
     if (part->status() == 1 &&
         sqrt(part->momentum().px() * part->momentum().px() + part->momentum().py() * part->momentum().py()) == 0)
       continue;
 
     if (mid <= 0) {
       vertice->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[ipart];
         prod_vertex->add_particle_in(mother);
 
         evt->add_vertex(prod_vertex);
         prods[ipart] = 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;
 }

◆ hadronize()

bool HijingHadronizer::hadronize ( )

Definition at line 278 of file HijingHadronizer.cc.

278 { return false; }

◆ initializeForExternalPartons()

bool gen::HijingHadronizer::initializeForExternalPartons ( )

◆ initializeForInternalPartons()

bool HijingHadronizer::initializeForInternalPartons ( )

Definition at line 254 of file HijingHadronizer.cc.

References call_hijset(), efrm_, frame_, iap_, iat_, izp_, izt_, proj_, and targ_.

                                                     {
   //initialize pythia5
 
   // std::string dumstr = "";
   // call_pygive(dumstr);
 
   // initialize hijing
   LogInfo("HIJINGinAction") << "##### Calling HIJSET(" << efrm_ << "," << frame_ << "," << proj_ << "," << targ_ << ","
                             << iap_ << "," << izp_ << "," << iat_ << "," << izt_ << ") ####";
   call_hijset(efrm_, frame_, proj_, targ_, iap_, izp_, iat_, izt_);
 
   return true;
 }

◆ readSettings()

bool gen::HijingHadronizer::readSettings ( int )

inline

Definition at line 36 of file HijingHadronizer.h.

36 { return true; }

◆ residualDecay()

bool HijingHadronizer::residualDecay ( )

Definition at line 282 of file HijingHadronizer.cc.

282 { return true; }

◆ rotateEvtPlane()

void HijingHadronizer::rotateEvtPlane ( )

private

Definition at line 271 of file HijingHadronizer.cc.

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

Referenced by generatePartonsAndHadronize().

                                       {
   phi0_ = 2. * pi * gen::hijran_(nullptr) - pi;
   sinphi0_ = sin(phi0_);
   cosphi0_ = cos(phi0_);
 }

◆ statistics()

void HijingHadronizer::statistics ( )

Definition at line 286 of file HijingHadronizer.cc.

286 { return; }

Member Data Documentation

◆ bmax_

double gen::HijingHadronizer::bmax_

private

Definition at line 61 of file HijingHadronizer.h.

Referenced by generatePartonsAndHadronize().

◆ bmin_

double gen::HijingHadronizer::bmin_

private

Definition at line 63 of file HijingHadronizer.h.

Referenced by generatePartonsAndHadronize().

◆ cosphi0_

double gen::HijingHadronizer::cosphi0_

private

Definition at line 79 of file HijingHadronizer.h.

Referenced by build_hijing(), build_hijing_vertex(), and rotateEvtPlane().

◆ efrm_

double gen::HijingHadronizer::efrm_

private

Definition at line 65 of file HijingHadronizer.h.

Referenced by initializeForInternalPartons().

◆ evt

HepMC::GenEvent* gen::HijingHadronizer::evt

private

Definition at line 59 of file HijingHadronizer.h.

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

◆ frame_

std::string gen::HijingHadronizer::frame_

private

Definition at line 66 of file HijingHadronizer.h.

Referenced by generatePartonsAndHadronize(), and initializeForInternalPartons().

◆ iap_

int gen::HijingHadronizer::iap_

private

Definition at line 69 of file HijingHadronizer.h.

Referenced by initializeForInternalPartons().

◆ iat_

int gen::HijingHadronizer::iat_

private

Definition at line 71 of file HijingHadronizer.h.

Referenced by initializeForInternalPartons().

◆ izp_

int gen::HijingHadronizer::izp_

private

Definition at line 70 of file HijingHadronizer.h.

Referenced by initializeForInternalPartons().

◆ izt_

int gen::HijingHadronizer::izt_

private

Definition at line 72 of file HijingHadronizer.h.

Referenced by initializeForInternalPartons().

◆ phi0_

double gen::HijingHadronizer::phi0_

private

Definition at line 77 of file HijingHadronizer.h.

Referenced by add_heavy_ion_rec(), and rotateEvtPlane().

◆ proj_

std::string gen::HijingHadronizer::proj_

private

Definition at line 67 of file HijingHadronizer.h.

Referenced by initializeForInternalPartons().

◆ pset_

edm::ParameterSet gen::HijingHadronizer::pset_

private

Definition at line 60 of file HijingHadronizer.h.

◆ rotate_

bool gen::HijingHadronizer::rotate_

private

Definition at line 80 of file HijingHadronizer.h.

Referenced by generatePartonsAndHadronize().

◆ sinphi0_

double gen::HijingHadronizer::sinphi0_

private

Definition at line 78 of file HijingHadronizer.h.

Referenced by build_hijing(), build_hijing_vertex(), and rotateEvtPlane().

◆ targ_

std::string gen::HijingHadronizer::targ_

private

Definition at line 68 of file HijingHadronizer.h.

Referenced by initializeForInternalPartons().

◆ theSharedResources

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

staticprivate

Definition at line 82 of file HijingHadronizer.h.

Referenced by doSharedResources().

Public Member Functions

Private Member Functions

Private Attributes

Static Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ HijingHadronizer()

◆ ~HijingHadronizer()

Member Function Documentation

◆ add_heavy_ion_rec()

◆ build_hijing()

◆ build_hijing_vertex()

◆ call_hijset()

◆ classname()

◆ decay()

◆ declareSpecialSettings()

◆ declareStableParticles()

◆ doSetRandomEngine()

◆ doSharedResources()

◆ finalizeEvent()

◆ generatePartonsAndHadronize()

◆ get_particles()

◆ hadronize()

◆ initializeForExternalPartons()

◆ initializeForInternalPartons()

◆ readSettings()

◆ residualDecay()

◆ rotateEvtPlane()

◆ statistics()

Member Data Documentation

◆ bmax_

◆ bmin_

◆ cosphi0_

◆ efrm_

◆ evt

◆ frame_

◆ iap_

◆ iat_

◆ izp_

◆ izt_

◆ phi0_

◆ proj_

◆ pset_

◆ rotate_

◆ sinphi0_

◆ targ_

◆ theSharedResources