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

Private Attributes
TPythia6 *	fPythia
	PYTHIA6 wrapper class. More...

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

Created:: August 17, 2004

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

Definition at line 36 of file Pythia6Hadro2019.h.

Constructor & Destructor Documentation

Pythia6Hadro2019::Pythia6Hadro2019 ( )

Definition at line 43 of file Pythia6Hadro2019.cxx.

References Initialize().

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

Pythia6Hadro2019::Pythia6Hadro2019 ( string config )

Definition at line 49 of file Pythia6Hadro2019.cxx.

References Initialize().

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

Pythia6Hadro2019::~Pythia6Hadro2019 ( )

virtual

Definition at line 55 of file Pythia6Hadro2019.cxx.

 {
 
 }

Member Function Documentation

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 296 of file Pythia6Hadro2019.cxx.

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

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

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 302 of file Pythia6Hadro2019.cxx.

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

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

void Pythia6Hadro2019::CopyOriginalDecayFlags ( void ) const

privatevirtual

Implements genie::PythiaBaseHadro2019.

Definition at line 214 of file Pythia6Hadro2019.cxx.

 {
 #ifdef __GENIE_PYTHIA6_ENABLED__
   fOriDecayFlag_pi0 = (fPythia->GetMDCY(fPythia->Pycomp(kPdgPi0),              1) == 1);
   fOriDecayFlag_K0  = (fPythia->GetMDCY(fPythia->Pycomp(kPdgK0),               1) == 1);
   fOriDecayFlag_K0b = (fPythia->GetMDCY(fPythia->Pycomp(kPdgAntiK0),           1) == 1);
   fOriDecayFlag_L0  = (fPythia->GetMDCY(fPythia->Pycomp(kPdgLambda),           1) == 1);
   fOriDecayFlag_L0b = (fPythia->GetMDCY(fPythia->Pycomp(kPdgAntiLambda),       1) == 1);
   fOriDecayFlag_Dm  = (fPythia->GetMDCY(fPythia->Pycomp(kPdgP33m1232_DeltaM),  1) == 1);
   fOriDecayFlag_D0  = (fPythia->GetMDCY(fPythia->Pycomp(kPdgP33m1232_Delta0),  1) == 1);
   fOriDecayFlag_Dp  = (fPythia->GetMDCY(fPythia->Pycomp(kPdgP33m1232_DeltaP),  1) == 1);
   fOriDecayFlag_Dpp = (fPythia->GetMDCY(fPythia->Pycomp(kPdgP33m1232_DeltaPP), 1) == 1);
 
 #ifdef __GENIE_LOW_LEVEL_MESG_ENABLED__
   LOG("Pythia6Had", 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 // __GENIE_PYTHIA6_ENABLED__
 }

bool Pythia6Hadro2019::Hadronize ( GHepRecord * event ) const

privatevirtual

Implements genie::PythiaBaseHadro2019.

Definition at line 76 of file Pythia6Hadro2019.cxx.

References genie::PythiaBaseHadro2019::fLeadingQuark, fPythia, 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(), genie::utils::kinematics::W(), and genie::GHepParticle::X4().

 {
 #ifdef __GENIE_PYTHIA6_ENABLED__
 
   LOG("Pythia6Had", pNOTICE) << "Running PYTHIA6 hadronizer";
 
   const Interaction * interaction = event->Summary();
   const Kinematics & kinematics = interaction->Kine();
   double W = kinematics.W();
 
   LOG("Pythia6Had", pNOTICE)
     << "Fragmentation: "
     << "q = " << fLeadingQuark << ", qq = " << fRemnantDiquark
     << ", W = " << W;
 
   // Hadronize
   int ip = 0;
   py2ent_(&ip, &fLeadingQuark, &fRemnantDiquark, &W); // hadronizer
 
   // Get LUJETS record
   fPythia->GetPrimaries();
   TClonesArray * pythia_particles =
        (TClonesArray *) fPythia->ImportParticles("All");
 
   // Copy PYTHIA container to a new TClonesArray so as to transfer ownership
   // of the container and of its elements to the calling method
 
   int np = pythia_particles->GetEntries();
   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
   unsigned int i = 0;
   TMCParticle * p = 0;
   TIter particle_iter(pythia_particles);
   while( (p = (TMCParticle *) particle_iter.Next()) ) {
 
      int particle_pdg_code      = p->GetKF();
      int pythia_particle_status = p->GetKS();
 
      // Sanity check
      if(pythia_particle_status == 1) {
        if( pdg::IsQuark  (particle_pdg_code) ||
            pdg::IsDiQuark(particle_pdg_code) )
        {
          LOG("Pythia6Had", pERROR)
             << "Hadronization failed! Bare quarks appear in final state!";
          return false;
        }
      }
 
      // 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(p->GetPx(), p->GetPy(), p->GetPz(), p->GetEnergy());
      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 == 1) ?
          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 + p->GetParent();
      int mother2   = -1;
      int daughter1 = (p->GetFirstChild() <= 0 ) ? -1 : mom  + p->GetFirstChild();
      int daughter2 = (p->GetLastChild()  <= 0 ) ? -1 : mom  + p->GetLastChild();
 
      // 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("Pythia6Had", pDEBUG)
           << "Adding final state particle pdgc = " << particle.Pdg()
           << " with status = " << particle.Status();
 
      // Insert the particle in the list
      event->AddParticle(particle);
   }
   return true;
 
 #else
   return false;
 #endif // __GENIE_PYTHIA6_ENABLED__
 }

void Pythia6Hadro2019::Initialize ( void )

privatevirtual

Reimplemented from genie::PythiaBaseHadro2019.

Definition at line 328 of file Pythia6Hadro2019.cxx.

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

Referenced by Pythia6Hadro2019().

 {
   PythiaBaseHadro2019::Initialize();
 #ifdef __GENIE_PYTHIA6_ENABLED__
   fPythia = TPythia6::Instance();
   // sync GENIE/PYTHIA6 seed number
   // PYTHIA6 is a singleton, so do this from RandomGen for all
   // GENIE algorithms that use PYTHIA6
   RandomGen::Instance();
 #endif
 }

void Pythia6Hadro2019::LoadConfig ( void )

privatevirtual

Reimplemented from genie::PythiaBaseHadro2019.

Definition at line 308 of file Pythia6Hadro2019.cxx.

Referenced by Configure().

 {
   PythiaBaseHadro2019::LoadConfig();
 
 #ifdef __GENIE_PYTHIA6_ENABLED__
   fPythia->SetPARJ(2,  fSSBarSuppression       );
   fPythia->SetPARJ(21, fGaussianPt2            );
   fPythia->SetPARJ(23, fNonGaussianPt2Tail     );
   fPythia->SetPARJ(33, fRemainingECutoff       );
   fPythia->SetPARJ(1,  fDiQuarkSuppression     );
   fPythia->SetPARJ(11, fLightVMesonSuppression );
   fPythia->SetPARJ(12, fSVMesonSuppression     );
   fPythia->SetPARJ(41, fLunda                  );
   fPythia->SetPARJ(42, fLundb                  );
   fPythia->SetPARJ(45, fLundaDiq               );
 #endif
 
   LOG("Pythia6Had", pDEBUG) << this->GetConfig() ;
 }

void Pythia6Hadro2019::ProcessEventRecord ( GHepRecord * event ) const

virtual

Reimplemented from genie::PythiaBaseHadro2019.

Definition at line 60 of file Pythia6Hadro2019.cxx.

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

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

void Pythia6Hadro2019::RestoreOriginalDecayFlags ( void ) const

privatevirtual

Implements genie::PythiaBaseHadro2019.

Definition at line 270 of file Pythia6Hadro2019.cxx.

 {
 #ifdef __GENIE_PYTHIA6_ENABLED__
 
 fPythia->SetMDCY(fPythia->Pycomp(kPdgPi0),
    1, (fOriDecayFlag_pi0) ? 1 : 0);
 fPythia->SetMDCY(fPythia->Pycomp(kPdgK0),
    1, (fOriDecayFlag_K0)  ? 1 : 0);
 fPythia->SetMDCY(fPythia->Pycomp(kPdgAntiK0),
    1, (fOriDecayFlag_K0b) ? 1 : 0);
 fPythia->SetMDCY(fPythia->Pycomp(kPdgLambda),
    1, (fOriDecayFlag_L0)  ? 1 : 0);
 fPythia->SetMDCY(fPythia->Pycomp(kPdgAntiLambda),
    1, (fOriDecayFlag_L0b) ? 1 : 0);
 fPythia->SetMDCY(fPythia->Pycomp(kPdgP33m1232_DeltaM),
    1, (fOriDecayFlag_Dm)  ? 1 : 0);
 fPythia->SetMDCY(fPythia->Pycomp(kPdgP33m1232_Delta0),
    1, (fOriDecayFlag_D0)  ? 1 : 0);
 fPythia->SetMDCY(fPythia->Pycomp(kPdgP33m1232_DeltaP),
    1, (fOriDecayFlag_Dp)  ? 1 : 0);
 fPythia->SetMDCY(fPythia->Pycomp(kPdgP33m1232_DeltaPP),
    1, (fOriDecayFlag_Dpp) ? 1 : 0);
 
 #endif  // __GENIE_PYTHIA6_ENABLED__
 }

void Pythia6Hadro2019::SetDesiredDecayFlags ( void ) const

privatevirtual

Implements genie::PythiaBaseHadro2019.

Definition at line 244 of file Pythia6Hadro2019.cxx.

 {
 #ifdef __GENIE_PYTHIA6_ENABLED__
 
   fPythia->SetMDCY(fPythia->Pycomp(kPdgPi0),
       1, (fReqDecayFlag_pi0) ? 1 : 0);
   fPythia->SetMDCY(fPythia->Pycomp(kPdgK0),
      1, (fReqDecayFlag_K0)  ? 1 : 0);
   fPythia->SetMDCY(fPythia->Pycomp(kPdgAntiK0),
      1, (fReqDecayFlag_K0b) ? 1 : 0);
   fPythia->SetMDCY(fPythia->Pycomp(kPdgLambda),
      1, (fReqDecayFlag_L0)  ? 1 : 0);
   fPythia->SetMDCY(fPythia->Pycomp(kPdgAntiLambda),
      1, (fReqDecayFlag_L0b) ? 1 : 0);
   fPythia->SetMDCY(fPythia->Pycomp(kPdgP33m1232_DeltaM),
      1, (fReqDecayFlag_Dm)  ? 1 : 0);
   fPythia->SetMDCY(fPythia->Pycomp(kPdgP33m1232_Delta0),
      1, (fReqDecayFlag_D0)  ? 1 : 0);
   fPythia->SetMDCY(fPythia->Pycomp(kPdgP33m1232_DeltaP),
      1, (fReqDecayFlag_Dp)  ? 1 : 0);
   fPythia->SetMDCY(fPythia->Pycomp(kPdgP33m1232_DeltaPP),
      1, (fReqDecayFlag_Dpp) ? 1 : 0);
 
 #endif // __GENIE_PYTHIA6_ENABLED__
 }

Member Data Documentation

TPythia6* genie::Pythia6Hadro2019::fPythia

mutableprivate

PYTHIA6 wrapper class.

Definition at line 63 of file Pythia6Hadro2019.h.

Referenced by CopyOriginalDecayFlags(), Hadronize(), Initialize(), LoadConfig(), RestoreOriginalDecayFlags(), and SetDesiredDecayFlags().

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

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation