gen::HijingHadronizer Class Reference

#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	residualDecay ()
void	statistics ()
virtual	~HijingHadronizer ()
Public Member Functions inherited from gen::BaseHadronizer
	BaseHadronizer (edm::ParameterSet const &ps)
edm::Event &	getEDMEvent () const
HepMC::GenEvent *	getGenEvent ()
GenEventInfoProduct *	getGenEventInfo ()
GenRunInfoProduct &	getGenRunInfo ()
const boost::shared_ptr < lhef::LHERunInfo > &	getLHERunInfo () const
void	resetEvent (HepMC::GenEvent *event)
void	resetEventInfo (GenEventInfoProduct *eventInfo)
void	setEDMEvent (edm::Event &event)
void	setLHEEvent (lhef::LHEEvent *event)
void	setLHERunInfo (lhef::LHERunInfo *runInfo)
	~BaseHadronizer ()

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

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

Detailed Description

Definition at line 27 of file HijingHadronizer.h.

Constructor & Destructor Documentation

HijingHadronizer::HijingHadronizer

(

const edm::ParameterSet &

pset

)

Definition at line 43 of file HijingHadronizer.cc.

References edm::RandomNumberGenerator::getEngine(), and hijRandomEngine.

                                                           :
    BaseHadronizer(pset),
    evt(0), 
    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
  Service<RandomNumberGenerator> rng;
  hijRandomEngine = &(rng->getEngine());

}

HijingHadronizer::~HijingHadronizer ( )

virtual

Definition at line 70 of file HijingHadronizer.cc.

{
  // destructor
}

Member Function Documentation

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

private

Definition at line 76 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 101 of file HijingHadronizer.cc.

References cosphi0_, configurableAnalysis::GenParticle, himain2, getHLTprescales::index, gen::p, sinphi0_, ExpressReco_HICollisions_FallBack::x, and ExpressReco_HICollisions_FallBack::y.

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

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

Referenced by get_particles().

{
   // build verteces for the hijing stored events                        
   HepMC::GenVertex* vertex = new HepMC::GenVertex(HepMC::FourVector(0,0,0,0),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 217 of file HijingHadronizer.cc.

References HIJSET.

Referenced by initializeForInternalPartons().

{
  // initialize hydjet  
   HIJSET(efrm,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 279 of file HijingHadronizer.cc.

{  
   return "gen::HijingHadronizer";
}

bool HijingHadronizer::decay

(

)

Definition at line 261 of file HijingHadronizer.cc.

{
   return true;
}

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

inline

Definition at line 39 of file HijingHadronizer.h.

{ return true; }

bool HijingHadronizer::declareStableParticles

(

const std::vector< int >

pdg

)

Definition at line 242 of file HijingHadronizer.cc.

{
   return true;
}

void HijingHadronizer::finalizeEvent

(

)

Definition at line 271 of file HijingHadronizer.cc.

                                    {
    return;
}

bool HijingHadronizer::generatePartonsAndHadronize

(

)

Definition at line 132 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
   HIJING(frame_.data(), bmin_, 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 155 of file HijingHadronizer.cc.

References build_hijing(), build_hijing_vertex(), configurableAnalysis::GenParticle, himain1, himain2, i, and LogDebug.

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

         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]=0; // 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 256 of file HijingHadronizer.cc.

{
   return false;
}

bool gen::HijingHadronizer::initializeForExternalPartons

(

)

bool HijingHadronizer::initializeForInternalPartons

(

)

Definition at line 225 of file HijingHadronizer.cc.

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

                                                   {

  //initialize pythia5

  if(0){
    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 HijingHadronizer::residualDecay

(

)

Definition at line 266 of file HijingHadronizer.cc.

{  
   return true;
}

void HijingHadronizer::rotateEvtPlane ( )

private

Definition at line 248 of file HijingHadronizer.cc.

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

Referenced by generatePartonsAndHadronize().

                                     {

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

void HijingHadronizer::statistics

(

)

Definition at line 275 of file HijingHadronizer.cc.

                                 {
    return;
}

Member Data Documentation

double gen::HijingHadronizer::bmax_

private

Definition at line 58 of file HijingHadronizer.h.

Referenced by generatePartonsAndHadronize().

double gen::HijingHadronizer::bmin_

private

Definition at line 60 of file HijingHadronizer.h.

Referenced by generatePartonsAndHadronize().

double gen::HijingHadronizer::cosphi0_

private

Definition at line 76 of file HijingHadronizer.h.

Referenced by build_hijing(), and rotateEvtPlane().

double gen::HijingHadronizer::efrm_

private

Definition at line 62 of file HijingHadronizer.h.

Referenced by initializeForInternalPartons().

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

private

Definition at line 56 of file HijingHadronizer.h.

Referenced by generatePartonsAndHadronize().

std::string gen::HijingHadronizer::frame_

private

Definition at line 63 of file HijingHadronizer.h.

Referenced by generatePartonsAndHadronize(), and initializeForInternalPartons().

int gen::HijingHadronizer::iap_

private

Definition at line 66 of file HijingHadronizer.h.

Referenced by initializeForInternalPartons().

int gen::HijingHadronizer::iat_

private

Definition at line 68 of file HijingHadronizer.h.

Referenced by initializeForInternalPartons().

int gen::HijingHadronizer::izp_

private

Definition at line 67 of file HijingHadronizer.h.

Referenced by initializeForInternalPartons().

int gen::HijingHadronizer::izt_

private

Definition at line 69 of file HijingHadronizer.h.

Referenced by initializeForInternalPartons().

double gen::HijingHadronizer::phi0_

private

Definition at line 74 of file HijingHadronizer.h.

Referenced by add_heavy_ion_rec(), and rotateEvtPlane().

std::string gen::HijingHadronizer::proj_

private

Definition at line 64 of file HijingHadronizer.h.

Referenced by initializeForInternalPartons().

edm::ParameterSet gen::HijingHadronizer::pset_

private

Definition at line 57 of file HijingHadronizer.h.

bool gen::HijingHadronizer::rotate_

private

Definition at line 77 of file HijingHadronizer.h.

Referenced by generatePartonsAndHadronize().

double gen::HijingHadronizer::sinphi0_

private

Definition at line 75 of file HijingHadronizer.h.

Referenced by build_hijing(), and rotateEvtPlane().

std::string gen::HijingHadronizer::targ_

private

Definition at line 65 of file HijingHadronizer.h.

Referenced by initializeForInternalPartons().