4 //+SEQ,CopyRight,T=NOINCLUDE.
6 //////////////////////////////////////////////////////////////////////////
10 // This class implements an interface to the DPMJET 3.0 event generator.//
12 //////////////////////////////////////////////////////////////////////////
14 #ifndef ROOT_TGenerator
16 #include "TGenerator.h"
19 typedef enum {kDpmMb, kDpmMbNonDiffr} DpmProcess_t;
21 class TDPMjet : public TGenerator {
26 TDPMjet(DpmProcess_t ip, Int_t Ip, Int_t Ipz, Int_t It, Int_t Itz, Double_t Epn, Double_t CMEn);
28 virtual ~TDPMjet() {;}
30 virtual void Initialize();
32 virtual void GenerateEvent();
34 virtual Int_t ImportParticles(TClonesArray *particles, Option_t *option="");
35 virtual TObjArray *ImportParticles(Option_t * /*option*/) {return 0;}
37 // Parameters for the generation:
38 virtual void SetNEvent(Int_t iev) {fNEvent = iev;}
39 virtual Int_t GetNEvent() {return fNEvent;}
41 virtual void SetfIp(Int_t Ip) {fIp = Ip;}
42 virtual Int_t GetfIp() const {return fIp;}
44 virtual void SetfIpz(Int_t Ipz) {fIpz = Ipz;}
45 virtual Int_t GetfIpz() const {return fIpz;}
47 virtual void SetfIt(Int_t It) {fIt = It;}
48 virtual Int_t GetfIt() const {return fIt;}
50 virtual void SetfItz(Int_t Itz) {fItz = Itz;}
51 virtual Int_t GetfItz() const {return fItz;}
53 virtual void SetfEpn(Double_t Epn) {fEpn = Epn;}
54 virtual Double_t GetfEpn() const {return fEpn;}
56 virtual void SetfCMEn(Double_t CMEn) {fCMEn = CMEn;}
57 virtual Double_t GetfCMEn() const {return fCMEn;}
59 virtual void SetfIdp(Int_t idp) {fIdp = idp;}
60 virtual Int_t GetfIdp() const {return fIdp;}
62 virtual void SetbRange(Double_t bmin, Double_t bmax)
63 {fBmin = bmin; fBmax = bmax;}
64 virtual Double_t GetMinImpactParameter() const {return fBmin;}
65 virtual Double_t GetMaxImpactParameter() const {return fBmax;}
67 virtual void SetfFCentr(Int_t icentr) {fFCentr = icentr;}
68 virtual Int_t GetfFCentr() const {return fFCentr;}
70 virtual void SetPi0Decay(Int_t iPi0) {fPi0Decay = iPi0;}
73 // Access to DPMJET common blocks:
74 virtual Int_t GetEvNum() const;
75 virtual Int_t GetEntriesNum() const;
76 virtual Int_t GetNumStablePc() const;
77 virtual Float_t GetTotEnergy() const;
78 virtual Int_t GetStatusCode(Int_t evnum) const;
79 virtual Int_t GetPDGCode(Int_t evnum) const;
80 virtual Double_t Getpx(Int_t evnum) const;
81 virtual Double_t Getpy(Int_t evnum) const;
82 virtual Double_t Getpz(Int_t evnum) const;
83 virtual Double_t GetEnergy(Int_t evnum) const;
84 virtual Double_t GetMass(Int_t evnum) const;
86 virtual Int_t GetFragmentA(Int_t evnum) const;
87 virtual Int_t GetFragmentZ(Int_t evnum) const;
89 virtual Double_t GetXSFrac() const;
90 virtual Double_t GetBImpac() const;
91 virtual Double_t GetProjRadius() const;
92 virtual Double_t GetTargRadius() const;
93 virtual Int_t GetProjWounded() const;
94 virtual Int_t GetTargWounded() const;
95 virtual Int_t GetProjSpectators() const;
96 virtual Int_t GetTargSpectators() const;
99 // Access to DPMJET routines:
100 virtual void Dt_Dtuini(int nevts, double epn, int npmass, int npchar,
101 int ntmass, int ntchar, int idp, int iemu);
103 virtual void Dt_Kkinc(int npmass, int npchar, int ntmass, int ntchar,
104 int idp, double elab, int kkmat, int irej);
106 virtual void Pho_Phist(int imode, double weight);
108 virtual void Dt_Dtuout();
110 virtual void Dt_Rndm(int idummy);
111 virtual void Dt_Rndmst(int na1, int na2, int na3, int nb1);
112 virtual void Dt_Rndmin(int u, int c, int cd, int cm, int i, int j);
113 virtual void Dt_Rndmou(int u, int c, int cd, int cm, int i, int j);
117 Int_t fNEvent; // Event number to be generated
118 Int_t fIp; // Projectile mass
119 Int_t fIpz; // Projectile charge
120 Int_t fIt; // Target mass
121 Int_t fItz; // Target charge
122 Float_t fEpn; // Beam energy
123 Float_t fPpn; // Beam momentum
124 Float_t fCMEn; // Energy in CM
125 Int_t fIdp; // Internal particle code
126 Float_t fBmin; // Minimum impact parameter
127 Float_t fBmax; // Maximum impact parameter
128 Int_t fFCentr; // Flag to force central collisions
129 Int_t fPi0Decay;// Flag for pi0 decays
130 DpmProcess_t fProcess; // Process type
132 ClassDef(TDPMjet,2) //Interface to DPMJET Event Generator