]>
Commit | Line | Data |
---|---|---|
a0134e39 | 1 | #ifndef ALIGENPYTHIA_H |
2 | #define ALIGENPYTHIA_H | |
3da30618 | 3 | /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. * |
4 | * See cxx source for full Copyright notice */ | |
5 | ||
6 | /* $Id$ */ | |
7 | ||
f87cfe57 | 8 | |
fe4da5cc | 9 | #include "AliGenerator.h" |
f1a48a38 | 10 | #include "AliPythia.h" |
11 | #include "AliDecayer.h" | |
f87cfe57 | 12 | #include <TArrayI.h> |
13 | ||
14 | class AliPythia; | |
15 | class TParticle; | |
fe4da5cc | 16 | |
17 | class AliGenPythia : public AliGenerator | |
18 | { | |
fe4da5cc | 19 | public: |
20 | AliGenPythia(); | |
21 | AliGenPythia(Int_t npart); | |
f87cfe57 | 22 | AliGenPythia(const AliGenPythia &Pythia); |
fe4da5cc | 23 | virtual ~AliGenPythia(); |
24 | virtual void Generate(); | |
25 | virtual void Init(); | |
26 | // select process type | |
f1a48a38 | 27 | virtual void SetProcess(Process_t proc = kPyCharm) {fProcess = proc;} |
fe4da5cc | 28 | // select structure function |
f1a48a38 | 29 | virtual void SetStrucFunc(StrucFunc_t func = kGRV_HO) {fStrucFunc = func;} |
fe4da5cc | 30 | // select pt of hard scattering |
f1a48a38 | 31 | virtual void SetPtHard(Float_t ptmin = 0, Float_t ptmax = 1.e10) |
32 | {fPtHardMin = ptmin; fPtHardMax = ptmax; } | |
fe4da5cc | 33 | // set centre of mass energy |
f1a48a38 | 34 | virtual void SetEnergyCMS(Float_t energy = 5500) {fEnergyCMS = energy;} |
fe4da5cc | 35 | // force decay type |
f1a48a38 | 36 | virtual void SetForceDecay(Decay_t decay = kSemiMuonic) {fForceDecay = decay;} |
811826d8 | 37 | // treat protons as inside nuclei |
38 | virtual void SetNuclei(Int_t a1, Int_t a2); | |
fe4da5cc | 39 | // get cross section of process |
40 | virtual Float_t GetXsection() {return fXsection;} | |
a8228d85 | 41 | // Check PDG code |
f1a48a38 | 42 | virtual Int_t CheckPDGCode(Int_t pdgcode); |
9faa0850 | 43 | virtual void FinishRun(); |
44 | ||
f87cfe57 | 45 | // Assignment Operator |
46 | AliGenPythia & operator=(const AliGenPythia & rhs); | |
a0134e39 | 47 | protected: |
48 | Process_t fProcess; // Process type | |
49 | StrucFunc_t fStrucFunc; // Structure Function | |
50 | Decay_t fForceDecay; // Decay channel are forced | |
51 | Float_t fEnergyCMS; // Centre of mass energy | |
52 | Float_t fKineBias; // Bias from kinematic selection | |
53 | Int_t fTrials; // Number of trials | |
54 | TArrayI fParentSelect; // Parent particles to be selected | |
55 | TArrayI fChildSelect; // Decay products to be selected | |
56 | Float_t fXsection; // Cross-section | |
a8a6107b | 57 | AliPythia *fPythia; //! Pythia |
a0134e39 | 58 | Float_t fPtHardMin; // lower pT-hard cut |
59 | Float_t fPtHardMax; // higher pT-hard cut | |
811826d8 | 60 | Int_t fNucA1; // mass number nucleus side 1 |
61 | Int_t fNucA2; // mass number nucleus side 2 | |
f1a48a38 | 62 | Bool_t fFullEvent; // Write Full event if true |
c78d7f8c | 63 | AliDecayer *fDecayer; // ! pointer to the decayer instance |
a0134e39 | 64 | private: |
65 | // check if particle is selected as parent particle | |
66 | Bool_t ParentSelected(Int_t ip); | |
67 | // check if particle is selected as child particle | |
68 | Bool_t ChildSelected(Int_t ip); | |
69 | // all kinematic selection cuts go here | |
70 | Bool_t KinematicSelection(TParticle *particle); | |
71 | // adjust the weight from kinematic cuts | |
72 | void AdjustWeights(); | |
73 | ||
f87cfe57 | 74 | ClassDef(AliGenPythia,1) // AliGenerator interface to Pythia |
fe4da5cc | 75 | }; |
76 | #endif | |
77 | ||
78 | ||
79 | ||
80 | ||
81 |