|
|
|
#include <AMPTHadronizer.h>
Public Member Functions | |
| AMPTHadronizer (const edm::ParameterSet &) | |
| 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 () |
| bool | initializeForExternalPartons () |
| bool | initializeForInternalPartons () |
| bool | readSettings (int) |
| bool | residualDecay () |
| void | statistics () |
| virtual | ~AMPTHadronizer () |
Private Member Functions | |
| void | add_heavy_ion_rec (HepMC::GenEvent *evt) |
| bool | ampt_init (const edm::ParameterSet &pset) |
| HepMC::GenParticle * | build_ampt (int index, int barcode) |
| HepMC::GenVertex * | build_ampt_vertex (int i, int id) |
| bool | call_amptset (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 | alpha_ |
| int | amptmode_ |
| double | bmax_ |
| double | bmin_ |
| double | cosphi0_ |
| int | deuteronfactor_ |
| int | deuteronmode_ |
| int | deuteronxsec_ |
| int | diquarkembedding_ |
| double | diquarkpx_ |
| double | diquarkpy_ |
| double | diquarkx_ |
| double | diquarky_ |
| int | doInitialAndFinalRadiation_ |
| double | dpcoal_ |
| double | drcoal_ |
| double | dt_ |
| double | efrm_ |
| HepMC::GenEvent * | evt |
| std::string | frame_ |
| int | iap_ |
| int | iat_ |
| int | izp_ |
| int | izpc_ |
| int | izt_ |
| bool | ks0decay_ |
| int | ktkick_ |
| int | maxmiss_ |
| double | minijetpt_ |
| double | mu_ |
| int | nsembd_ |
| int | ntmax_ |
| double | phi0_ |
| bool | phidecay_ |
| bool | popcornmode_ |
| double | popcornpar_ |
| std::string | proj_ |
| double | psembd_ |
| edm::ParameterSet | pset_ |
| double | pthard_ |
| bool | quenchingmode_ |
| double | quenchingpar_ |
| bool | rotate_ |
| double | shadowingfactor_ |
| bool | shadowingmode_ |
| bool | shadowingmodflag_ |
| double | sinphi0_ |
| double | stringFragA_ |
| double | stringFragB_ |
| std::string | targ_ |
| double | tmaxembd_ |
Definition at line 32 of file AMPTHadronizer.h.
| AMPTHadronizer::AMPTHadronizer | ( | const edm::ParameterSet & | pset | ) |
Definition at line 54 of file AMPTHadronizer.cc.
References _amptRandomEngine.
:
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")),
amptmode_(pset.getParameter<int>("amptmode")),
ntmax_(pset.getParameter<int>("ntmax")),
dt_(pset.getParameter<double>("dt")),
stringFragA_(pset.getParameter<double>("stringFragA")),
stringFragB_(pset.getParameter<double>("stringFragB")),
popcornmode_(pset.getParameter<bool>("popcornmode")),
popcornpar_(pset.getParameter<double>("popcornpar")),
shadowingmode_(pset.getParameter<bool>("shadowingmode")),
quenchingmode_(pset.getParameter<bool>("quenchingmode")),
quenchingpar_(pset.getParameter<double>("quenchingpar")),
pthard_(pset.getParameter<double>("pthard")),
mu_(pset.getParameter<double>("mu")),
izpc_(pset.getParameter<int>("izpc")),
alpha_(pset.getParameter<double>("alpha")),
dpcoal_(pset.getParameter<double>("dpcoal")),
drcoal_(pset.getParameter<double>("drcoal")),
ks0decay_(pset.getParameter<bool>("ks0decay")),
phidecay_(pset.getParameter<bool>("phidecay")),
deuteronmode_(pset.getParameter<int>("deuteronmode")),
deuteronfactor_(pset.getParameter<int>("deuteronfactor")),
deuteronxsec_(pset.getParameter<int>("deuteronxsec")),
minijetpt_(pset.getParameter<double>("minijetpt")),
maxmiss_(pset.getParameter<int>("maxmiss")),
doInitialAndFinalRadiation_(pset.getParameter<int>("doInitialAndFinalRadiation")),
ktkick_(pset.getParameter<int>("ktkick")),
diquarkembedding_(pset.getParameter<int>("diquarkembedding")),
diquarkpx_(pset.getParameter<double>("diquarkpx")),
diquarkpy_(pset.getParameter<double>("diquarkpy")),
diquarkx_(pset.getParameter<double>("diquarkx")),
diquarky_(pset.getParameter<double>("diquarky")),
phi0_(0.),
sinphi0_(0.),
cosphi0_(1.),
rotate_(pset.getParameter<bool>("rotateEventPlane"))
{
// Default constructor
edm::Service<RandomNumberGenerator> rng;
_amptRandomEngine = &(rng->getEngine());
}
| AMPTHadronizer::~AMPTHadronizer | ( | ) | [virtual] |
Definition at line 110 of file AMPTHadronizer.cc.
{
}
| void AMPTHadronizer::add_heavy_ion_rec | ( | HepMC::GenEvent * | evt | ) | [private] |
Definition at line 115 of file AMPTHadronizer.cc.
References hmain1, hparnt, and phi0_.
Referenced by generatePartonsAndHadronize().
{
// heavy ion record in the final CMSSW Event
HepMC::HeavyIon* hi = new HepMC::HeavyIon(
hmain1.jatt, // Ncoll_hard/N of SubEvents
hmain1.np, // Npart_proj
hmain1.nt, // Npart_targ
hmain1.n0+hmain1.n01+hmain1.n10+hmain1.n11, // Ncoll
0, // spectator_neutrons
0, // spectator_protons
hmain1.n01, // N_Nwounded_collisions
hmain1.n10, // Nwounded_N_collisions
hmain1.n11, // Nwounded_Nwounded_collisions
hparnt.hint1[18], // impact_parameter in [fm]
phi0_, // event_plane_angle
0, // eccentricity
hparnt.hint1[11] // sigma_inel_NN
);
evt->set_heavy_ion(*hi);
delete hi;
}
| bool AMPTHadronizer::ampt_init | ( | const edm::ParameterSet & | pset | ) | [private] |
Definition at line 263 of file AMPTHadronizer.cc.
References alpha_, amptmode_, anim, coal, deuteronfactor_, deuteronmode_, deuteronxsec_, diquarkembedding_, diquarkpx_, diquarkpy_, diquarkx_, diquarky_, doInitialAndFinalRadiation_, dpcoal_, drcoal_, dt_, embed, hparnt, input1, input2, izpc_, ks0decay_, ktkick_, ludat1, maxmiss_, minijetpt_, mu_, ntmax_, para2, para8, phidcy, phidecay_, popcorn, popcornmode_, popcornpar_, pthard_, quenchingmode_, quenchingpar_, resdcy, shadowingmode_, stringFragA_, and stringFragB_.
Referenced by initializeForInternalPartons().
{
anim.isoft=amptmode_;
input2.ntmax=ntmax_;
input1.dt=dt_;
ludat1.parj[40]=stringFragA_;
ludat1.parj[41]=stringFragB_;
popcorn.ipop=popcornmode_;
ludat1.parj[4]=popcornpar_;
hparnt.ihpr2[5]=shadowingmode_;
hparnt.ihpr2[3]=quenchingmode_;
hparnt.hipr1[13]=quenchingpar_;
hparnt.hipr1[7]=pthard_;
para2.xmu=mu_;
anim.izpc=izpc_;
para2.alpha=alpha_;
coal.dpcoal=dpcoal_;
coal.drcoal=drcoal_;
resdcy.iksdcy=ks0decay_;
phidcy.iphidcy=phidecay_;
para8.idpert=deuteronmode_;
para8.npertd=deuteronfactor_;
para8.idxsec=deuteronxsec_;
phidcy.pttrig=minijetpt_;
phidcy.maxmiss=maxmiss_;
hparnt.ihpr2[1]=doInitialAndFinalRadiation_;
hparnt.ihpr2[4]=ktkick_;
embed.iembed=diquarkembedding_;
embed.pxqembd=diquarkpx_;
embed.pyqembd=diquarkpy_;
embed.xembd=diquarkx_;
embed.yembd=diquarky_;
return true;
}
| HepMC::GenParticle * AMPTHadronizer::build_ampt | ( | int | index, |
| int | barcode | ||
| ) | [private] |
Definition at line 138 of file AMPTHadronizer.cc.
References cosphi0_, alignCSCRings::e, configurableAnalysis::GenParticle, hbt, getHLTprescales::index, INVFLV, m, gen::p, sinphi0_, mathSSE::sqrt(), and ntuplemaker::status.
Referenced by get_particles().
{
// Build particle object corresponding to index in ampt
float px0 = hbt.plast[index][0];
float py0 = hbt.plast[index][1];
float pz0 = hbt.plast[index][2];
float m = hbt.plast[index][3];
float px = px0*cosphi0_-py0*sinphi0_;
float py = py0*cosphi0_+px0*sinphi0_;
float pz = pz0;
float e = sqrt(px*px+py*py+pz*pz+m*m);
int status = 1;
HepMC::GenParticle* p = new HepMC::GenParticle(
HepMC::FourVector(px,
py,
pz,
e),
INVFLV(hbt.lblast[index]),// id
status // status
);
p->suggest_barcode(barcode);
return p;
}
| HepMC::GenVertex * AMPTHadronizer::build_ampt_vertex | ( | int | i, |
| int | id | ||
| ) | [private] |
Definition at line 167 of file AMPTHadronizer.cc.
Referenced by get_particles().
{
// build verteces for the ampt stored events
HepMC::GenVertex* vertex = new HepMC::GenVertex(HepMC::FourVector(0,0,0,0),id);
return vertex;
}
| bool AMPTHadronizer::call_amptset | ( | double | efrm, |
| std::string | frame, | ||
| std::string | proj, | ||
| std::string | targ, | ||
| int | iap, | ||
| int | izp, | ||
| int | iat, | ||
| int | izt | ||
| ) | [private] |
Definition at line 256 of file AMPTHadronizer.cc.
References AMPTSET.
Referenced by initializeForInternalPartons().
| const char * AMPTHadronizer::classname | ( | ) | const |
Definition at line 351 of file AMPTHadronizer.cc.
{
return "gen::AMPTHadronizer";
}
| bool AMPTHadronizer::decay | ( | ) |
Definition at line 333 of file AMPTHadronizer.cc.
{
return true;
}
| bool gen::AMPTHadronizer::declareSpecialSettings | ( | const std::vector< std::string > & | ) | [inline] |
Definition at line 45 of file AMPTHadronizer.h.
{ return true; }
| bool AMPTHadronizer::declareStableParticles | ( | const std::vector< int > & | pdg | ) |
Definition at line 313 of file AMPTHadronizer.cc.
{
return true;
}
| void AMPTHadronizer::finalizeEvent | ( | ) |
Definition at line 343 of file AMPTHadronizer.cc.
{
return;
}
| bool AMPTHadronizer::generatePartonsAndHadronize | ( | ) |
Definition at line 174 of file AMPTHadronizer.cc.
References add_heavy_ion_rec(), AMPT, bmax_, bmin_, gen::BaseHadronizer::event(), evt, frame_, get_particles(), rotate_, and rotateEvtPlane().
{
// generate single event
if(rotate_) rotateEvtPlane();
// generate a AMPT event
AMPT(frame_.data(), bmin_, bmax_, strlen(frame_.data()));
// event information
HepMC::GenEvent *evt = new HepMC::GenEvent();
get_particles(evt);
add_heavy_ion_rec(evt);
event().reset(evt);
return true;
}
| bool AMPTHadronizer::get_particles | ( | HepMC::GenEvent * | evt | ) | [private] |
Definition at line 195 of file AMPTHadronizer.cc.
References build_ampt(), build_ampt_vertex(), configurableAnalysis::GenParticle, hbt, 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 = hbt.nlast;
for (unsigned int ipart = 0; ipart<knumpart; ipart++) {
int mid = 0;
particles.push_back(build_ampt(ipart,ipart+1));
prods.push_back(build_ampt_vertex(ipart,0));
mother_ids.push_back(mid);
LogDebug("DecayChain")<<"Mother index : "<<mid;
}
LogDebug("AMPT")<<"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 AMPTHadronizer::hadronize | ( | ) |
Definition at line 328 of file AMPTHadronizer.cc.
{
return false;
}
| bool gen::AMPTHadronizer::initializeForExternalPartons | ( | ) |
| bool AMPTHadronizer::initializeForInternalPartons | ( | ) |
Definition at line 300 of file AMPTHadronizer.cc.
References ampt_init(), call_amptset(), efrm_, frame_, iap_, iat_, izp_, izt_, proj_, pset_, and targ_.
| bool gen::AMPTHadronizer::readSettings | ( | int | ) | [inline] |
Definition at line 41 of file AMPTHadronizer.h.
{ return true; }
| bool AMPTHadronizer::residualDecay | ( | ) |
Definition at line 338 of file AMPTHadronizer.cc.
{
return true;
}
| void AMPTHadronizer::rotateEvtPlane | ( | ) | [private] |
Definition at line 319 of file AMPTHadronizer.cc.
References funct::cos(), cosphi0_, phi0_, pi, gen::ranart_(), funct::sin(), sinphi0_, run_regression::test, and zero.
Referenced by generatePartonsAndHadronize().
| void AMPTHadronizer::statistics | ( | ) |
Definition at line 347 of file AMPTHadronizer.cc.
{
return;
}
double gen::AMPTHadronizer::alpha_ [private] |
Definition at line 89 of file AMPTHadronizer.h.
Referenced by ampt_init().
int gen::AMPTHadronizer::amptmode_ [private] |
Definition at line 76 of file AMPTHadronizer.h.
Referenced by ampt_init().
double gen::AMPTHadronizer::bmax_ [private] |
Definition at line 64 of file AMPTHadronizer.h.
Referenced by generatePartonsAndHadronize().
double gen::AMPTHadronizer::bmin_ [private] |
Definition at line 66 of file AMPTHadronizer.h.
Referenced by generatePartonsAndHadronize().
double gen::AMPTHadronizer::cosphi0_ [private] |
Definition at line 113 of file AMPTHadronizer.h.
Referenced by build_ampt(), and rotateEvtPlane().
int gen::AMPTHadronizer::deuteronfactor_ [private] |
Definition at line 95 of file AMPTHadronizer.h.
Referenced by ampt_init().
int gen::AMPTHadronizer::deuteronmode_ [private] |
Definition at line 94 of file AMPTHadronizer.h.
Referenced by ampt_init().
int gen::AMPTHadronizer::deuteronxsec_ [private] |
Definition at line 96 of file AMPTHadronizer.h.
Referenced by ampt_init().
int gen::AMPTHadronizer::diquarkembedding_ [private] |
Definition at line 101 of file AMPTHadronizer.h.
Referenced by ampt_init().
double gen::AMPTHadronizer::diquarkpx_ [private] |
Definition at line 102 of file AMPTHadronizer.h.
Referenced by ampt_init().
double gen::AMPTHadronizer::diquarkpy_ [private] |
Definition at line 103 of file AMPTHadronizer.h.
Referenced by ampt_init().
double gen::AMPTHadronizer::diquarkx_ [private] |
Definition at line 104 of file AMPTHadronizer.h.
Referenced by ampt_init().
double gen::AMPTHadronizer::diquarky_ [private] |
Definition at line 105 of file AMPTHadronizer.h.
Referenced by ampt_init().
int gen::AMPTHadronizer::doInitialAndFinalRadiation_ [private] |
Definition at line 99 of file AMPTHadronizer.h.
Referenced by ampt_init().
double gen::AMPTHadronizer::dpcoal_ [private] |
Definition at line 90 of file AMPTHadronizer.h.
Referenced by ampt_init().
double gen::AMPTHadronizer::drcoal_ [private] |
Definition at line 91 of file AMPTHadronizer.h.
Referenced by ampt_init().
double gen::AMPTHadronizer::dt_ [private] |
Definition at line 78 of file AMPTHadronizer.h.
Referenced by ampt_init().
double gen::AMPTHadronizer::efrm_ [private] |
Definition at line 68 of file AMPTHadronizer.h.
Referenced by initializeForInternalPartons().
HepMC::GenEvent* gen::AMPTHadronizer::evt [private] |
Definition at line 62 of file AMPTHadronizer.h.
Referenced by generatePartonsAndHadronize().
std::string gen::AMPTHadronizer::frame_ [private] |
Definition at line 69 of file AMPTHadronizer.h.
Referenced by generatePartonsAndHadronize(), and initializeForInternalPartons().
int gen::AMPTHadronizer::iap_ [private] |
Definition at line 72 of file AMPTHadronizer.h.
Referenced by initializeForInternalPartons().
int gen::AMPTHadronizer::iat_ [private] |
Definition at line 74 of file AMPTHadronizer.h.
Referenced by initializeForInternalPartons().
int gen::AMPTHadronizer::izp_ [private] |
Definition at line 73 of file AMPTHadronizer.h.
Referenced by initializeForInternalPartons().
int gen::AMPTHadronizer::izpc_ [private] |
Definition at line 88 of file AMPTHadronizer.h.
Referenced by ampt_init().
int gen::AMPTHadronizer::izt_ [private] |
Definition at line 75 of file AMPTHadronizer.h.
Referenced by initializeForInternalPartons().
bool gen::AMPTHadronizer::ks0decay_ [private] |
Definition at line 92 of file AMPTHadronizer.h.
Referenced by ampt_init().
int gen::AMPTHadronizer::ktkick_ [private] |
Definition at line 100 of file AMPTHadronizer.h.
Referenced by ampt_init().
int gen::AMPTHadronizer::maxmiss_ [private] |
Definition at line 98 of file AMPTHadronizer.h.
Referenced by ampt_init().
double gen::AMPTHadronizer::minijetpt_ [private] |
Definition at line 97 of file AMPTHadronizer.h.
Referenced by ampt_init().
double gen::AMPTHadronizer::mu_ [private] |
Definition at line 87 of file AMPTHadronizer.h.
Referenced by ampt_init().
int gen::AMPTHadronizer::nsembd_ [private] |
Definition at line 106 of file AMPTHadronizer.h.
int gen::AMPTHadronizer::ntmax_ [private] |
Definition at line 77 of file AMPTHadronizer.h.
Referenced by ampt_init().
double gen::AMPTHadronizer::phi0_ [private] |
Definition at line 111 of file AMPTHadronizer.h.
Referenced by add_heavy_ion_rec(), and rotateEvtPlane().
bool gen::AMPTHadronizer::phidecay_ [private] |
Definition at line 93 of file AMPTHadronizer.h.
Referenced by ampt_init().
bool gen::AMPTHadronizer::popcornmode_ [private] |
Definition at line 81 of file AMPTHadronizer.h.
Referenced by ampt_init().
double gen::AMPTHadronizer::popcornpar_ [private] |
Definition at line 82 of file AMPTHadronizer.h.
Referenced by ampt_init().
std::string gen::AMPTHadronizer::proj_ [private] |
Definition at line 70 of file AMPTHadronizer.h.
Referenced by initializeForInternalPartons().
double gen::AMPTHadronizer::psembd_ [private] |
Definition at line 107 of file AMPTHadronizer.h.
edm::ParameterSet gen::AMPTHadronizer::pset_ [private] |
Definition at line 63 of file AMPTHadronizer.h.
Referenced by initializeForInternalPartons().
double gen::AMPTHadronizer::pthard_ [private] |
Definition at line 86 of file AMPTHadronizer.h.
Referenced by ampt_init().
bool gen::AMPTHadronizer::quenchingmode_ [private] |
Definition at line 84 of file AMPTHadronizer.h.
Referenced by ampt_init().
double gen::AMPTHadronizer::quenchingpar_ [private] |
Definition at line 85 of file AMPTHadronizer.h.
Referenced by ampt_init().
bool gen::AMPTHadronizer::rotate_ [private] |
Definition at line 114 of file AMPTHadronizer.h.
Referenced by generatePartonsAndHadronize().
double gen::AMPTHadronizer::shadowingfactor_ [private] |
Definition at line 110 of file AMPTHadronizer.h.
bool gen::AMPTHadronizer::shadowingmode_ [private] |
Definition at line 83 of file AMPTHadronizer.h.
Referenced by ampt_init().
bool gen::AMPTHadronizer::shadowingmodflag_ [private] |
Definition at line 109 of file AMPTHadronizer.h.
double gen::AMPTHadronizer::sinphi0_ [private] |
Definition at line 112 of file AMPTHadronizer.h.
Referenced by build_ampt(), and rotateEvtPlane().
double gen::AMPTHadronizer::stringFragA_ [private] |
Definition at line 79 of file AMPTHadronizer.h.
Referenced by ampt_init().
double gen::AMPTHadronizer::stringFragB_ [private] |
Definition at line 80 of file AMPTHadronizer.h.
Referenced by ampt_init().
std::string gen::AMPTHadronizer::targ_ [private] |
Definition at line 71 of file AMPTHadronizer.h.
Referenced by initializeForInternalPartons().
double gen::AMPTHadronizer::tmaxembd_ [private] |
Definition at line 108 of file AMPTHadronizer.h.
1.7.3