]>
Commit | Line | Data |
---|---|---|
110410f9 | 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 | ||
9 | #include "AliGenerator.h" | |
10 | #include "GenTypeDefs.h" | |
3b8bd63f | 11 | #include <TString.h> |
110410f9 | 12 | #include <TArrayI.h> |
13 | ||
65fb704d | 14 | class THijing; |
110410f9 | 15 | class TArrayI; |
16 | class TParticle; | |
b23a502f | 17 | class TClonesArray; |
110410f9 | 18 | |
19 | class AliGenHijing : public AliGenerator | |
20 | { | |
21 | enum {kNoTrigger, kHardProcesses, kDirectPhotons}; | |
22 | ||
23 | public: | |
24 | AliGenHijing(); | |
25 | AliGenHijing(Int_t npart); | |
26 | AliGenHijing(const AliGenHijing &Hijing); | |
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;} | |
59d7832e | 32 | virtual void SetReferenceFrame(TString frame="CMS") |
3b8bd63f | 33 | {fFrame=frame;} |
59d7832e | 34 | virtual void SetProjectile(TString proj="A", |
3b8bd63f | 35 | Int_t a=208, Int_t z=82) |
110410f9 | 36 | {fProjectile = proj; |
37 | fAProjectile = a; | |
38 | fZProjectile = z;} | |
59d7832e | 39 | virtual void SetTarget(TString tar="A", |
3b8bd63f | 40 | Int_t a=208, Int_t z=82) |
110410f9 | 41 | {fTarget = tar; |
42 | fATarget = a; | |
43 | fZTarget =z;} | |
44 | virtual void SetImpactParameterRange(Float_t bmin = 0, Float_t bmax =15.) | |
45 | { fMinImpactParam=bmin; | |
46 | fMaxImpactParam=bmax; | |
47 | } | |
48 | virtual void KeepFullEvent(); | |
49 | virtual void SetJetQuenching(Int_t flag=1) {fQuench = flag;} | |
50 | virtual void SetShadowing(Int_t flag=1) {fShadowing = flag;} | |
51 | virtual void SetDecaysOff(Int_t flag=1) {fDecaysOff = flag;} | |
52 | virtual void SetTrigger(Int_t flag=kNoTrigger) {fTrigger = flag;} | |
eb342d2d | 53 | virtual void SetFlavor(Int_t flag=0) {fFlavor = flag;} |
110410f9 | 54 | virtual void SetEvaluate(Int_t flag=0) {fEvaluate = flag;} |
55 | virtual void SetSelectAll(Int_t flag=0) {fSelectAll = flag;} | |
2705454f | 56 | virtual void SetSpectators(Int_t spects=1) {fSpectators = spects;} |
eb342d2d | 57 | AliGenHijing & operator=(const AliGenHijing & rhs); |
110410f9 | 58 | // Physics Routines |
59 | virtual void EvaluateCrossSections(); | |
eb342d2d | 60 | protected: |
61 | Bool_t SelectFlavor(Int_t pid); | |
3b8bd63f | 62 | void MakeHeader(); |
63 | ||
110410f9 | 64 | protected: |
3b8bd63f | 65 | TString fFrame; // Reference frame |
66 | TString fProjectile; // Projectile | |
67 | TString fTarget; // Target | |
68 | Int_t fAProjectile; // Projectile A | |
69 | Int_t fZProjectile; // Projectile Z | |
70 | Int_t fATarget; // Target A | |
71 | Int_t fZTarget; // Target Z | |
110410f9 | 72 | Float_t fMinImpactParam; // minimum impact parameter |
73 | Float_t fMaxImpactParam; // maximum impact parameter | |
3b8bd63f | 74 | Int_t fKeep; // Flag to keep full event information |
75 | Int_t fQuench; // Flag to switch on jet quenching | |
76 | Int_t fShadowing; // Flag to switch on voclear effects on parton distribution function | |
77 | Int_t fDecaysOff; // Flag to turn off decays of pi0, K_s, D, Lambda, sigma | |
78 | Int_t fTrigger; // Trigger type | |
79 | Int_t fEvaluate; // Evaluate total and partial cross-sections | |
80 | Int_t fSelectAll; // Flag to write the full event | |
81 | Int_t fFlavor; // Selected particle flavor 4: charm+beauty 5: beauty | |
82 | Decay_t fForceDecay; // Decay channel are forced | |
83 | Float_t fEnergyCMS; // Centre of mass energy | |
84 | Float_t fKineBias; // Bias from kinematic selection | |
85 | Int_t fTrials; // Number of trials | |
86 | TArrayI fParentSelect; // Parent particles to be selected | |
87 | TArrayI fChildSelect; // Decay products to be selected | |
88 | Float_t fXsection; // Cross-section | |
89 | THijing *fHijing; // Hijing | |
90 | Float_t fPtHardMin; // lower pT-hard cut | |
91 | Float_t fPtHardMax; // higher pT-hard cut | |
2705454f | 92 | Int_t fSpectators; // put spectators on stack |
110410f9 | 93 | |
94 | private: | |
95 | // check if particle is selected as parent particle | |
96 | Bool_t ParentSelected(Int_t ip); | |
97 | // check if particle is selected as child particle | |
98 | Bool_t ChildSelected(Int_t ip); | |
99 | // all kinematic selection cuts go here | |
100 | Bool_t KinematicSelection(TParticle *particle); | |
101 | // adjust the weight from kinematic cuts | |
102 | void AdjustWeights(); | |
103 | // check seleted daughters | |
104 | Bool_t DaughtersSelection(TParticle* iparticle, TClonesArray* particles); | |
608708ba | 105 | // check if stable |
106 | Bool_t Stable(TParticle* particle); | |
107 | ||
110410f9 | 108 | ClassDef(AliGenHijing,1) // AliGenerator interface to Hijing |
109 | }; | |
110 | #endif | |
111 | ||
112 | ||
113 | ||
114 | ||
115 |