Andreopoulos - Gallagher (AG) GENIE Charm Hadronization model. More...

#include <AGCharmPythia8Hadro2023.h>

Inheritance diagram for genie::AGCharmPythia8Hadro2023:

Collaboration diagram for genie::AGCharmPythia8Hadro2023:

Public Member Functions
	AGCharmPythia8Hadro2023 ()
	AGCharmPythia8Hadro2023 (string config)
virtual	~AGCharmPythia8Hadro2023 ()
Public Member Functions inherited from genie::AGCharmPythiaBaseHadro2023
void	ProcessEventRecord (GHepRecord *event) const
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
void	Initialize (void) const
bool	HadronizeRemnant (int qrkSyst1, int qrkSyst2, double WR, TLorentzVector p4R, unsigned int &rpos, TClonesArray *particle_list) const
void	LoadConfig (void)
void	CopyOriginalDecayFlags (void) const
void	SetDesiredDecayFlags (void) const
void	RestoreOriginalDecayFlags (void) const

Private Attributes
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

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::AGCharmPythiaBaseHadro2023
	AGCharmPythiaBaseHadro2023 ()
	AGCharmPythiaBaseHadro2023 (string config)
	AGCharmPythiaBaseHadro2023 (string name, string config)
virtual	~AGCharmPythiaBaseHadro2023 ()
void	LoadConfig (void)
void	Initialize (void) const
void	Configure (const Registry &config)
void	Configure (string config)
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::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

Andreopoulos - Gallagher (AG) GENIE Charm Hadronization model.

     The model relies on empirical charm fragmentation and pT functions,
     as well as on experimentally-determined charm fractions, to produce
     the ID and 4-momentum of charmed hadron in charm production events.

     The remnant (non-charm) system is hadronised by a call to PYTHIA.

     Is a concrete implementation of the EventRecordVisitorI interface.

Author: Costas Andreopoulos <constantinos.andreopoulos \at cern.ch> University of Liverpool & STFC Rutherford Appleton Laboratory

Hugh Gallagher galla.nosp@m.g@mi.nosp@m.nos.p.nosp@m.hy.t.nosp@m.ufts..nosp@m.edu Tufts University

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 43 of file AGCharmPythia8Hadro2023.h.

Constructor & Destructor Documentation

◆ AGCharmPythia8Hadro2023() [1/2]

AGCharmPythia8Hadro2023::AGCharmPythia8Hadro2023 ( )

Definition at line 45 of file AGCharmPythia8Hadro2023.cxx.

                                                 :
AGCharmPythiaBaseHadro2023("genie::AGCharmPythia8Hadro2023")
{
  this->Initialize();
}

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

◆ AGCharmPythia8Hadro2023() [2/2]

AGCharmPythia8Hadro2023::AGCharmPythia8Hadro2023 ( string config )

Definition at line 51 of file AGCharmPythia8Hadro2023.cxx.

                                                              :
AGCharmPythiaBaseHadro2023("genie::AGCharmPythia8Hadro2023", config)
{
  this->Initialize();
}

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

◆ ~AGCharmPythia8Hadro2023()

AGCharmPythia8Hadro2023::~AGCharmPythia8Hadro2023 ( )

virtual

Definition at line 57 of file AGCharmPythia8Hadro2023.cxx.

{
#ifdef __GENIE_PYTHIA8_ENABLED__
 
#endif
}

Member Function Documentation

◆ CopyOriginalDecayFlags()

void AGCharmPythia8Hadro2023::CopyOriginalDecayFlags ( void ) const

private

Definition at line 283 of file AGCharmPythia8Hadro2023.cxx.

{
#ifdef __GENIE_PYTHIA8_ENABLED__
  Pythia8::Pythia* gPythia = Pythia8Singleton::Instance()->Pythia8();
 
  fOriDecayFlag_pi0 = gPythia->particleData.canDecay(kPdgPi0);
  fOriDecayFlag_K0  = gPythia->particleData.canDecay(kPdgK0);
  fOriDecayFlag_K0b = gPythia->particleData.canDecay(kPdgAntiK0);
  fOriDecayFlag_L0  = gPythia->particleData.canDecay(kPdgLambda);
  fOriDecayFlag_L0b = gPythia->particleData.canDecay(kPdgAntiLambda);
  fOriDecayFlag_Dm  = gPythia->particleData.canDecay(kPdgP33m1232_DeltaM);
  fOriDecayFlag_D0  = gPythia->particleData.canDecay(kPdgP33m1232_Delta0);
  fOriDecayFlag_Dp  = gPythia->particleData.canDecay(kPdgP33m1232_DeltaP);
  fOriDecayFlag_Dpp = gPythia->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
}

References fOriDecayFlag_D0, fOriDecayFlag_Dm, fOriDecayFlag_Dp, fOriDecayFlag_Dpp, fOriDecayFlag_K0, fOriDecayFlag_K0b, fOriDecayFlag_L0, fOriDecayFlag_L0b, fOriDecayFlag_pi0, genie::Pythia8Singleton::Instance(), 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.

Referenced by HadronizeRemnant().

◆ HadronizeRemnant()

bool AGCharmPythia8Hadro2023::HadronizeRemnant	(	int	qrkSyst1,
		int	qrkSyst2,
		double	WR,
		TLorentzVector	p4R,
		unsigned int &	rpos,
		TClonesArray *	particle_list ) const

privatevirtual

Implements genie::AGCharmPythiaBaseHadro2023.

Definition at line 126 of file AGCharmPythia8Hadro2023.cxx.

{
#ifdef __GENIE_PYTHIA8_ENABLED__
  Pythia8::Pythia* gPythia = Pythia8Singleton::Instance()->Pythia8();
 
  CopyOriginalDecayFlags();
  SetDesiredDecayFlags();
 
  gPythia->event.reset();
 
  double m1 = gPythia->particleData.m0(qrkSyst1);
  double m2 = gPythia->particleData.m0(qrkSyst2);
 
  LOG("AGCharmPythia8Hadro2023", pINFO) // debug
    << " qrkSyst1 " << qrkSyst1 << " m1 " << m1
    << " qrkSyst2 " << qrkSyst2 << " m2 " << m2
    << " WR " << WR << " p4R " << genie::utils::print::P4AsString(&p4R);
 
  double pz1cm = 0.5 * Pythia8::sqrtpos(
    (WR + m1 + m2)*(WR - m1 -m2)*(WR - m1 + m2)*(WR + m1 - m2) ) / WR;
 
  double pz2cm = - pz1cm;
  double e1    = sqrt(m1*m1 + pz1cm*pz1cm);
  double e2    = sqrt(m2*m2 + pz2cm*pz2cm);
 
  LOG("AGCharmPythia8Hadro2023", pINFO) // debug
    << " qrkSyst1 pz=" << pz1cm << ", E=" << e1
    << " qrkSyst2 pz=" << pz2cm << ", E=" << e2;
 
  // status codes and anti/collor tags must complement each other
  // id, status, int col, int acol, px,py,pz,E,m,scale=0,pol=9
  gPythia->event.append(qrkSyst1,23,101,  0,0.,0.,pz1cm,e1,m1);
  gPythia->event.append(qrkSyst2,23,  0,101,0.,0.,pz2cm,e2,m2);
  //gPythia->event.list();
 
  LOG("AGCharmPythia8Hadro2023", pDEBUG)
    << "Generating next PYTHIA8 event";
  // Generating next pythia8 event
  gPythia->next();
 
  //gPythia->event.list();
  //gPythia->stat();
 
  // get the LUJETss record
  Pythia8::Event &fEvent = gPythia->event;
  int np = fEvent.size();
  assert(np>0);
 
  // Vector defining rotation from LAB to LAB' (z:= \vec{p4R}J)
  TVector3 unitvq = p4R.Vect().Unit();
 
  // Boost velocity LAB' -> HCM
  TVector3 beta(0,0,p4R.P()/p4R.Energy());
 
  GHepStatus_t istfin = kIStHadronInTheNucleus; /// should be testing on is_nucleus
  // from interaction->InitState().Tgt().IsNucleus();
 
  int mom = -1; // event->FinalStateHadronicSystemPosition();
  const TLorentzVector vtx(0,0,0,0);
 
  for (int i=0; i<np; ++i) {
    if (fEvent[i].id() == 90) continue;   // ignore (system) pseudoparticle
 
    int particle_pdg_code      = fEvent[i].id();
    int pythia_particle_status = fEvent[i].status();
 
    GHepParticle* bremn = 0; // boosted remnant
 
    // sanity check
    if (pythia_particle_status > 0 ) {
      if ( pdg::IsQuark  (particle_pdg_code ) ||
           pdg::IsDiQuark(particle_pdg_code )    )
        {
          LOG("AGCharmPythia8Hadro2023", pERROR)
            << "Hadronization failed! Bare quarks appear in final state!";
          return false;
        }
    }
    // copy the initial quark systemcs (satus = -23) and all undecayed particles
    // ignore particle we asked PYTHIA to decay but record their decay products
    bool copy = (pythia_particle_status==-23) ||
      (pythia_particle_status > 0);
    if (copy) {
      // the fragmenation products are generated in the remnant hadronic CM frame
      // where the z>0 asix is the \vec{p4R} direction.
      // for each particle returned by the hadronizer:
      //   -- boost it back to LAB' frame {z:\vec{p4T}} / 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());
 
      //!!!!
      istfin = GHepStatus_t(1);
 
      // set the proper GENIE status according to a number of things:
      GHepStatus_t ist = (pythia_particle_status > 0) ?
        istfin : kIStDISPreFragmHadronicState;
 
      // handle gammas, and leptons that might come from internal pythia decay
      // 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; }
 
      int mother1 = mom;
      int mother2 = -1;
 
      int daughter1 = -1;
      int daughter2 = -1;
 
      // create GHepParticle
 
      bremn = new ((*particle_list)[rpos++]) 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
                                                          );
 
      //std::cout << "bremn " << rpos-1 << "||" << *bremn << std::endl;
 
    } // if (copy)
  } // loop over particles in pythia event record
 
  RestoreOriginalDecayFlags();
 
  return true;
 
#else
  LOG("AGCharmPythia8Hadro2023", pFATAL)
    << "calling GENIE/PYTHIA8 charm hadronization without enabling PYTHIA8"
    << " qrkSyst " << qrkSyst1 << "," << qrkSyst2 << " WR " << WR;
  gAbortingInErr = true;
  std::exit(1);
#endif
 
 
}

References CopyOriginalDecayFlags(), e, genie::gAbortingInErr, genie::Pythia8Singleton::Instance(), genie::pdg::IsChargedLepton(), genie::pdg::IsDiQuark(), genie::pdg::IsNeutralLepton(), genie::pdg::IsQuark(), genie::kIStDISPreFragmHadronicState, genie::kIStHadronInTheNucleus, genie::kIStStableFinalState, genie::kPdgGamma, LOG, genie::utils::print::P4AsString(), pDEBUG, pERROR, pFATAL, pINFO, RestoreOriginalDecayFlags(), and SetDesiredDecayFlags().

◆ Initialize()

void AGCharmPythia8Hadro2023::Initialize ( void ) const

private

Definition at line 64 of file AGCharmPythia8Hadro2023.cxx.

{
  AGCharmPythiaBaseHadro2023::Initialize();
#ifdef __GENIE_PYTHIA8_ENABLED__
  Pythia8::Pythia* gPythia = Pythia8Singleton::Instance()->Pythia8();
 
  gPythia->readString("ProcessLevel:all = off");
  gPythia->readString("Print:quiet = on");
 
  // sync GENIE and PYTHIA8 seeds
  RandomGen * rnd = RandomGen::Instance();
  long int seed = rnd->GetSeed();
  gPythia->readString("Random:setSeed = on");
  gPythia->settings.mode("Random:seed", seed);
  LOG("AGCharmPythia8Hadro2023", pINFO)
    << "PYTHIA8 seed = " << gPythia->settings.mode("Random:seed");
  gPythia->init();
 
  fOriDecayFlag_pi0       = false;
  fOriDecayFlag_K0        = false;
  fOriDecayFlag_K0b       = false;
  fOriDecayFlag_L0        = false;
  fOriDecayFlag_L0b       = false;
  fOriDecayFlag_Dm        = false;
  fOriDecayFlag_D0        = false;
  fOriDecayFlag_Dp        = false;
  fOriDecayFlag_Dpp       = false;
  fReqDecayFlag_pi0       = false;
  fReqDecayFlag_K0        = false;
  fReqDecayFlag_K0b       = false;
  fReqDecayFlag_L0        = false;
  fReqDecayFlag_L0b       = false;
  fReqDecayFlag_Dm        = false;
  fReqDecayFlag_D0        = false;
  fReqDecayFlag_Dp        = false;
  fReqDecayFlag_Dpp       = false;
 
#else
  LOG("AGCharmPythia8Hadro2023", pFATAL)
    << "calling GENIE/PYTHIA8 charm hadronization without enabling PYTHIA8";
  gAbortingInErr = true;
  std::exit(1);
#endif
 
}

References fOriDecayFlag_D0, fOriDecayFlag_Dm, fOriDecayFlag_Dp, fOriDecayFlag_Dpp, fOriDecayFlag_K0, fOriDecayFlag_K0b, fOriDecayFlag_L0, fOriDecayFlag_L0b, fOriDecayFlag_pi0, fReqDecayFlag_D0, fReqDecayFlag_Dm, fReqDecayFlag_Dp, fReqDecayFlag_Dpp, fReqDecayFlag_K0, fReqDecayFlag_K0b, fReqDecayFlag_L0, fReqDecayFlag_L0b, fReqDecayFlag_pi0, genie::gAbortingInErr, genie::RandomGen::GetSeed(), genie::AGCharmPythiaBaseHadro2023::Initialize(), genie::Pythia8Singleton::Instance(), genie::RandomGen::Instance(), LOG, pFATAL, and pINFO.

Referenced by AGCharmPythia8Hadro2023(), and AGCharmPythia8Hadro2023().

◆ LoadConfig()

void AGCharmPythia8Hadro2023::LoadConfig ( void )

private

Definition at line 110 of file AGCharmPythia8Hadro2023.cxx.

{
  AGCharmPythiaBaseHadro2023::LoadConfig();
  // Set required PYTHIA decay flags
  fReqDecayFlag_pi0       = false; // don't decay pi0
  fReqDecayFlag_K0        = false; // don't decay K0
  fReqDecayFlag_K0b       = false; // don't decay \bar{K0}
  fReqDecayFlag_L0        = false; // don't decay Lambda0
  fReqDecayFlag_L0b       = false; // don't decay \bar{Lambda0}
  fReqDecayFlag_Dm        = true;  // decay Delta-
  fReqDecayFlag_D0        = true;  // decay Delta0
  fReqDecayFlag_Dp        = true;  // decay Delta+
  fReqDecayFlag_Dpp       = true;  // decay Delta++
}

References fReqDecayFlag_D0, fReqDecayFlag_Dm, fReqDecayFlag_Dp, fReqDecayFlag_Dpp, fReqDecayFlag_K0, fReqDecayFlag_K0b, fReqDecayFlag_L0, fReqDecayFlag_L0b, fReqDecayFlag_pi0, and genie::AGCharmPythiaBaseHadro2023::LoadConfig().

◆ RestoreOriginalDecayFlags()

void AGCharmPythia8Hadro2023::RestoreOriginalDecayFlags ( void ) const

private

Definition at line 332 of file AGCharmPythia8Hadro2023.cxx.

{
#ifdef __GENIE_PYTHIA8_ENABLED__
  Pythia8::Pythia* gPythia = Pythia8Singleton::Instance()->Pythia8();
 
  gPythia->particleData.mayDecay(kPdgPi0,              fOriDecayFlag_pi0 );
  gPythia->particleData.mayDecay(kPdgK0,               fOriDecayFlag_K0  );
  gPythia->particleData.mayDecay(kPdgAntiK0,           fOriDecayFlag_K0b );
  gPythia->particleData.mayDecay(kPdgLambda,           fOriDecayFlag_L0  );
  gPythia->particleData.mayDecay(kPdgAntiLambda,       fOriDecayFlag_L0b );
  gPythia->particleData.mayDecay(kPdgP33m1232_DeltaM,  fOriDecayFlag_Dm  );
  gPythia->particleData.mayDecay(kPdgP33m1232_Delta0,  fOriDecayFlag_D0  );
  gPythia->particleData.mayDecay(kPdgP33m1232_DeltaP,  fOriDecayFlag_Dp  );
  gPythia->particleData.mayDecay(kPdgP33m1232_DeltaPP, fOriDecayFlag_Dpp );
#endif
}

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

Referenced by HadronizeRemnant().

◆ SetDesiredDecayFlags()

void AGCharmPythia8Hadro2023::SetDesiredDecayFlags ( void ) const

private

Definition at line 315 of file AGCharmPythia8Hadro2023.cxx.

{
#ifdef __GENIE_PYTHIA8_ENABLED__
  Pythia8::Pythia* gPythia = Pythia8Singleton::Instance()->Pythia8();
 
  gPythia->particleData.mayDecay(kPdgPi0,              fReqDecayFlag_pi0 );
  gPythia->particleData.mayDecay(kPdgK0,               fReqDecayFlag_K0  );
  gPythia->particleData.mayDecay(kPdgAntiK0,           fReqDecayFlag_K0b );
  gPythia->particleData.mayDecay(kPdgLambda,           fReqDecayFlag_L0  );
  gPythia->particleData.mayDecay(kPdgAntiLambda,       fReqDecayFlag_L0b );
  gPythia->particleData.mayDecay(kPdgP33m1232_DeltaM,  fReqDecayFlag_Dm  );
  gPythia->particleData.mayDecay(kPdgP33m1232_Delta0,  fReqDecayFlag_D0  );
  gPythia->particleData.mayDecay(kPdgP33m1232_DeltaP,  fReqDecayFlag_Dp  );
  gPythia->particleData.mayDecay(kPdgP33m1232_DeltaPP, fReqDecayFlag_Dpp );
#endif
}