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"
20 class AliGenEventHeader;
23 class AliGenDPMjetEventHeader;
27 class AliGenDPMjet : public AliGenMC
32 AliGenDPMjet(Int_t npart);
33 AliGenDPMjet(const AliGenDPMjet &Dpmjet);
34 virtual ~AliGenDPMjet();
35 virtual void Generate();
37 virtual void FinishRun();
38 virtual void SetEnergyCMS(Float_t energy = 14000.) {fEnergyCMS = energy;}
39 virtual void SetProjectileBeamEnergy(Float_t benergy = 7000.) {fBeamEn = benergy;}
40 virtual void SetImpactParameterRange(Float_t bmin=0., Float_t bmax=1.)
41 {fMinImpactParam=bmin; fMaxImpactParam=bmax;}
42 virtual void SetProcess(DpmProcess_t iproc) {fProcess = iproc;}
43 virtual void SetCentral(Int_t icentr=-2) {fICentr = icentr;}
44 virtual void SetFlavor(Int_t flag=0) {fFlavor = flag;}
45 virtual void SetSelectAll(Int_t flag=0) {fSelectAll = flag;}
46 virtual void SetSpectators(Int_t spects=1) {fSpectators = spects;}
47 virtual void SetBoostLHC(Int_t flag=0) {fLHC = flag;}
48 virtual void SetPi0Decay(Int_t iPi0) {fPi0Decay = iPi0;}
49 virtual void SetDecayAll(Int_t iDec) {fDecayAll = iDec;}
50 virtual void GetImpactParameterRange(Float_t& bmin, Float_t& bmax)
51 {bmin = fMinImpactParam; bmax = fMaxImpactParam;}
52 virtual Int_t GetSpectators() {return fSpectators;}
53 virtual Int_t GetFlavor() {return fFlavor;}
55 virtual void SetGenImpPar(Float_t bValue) {fGenImpPar=bValue;}
56 virtual Float_t GetGenImpPar() {return fGenImpPar;}
57 virtual void SetTriggerChargedMultiplicity(Int_t multiplicity,
58 Float_t etamax = 0, Float_t ptmin = -1.)
59 {fTriggerMultiplicity = multiplicity; fTriggerMultiplicityEta = etamax;
60 fTriggerMultiplicityPtMin = ptmin;}
62 AliGenDPMjet & operator=(const AliGenDPMjet & rhs);
63 //void AddHeader(AliGenEventHeader* header);
65 void SetTuneForDiff(Bool_t a=kTRUE) {fkTuneForDiff=a;}
67 virtual void SetFragmentProd(Bool_t val) {fFragmentation = val;}
68 virtual Bool_t ProvidesCollisionGeometry() const {return kTRUE;}
71 Bool_t SelectFlavor(Int_t pid);
75 Float_t fBeamEn; // beam energy
76 Float_t fMinImpactParam; // minimum impact parameter
77 Float_t fMaxImpactParam; // maximum impact parameter
78 Int_t fICentr; // Flag to force central production
79 Int_t fSelectAll; // Flag to write the full event
80 Int_t fFlavor; // Selected particle flavor 4: charm+beauty 5: beauty
81 Int_t fTrials; // Number of trials
82 Int_t fSpectators; // put spectators on stack
83 Int_t fSpecn; // Num. of spectator neutrons
84 Int_t fSpecp; // Num. of spectator protons
85 TDPMjet *fDPMjet; // DPMjet
86 Int_t fNoGammas; // Don't write gammas if flag "on"
87 Int_t fLHC; // Assume LHC as lab frame
88 Int_t fPi0Decay; // Flag for pi0 decays
89 Int_t fDecayAll; // Flag to switch on long-lived particle decays
90 Float_t fGenImpPar; // GeneratedImpactParameter
91 DpmProcess_t fProcess; // Process type
92 // Multiplicity Trigger
93 Int_t fTriggerMultiplicity; // Triggered multiplicity
94 Float_t fTriggerMultiplicityEta; // Triggered multiplicity eta cut
95 Float_t fTriggerMultiplicityPtMin; // Triggered multiplicity min pt
97 Bool_t fkTuneForDiff; // Phojet tune
100 Bool_t fFragmentation; // Allows evaporation and fragments production
101 AliGenDPMjetEventHeader fHeader; // MC header
104 // adjust the weight from kinematic cuts
105 void AdjustWeights();
106 // check seleted daughters
107 Bool_t DaughtersSelection(TParticle* iparticle);
109 Bool_t Stable(TParticle* particle);
111 Bool_t CheckDiffraction();
112 Bool_t GetWeightsDiffraction(Double_t M, Double_t &Mmin, Double_t &Mmax,
113 Double_t &wSD, Double_t &wDD, Double_t &wND);
115 ClassDef(AliGenDPMjet,6) // AliGenerator interface to DPMJET