]>
Commit | Line | Data |
---|---|---|
762e8424 | 1 | #ifndef ALIJETMODELBASETASK_H |
2 | #define ALIJETMODELBASETASK_H | |
3 | ||
4 | // $Id$ | |
5 | ||
6 | class TClonesArray; | |
7 | class AliEMCALGeometry; | |
8 | class AliVCluster; | |
9 | class AliPicoTrack; | |
b16bb001 | 10 | class AliVCaloCells; |
787a3c4f | 11 | class AliAODMCParticle; |
5be3857d | 12 | class AliNamedArrayI; |
43032ce2 | 13 | class TF2; |
65bb5510 | 14 | |
15 | #include <TH1F.h> | |
16 | #include <TF1.h> | |
762e8424 | 17 | |
18 | #include "AliAnalysisTaskSE.h" | |
19 | ||
20 | class AliJetModelBaseTask : public AliAnalysisTaskSE { | |
21 | public: | |
22 | AliJetModelBaseTask(); | |
b16bb001 | 23 | AliJetModelBaseTask(const char *name, Bool_t drawqa=kFALSE); |
762e8424 | 24 | virtual ~AliJetModelBaseTask(); |
25 | ||
a825589f | 26 | void SetEtaRange(Float_t min, Float_t max) { fEtaMin = min; fEtaMax = max; } |
a825589f | 27 | void SetPhiRange(Float_t min, Float_t max) { fPhiMin = min; fPhiMax = max; } |
43032ce2 | 28 | void SetPtRange(Float_t min, Float_t max) { fPtMin = min; fPtMax = max; } |
65bb5510 | 29 | void SetPtSpectrum(TH1 *f) { fPtSpectrum = f; } |
27f77648 | 30 | void SetPtSpectrum(TF1 *f) { fPtSpectrum = new TH1F("ptSpectrum","ptSpectrum",1000,f->GetXmin(),f->GetXmax()); |
65bb5510 | 31 | fPtSpectrum->Add(f); } |
43032ce2 | 32 | void SetPtPhiEvPlDistribution(TF2 *f) { fPtPhiEvPlDistribution = f; } |
27f77648 | 33 | void SetDensitySpectrum(TH1 *f) { fDensitySpectrum = f; } |
34 | void SetDensitySpectrum(TF1 *f) { fDensitySpectrum = new TH1F("densitypectrum","densitypectrum",1000,f->GetXmin(),f->GetXmax()); | |
35 | fDensitySpectrum->Add(f); } | |
50da5cf2 | 36 | void SetDifferentialV2(TF1* f) { fDifferentialV2 = f; } |
37 | void SetAddV2(Bool_t b) { fAddV2 = b; } | |
38 | void SetAddFlowFluctuations(Bool_t b) { fFlowFluctuations = b;} | |
4358e58a | 39 | void SetMC(Bool_t a) { fIsMC = a ; } |
b16bb001 | 40 | void SetCopyArray(Bool_t copy) { fCopyArray = copy; } |
a825589f | 41 | void SetTracksName(const char *n) { fTracksName = n; } |
b16bb001 | 42 | void SetClusName(const char *n) { fCaloName = n; } |
43 | void SetCellsName(const char *n) { fCellsName = n; } | |
787a3c4f | 44 | void SetMCParticlesName(const char *n) { fMCParticlesName = n; } |
b16bb001 | 45 | void SetSuffix(const char *s) { fSuffix = s; } |
b16bb001 | 46 | void SetGeometryName(const char *n) { fGeomName = n; } |
f660c2d6 | 47 | void SetMarkMC(Int_t m) { fMarkMC = m; } |
b16bb001 | 48 | virtual void SetNClusters(Int_t n) { fNClusters = n; } |
49 | virtual void SetNCells(Int_t n) { fNCells = n; } | |
50 | virtual void SetNTracks(Int_t n) { fNTracks = n; } | |
51 | ||
762e8424 | 52 | protected: |
50da5cf2 | 53 | void UserExec(Option_t* /*option*/); |
54 | void UserCreateOutputObjects(); | |
b16bb001 | 55 | Int_t SetNumberOfOutCells(Int_t n); // set the number of cells |
56 | Int_t AddCell(Double_t e = -1, Double_t eta = -999, Double_t phi = -1); // add a cell; if values are -1 generate random parameters | |
50da5cf2 | 57 | Int_t AddCell(Double_t e, Int_t absId, Double_t time = 0, Int_t label=0); // add a cell with given energy, position and times |
787a3c4f | 58 | AliVCluster *AddCluster(Double_t e = -1, Double_t eta = -999, Double_t phi = -1, Int_t label=0); // add a cluster; if values are -1 generate random parameters |
50da5cf2 | 59 | AliVCluster *AddCluster(Double_t e, Int_t absId, Int_t label=0); // add a cluster with given energy and position |
60 | AliVCluster *AddCluster(AliVCluster *oc); // add a cluster (copy) | |
6fd2f1a9 | 61 | AliPicoTrack *AddTrack(Double_t pt = -1, Double_t eta = -999, Double_t phi = -1, Byte_t type=0, |
56bd3193 | 62 | Double_t etaemc=0, Double_t phiemc=0, Double_t ptemc=0, Bool_t ise=kFALSE, |
4d3b366f | 63 | Int_t label=0, Short_t charge=1, Double_t mass = 0.1396); // add a track; if values are -1 generate random parameters |
787a3c4f | 64 | AliAODMCParticle *AddMCParticle(AliAODMCParticle *part, Int_t origIndex); // add a MC particle |
50da5cf2 | 65 | void AddV2(Double_t &phi, Double_t &pt) const; |
b16bb001 | 66 | void CopyCells(); |
e44e8726 | 67 | void CopyClusters(); |
68 | void CopyTracks(); | |
787a3c4f | 69 | void CopyMCParticles(); |
ffe32451 | 70 | void GetRandomCell(Double_t &eta, Double_t &phi, Int_t &absId); // generate a random cell in the calorimeter |
2103dc6a | 71 | Double_t GetRandomEta(Bool_t emcal=kFALSE); // generate a random eta value in the given range |
72 | Double_t GetRandomPhi(Bool_t emcal=kFALSE); // generate a random phi value in the given range | |
ffe32451 | 73 | Double_t GetRandomPt(); // generate a random pt value in the given range |
43032ce2 | 74 | void GetRandomParticle(Double_t &pt, Double_t &eta, Double_t &phi, Bool_t emcal=kFALSE); // generate a particle with random eta,phi,pt values |
b16bb001 | 75 | virtual Bool_t ExecOnce(); // intialize task |
ffe32451 | 76 | virtual void Run(); // do jet model action |
762e8424 | 77 | |
78 | TString fGeomName; // EMCal geometry name | |
79 | TString fTracksName; // name of track collection | |
80 | TString fOutTracksName; // name of output track collection | |
81 | TString fCaloName; // name of calo cluster collection | |
82 | TString fOutCaloName; // name of output cluster collection | |
b16bb001 | 83 | TString fCellsName; // name of calo cells collection |
84 | TString fOutCellsName; // name of output cells collection | |
787a3c4f | 85 | TString fMCParticlesName; // name of MC particle collection |
86 | TString fOutMCParticlesName; // name of output MC particle collection | |
4358e58a | 87 | Bool_t fIsMC; // whether the current event is MC or not |
762e8424 | 88 | TString fSuffix; // suffix to add in the name of new collections |
89 | Float_t fEtaMin; // eta minimum value | |
90 | Float_t fEtaMax; // eta maximum value | |
91 | Float_t fPhiMin; // phi minimum value | |
92 | Float_t fPhiMax; // phi maximum value | |
93 | Float_t fPtMin; // pt minimum value | |
94 | Float_t fPtMax; // pt maximum value | |
95 | Bool_t fCopyArray; // whether or not the array will be copied to a new one before modelling | |
96 | Int_t fNClusters; // how many clusters are being processed | |
b16bb001 | 97 | Int_t fNCells; // how many cells are being processed |
762e8424 | 98 | Int_t fNTracks; // how many tracks are being processed |
f660c2d6 | 99 | Int_t fMarkMC; // which MC label is to be used (default=100) |
27f77648 | 100 | TH1 *fPtSpectrum; // pt spectrum to extract random pt values |
43032ce2 | 101 | TF2 *fPtPhiEvPlDistribution; // pt vs. (phi-psi) distribution to extract random pt/phi values |
27f77648 | 102 | TH1 *fDensitySpectrum; // particle density spectrum to extract random density values |
50da5cf2 | 103 | TF1 *fDifferentialV2; // v2 as function of pt |
104 | Bool_t fAddV2; // add v2 sampled from a tf1 | |
105 | Bool_t fFlowFluctuations; // introduce gaussian flow fluctuation | |
b16bb001 | 106 | Bool_t fQAhistos; // draw QA histograms |
43032ce2 | 107 | Double_t fPsi; //!simmetry plane for the elliptic flow |
4358e58a | 108 | Bool_t fIsInit; //!=true if initialized |
762e8424 | 109 | AliEMCALGeometry *fGeom; //!pointer to EMCal geometry |
6fd2f1a9 | 110 | Double_t fVertex[3]; //!event vertex |
762e8424 | 111 | TClonesArray *fClusters; //!cluster collection |
112 | TClonesArray *fOutClusters; //!output cluster collection | |
113 | TClonesArray *fTracks; //!track collection | |
114 | TClonesArray *fOutTracks; //!output track collection | |
b16bb001 | 115 | AliVCaloCells *fCaloCells; //!cells collection |
116 | AliVCaloCells *fOutCaloCells; //!output cells collection | |
117 | Int_t fAddedCells; //!number of added cells | |
787a3c4f | 118 | TClonesArray *fMCParticles; //!MC particles collection |
5be3857d | 119 | AliNamedArrayI *fMCParticlesMap; //!MC particles mapping |
787a3c4f | 120 | TClonesArray *fOutMCParticles; //!output MC particles collection |
5be3857d | 121 | AliNamedArrayI *fOutMCParticlesMap; //!MC particles mapping |
787a3c4f | 122 | Int_t fMCLabelShift; //!MC label shift |
b16bb001 | 123 | Bool_t fEsdMode; //!ESD/AOD mode |
124 | TList *fOutput; //!output list for QA histograms | |
762e8424 | 125 | |
126 | private: | |
127 | AliJetModelBaseTask(const AliJetModelBaseTask&); // not implemented | |
128 | AliJetModelBaseTask &operator=(const AliJetModelBaseTask&); // not implemented | |
129 | ||
50da5cf2 | 130 | ClassDef(AliJetModelBaseTask, 11) // Jet modelling task |
762e8424 | 131 | }; |
132 | #endif |