| | 1 | #ifndef ALIGENHIJING_H |
| | 2 | #define ALIGENHIJING_H |
| | 3 | /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. * |
| | 4 | * See cxx source for full Copyright notice */ |
| | 5 | |
| | 6 | /* $Id$ */ |
| | 7 | |
| | 8 | // Generator using HIJING as an external generator |
| | 9 | // The main HIJING options are accessable for the user through this interface. |
| | 10 | // andreas.morsch@cern.ch |
| | 11 | |
| | 12 | #include "AliGenMC.h" |
| | 13 | #include <TString.h> |
| | 14 | |
| | 15 | class THijing; |
| | 16 | class TParticle; |
| | 17 | class TClonesArray; |
| | 18 | class TGraph; |
| | 19 | |
| | 20 | class AliGenHijing : public AliGenMC |
| | 21 | { |
| | 22 | enum {kNoTrigger, kHardProcesses, kDirectPhotons}; |
| | 23 | |
| | 24 | public: |
| | 25 | AliGenHijing(); |
| | 26 | AliGenHijing(Int_t npart); |
| | 27 | virtual ~AliGenHijing(); |
| | 28 | virtual void Generate(); |
| | 29 | virtual void Init(); |
| | 30 | // set centre of mass energy |
| | 31 | virtual void SetEnergyCMS(Float_t energy=5500) {fEnergyCMS=energy;} |
| | 32 | virtual void SetReferenceFrame(TString frame="CMS") |
| | 33 | {fFrame=frame;} |
| | 34 | |
| | 35 | virtual void SetImpactParameterRange(Float_t bmin = 0, Float_t bmax = 15.) |
| | 36 | {fMinImpactParam=bmin; fMaxImpactParam=bmax;} |
| | 37 | virtual void KeepFullEvent(); |
| | 38 | virtual void SetJetQuenching(Int_t flag=1) {fQuench = flag;} |
| | 39 | virtual void SetShadowing(Int_t flag=1) {fShadowing = flag;} |
| | 40 | virtual void SetDecaysOff(Int_t flag=1) {fDecaysOff = flag;} |
| | 41 | virtual void SetTrigger(Int_t flag=kNoTrigger) {fTrigger = flag;} |
| | 42 | virtual void SetFlavor(Int_t flag=0) {fFlavor = flag;} |
| | 43 | virtual void SetEvaluate(Int_t flag=0) {fEvaluate = flag;} |
| | 44 | virtual void SetSelectAll(Int_t flag=0) {fSelectAll = flag;} |
| | 45 | virtual void SetRadiation(Int_t flag=3) {fRadiation = flag;} |
| | 46 | virtual void SetSpectators(Int_t spects=1) {fSpectators = spects;} |
| | 47 | virtual void SetPtJet(Float_t ptmin) {fPtMinJet = ptmin;} |
| | 48 | virtual void SetSimpleJets(Int_t flag=0) {fSimpleJet = flag;} |
| | 49 | virtual void SetNoGammas(Int_t flag=0) {fNoGammas = flag;} |
| | 50 | |
| | 51 | virtual void SetJetEtaRange(Float_t etamin = -20., Float_t etamax = 20.) |
| | 52 | {fEtaMinJet = etamin; fEtaMaxJet = etamax;} |
| | 53 | virtual void SetJetPhiRange(Float_t phimin = -180., Float_t phimax = 180.) |
| | 54 | {fPhiMinJet = TMath::Pi()*phimin/180.; fPhiMaxJet = TMath::Pi()*phimax/180.;} |
| | 55 | virtual void SetBoostLHC(Int_t flag = 0) {fLHC = flag;} |
| | 56 | virtual void SetRandomPz(Bool_t flag = 0) {fRandomPz = flag;} |
| | 57 | virtual void SwitchOffHeavyQuarks(Bool_t flag = kTRUE) {fNoHeavyQuarks = flag;} |
| | 58 | |
| | 59 | |
| | 60 | // Getters |
| | 61 | virtual Float_t GetEnergyCMS() const {return fEnergyCMS;} |
| | 62 | virtual TString GetReferenceFrame() const {return fFrame;} |
| | 63 | virtual void GetProjectile(TString& tar, Int_t& a, Int_t& z) const |
| | 64 | {tar = fProjectile; a = fAProjectile; z = fZProjectile;} |
| | 65 | virtual void GetTarget(TString& tar, Int_t& a, Int_t& z) const |
| | 66 | {tar = fTarget; a = fATarget; z = fZTarget;} |
| | 67 | virtual void GetImpactParameterRange(Float_t& bmin, Float_t& bmax) const |
| | 68 | {bmin = fMinImpactParam; bmax = fMaxImpactParam;} |
| | 69 | virtual Int_t GetJetQuenching() const {return fQuench;} |
| | 70 | virtual Int_t GetShadowing() const {return fShadowing;} |
| | 71 | virtual Int_t GetTrigger() const {return fTrigger;} |
| | 72 | virtual Int_t GetFlavor() const {return fFlavor;} |
| | 73 | virtual Int_t GetRadiation() const {return fRadiation;} |
| | 74 | virtual Int_t GetSpectators() const {return fSpectators;} |
| | 75 | virtual Float_t GetPtJet() const {return fPtMinJet;} |
| | 76 | virtual void GetJetEtaRange(Float_t& etamin, Float_t& etamax) const |
| | 77 | {etamin = fEtaMinJet; etamax = fEtaMaxJet;} |
| | 78 | virtual void GetJetPhiRange(Float_t& phimin, Float_t& phimax) const |
| | 79 | {phimin = fPhiMinJet*180./TMath::Pi(); phimax = fPhiMaxJet*180./TMath::Pi();} |
| | 80 | |
| | 81 | |
| | 82 | // Physics Routines |
| | 83 | virtual Bool_t ProvidesCollisionGeometry() const {return kTRUE;} |
| | 84 | virtual void EvaluateCrossSections(); |
| | 85 | virtual TGraph* CrossSection() {return fDsigmaDb;} |
| | 86 | virtual TGraph* BinaryCollisions() {return fDnDb;} |
| | 87 | virtual Bool_t CheckTrigger(); |
| | 88 | // |
| | 89 | protected: |
| | 90 | Bool_t SelectFlavor(Int_t pid); |
| | 91 | void MakeHeader(); |
| | 92 | protected: |
| | 93 | TString fFrame; // Reference frame |
| | 94 | Float_t fMinImpactParam; // minimum impact parameter |
| | 95 | Float_t fMaxImpactParam; // maximum impact parameter |
| | 96 | Int_t fKeep; // Flag to keep full event information |
| | 97 | Int_t fQuench; // Flag to switch on jet quenching |
| | 98 | Int_t fShadowing; // Flag to switch on nuclear effects on parton distribution function |
| | 99 | Int_t fDecaysOff; // Flag to turn off decays of pi0, K_s, D, Lambda, sigma |
| | 100 | Int_t fTrigger; // Trigger type |
| | 101 | Int_t fEvaluate; // Evaluate total and partial cross-sections |
| | 102 | Int_t fSelectAll; // Flag to write the full event |
| | 103 | Int_t fFlavor; // Selected particle flavor 4: charm+beauty 5: beauty |
| | 104 | Float_t fEnergyCMS; // Centre of mass energy |
| | 105 | Float_t fKineBias; // Bias from kinematic selection |
| | 106 | Int_t fTrials; // Number of trials |
| | 107 | Float_t fXsection; // Cross-section |
| | 108 | THijing *fHijing; //!Hijing |
| | 109 | Float_t fPtHardMin; // lower pT-hard cut |
| | 110 | Float_t fPtHardMax; // higher pT-hard cut |
| | 111 | Int_t fSpectators; // put spectators on stack |
| | 112 | TGraph* fDsigmaDb; // dSigma/db for the system |
| | 113 | TGraph* fDnDb; // dNBinaryCollisions/db |
| | 114 | Float_t fPtMinJet; // Minimum Pt of triggered Jet |
| | 115 | Float_t fEtaMinJet; // Minimum eta of triggered Jet |
| | 116 | Float_t fEtaMaxJet; // Maximum eta of triggered Jet |
| | 117 | Float_t fPhiMinJet; // At least one of triggered Jets must be in this |
| | 118 | Float_t fPhiMaxJet; // phi range |
| | 119 | Int_t fRadiation; // Flag to switch on/off initial and final state radiation |
| | 120 | Int_t fSimpleJet; // Flag to produce simple tiggered jet topology |
| | 121 | Int_t fNoGammas; // Don't write gammas if flag "on" |
| | 122 | // ZDC proposal (by Chiara) to store num. of SPECTATORS protons and neutrons |
| | 123 | Int_t fProjectileSpecn;// Num. of spectator neutrons from projectile nucleus |
| | 124 | Int_t fProjectileSpecp;// Num. of spectator protons from projectile nucleus |
| | 125 | Int_t fTargetSpecn; // Num. of spectator neutrons from target nucleus |
| | 126 | Int_t fTargetSpecp; // Num. of spectator protons from target nucleus |
| | 127 | Int_t fLHC; // Assume LHC as lab frame |
| | 128 | Bool_t fRandomPz; // Randomise sign of pz event by event |
| | 129 | Bool_t fNoHeavyQuarks; // If true no heavy quarks are produced |
| | 130 | |
| | 131 | private: |
| | 132 | AliGenHijing(const AliGenHijing &Hijing); |
| | 133 | AliGenHijing & operator=(const AliGenHijing & rhs); |
| | 134 | |
| | 135 | // adjust the weight from kinematic cuts |
| | 136 | void AdjustWeights(); |
| | 137 | // check seleted daughters |
| | 138 | Bool_t DaughtersSelection(TParticle* iparticle); |
| | 139 | // check if stable |
| | 140 | Bool_t Stable(TParticle* particle) const; |
| | 141 | |
| | 142 | ClassDef(AliGenHijing, 6) // AliGenerator interface to Hijing |
| | 143 | }; |
| | 144 | #endif |
| | 145 | |
| | 146 | |
| | 147 | |
| | 148 | |
| | 149 | |
git://git.uio.no / u / mrichter / AliRoot.git / blame_incremental