LeptoHadPythia8.h

Go to the documentation of this file.

//____________________________________________________________________________
/*!
 
\class    genie::LeptoHadPythia8
 
\brief    Provides access to the PYTHIA8 hadronization. \n
 
\author   Alfonso Garcia <aagarciasoto \at km3net.de>
          IFIC (Valencia)
 
\created  December 12, 2024
 
\cpright  Copyright (c) 2003-2025, The GENIE Collaboration
          For the full text of the license visit http://copyright.genie-mc.org
          or see $GENIE/LICENSE
*/
//____________________________________________________________________________
 
#ifndef _LEPTO_HAD_PYTHIA8__H_
#define _LEPTO_HAD_PYTHIA8__H_
 
#include "Framework/EventGen/EventRecordVisitorI.h"
#include "Framework/Conventions/Constants.h"
#include "Framework/Numerical/RandomGen.h"
#include "Framework/Numerical/MathUtils.h"
#include "Framework/GHEP/GHepRecord.h"
#include "Framework/GHEP/GHepParticle.h"
#include "Framework/GHEP/GHepFlags.h"
#include "Framework/ParticleData/PDGCodes.h"
#include "Framework/ParticleData/PDGUtils.h"
#include "Framework/ParticleData/PDGLibrary.h"
#include "Framework/EventGen/EVGThreadException.h"
 
#ifdef __GENIE_PYTHIA8_ENABLED__
#include "Framework/Utils/Pythia8Singleton.h"
#endif
 
using namespace genie;
using namespace genie::constants;
using namespace genie::utils::math;
 
namespace genie {
 
class GHepParticle;
 
class LeptoHadPythia8 : public EventRecordVisitorI {
 
public:
  LeptoHadPythia8();
  LeptoHadPythia8(string config);
  virtual ~LeptoHadPythia8();
 
  //-- implement the HadronizationModelI interface
  void ProcessEventRecord(GHepRecord * event) const;
 
  //-- overload the Algorithm::Configure() methods to load private data
  //   members from configuration options
  void Configure(const Registry & config);
  void Configure(string config);
 
private:
 
  bool           Hadronize        (GHepRecord * event) const;
  void           Initialize       (void)               const;
  void           LoadConfig       (void);
 
  int            getMeson         (int,int,double)     const;    ///< create meson
  int            getBaryon        (int,int,double)     const;    ///< create baryon
  double         getRandomZ       (double,double)      const;    ///< fragmentation
  double Afrag;                                                ///< fragmentation parameter
  double Bfrag;                                                ///< fragmentation parameter
  double mesonRateSum[3];                                      ///<meson parameter
  double CGOct[6] = { 0.75, 0.5, 0., 0.1667, 0.0833, 0.1667};  ///<baryon parameter
  double CGDec[6] = { 0.00, 0.0, 1., 0.3333, 0.6667, 0.3333};  ///<baryon parameter
  double CGSum[6];                                             ///< baryon parameter
 
  //-- configuration parameters
  int    fMaxIterHad;         // Maxmium number of iterations to look for a combination of hadrons
  double fPrimordialKT;       // Width of Gaussian distribution for the primordial transverse momentum kT of partons in the nucleon.
  double fRemnantPT;          // Width of Gaussian distribution in transverse momentum when a non-trivial target remnant is split into two particles
  double fMinESinglet;        // It is, with quark masses added, used to define the minimum allowable energy of a colour-singlet parton system.
  double fWmin;               // Minimum value of W
 
  // PYTHIA physics configuration parameters used
  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
 
#ifdef __GENIE_PYTHIA8_ENABLED__
  mutable Pythia8::Pythia * fPythia;         ///< PYTHIA8 object
#endif
 
};
 
}         // genie namespace
 
#endif    // _LEPTO_HAD_PYTHIA8__H_