#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

HepMC::GenEvent *	getGenEvent ()

GenEventInfoProduct *	getGenEventInfo ()

virtual GenLumiInfoHeader *	getGenLumiInfoHeader () const

GenRunInfoProduct &	getGenRunInfo ()

const boost::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 (HepMC::GenEvent *event)

void	resetEventInfo (GenEventInfoProduct *eventInfo)

virtual bool	select (HepMC::GenEvent *) const

void	setEDMEvent (edm::Event &event)

void	setLHEEvent (lhef::LHEEvent *event)

void	setLHERunInfo (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_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::auto_ptr< HepMC::GenEvent > &	event ()

std::auto_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 28 of file HijingHadronizer.h.

Constructor & Destructor Documentation

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

Definition at line 60 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 ( )

override

Definition at line 84 of file HijingHadronizer.cc.

 {
   // destructor
 }

Member Function Documentation

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

private

Definition at line 98 of file HijingHadronizer.cc.

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

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

private

Definition at line 123 of file HijingHadronizer.cc.

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

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

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

private

Definition at line 152 of file HijingHadronizer.cc.

References cosphi0_, himain2, mps_fire::i, sinphi0_, and lumiQTWidget::t.

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

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 259 of file HijingHadronizer.cc.

References HIJSET.

Referenced by initializeForInternalPartons().

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

const char * HijingHadronizer::classname ( ) const

Definition at line 321 of file HijingHadronizer.cc.

 {  
    return "gen::HijingHadronizer";
 }

bool HijingHadronizer::decay ( )

Definition at line 303 of file HijingHadronizer.cc.

 {
    return true;
 }

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

inline

Definition at line 41 of file HijingHadronizer.h.

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

41 { return true; }

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

Definition at line 284 of file HijingHadronizer.cc.

 {
    return true;
 }

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

overrideprivatevirtual

Reimplemented from gen::BaseHadronizer.

Definition at line 91 of file HijingHadronizer.cc.

References hijRandomEngine, and gen::v.

 {
   hijRandomEngine = v;
 }

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

inlineoverrideprivatevirtual

Reimplemented from gen::BaseHadronizer.

Definition at line 50 of file HijingHadronizer.h.

References GenParticle::GenParticle, mps_fire::i, and AlCaHLTBitMon_QueryRunRegistry::string.

50 { return theSharedResources; }

gen::HijingHadronizer::theSharedResources

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

Definition: HijingHadronizer.h:84

void HijingHadronizer::finalizeEvent ( )

Definition at line 313 of file HijingHadronizer.cc.

                                     {
     return;
 }

bool HijingHadronizer::generatePartonsAndHadronize ( )

Definition at line 166 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_.data(), f_bmin, f_bmax, strlen(frame_.data()));
 
    // 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;
 }

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

private

Definition at line 191 of file HijingHadronizer.cc.

References build_hijing(), build_hijing_vertex(), GenParticle::GenParticle, himain1, himain2, mps_fire::i, LogDebug, HadronAndPartonSelector_cfi::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;
 }

bool HijingHadronizer::hadronize ( )

Definition at line 298 of file HijingHadronizer.cc.

 {
    return false;
 }

bool gen::HijingHadronizer::initializeForExternalPartons ( )

bool HijingHadronizer::initializeForInternalPartons ( )

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

bool gen::HijingHadronizer::readSettings ( int )

inline

Definition at line 37 of file HijingHadronizer.h.

37 { return true; }

bool HijingHadronizer::residualDecay ( )

Definition at line 308 of file HijingHadronizer.cc.

 {  
    return true;
 }

void HijingHadronizer::rotateEvtPlane ( )

private

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

void HijingHadronizer::statistics ( )

Definition at line 317 of file HijingHadronizer.cc.

                                  {
     return;
 }

Member Data Documentation

double gen::HijingHadronizer::bmax_

private

Definition at line 63 of file HijingHadronizer.h.

Referenced by generatePartonsAndHadronize().

double gen::HijingHadronizer::bmin_

private

Definition at line 65 of file HijingHadronizer.h.

Referenced by generatePartonsAndHadronize().

double gen::HijingHadronizer::cosphi0_

private

Definition at line 81 of file HijingHadronizer.h.

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

double gen::HijingHadronizer::efrm_

private

Definition at line 67 of file HijingHadronizer.h.

Referenced by initializeForInternalPartons().

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

private

Definition at line 61 of file HijingHadronizer.h.

Referenced by generatePartonsAndHadronize().

std::string gen::HijingHadronizer::frame_

private

Definition at line 68 of file HijingHadronizer.h.

Referenced by generatePartonsAndHadronize(), and initializeForInternalPartons().

int gen::HijingHadronizer::iap_

private

Definition at line 71 of file HijingHadronizer.h.

Referenced by initializeForInternalPartons().

int gen::HijingHadronizer::iat_

private

Definition at line 73 of file HijingHadronizer.h.

Referenced by initializeForInternalPartons().

int gen::HijingHadronizer::izp_

private

Definition at line 72 of file HijingHadronizer.h.

Referenced by initializeForInternalPartons().

int gen::HijingHadronizer::izt_

private

Definition at line 74 of file HijingHadronizer.h.

Referenced by initializeForInternalPartons().

double gen::HijingHadronizer::phi0_

private

Definition at line 79 of file HijingHadronizer.h.

Referenced by add_heavy_ion_rec(), and rotateEvtPlane().

std::string gen::HijingHadronizer::proj_

private

Definition at line 69 of file HijingHadronizer.h.

Referenced by initializeForInternalPartons().

edm::ParameterSet gen::HijingHadronizer::pset_

private

Definition at line 62 of file HijingHadronizer.h.

bool gen::HijingHadronizer::rotate_

private

Definition at line 82 of file HijingHadronizer.h.

Referenced by generatePartonsAndHadronize().

double gen::HijingHadronizer::sinphi0_

private

Definition at line 80 of file HijingHadronizer.h.

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

std::string gen::HijingHadronizer::targ_

private

Definition at line 70 of file HijingHadronizer.h.

Referenced by initializeForInternalPartons().

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

staticprivate

Definition at line 84 of file HijingHadronizer.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