GENIEGenerator
Loading...
Searching...
No Matches
Public Member Functions | Private Member Functions | List of all members
genie::Pythia6Hadro2019 Class Reference

Provides access to the PYTHIA hadronization models.
Is a concrete implementation of the EventRecordVisitorI interface. More...

#include <Pythia6Hadro2019.h>

Inheritance diagram for genie::Pythia6Hadro2019:
[legend]
Collaboration diagram for genie::Pythia6Hadro2019:
[legend]

Public Member Functions

 Pythia6Hadro2019 ()
 Pythia6Hadro2019 (string config)
virtual ~Pythia6Hadro2019 ()
void ProcessEventRecord (GHepRecord *event) const
void Configure (const Registry &config)
void Configure (string config)
Public Member Functions inherited from genie::EventRecordVisitorI
virtual ~EventRecordVisitorI ()
Public Member Functions inherited from genie::Algorithm
virtual ~Algorithm ()
virtual void FindConfig (void)
virtual const Registry & GetConfig (void) const
Registry * GetOwnedConfig (void)
virtual const AlgId & Id (void) const
 Get algorithm ID.
virtual AlgStatus_t GetStatus (void) const
 Get algorithm status.
virtual bool AllowReconfig (void) const
virtual AlgCmp_t Compare (const Algorithm *alg) const
 Compare with input algorithm.
virtual void SetId (const AlgId &id)
 Set algorithm ID.
virtual void SetId (string name, string config)
const Algorithm * SubAlg (const RgKey &registry_key) const
void AdoptConfig (void)
void AdoptSubstructure (void)
virtual void Print (ostream &stream) const
 Print algorithm info.

Private Member Functions

bool Hadronize (GHepRecord *event) const
void CopyOriginalDecayFlags (void) const
void SetDesiredDecayFlags (void) const
void RestoreOriginalDecayFlags (void) const
void LoadConfig (void)
void Initialize (void)

Additional Inherited Members

Static Public Member Functions inherited from genie::Algorithm
static string BuildParamVectKey (const std::string &comm_name, unsigned int i)
static string BuildParamVectSizeKey (const std::string &comm_name)
static string BuildParamMatKey (const std::string &comm_name, unsigned int i, unsigned int j)
static string BuildParamMatRowSizeKey (const std::string &comm_name)
static string BuildParamMatColSizeKey (const std::string &comm_name)
Protected Member Functions inherited from genie::PythiaBaseHadro2019
 PythiaBaseHadro2019 ()
 PythiaBaseHadro2019 (string name)
 PythiaBaseHadro2019 (string name, string config)
virtual ~PythiaBaseHadro2019 ()
virtual void MakeQuarkDiquarkAssignments (const Interaction *in) const
virtual bool AssertValidity (const Interaction *in) const
Protected Member Functions inherited from genie::EventRecordVisitorI
 EventRecordVisitorI ()
 EventRecordVisitorI (string name)
 EventRecordVisitorI (string name, string config)
Protected Member Functions inherited from genie::Algorithm
 Algorithm ()
 Algorithm (string name)
 Algorithm (string name, string config)
void Initialize (void)
void DeleteConfig (void)
void DeleteSubstructure (void)
Registry * ExtractLocalConfig (const Registry &in) const
Registry * ExtractLowerConfig (const Registry &in, const string &alg_key) const
 Split an incoming configuration Registry into a block valid for the sub-algo identified by alg_key.
template<class T>
bool GetParam (const RgKey &name, T &p, bool is_top_call=true) const
template<class T>
bool GetParamDef (const RgKey &name, T &p, const T &def) const
template<class T>
int GetParamVect (const std::string &comm_name, std::vector< T > &v, bool is_top_call=true) const
 Handle to load vectors of parameters.
int GetParamVectKeys (const std::string &comm_name, std::vector< RgKey > &k, bool is_top_call=true) const
template<class T>
int GetParamMat (const std::string &comm_name, TMatrixT< T > &mat, bool is_top_call=true) const
 Handle to load matrix of parameters.
template<class T>
int GetParamMatSym (const std::string &comm_name, TMatrixTSym< T > &mat, bool is_top_call=true) const
int GetParamMatKeys (const std::string &comm_name, std::vector< RgKey > &k, bool is_top_call=true) const
int AddTopRegistry (Registry *rp, bool owns=true)
 add registry with top priority, also update ownership
int AddLowRegistry (Registry *rp, bool owns=true)
 add registry with lowest priority, also update ownership
int MergeTopRegistry (const Registry &r)
int AddTopRegisties (const vector< Registry * > &rs, bool owns=false)
 Add registries with top priority, also udated Ownerships.
Protected Attributes inherited from genie::PythiaBaseHadro2019
int fLeadingQuark
int fRemnantDiquark
double fSSBarSuppression
 ssbar suppression
double fGaussianPt2
 gaussian pt2 distribution width
double fNonGaussianPt2Tail
 non gaussian pt2 tail parameterization
double fRemainingECutoff
 remaining E cutoff stopping fragmentation
double fDiQuarkSuppression
 di-quark suppression parameter
double fLightVMesonSuppression
 light vector meson suppression
double fSVMesonSuppression
 strange vector meson suppression
double fLunda
 Lund a parameter.
double fLundb
 Lund b parameter.
double fLundaDiq
 adjustment of Lund a for di-quark
bool fOriDecayFlag_pi0
bool fOriDecayFlag_K0
bool fOriDecayFlag_K0b
bool fOriDecayFlag_L0
bool fOriDecayFlag_L0b
bool fOriDecayFlag_Dm
bool fOriDecayFlag_D0
bool fOriDecayFlag_Dp
bool fOriDecayFlag_Dpp
bool fReqDecayFlag_pi0
bool fReqDecayFlag_K0
bool fReqDecayFlag_K0b
bool fReqDecayFlag_L0
bool fReqDecayFlag_L0b
bool fReqDecayFlag_Dm
bool fReqDecayFlag_D0
bool fReqDecayFlag_Dp
bool fReqDecayFlag_Dpp
Protected Attributes inherited from genie::Algorithm
bool fAllowReconfig
bool fOwnsSubstruc
 true if it owns its substructure (sub-algs,...)
AlgId fID
 algorithm name and configuration set
vector< Registry * > fConfVect
vector< bool > fOwnerships
 ownership for every registry in fConfVect
AlgStatus_t fStatus
 algorithm execution status
AlgMap * fOwnedSubAlgMp
 local pool for owned sub-algs (taken out of the factory pool)

Detailed Description

Provides access to the PYTHIA hadronization models.
Is a concrete implementation of the EventRecordVisitorI interface.

Author
Costas Andreopoulos <c.andreopoulos \at cern.ch> University of Liverpool
Created:\n August 17, 2004
License:\n Copyright (c) 2003-2025, The GENIE Collaboration
For the full text of the license visit http://copyright.genie-mc.org

Definition at line 34 of file Pythia6Hadro2019.h.

Constructor & Destructor Documentation

◆ Pythia6Hadro2019() [1/2]

Pythia6Hadro2019::Pythia6Hadro2019 ( )

Definition at line 43 of file Pythia6Hadro2019.cxx.

43 :
44PythiaBaseHadro2019("genie::Pythia6Hadro2019")
45{
46 this->Initialize();
47}

References Initialize(), and genie::PythiaBaseHadro2019::PythiaBaseHadro2019().

◆ Pythia6Hadro2019() [2/2]

Pythia6Hadro2019::Pythia6Hadro2019 ( string config)

Definition at line 49 of file Pythia6Hadro2019.cxx.

49 :
50PythiaBaseHadro2019("genie::Pythia6Hadro2019", config)
51{
52 this->Initialize();
53}

References Initialize(), and genie::PythiaBaseHadro2019::PythiaBaseHadro2019().

◆ ~Pythia6Hadro2019()

Pythia6Hadro2019::~Pythia6Hadro2019 ( )
virtual

Definition at line 55 of file Pythia6Hadro2019.cxx.

56{
57
58}

Member Function Documentation

◆ Configure() [1/2]

void Pythia6Hadro2019::Configure ( const Registry & config)
virtual

Configure the algorithm with an external registry The registry is merged with the top level registry if it is owned, Otherwise a copy of it is added with the highest priority

Reimplemented from genie::Algorithm.

Definition at line 294 of file Pythia6Hadro2019.cxx.

295{
296 Algorithm::Configure(config);
297 this->LoadConfig();
298}
genie::Algorithm::Configure
virtual void Configure(const Registry &config)
Definition Algorithm.cxx:62

References genie::Algorithm::Configure(), and LoadConfig().

◆ Configure() [2/2]

void Pythia6Hadro2019::Configure ( string config)
virtual

Configure the algorithm from the AlgoConfigPool based on param_set string given in input An algorithm contains a vector of registries coming from different xml configuration files, which are loaded according a very precise prioriy This methods will load a number registries in order of priority: 1) "Tunable" parameter set from CommonParametes. This is loaded with the highest prioriry and it is designed to be used for tuning procedure Usage not expected from the user. 2) For every string defined in "CommonParame" the corresponding parameter set will be loaded from CommonParameter.xml 3) parameter set specified by the config string and defined in the xml file of the algorithm 4) if config is not "Default" also the Default parameter set from the same xml file will be loaded Effectively this avoids the repetion of a parameter when it is not changed in the requested configuration

Reimplemented from genie::Algorithm.

Definition at line 300 of file Pythia6Hadro2019.cxx.

301{
302 Algorithm::Configure(config);
303 this->LoadConfig();
304}

References genie::Algorithm::Configure(), and LoadConfig().

◆ CopyOriginalDecayFlags()

void Pythia6Hadro2019::CopyOriginalDecayFlags ( void ) const
privatevirtual

Implements genie::PythiaBaseHadro2019.

Definition at line 212 of file Pythia6Hadro2019.cxx.

213{
214#ifdef __GENIE_PYTHIA6_ENABLED__
215 fOriDecayFlag_pi0 = (fPythia->GetMDCY(fPythia->Pycomp(kPdgPi0), 1) == 1);
216 fOriDecayFlag_K0 = (fPythia->GetMDCY(fPythia->Pycomp(kPdgK0), 1) == 1);
217 fOriDecayFlag_K0b = (fPythia->GetMDCY(fPythia->Pycomp(kPdgAntiK0), 1) == 1);
218 fOriDecayFlag_L0 = (fPythia->GetMDCY(fPythia->Pycomp(kPdgLambda), 1) == 1);
219 fOriDecayFlag_L0b = (fPythia->GetMDCY(fPythia->Pycomp(kPdgAntiLambda), 1) == 1);
220 fOriDecayFlag_Dm = (fPythia->GetMDCY(fPythia->Pycomp(kPdgP33m1232_DeltaM), 1) == 1);
221 fOriDecayFlag_D0 = (fPythia->GetMDCY(fPythia->Pycomp(kPdgP33m1232_Delta0), 1) == 1);
222 fOriDecayFlag_Dp = (fPythia->GetMDCY(fPythia->Pycomp(kPdgP33m1232_DeltaP), 1) == 1);
223 fOriDecayFlag_Dpp = (fPythia->GetMDCY(fPythia->Pycomp(kPdgP33m1232_DeltaPP), 1) == 1);
224
225#ifdef __GENIE_LOW_LEVEL_MESG_ENABLED__
226 LOG("Pythia6Had", pDEBUG)
227 << "Original PYTHIA6 decay flags:"
228 << "\n pi0 = " << fOriDecayFlag_pi0
229 << "\n K0 = " << fOriDecayFlag_K0
230 << "\n \bar{K0} = " << fOriDecayFlag_K0b
231 << "\n Lambda = " << fOriDecayFlag_L0
232 << "\n \bar{Lambda0} = " << fOriDecayFlag_L0b
233 << "\n D- = " << fOriDecayFlag_Dm
234 << "\n D0 = " << fOriDecayFlag_D0
235 << "\n D+ = " << fOriDecayFlag_Dp
236 << "\n D++ = " << fOriDecayFlag_Dpp;
237#endif
238
239#endif // __GENIE_PYTHIA6_ENABLED__
240}
pDEBUG
#define pDEBUG
Definition Messenger.h:63
LOG
#define LOG(stream, priority)
A macro that returns the requested log4cpp::Category appending a string (using the FILE,...
Definition Messenger.h:96
genie::kPdgP33m1232_Delta0
const int kPdgP33m1232_Delta0
Definition PDGCodes.h:105
genie::kPdgPi0
const int kPdgPi0
Definition PDGCodes.h:160
genie::kPdgP33m1232_DeltaPP
const int kPdgP33m1232_DeltaPP
Definition PDGCodes.h:107
genie::kPdgAntiLambda
const int kPdgAntiLambda
Definition PDGCodes.h:86
genie::kPdgP33m1232_DeltaM
const int kPdgP33m1232_DeltaM
Definition PDGCodes.h:104
genie::kPdgLambda
const int kPdgLambda
Definition PDGCodes.h:85
genie::kPdgP33m1232_DeltaP
const int kPdgP33m1232_DeltaP
Definition PDGCodes.h:106
genie::kPdgAntiK0
const int kPdgAntiK0
Definition PDGCodes.h:175
genie::kPdgK0
const int kPdgK0
Definition PDGCodes.h:174

References genie::PythiaBaseHadro2019::fOriDecayFlag_D0, genie::PythiaBaseHadro2019::fOriDecayFlag_Dm, genie::PythiaBaseHadro2019::fOriDecayFlag_Dp, genie::PythiaBaseHadro2019::fOriDecayFlag_Dpp, genie::PythiaBaseHadro2019::fOriDecayFlag_K0, genie::PythiaBaseHadro2019::fOriDecayFlag_K0b, genie::PythiaBaseHadro2019::fOriDecayFlag_L0, genie::PythiaBaseHadro2019::fOriDecayFlag_L0b, genie::PythiaBaseHadro2019::fOriDecayFlag_pi0, genie::kPdgAntiK0, genie::kPdgAntiLambda, genie::kPdgK0, genie::kPdgLambda, genie::kPdgP33m1232_Delta0, genie::kPdgP33m1232_DeltaM, genie::kPdgP33m1232_DeltaP, genie::kPdgP33m1232_DeltaPP, genie::kPdgPi0, LOG, and pDEBUG.

◆ Hadronize()

bool Pythia6Hadro2019::Hadronize ( GHepRecord * event) const
privatevirtual

Implements genie::PythiaBaseHadro2019.

Definition at line 76 of file Pythia6Hadro2019.cxx.

81{
82#ifdef __GENIE_PYTHIA6_ENABLED__
83
84 LOG("Pythia6Had", pNOTICE) << "Running PYTHIA6 hadronizer";
85
86 const Interaction * interaction = event->Summary();
87 const Kinematics & kinematics = interaction->Kine();
88 double W = kinematics.W();
89
90 LOG("Pythia6Had", pNOTICE)
91 << "Fragmentation: "
92 << "q = " << fLeadingQuark << ", qq = " << fRemnantDiquark
93 << ", W = " << W;
94
95 // Hadronize
96 int ip = 0;
97 py2ent_(&ip, &fLeadingQuark, &fRemnantDiquark, &W); // hadronizer
98
99 // Get LUJETS record
100 fPythia->GetPrimaries();
101 TClonesArray * pythia_particles =
102 (TClonesArray *) fPythia->ImportParticles("All");
103
104 // Copy PYTHIA container to a new TClonesArray so as to transfer ownership
105 // of the container and of its elements to the calling method
106
107 int np = pythia_particles->GetEntries();
108 assert(np>0);
109
110 // Hadronic 4vec
111 TLorentzVector p4Had = kinematics.HadSystP4();
112
113 // Vector defining rotation from LAB to LAB' (z:= \vec{phad})
114 TVector3 unitvq = p4Had.Vect().Unit();
115
116 // Boost velocity LAB' -> HCM
117 TVector3 beta(0,0,p4Had.P()/p4Had.Energy());
118
119 // Check target and decide appropriate status code for f/s particles
120 bool is_nucleus = interaction->InitState().Tgt().IsNucleus();
121 GHepStatus_t istfin = is_nucleus ? kIStHadronInTheNucleus : kIStStableFinalState ;
122
123 // Get the index of the mother of the hadronic system
124 int mom = event->FinalStateHadronicSystemPosition();
125 assert(mom!=-1);
126
127 // Get the neutrino vertex position (all hadrons positions set to this point)
128 GHepParticle * neutrino = event->Probe();
129 const TLorentzVector & vtx = *(neutrino->X4());
130
131 // Loop over PYTHIA6 event particles and copy relevant entries
132 unsigned int i = 0;
133 TMCParticle * p = 0;
134 TIter particle_iter(pythia_particles);
135 while( (p = (TMCParticle *) particle_iter.Next()) ) {
136
137 int particle_pdg_code = p->GetKF();
138 int pythia_particle_status = p->GetKS();
139
140 // Sanity check
141 if(pythia_particle_status == 1) {
142 if( pdg::IsQuark (particle_pdg_code) ||
143 pdg::IsDiQuark(particle_pdg_code) )
144 {
145 LOG("Pythia6Had", pERROR)
146 << "Hadronization failed! Bare quarks appear in final state!";
147 return false;
148 }
149 }
150
151 // The fragmentation products are generated in the hadronic CM frame
152 // where the z>0 axis is the \vec{phad} direction. For each particle
153 // returned by the hadronizer:
154 // - boost it back to LAB' frame {z:=\vec{phad}} / doesn't affect pT
155 // - rotate its 3-momentum from LAB' to LAB
156 TLorentzVector p4o(p->GetPx(), p->GetPy(), p->GetPz(), p->GetEnergy());
157 p4o.Boost(beta);
158 TVector3 p3 = p4o.Vect();
159 p3.RotateUz(unitvq);
160 TLorentzVector p4(p3,p4o.Energy());
161
162 // Set the proper GENIE status according to a number of things:
163 // interaction on a nucleus or nucleon, particle type
164 GHepStatus_t ist = (pythia_particle_status == 1) ?
165 istfin : kIStDISPreFragmHadronicState;
166 // Handle gammas, and leptons that might come from internal pythia decays
167 // mark them as final state particles
168 bool is_gamma = (particle_pdg_code == kPdgGamma);
169 bool is_nu = pdg::IsNeutralLepton(particle_pdg_code);
170 bool is_lchg = pdg::IsChargedLepton(particle_pdg_code);
171 bool not_hadr = is_gamma || is_nu || is_lchg;
172 if(not_hadr) { ist = kIStStableFinalState; }
173
174 // Set mother/daugher indices
175 int mother1 = mom + p->GetParent();
176 int mother2 = -1;
177 int daughter1 = (p->GetFirstChild() <= 0 ) ? -1 : mom + p->GetFirstChild();
178 int daughter2 = (p->GetLastChild() <= 0 ) ? -1 : mom + p->GetLastChild();
179
180 // Create GHepParticle
181 GHepParticle particle = GHepParticle(
182 particle_pdg_code, // pdg
183 ist, // status
184 mother1, // first parent
185 mother2, // second parent
186 daughter1, // first daughter
187 daughter2, // second daughter
188 p4.Px(), // px
189 p4.Py(), // py
190 p4.Pz(), // pz
191 p4.Energy(), // e
192 vtx.X(), // x
193 vtx.Y(), // y
194 vtx.Z(), // z
195 vtx.T() // t
196 );
197
198 LOG("Pythia6Had", pDEBUG)
199 << "Adding final state particle pdgc = " << particle.Pdg()
200 << " with status = " << particle.Status();
201
202 // Insert the particle in the list
203 event->AddParticle(particle);
204 }
205 return true;
206
207#else
208 return false;
209#endif // __GENIE_PYTHIA6_ENABLED__
210}
pNOTICE
#define pNOTICE
Definition Messenger.h:61
pERROR
#define pERROR
Definition Messenger.h:59
py2ent_
void py2ent_(int *, int *, int *, double *)
genie::GHepParticle::Pdg
int Pdg(void) const
Definition GHepParticle.h:63
genie::GHepParticle::X4
const TLorentzVector * X4(void) const
Definition GHepParticle.h:79
genie::GHepParticle::Status
GHepStatus_t Status(void) const
Definition GHepParticle.h:64
genie::InitialState::Tgt
const Target & Tgt(void) const
Definition InitialState.h:66
genie::Interaction::Kine
const Kinematics & Kine(void) const
Definition Interaction.h:71
genie::Interaction::InitState
const InitialState & InitState(void) const
Definition Interaction.h:69
genie::Kinematics::HadSystP4
const TLorentzVector & HadSystP4(void) const
Definition Kinematics.h:66
genie::Kinematics::W
double W(bool selected=false) const
Definition Kinematics.cxx:157
genie::Target::IsNucleus
bool IsNucleus(void) const
Definition Target.cxx:272
genie::pdg::IsQuark
bool IsQuark(int pdgc)
Definition PDGUtils.cxx:250
genie::pdg::IsChargedLepton
bool IsChargedLepton(int pdgc)
Definition PDGUtils.cxx:101
genie::pdg::IsNeutralLepton
bool IsNeutralLepton(int pdgc)
Definition PDGUtils.cxx:95
genie::pdg::IsDiQuark
bool IsDiQuark(int pdgc)
Definition PDGUtils.cxx:231
genie::utils::kinematics::W
double W(const Interaction *const i)
Definition KineUtils.cxx:1101
genie::kIStHadronInTheNucleus
@ kIStHadronInTheNucleus
Definition GHepStatus.h:37
genie::kIStDISPreFragmHadronicState
@ kIStDISPreFragmHadronicState
Definition GHepStatus.h:35
genie::kIStStableFinalState
@ kIStStableFinalState
Definition GHepStatus.h:30
genie::GHepStatus_t
enum genie::EGHepStatus GHepStatus_t
genie::kPdgGamma
const int kPdgGamma
Definition PDGCodes.h:189

References genie::PythiaBaseHadro2019::fLeadingQuark, genie::PythiaBaseHadro2019::fRemnantDiquark, genie::Kinematics::HadSystP4(), genie::Interaction::InitState(), genie::pdg::IsChargedLepton(), genie::pdg::IsDiQuark(), genie::pdg::IsNeutralLepton(), genie::Target::IsNucleus(), genie::pdg::IsQuark(), genie::Interaction::Kine(), genie::kIStDISPreFragmHadronicState, genie::kIStHadronInTheNucleus, genie::kIStStableFinalState, genie::kPdgGamma, LOG, pDEBUG, genie::GHepParticle::Pdg(), pERROR, pNOTICE, py2ent_(), genie::GHepParticle::Status(), genie::InitialState::Tgt(), genie::Kinematics::W(), and genie::GHepParticle::X4().

◆ Initialize()

void Pythia6Hadro2019::Initialize ( void )
privatevirtual

Reimplemented from genie::PythiaBaseHadro2019.

Definition at line 326 of file Pythia6Hadro2019.cxx.

327{
328 PythiaBaseHadro2019::Initialize();
329#ifdef __GENIE_PYTHIA6_ENABLED__
330 fPythia = TPythia6::Instance();
331 // sync GENIE/PYTHIA6 seed number
332 // PYTHIA6 is a singleton, so do this from RandomGen for all
333 // GENIE algorithms that use PYTHIA6
334 RandomGen::Instance();
335#endif
336}
genie::RandomGen::Instance
static RandomGen * Instance()
Access instance.
Definition RandomGen.cxx:74

References genie::PythiaBaseHadro2019::Initialize(), and genie::RandomGen::Instance().

Referenced by Pythia6Hadro2019(), and Pythia6Hadro2019().

◆ LoadConfig()

void Pythia6Hadro2019::LoadConfig ( void )
privatevirtual

Reimplemented from genie::PythiaBaseHadro2019.

Definition at line 306 of file Pythia6Hadro2019.cxx.

307{
308 PythiaBaseHadro2019::LoadConfig();
309
310#ifdef __GENIE_PYTHIA6_ENABLED__
311 fPythia->SetPARJ(2, fSSBarSuppression );
312 fPythia->SetPARJ(21, fGaussianPt2 );
313 fPythia->SetPARJ(23, fNonGaussianPt2Tail );
314 fPythia->SetPARJ(33, fRemainingECutoff );
315 fPythia->SetPARJ(1, fDiQuarkSuppression );
316 fPythia->SetPARJ(11, fLightVMesonSuppression );
317 fPythia->SetPARJ(12, fSVMesonSuppression );
318 fPythia->SetPARJ(41, fLunda );
319 fPythia->SetPARJ(42, fLundb );
320 fPythia->SetPARJ(45, fLundaDiq );
321#endif
322
323 LOG("Pythia6Had", pDEBUG) << this->GetConfig() ;
324}
genie::Algorithm::GetConfig
virtual const Registry & GetConfig(void) const
Definition Algorithm.cxx:246
genie::PythiaBaseHadro2019::fLundb
double fLundb
Lund b parameter.
Definition PythiaBaseHadro2019.h:63
genie::PythiaBaseHadro2019::fSSBarSuppression
double fSSBarSuppression
ssbar suppression
Definition PythiaBaseHadro2019.h:55
genie::PythiaBaseHadro2019::fNonGaussianPt2Tail
double fNonGaussianPt2Tail
non gaussian pt2 tail parameterization
Definition PythiaBaseHadro2019.h:57
genie::PythiaBaseHadro2019::fDiQuarkSuppression
double fDiQuarkSuppression
di-quark suppression parameter
Definition PythiaBaseHadro2019.h:59
genie::PythiaBaseHadro2019::fLundaDiq
double fLundaDiq
adjustment of Lund a for di-quark
Definition PythiaBaseHadro2019.h:64
genie::PythiaBaseHadro2019::fGaussianPt2
double fGaussianPt2
gaussian pt2 distribution width
Definition PythiaBaseHadro2019.h:56
genie::PythiaBaseHadro2019::fRemainingECutoff
double fRemainingECutoff
remaining E cutoff stopping fragmentation
Definition PythiaBaseHadro2019.h:58
genie::PythiaBaseHadro2019::fLunda
double fLunda
Lund a parameter.
Definition PythiaBaseHadro2019.h:62
genie::PythiaBaseHadro2019::fSVMesonSuppression
double fSVMesonSuppression
strange vector meson suppression
Definition PythiaBaseHadro2019.h:61
genie::PythiaBaseHadro2019::fLightVMesonSuppression
double fLightVMesonSuppression
light vector meson suppression
Definition PythiaBaseHadro2019.h:60

References genie::PythiaBaseHadro2019::fDiQuarkSuppression, genie::PythiaBaseHadro2019::fGaussianPt2, genie::PythiaBaseHadro2019::fLightVMesonSuppression, genie::PythiaBaseHadro2019::fLunda, genie::PythiaBaseHadro2019::fLundaDiq, genie::PythiaBaseHadro2019::fLundb, genie::PythiaBaseHadro2019::fNonGaussianPt2Tail, genie::PythiaBaseHadro2019::fRemainingECutoff, genie::PythiaBaseHadro2019::fSSBarSuppression, genie::PythiaBaseHadro2019::fSVMesonSuppression, genie::Algorithm::GetConfig(), genie::PythiaBaseHadro2019::LoadConfig(), LOG, and pDEBUG.

Referenced by Configure(), and Configure().

◆ ProcessEventRecord()

void Pythia6Hadro2019::ProcessEventRecord ( GHepRecord * event) const
virtual

Reimplemented from genie::PythiaBaseHadro2019.

Definition at line 60 of file Pythia6Hadro2019.cxx.

65{
66#ifdef __GENIE_PYTHIA6_ENABLED__
67 PythiaBaseHadro2019::ProcessEventRecord(event);
68#else
69 LOG("Pythia6Had", pFATAL)
70 << "Calling GENIE/PYTHIA6 hadronization modules without enabling PYTHIA6";
71 gAbortingInErr = true;
72 std::exit(1);
73#endif
74}
pFATAL
#define pFATAL
Definition Messenger.h:56
genie::PythiaBaseHadro2019::ProcessEventRecord
virtual void ProcessEventRecord(GHepRecord *event) const
Definition PythiaBaseHadro2019.cxx:55
genie::gAbortingInErr
bool gAbortingInErr
Definition Messenger.cxx:34

References genie::gAbortingInErr, LOG, pFATAL, and genie::PythiaBaseHadro2019::ProcessEventRecord().

◆ RestoreOriginalDecayFlags()

void Pythia6Hadro2019::RestoreOriginalDecayFlags ( void ) const
privatevirtual

Implements genie::PythiaBaseHadro2019.

Definition at line 268 of file Pythia6Hadro2019.cxx.

269{
270#ifdef __GENIE_PYTHIA6_ENABLED__
271
272fPythia->SetMDCY(fPythia->Pycomp(kPdgPi0),
273 1, (fOriDecayFlag_pi0) ? 1 : 0);
274fPythia->SetMDCY(fPythia->Pycomp(kPdgK0),
275 1, (fOriDecayFlag_K0) ? 1 : 0);
276fPythia->SetMDCY(fPythia->Pycomp(kPdgAntiK0),
277 1, (fOriDecayFlag_K0b) ? 1 : 0);
278fPythia->SetMDCY(fPythia->Pycomp(kPdgLambda),
279 1, (fOriDecayFlag_L0) ? 1 : 0);
280fPythia->SetMDCY(fPythia->Pycomp(kPdgAntiLambda),
281 1, (fOriDecayFlag_L0b) ? 1 : 0);
282fPythia->SetMDCY(fPythia->Pycomp(kPdgP33m1232_DeltaM),
283 1, (fOriDecayFlag_Dm) ? 1 : 0);
284fPythia->SetMDCY(fPythia->Pycomp(kPdgP33m1232_Delta0),
285 1, (fOriDecayFlag_D0) ? 1 : 0);
286fPythia->SetMDCY(fPythia->Pycomp(kPdgP33m1232_DeltaP),
287 1, (fOriDecayFlag_Dp) ? 1 : 0);
288fPythia->SetMDCY(fPythia->Pycomp(kPdgP33m1232_DeltaPP),
289 1, (fOriDecayFlag_Dpp) ? 1 : 0);
290
291#endif // __GENIE_PYTHIA6_ENABLED__
292}

References genie::PythiaBaseHadro2019::fOriDecayFlag_D0, genie::PythiaBaseHadro2019::fOriDecayFlag_Dm, genie::PythiaBaseHadro2019::fOriDecayFlag_Dp, genie::PythiaBaseHadro2019::fOriDecayFlag_Dpp, genie::PythiaBaseHadro2019::fOriDecayFlag_K0, genie::PythiaBaseHadro2019::fOriDecayFlag_K0b, genie::PythiaBaseHadro2019::fOriDecayFlag_L0, genie::PythiaBaseHadro2019::fOriDecayFlag_L0b, genie::PythiaBaseHadro2019::fOriDecayFlag_pi0, genie::kPdgAntiK0, genie::kPdgAntiLambda, genie::kPdgK0, genie::kPdgLambda, genie::kPdgP33m1232_Delta0, genie::kPdgP33m1232_DeltaM, genie::kPdgP33m1232_DeltaP, genie::kPdgP33m1232_DeltaPP, and genie::kPdgPi0.

◆ SetDesiredDecayFlags()

void Pythia6Hadro2019::SetDesiredDecayFlags ( void ) const
privatevirtual

Implements genie::PythiaBaseHadro2019.

Definition at line 242 of file Pythia6Hadro2019.cxx.

243{
244#ifdef __GENIE_PYTHIA6_ENABLED__
245
246 fPythia->SetMDCY(fPythia->Pycomp(kPdgPi0),
247 1, (fReqDecayFlag_pi0) ? 1 : 0);
248 fPythia->SetMDCY(fPythia->Pycomp(kPdgK0),
249 1, (fReqDecayFlag_K0) ? 1 : 0);
250 fPythia->SetMDCY(fPythia->Pycomp(kPdgAntiK0),
251 1, (fReqDecayFlag_K0b) ? 1 : 0);
252 fPythia->SetMDCY(fPythia->Pycomp(kPdgLambda),
253 1, (fReqDecayFlag_L0) ? 1 : 0);
254 fPythia->SetMDCY(fPythia->Pycomp(kPdgAntiLambda),
255 1, (fReqDecayFlag_L0b) ? 1 : 0);
256 fPythia->SetMDCY(fPythia->Pycomp(kPdgP33m1232_DeltaM),
257 1, (fReqDecayFlag_Dm) ? 1 : 0);
258 fPythia->SetMDCY(fPythia->Pycomp(kPdgP33m1232_Delta0),
259 1, (fReqDecayFlag_D0) ? 1 : 0);
260 fPythia->SetMDCY(fPythia->Pycomp(kPdgP33m1232_DeltaP),
261 1, (fReqDecayFlag_Dp) ? 1 : 0);
262 fPythia->SetMDCY(fPythia->Pycomp(kPdgP33m1232_DeltaPP),
263 1, (fReqDecayFlag_Dpp) ? 1 : 0);
264
265#endif // __GENIE_PYTHIA6_ENABLED__
266}

References genie::PythiaBaseHadro2019::fReqDecayFlag_D0, genie::PythiaBaseHadro2019::fReqDecayFlag_Dm, genie::PythiaBaseHadro2019::fReqDecayFlag_Dp, genie::PythiaBaseHadro2019::fReqDecayFlag_Dpp, genie::PythiaBaseHadro2019::fReqDecayFlag_K0, genie::PythiaBaseHadro2019::fReqDecayFlag_K0b, genie::PythiaBaseHadro2019::fReqDecayFlag_L0, genie::PythiaBaseHadro2019::fReqDecayFlag_L0b, genie::PythiaBaseHadro2019::fReqDecayFlag_pi0, genie::kPdgAntiK0, genie::kPdgAntiLambda, genie::kPdgK0, genie::kPdgLambda, genie::kPdgP33m1232_Delta0, genie::kPdgP33m1232_DeltaM, genie::kPdgP33m1232_DeltaP, genie::kPdgP33m1232_DeltaPP, and genie::kPdgPi0.

The documentation for this class was generated from the following files:
Generated on for GENIEGenerator by doxygen 1.14.0