genie::Pythia8Hadro2019 Class Reference

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

#include <Pythia8Hadro2019.h>

Inheritance diagram for genie::Pythia8Hadro2019:

Collaboration diagram for genie::Pythia8Hadro2019:

Public Member Functions
	Pythia8Hadro2019 ()
	Pythia8Hadro2019 (string config)
virtual	~Pythia8Hadro2019 ()
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. More...
virtual AlgStatus_t	GetStatus (void) const
	Get algorithm status. More...
virtual bool	AllowReconfig (void) const
virtual AlgCmp_t	Compare (const Algorithm *alg) const
	Compare with input algorithm. More...
virtual void	SetId (const AlgId &id)
	Set algorithm ID. More...
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. More...

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. More...
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. More...
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. More...
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 More...
int	AddLowRegistry (Registry *rp, bool owns=true)
	add registry with lowest priority, also update ownership More...
int	MergeTopRegistry (const Registry &r)
int	AddTopRegisties (const vector< Registry * > &rs, bool owns=false)
	Add registries with top priority, also udated Ownerships. More...
Protected Attributes inherited from genie::PythiaBaseHadro2019
int	fLeadingQuark
int	fRemnantDiquark
double	fSSBarSuppression
	ssbar suppression More...
double	fGaussianPt2
	gaussian pt2 distribution width More...
double	fNonGaussianPt2Tail
	non gaussian pt2 tail parameterization More...
double	fRemainingECutoff
	remaining E cutoff stopping fragmentation More...
double	fDiQuarkSuppression
	di-quark suppression parameter More...
double	fLightVMesonSuppression
	light vector meson suppression More...
double	fSVMesonSuppression
	strange vector meson suppression More...
double	fLunda
	Lund a parameter. More...
double	fLundb
	Lund b parameter. More...
double	fLundaDiq
	adjustment of Lund a for di-quark More...
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,...) More...
AlgId	fID
	algorithm name and configuration set More...
vector< Registry * >	fConfVect
vector< bool >	fOwnerships
	ownership for every registry in fConfVect More...
AlgStatus_t	fStatus
	algorithm execution status More...
AlgMap *	fOwnedSubAlgMp
	local pool for owned sub-algs (taken out of the factory pool) More...

Detailed Description

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

Author

Costas Andreopoulos <c.andreopoulos cern.ch> University of Liverpool

Shivesh Mandalia s.p.m.nosp@m.anda.nosp@m.lia@q.nosp@m.mul..nosp@m.ac.uk Queen Mary University of London

Created:

October 17, 2019

License:

Copyright (c) 2003-2024, The GENIE Collaboration For the full text of the license visit http://copyright.genie-mc.org

Definition at line 37 of file Pythia8Hadro2019.h.

Constructor & Destructor Documentation

Pythia8Hadro2019::Pythia8Hadro2019

(

)

Definition at line 44 of file Pythia8Hadro2019.cxx.

References Initialize().

                                   :
PythiaBaseHadro2019("genie::Pythia8Hadro2019")
{
  this->Initialize();
}

Pythia8Hadro2019::Pythia8Hadro2019

(

string

config

)

Definition at line 50 of file Pythia8Hadro2019.cxx.

References Initialize().

                                                :
PythiaBaseHadro2019("genie::Pythia8Hadro2019", config)
{
  this->Initialize();
}

Pythia8Hadro2019::~Pythia8Hadro2019 ( )

virtual

Definition at line 56 of file Pythia8Hadro2019.cxx.

{
#ifdef __GENIE_PYTHIA8_ENABLED__
  delete fPythia;
#endif
}

Member Function Documentation

void Pythia8Hadro2019::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 320 of file Pythia8Hadro2019.cxx.

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

{
  Algorithm::Configure(config);
  this->LoadConfig();
}

void Pythia8Hadro2019::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 326 of file Pythia8Hadro2019.cxx.

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

{
  Algorithm::Configure(config);
  this->LoadConfig();
}

void Pythia8Hadro2019::CopyOriginalDecayFlags ( void  ) const

privatevirtual

Implements genie::PythiaBaseHadro2019.

Definition at line 260 of file Pythia8Hadro2019.cxx.

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.

{
#ifdef __GENIE_PYTHIA8_ENABLED__
  fOriDecayFlag_pi0 = fPythia->particleData.canDecay(kPdgPi0);
  fOriDecayFlag_K0  = fPythia->particleData.canDecay(kPdgK0);
  fOriDecayFlag_K0b = fPythia->particleData.canDecay(kPdgAntiK0);
  fOriDecayFlag_L0  = fPythia->particleData.canDecay(kPdgLambda);
  fOriDecayFlag_L0b = fPythia->particleData.canDecay(kPdgAntiLambda);
  fOriDecayFlag_Dm  = fPythia->particleData.canDecay(kPdgP33m1232_DeltaM);
  fOriDecayFlag_D0  = fPythia->particleData.canDecay(kPdgP33m1232_Delta0);
  fOriDecayFlag_Dp  = fPythia->particleData.canDecay(kPdgP33m1232_DeltaP);
  fOriDecayFlag_Dpp = fPythia->particleData.canDecay(kPdgP33m1232_DeltaPP);

#ifdef __GENIE_LOW_LEVEL_MESG_ENABLED__
    LOG("Pythia8Had", pDEBUG)
       << "Original PYTHIA6 decay flags:"
       << "\n  pi0           =  " << fOriDecayFlag_pi0
       << "\n  K0            =  " << fOriDecayFlag_K0
       << "\n  \bar{K0}      =  " << fOriDecayFlag_K0b
       << "\n  Lambda        =  " << fOriDecayFlag_L0
       << "\n  \bar{Lambda0} =  " << fOriDecayFlag_L0b
       << "\n  D-            =  " << fOriDecayFlag_Dm
       << "\n  D0            =  " << fOriDecayFlag_D0
       << "\n  D+            =  " << fOriDecayFlag_Dp
       << "\n  D++           =  " << fOriDecayFlag_Dpp;
#endif

#endif
}

bool Pythia8Hadro2019::Hadronize ( GHepRecord *  event ) const

privatevirtual

Implements genie::PythiaBaseHadro2019.

Definition at line 79 of file Pythia8Hadro2019.cxx.

References e, 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, pINFO, pNOTICE, genie::GHepParticle::Status(), genie::InitialState::Tgt(), genie::Kinematics::W(), genie::utils::kinematics::W(), and genie::GHepParticle::X4().

{
#ifdef __GENIE_PYTHIA8_ENABLED__
  LOG("Pythia8Had", pNOTICE) << "Running PYTHIA8 hadronizer";

  const Interaction * interaction = event->Summary();
  const Kinematics & kinematics = interaction->Kine();
  double W = kinematics.W();

  LOG("Pythia8Had", pNOTICE)
    << "Fragmentation: "
    << "q = " << fLeadingQuark << ", qq = " << fRemnantDiquark
    << ", W = " << W << " GeV";

  // Hadronize

  LOG("Pythia8Had", pDEBUG) << "Reseting PYTHIA8 event";
  fPythia->event.reset();

  // Get quark/diquark masses
  double mA = fPythia->particleData.m0(fLeadingQuark);
  double mB = fPythia->particleData.m0(fRemnantDiquark);

  LOG("Pythia8Had", pINFO)
    << "Leading quark mass = " << mA
    << " GeV, remnant diqurak mass = " << mB << ", GeV";

  // Calculate quark/diquark energy/momentum
  double pzAcm = 0.5 * Pythia8::sqrtpos(
     (W + mA + mB) * (W - mA - mB) * (W - mA + mB) * (W + mA - mB) ) / W;
  double pzBcm = -pzAcm;
  double eA    = sqrt(mA*mA + pzAcm*pzAcm);
  double eB    = sqrt(mB*mB + pzBcm*pzBcm);

  LOG("Pythia8Had", pINFO)
   << "Quark: (pz = " << pzAcm << ", E = " << eA << ") GeV, "
   << "Diquark: (pz = " << pzBcm << ", E = " << eB << ") GeV";

  // Pythia8 status code for outgoing particles of the hardest subprocesses is 23
  // anti/colour tags for these 2 particles must complement each other
  LOG("Pythia8Had", pDEBUG) << "Appending quark/diquark into the PYTHIA8 event";
  fPythia->event.append(fLeadingQuark,   23, 101, 0, 0., 0., pzAcm, eA, mA);
  fPythia->event.append(fRemnantDiquark, 23, 0, 101, 0., 0., pzBcm, eB, mB);
  fPythia->event.list();

  LOG("Pythia8Had", pDEBUG) << "Generating next PYTHIA8 event";
  fPythia->next();

  // List the event information
  fPythia->event.list();
  fPythia->stat();

  // Get LUJETS record
  LOG("Pythia8Had", pDEBUG) << "Copying PYTHIA8 event record into GENIE's";
  Pythia8::Event &fEvent = fPythia->event;
  int np = fEvent.size();
  assert(np>0);

  TClonesArray * particle_list = new TClonesArray("genie::GHepParticle", np);
  particle_list->SetOwner(true);

  // Hadronic 4vec
  TLorentzVector p4Had = kinematics.HadSystP4();

  // Vector defining rotation from LAB to LAB' (z:= \vec{phad})
  TVector3 unitvq = p4Had.Vect().Unit();

  // Boost velocity LAB' -> HCM
  TVector3 beta(0,0,p4Had.P()/p4Had.Energy());

  // Check target and decide appropriate status code for f/s particles
  bool is_nucleus = interaction->InitState().Tgt().IsNucleus();
  GHepStatus_t istfin = is_nucleus ? kIStHadronInTheNucleus : kIStStableFinalState ;

  // Get the index of the mother of the hadronic system
  int mom = event->FinalStateHadronicSystemPosition();
  assert(mom!=-1);

  // Get the neutrino vertex position (all hadrons positions set to this point)
  GHepParticle * neutrino  = event->Probe();
  const TLorentzVector & vtx = *(neutrino->X4());

  // Loop over PYTHIA8 event particles and copy relevant entries
  for (int i = 0; i < np; i++) {

     if (fEvent[i].id() == 90) continue; // ignore (system) pseudoparticle

     // Get PYTHIA8 particle ID and status codes
     int particle_pdg_code      = fEvent[i].id();
     int pythia_particle_status = fEvent[i].status();

     // Sanity check
     if(pythia_particle_status > 0) {
       if( pdg::IsQuark  (particle_pdg_code) ||
           pdg::IsDiQuark(particle_pdg_code) )
       {
          LOG("Pythia8Had", pERROR)
             << "Hadronization failed! Bare quarks appear in final state!";
          return false;
       }
     }

     // Copy the initial q, qq (status = -23) and all undecayed particles
     // Ignore particles we asked PYTHIA to decay (eg omega, Delta) but
     // record their decay products
     bool copy = (pythia_particle_status==-23) || (pythia_particle_status > 0);
     if(copy) {
        // The fragmentation products are generated in the hadronic CM frame
        // where the z>0 axis is the \vec{phad} direction. For each particle
        // returned by the hadronizer:
        // - boost it back to LAB' frame {z:=\vec{phad}} / doesn't affect pT
        // - rotate its 3-momentum from LAB' to LAB
        TLorentzVector p4o(
          fEvent[i].px(), fEvent[i].py(), fEvent[i].pz(), fEvent[i].e());
        p4o.Boost(beta);
        TVector3 p3 = p4o.Vect();
        p3.RotateUz(unitvq);
        TLorentzVector p4(p3,p4o.Energy());

        // Set the proper GENIE status according to a number of things:
        // interaction on a nucleus or nucleon, particle type
        GHepStatus_t ist = (pythia_particle_status > 0) ?
           istfin : kIStDISPreFragmHadronicState;
        // Handle gammas, and leptons that might come from internal pythia decays
        // mark them as final state particles
        bool is_gamma = (particle_pdg_code == kPdgGamma);
        bool is_nu    = pdg::IsNeutralLepton(particle_pdg_code);
        bool is_lchg  = pdg::IsChargedLepton(particle_pdg_code);
        bool not_hadr = is_gamma || is_nu || is_lchg;
        if(not_hadr)  { ist = kIStStableFinalState; }

        // Set mother/daugher indices
        // int mother1 = mom+ fEvent[i].mother1();
        // int mother2 = (pythia_particle_status > 0) ? mom + fEvent[i].mother2() : -1;
        int mother1 = mom; // fEvent[i].mother1();
        int mother2 = -1; //(pythia_particle_status > 0) ? mom + fEvent[i].mother2() : -1;
        if(pythia_particle_status > 0) {
          mother1 = mom+1;
          mother2 = mom+2;
        }
        int daughter1 = -1;//(fEvent[i].daughter1() <= 0 ) ? -1 : mom  + fEvent[i].daughter1();
        int daughter2 = -1;//(fEvent[i].daughter1() <= 0 ) ? -1 : mom  + fEvent[i].daughter2();

        // Create GHepParticle
        GHepParticle particle = GHepParticle(
            particle_pdg_code, // pdg
            ist,               // status
            mother1,           // first parent
            mother2,           // second parent
            daughter1,         // first daughter
            daughter2,         // second daughter
            p4.Px(),           // px
            p4.Py(),           // py
            p4.Pz(),           // pz
            p4.Energy(),       // e
            vtx.X(),           // x
            vtx.Y(),           // y
            vtx.Z(),           // z
            vtx.T()            // t
        );

        LOG("Pythia8Had", pDEBUG)
             << "Adding final state particle pdgc = " << particle.Pdg()
             << " with status = " << particle.Status();

        // Insert the particle in the list
        event->AddParticle(particle);
     }// copy?
  }// loop over particles

  return true;

#else
  return false;
#endif
}

void Pythia8Hadro2019::Initialize ( void  )

privatevirtual

Reimplemented from genie::PythiaBaseHadro2019.

Definition at line 355 of file Pythia8Hadro2019.cxx.

References genie::RandomGen::GetSeed(), genie::RandomGen::Instance(), LOG, and pINFO.

Referenced by Pythia8Hadro2019().

{
#ifdef __GENIE_PYTHIA8_ENABLED__
  fPythia = new Pythia8::Pythia();

  fPythia->readString("ProcessLevel:all = off");
  fPythia->readString("Print:quiet      = on");

  // sync GENIE and PYTHIA8 seeds
  RandomGen * rnd = RandomGen::Instance();
  long int seed = rnd->GetSeed();
  fPythia->readString("Random:setSeed = on");
  fPythia->settings.mode("Random:seed", seed);
  LOG("Pythia8Had", pINFO)
    << "PYTHIA8  seed = " << fPythia->settings.mode("Random:seed");

  fPythia->init();

#endif
}

void Pythia8Hadro2019::LoadConfig ( void  )

privatevirtual

Reimplemented from genie::PythiaBaseHadro2019.

Definition at line 332 of file Pythia8Hadro2019.cxx.

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

{
  PythiaBaseHadro2019::LoadConfig();

#ifdef __GENIE_PYTHIA8_ENABLED__
  fPythia->settings.parm("StringFlav:probStoUD",         fSSBarSuppression);
  fPythia->settings.parm("Diffraction:primKTwidth",      fGaussianPt2);
  fPythia->settings.parm("StringPT:enhancedFraction",    fNonGaussianPt2Tail);
  fPythia->settings.parm("StringFragmentation:stopMass", fRemainingECutoff);
  fPythia->settings.parm("StringFlav:probQQtoQ",         fDiQuarkSuppression);
  fPythia->settings.parm("StringFlav:mesonUDvector",     fLightVMesonSuppression);
  fPythia->settings.parm("StringFlav:mesonSvector",      fSVMesonSuppression);
  fPythia->settings.parm("StringZ:aLund",                fLunda);
  fPythia->settings.parm("StringZ:bLund",                fLundb);
  fPythia->settings.parm("StringZ:aExtraDiquark",        fLundaDiq);

  fPythia->init(); // needed again to read the above?

#endif

  LOG("Pythia8Had", pDEBUG) << this->GetConfig();
}

void Pythia8Hadro2019::ProcessEventRecord ( GHepRecord *  event ) const

virtual

Reimplemented from genie::PythiaBaseHadro2019.

Definition at line 63 of file Pythia8Hadro2019.cxx.

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

{
#ifdef __GENIE_PYTHIA8_ENABLED__
  PythiaBaseHadro2019::ProcessEventRecord(event);
#else
  LOG("Pythia8Had", pFATAL)
    << "Calling GENIE/PYTHIA8 hadronization modules without enabling PYTHIA8";
  gAbortingInErr = true;
  std::exit(1);
#endif
}

void Pythia8Hadro2019::RestoreOriginalDecayFlags ( void  ) const

privatevirtual

Implements genie::PythiaBaseHadro2019.

Definition at line 305 of file Pythia8Hadro2019.cxx.

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.

{
#ifdef __GENIE_PYTHIA8_ENABLED__
  fPythia->particleData.mayDecay(kPdgPi0,              fOriDecayFlag_pi0 );
  fPythia->particleData.mayDecay(kPdgK0,               fOriDecayFlag_K0  );
  fPythia->particleData.mayDecay(kPdgAntiK0,           fOriDecayFlag_K0b );
  fPythia->particleData.mayDecay(kPdgLambda,           fOriDecayFlag_L0  );
  fPythia->particleData.mayDecay(kPdgAntiLambda,       fOriDecayFlag_L0b );
  fPythia->particleData.mayDecay(kPdgP33m1232_DeltaM,  fOriDecayFlag_Dm  );
  fPythia->particleData.mayDecay(kPdgP33m1232_Delta0,  fOriDecayFlag_D0  );
  fPythia->particleData.mayDecay(kPdgP33m1232_DeltaP,  fOriDecayFlag_Dp  );
  fPythia->particleData.mayDecay(kPdgP33m1232_DeltaPP, fOriDecayFlag_Dpp );
#endif
}

void Pythia8Hadro2019::SetDesiredDecayFlags ( void  ) const

privatevirtual

Implements genie::PythiaBaseHadro2019.

Definition at line 290 of file Pythia8Hadro2019.cxx.

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.

{
#ifdef __GENIE_PYTHIA8_ENABLED__
  fPythia->particleData.mayDecay(kPdgPi0,              fReqDecayFlag_pi0 );
  fPythia->particleData.mayDecay(kPdgK0,               fReqDecayFlag_K0  );
  fPythia->particleData.mayDecay(kPdgAntiK0,           fReqDecayFlag_K0b );
  fPythia->particleData.mayDecay(kPdgLambda,           fReqDecayFlag_L0  );
  fPythia->particleData.mayDecay(kPdgAntiLambda,       fReqDecayFlag_L0b );
  fPythia->particleData.mayDecay(kPdgP33m1232_DeltaM,  fReqDecayFlag_Dm  );
  fPythia->particleData.mayDecay(kPdgP33m1232_Delta0,  fReqDecayFlag_D0  );
  fPythia->particleData.mayDecay(kPdgP33m1232_DeltaP,  fReqDecayFlag_Dp  );
  fPythia->particleData.mayDecay(kPdgP33m1232_DeltaPP, fReqDecayFlag_Dpp );
#endif
}

---

The documentation for this class was generated from the following files:

• Pythia8Hadro2019.h
• Pythia8Hadro2019.cxx