#include <TPythia6.h>
#include <AliRndm.h>
-
+#include <AliStructFuncType.h>
typedef enum
-{kPyCharm, kPyBeauty, kPyCharmUnforced, kPyBeautyUnforced, kPyJpsi, kPyJpsiChi, kPyMb, kPyJets, kPyDirectGamma, kPyCharmPbMNR, kPyD0PbMNR}
+{kPyCharm, kPyBeauty, kPyCharmUnforced, kPyBeautyUnforced,
+ kPyJpsi, kPyJpsiChi, kPyMb, kPyMbNonDiffr, kPyJets, kPyDirectGamma,
+ kPyCharmPbMNR, kPyD0PbMNR, kPyBeautyPbMNR}
Process_t;
-
+/*
typedef enum
{
- kDO_Set_1=1006,
- kGRV_LO=5005,
- kGRV_HO=5006,
- kMRS_D_minus=3031,
- kMRS_D0=3030,
- kMRS_G=3041,
- kCTEQ_2pM=4024,
- kCTEQ_4L=4032,
- kCTEQ_4M=4034,
- kMRST_cg_LO=3072,
- kCTEQ_5L=4046,
- kGRV_LO_98=5012
+ kDOSet1 = 1006,
+ kGRVLO = 5005,
+ kGRVHO = 5006,
+ kMRSDminus = 3031,
+ kMRSD0 = 3030,
+ kMRSG = 3041,
+ kCTEQ2pM = 4024,
+ kCTEQ4L = 4032,
+ kCTEQ4M = 4034,
+ kMRSTcgLO = 3072,
+ kCTEQ5L = 4046,
+ kGRVLO98 = 5012
}
StrucFunc_t;
-
+*/
class AliPythia : public TPythia6, public AliRndm
{
virtual void ProcInit
(Process_t process, Float_t energy, StrucFunc_t strucfunc);
// treat protons as inside nuclei
- virtual void SetNuclei(Int_t a1, Int_t a2);
+ virtual void SetNuclei(Int_t a1, Int_t a2);
// Print particle properties
virtual void PrintParticles();
virtual void ResetDecayTable();
virtual void SetDecayTable();
+ virtual void Pycell(Int_t& nclus);
+ virtual void Pyclus(Int_t& nclus);
// return instance of the singleton
static AliPythia* Instance();
Float_t fEcms; // Centre of mass energy
StrucFunc_t fStrucFunc; // Structure function
Int_t fDefMDCY[501]; // ! Default decay switches per particle
- Int_t fDefMDME[2000]; // ! Default decay switches per mode
+ Int_t fDefMDME[2001]; // ! Default decay switches per mode
static AliPythia* fgAliPythia; // Pointer to single instance
private:
AliPythia();