3 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
4 * See cxx source for full Copyright notice */
6 // Generator using DPMJET as an external generator
7 // The main DPMJET options are accessable for the user through this interface.
11 #include "AliGenDPMjetEventHeader.h"
18 class AliGenEventHeader;
21 class AliGenDPMjetEventHeader;
25 class AliGenDPMjet : public AliGenMC
30 AliGenDPMjet(Int_t npart);
31 AliGenDPMjet(const AliGenDPMjet &Dpmjet);
32 virtual ~AliGenDPMjet();
33 virtual void Generate();
35 virtual void FinishRun();
36 virtual void SetEnergyCMS(Float_t energy = 14000.) {fEnergyCMS = energy;}
37 virtual void SetProjectileBeamEnergy(Float_t benergy = 7000.) {fBeamEn = benergy;}
38 virtual void SetImpactParameterRange(Float_t bmin=0., Float_t bmax=1.)
39 {fMinImpactParam=bmin; fMaxImpactParam=bmax;}
40 virtual void SetProcess(DpmProcess_t iproc) {fProcess = iproc;}
41 virtual void SetCentral(Int_t icentr=-2) {fICentr = icentr;}
42 virtual void SetFlavor(Int_t flag=0) {fFlavor = flag;}
43 virtual void SetSelectAll(Int_t flag=0) {fSelectAll = flag;}
44 virtual void SetSpectators(Int_t spects=1) {fSpectators = spects;}
45 virtual void SetBoostLHC(Int_t flag=0) {fLHC = flag;}
46 virtual void SetPi0Decay(Int_t iPi0) {fPi0Decay = iPi0;}
47 virtual void SetDecayAll(Int_t iDec) {fDecayAll = iDec;}
48 virtual void SetIonPDGCodes();
49 virtual void GetImpactParameterRange(Float_t& bmin, Float_t& bmax)
50 {bmin = fMinImpactParam; bmax = fMaxImpactParam;}
51 virtual Int_t GetSpectators() {return fSpectators;}
52 virtual Int_t GetFlavor() {return fFlavor;}
54 virtual void SetGenImpPar(Float_t bValue) {fGenImpPar=bValue;}
55 virtual Float_t GetGenImpPar() {return fGenImpPar;}
57 // Trigger on a single particle (not related to calorimeter trigger above)
58 virtual void SetTriggerParticle(Int_t particle = 0, Float_t etamax = 0.9, Float_t ptmin = -1, Float_t ptmax = 1000)
59 {fTriggerParticle = particle; fTriggerEta = etamax; fTriggerMinPt = ptmin; fTriggerMaxPt = ptmax;}
61 virtual void SetTriggerChargedMultiplicity(Int_t multiplicity,
62 Float_t etamax = 0, Float_t ptmin = -1.)
63 {fTriggerMultiplicity = multiplicity; fTriggerMultiplicityEta = etamax;
64 fTriggerMultiplicityPtMin = ptmin;}
66 AliGenDPMjet & operator=(const AliGenDPMjet & rhs);
67 //void AddHeader(AliGenEventHeader* header);
69 void SetTuneForDiff(Bool_t a=kTRUE) {fkTuneForDiff=a;}
71 virtual void SetFragmentProd(Bool_t val) {fFragmentation = val;}
72 virtual Bool_t ProvidesCollisionGeometry() const {return kTRUE;}
75 Bool_t SelectFlavor(Int_t pid);
79 Float_t fBeamEn; // beam energy
80 Float_t fMinImpactParam; // minimum impact parameter
81 Float_t fMaxImpactParam; // maximum impact parameter
82 Int_t fICentr; // Flag to force central production
83 Int_t fSelectAll; // Flag to write the full event
84 Int_t fFlavor; // Selected particle flavor 4: charm+beauty 5: beauty
85 Int_t fTrials; // Number of trials
86 Int_t fNprimaries; // Number of generated particles
87 Int_t fSpectators; // put spectators on stack
88 Int_t fSpecn; // Num. of spectator neutrons
89 Int_t fSpecp; // Num. of spectator protons
90 TDPMjet *fDPMjet; // DPMjet
91 Int_t fNoGammas; // Don't write gammas if flag "on"
92 Int_t fLHC; // Assume LHC as lab frame
93 Int_t fPi0Decay; // Flag for pi0 decays
94 Int_t fDecayAll; // Flag to switch on long-lived particle decays
95 Float_t fGenImpPar; // GeneratedImpactParameter
96 DpmProcess_t fProcess; // Process type
97 Int_t fTriggerParticle; // Trigger on this particle ...
98 Float_t fTriggerEta; // .. within |eta| < fTriggerEta
99 Float_t fTriggerMinPt; // .. within pt > fTriggerMinPt
100 Float_t fTriggerMaxPt; // .. within pt < fTriggerMaxPt
101 // Multiplicity Trigger
102 Int_t fTriggerMultiplicity; // Triggered multiplicity
103 Float_t fTriggerMultiplicityEta; // Triggered multiplicity eta cut
104 Float_t fTriggerMultiplicityPtMin; // Triggered multiplicity min pt
106 Bool_t fkTuneForDiff; // Phojet tune
109 Bool_t fFragmentation; // Allows evaporation and fragments production
110 AliGenDPMjetEventHeader fHeader; // MC header
113 // adjust the weight from kinematic cuts
114 void AdjustWeights();
115 // check seleted daughters
116 Bool_t DaughtersSelection(TParticle* iparticle);
118 Bool_t Stable(TParticle* particle);
120 Bool_t CheckDiffraction();
121 Bool_t GetWeightsDiffraction(Double_t M, Double_t &Mmin, Double_t &Mmax,
122 Double_t &wSD, Double_t &wDD, Double_t &wND);
124 ClassDef(AliGenDPMjet,6) // AliGenerator interface to DPMJET