]>
Commit | Line | Data |
---|---|---|
762e8424 | 1 | #ifndef ALIJETMODELBASETASK_H |
2 | #define ALIJETMODELBASETASK_H | |
3 | ||
4 | // $Id$ | |
5 | ||
6 | class TClonesArray; | |
787a3c4f | 7 | class TH1I; |
762e8424 | 8 | class AliEMCALGeometry; |
9 | class AliVCluster; | |
10 | class AliPicoTrack; | |
b16bb001 | 11 | class AliVCaloCells; |
787a3c4f | 12 | class AliAODMCParticle; |
5be3857d | 13 | class AliNamedArrayI; |
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 | ||
ffe32451 | 26 | void UserExec(Option_t* /*option*/); |
b16bb001 | 27 | void UserCreateOutputObjects(); |
762e8424 | 28 | |
a825589f | 29 | void SetEtaRange(Float_t min, Float_t max) { fEtaMin = min; fEtaMax = max; } |
a825589f | 30 | void SetPhiRange(Float_t min, Float_t max) { fPhiMin = min; fPhiMax = max; } |
31 | void SetPtRange(Float_t min, Float_t max) { fPtMin = min; fPtMax = max; } | |
65bb5510 | 32 | void SetPtSpectrum(TH1 *f) { fPtSpectrum = f; } |
33 | void SetPtSpectrum(TF1 *f) { fPtSpectrum = new TH1F("ptSpectrum","ptSpectrum",250,f->GetXmin(),f->GetXmax()); | |
34 | fPtSpectrum->Add(f); } | |
b16bb001 | 35 | |
36 | void SetCopyArray(Bool_t copy) { fCopyArray = copy; } | |
a825589f | 37 | void SetTracksName(const char *n) { fTracksName = n; } |
b16bb001 | 38 | void SetClusName(const char *n) { fCaloName = n; } |
39 | void SetCellsName(const char *n) { fCellsName = n; } | |
787a3c4f | 40 | void SetMCParticlesName(const char *n) { fMCParticlesName = n; } |
b16bb001 | 41 | void SetSuffix(const char *s) { fSuffix = s; } |
42 | ||
43 | void SetGeometryName(const char *n) { fGeomName = n; } | |
f660c2d6 | 44 | void SetMarkMC(Int_t m) { fMarkMC = m; } |
b16bb001 | 45 | virtual void SetNClusters(Int_t n) { fNClusters = n; } |
46 | virtual void SetNCells(Int_t n) { fNCells = n; } | |
47 | virtual void SetNTracks(Int_t n) { fNTracks = n; } | |
48 | ||
762e8424 | 49 | |
50 | protected: | |
b16bb001 | 51 | Int_t SetNumberOfOutCells(Int_t n); // set the number of cells |
52 | 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 | |
787a3c4f | 53 | 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 |
54 | 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 | |
55 | AliVCluster *AddCluster(Double_t e, Int_t absId, Int_t label=0); // add a cluster with given energy and position | |
b16bb001 | 56 | AliPicoTrack *AddTrack(Double_t pt = -1, Double_t eta = -999, Double_t phi = -1, |
787a3c4f | 57 | Byte_t type=0, Double_t etaemc=0, Double_t phiemc=0, Bool_t ise=kFALSE, Int_t label=0); // add a track; if values are -1 generate random parameters |
58 | AliAODMCParticle *AddMCParticle(AliAODMCParticle *part, Int_t origIndex); // add a MC particle | |
b16bb001 | 59 | void CopyCells(); |
e44e8726 | 60 | void CopyClusters(); |
61 | void CopyTracks(); | |
787a3c4f | 62 | void CopyMCParticles(); |
ffe32451 | 63 | void GetRandomCell(Double_t &eta, Double_t &phi, Int_t &absId); // generate a random cell in the calorimeter |
2103dc6a | 64 | Double_t GetRandomEta(Bool_t emcal=kFALSE); // generate a random eta value in the given range |
65 | Double_t GetRandomPhi(Bool_t emcal=kFALSE); // generate a random phi value in the given range | |
ffe32451 | 66 | Double_t GetRandomPt(); // generate a random pt value in the given range |
b16bb001 | 67 | virtual Bool_t ExecOnce(); // intialize task |
ffe32451 | 68 | virtual void Run(); // do jet model action |
762e8424 | 69 | |
70 | TString fGeomName; // EMCal geometry name | |
71 | TString fTracksName; // name of track collection | |
72 | TString fOutTracksName; // name of output track collection | |
73 | TString fCaloName; // name of calo cluster collection | |
74 | TString fOutCaloName; // name of output cluster collection | |
b16bb001 | 75 | TString fCellsName; // name of calo cells collection |
76 | TString fOutCellsName; // name of output cells collection | |
787a3c4f | 77 | TString fMCParticlesName; // name of MC particle collection |
78 | TString fOutMCParticlesName; // name of output MC particle collection | |
762e8424 | 79 | TString fSuffix; // suffix to add in the name of new collections |
80 | Float_t fEtaMin; // eta minimum value | |
81 | Float_t fEtaMax; // eta maximum value | |
82 | Float_t fPhiMin; // phi minimum value | |
83 | Float_t fPhiMax; // phi maximum value | |
84 | Float_t fPtMin; // pt minimum value | |
85 | Float_t fPtMax; // pt maximum value | |
86 | Bool_t fCopyArray; // whether or not the array will be copied to a new one before modelling | |
87 | Int_t fNClusters; // how many clusters are being processed | |
b16bb001 | 88 | Int_t fNCells; // how many cells are being processed |
762e8424 | 89 | Int_t fNTracks; // how many tracks are being processed |
f660c2d6 | 90 | Int_t fMarkMC; // which MC label is to be used (default=100) |
65bb5510 | 91 | TH1 *fPtSpectrum; // pt spectrum parametrization to extract random pt values |
b16bb001 | 92 | Bool_t fQAhistos; // draw QA histograms |
ffe32451 | 93 | Bool_t fIsInit; //=true if initialized |
762e8424 | 94 | AliEMCALGeometry *fGeom; //!pointer to EMCal geometry |
95 | TClonesArray *fClusters; //!cluster collection | |
96 | TClonesArray *fOutClusters; //!output cluster collection | |
97 | TClonesArray *fTracks; //!track collection | |
98 | TClonesArray *fOutTracks; //!output track collection | |
b16bb001 | 99 | AliVCaloCells *fCaloCells; //!cells collection |
100 | AliVCaloCells *fOutCaloCells; //!output cells collection | |
101 | Int_t fAddedCells; //!number of added cells | |
787a3c4f | 102 | TClonesArray *fMCParticles; //!MC particles collection |
5be3857d | 103 | AliNamedArrayI *fMCParticlesMap; //!MC particles mapping |
787a3c4f | 104 | TClonesArray *fOutMCParticles; //!output MC particles collection |
5be3857d | 105 | AliNamedArrayI *fOutMCParticlesMap; //!MC particles mapping |
787a3c4f | 106 | Int_t fMCLabelShift; //!MC label shift |
b16bb001 | 107 | Bool_t fEsdMode; //!ESD/AOD mode |
108 | TList *fOutput; //!output list for QA histograms | |
762e8424 | 109 | |
110 | private: | |
111 | AliJetModelBaseTask(const AliJetModelBaseTask&); // not implemented | |
112 | AliJetModelBaseTask &operator=(const AliJetModelBaseTask&); // not implemented | |
113 | ||
787a3c4f | 114 | ClassDef(AliJetModelBaseTask, 7) // Jet modelling task |
762e8424 | 115 | }; |
116 | #endif |