]>
Commit | Line | Data |
---|---|---|
e0331fd9 | 1 | #ifndef AliAnalysisTaskPhiCorrelations_H |
2 | #define AliAnalysisTaskPhiCorrelations_H | |
3 | ||
4 | /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. * | |
5 | * See cxx source for full Copyright notice */ | |
6 | ||
7 | //////////////////////////////////////////////////////////////////////// | |
8 | // | |
9 | // Analysis class for Underlying Event studies w.r.t. leading track | |
10 | // | |
11 | // Look for correlations on the tranverse regions w.r.t | |
12 | // the leading track in the event | |
13 | // | |
14 | // This class needs input AODs. | |
15 | // The output is a list of analysis-specific containers. | |
16 | // | |
17 | // The AOD can be either connected to the InputEventHandler | |
18 | // for a chain of AOD files | |
19 | // or | |
20 | // to the OutputEventHandler | |
21 | // for a chain of ESD files, | |
22 | // in this case the class should be in the train after the jet-finder | |
23 | // | |
24 | // Authors: | |
e0331fd9 | 25 | // Jan Fiete Grosse-Oetringhaus |
e0331fd9 | 26 | // |
27 | //////////////////////////////////////////////////////////////////////// | |
28 | ||
29 | #include "AliAnalysisTask.h" | |
30 | #include "AliUEHist.h" | |
2a910c25 | 31 | #include "TString.h" |
a1c31636 | 32 | #include "AliVParticle.h" |
33 | #include "AliLog.h" | |
408d1ac9 | 34 | #include "THn.h" // in cxx file causes .../THn.h:257: error: conflicting declaration ‘typedef class THnT<float> THnF’ |
e0331fd9 | 35 | |
36 | class AliAODEvent; | |
37 | class AliAnalyseLeadingTrackUE; | |
38 | class AliInputEventHandler; | |
39 | class AliMCEvent; | |
40 | class AliMCEventHandler; | |
41 | class AliUEHistograms; | |
42 | class AliVParticle; | |
43 | class TH1D; | |
44 | class TObjArray; | |
2a910c25 | 45 | class AliEventPoolManager; |
85bfac17 | 46 | class AliESDEvent; |
e0331fd9 | 47 | |
408d1ac9 | 48 | |
e0331fd9 | 49 | class AliAnalysisTaskPhiCorrelations : public AliAnalysisTask |
50 | { | |
51 | public: | |
52 | AliAnalysisTaskPhiCorrelations(const char* name="AliAnalysisTaskPhiCorrelations"); | |
53 | virtual ~AliAnalysisTaskPhiCorrelations(); | |
54 | ||
55 | ||
56 | // Implementation of interace methods | |
57 | virtual void ConnectInputData(Option_t *); | |
58 | virtual void CreateOutputObjects(); | |
59 | virtual void Exec(Option_t *option); | |
60 | ||
61 | // Setters/Getters | |
62 | // general configuration | |
63 | virtual void SetDebugLevel( Int_t level ) { fDebug = level; } | |
64 | virtual void SetMode(Int_t mode) { fMode = mode; } | |
65 | virtual void SetReduceMemoryFootprint(Bool_t flag) { fReduceMemoryFootprint = flag; } | |
eed401dc | 66 | virtual void SetEventMixing(Bool_t flag) { fFillMixed = flag; } |
ac647b0f | 67 | virtual void SetMixingTracks(Int_t tracks) { fMixingTracks = tracks; } |
c3294f09 | 68 | virtual void SetCompareCentralities(Bool_t flag) { fCompareCentralities = flag; } |
1bba939a | 69 | virtual void SetTwoTrackEfficiencyStudy(Bool_t flag) { fTwoTrackEfficiencyStudy = flag; } |
d4b3dbfc | 70 | virtual void SetTwoTrackEfficiencyCut(Float_t value = 0.02) { fTwoTrackEfficiencyCut = value; } |
3bbad7c1 | 71 | virtual void SetUseVtxAxis(Int_t flag) { fUseVtxAxis = flag; } |
9894bedd | 72 | virtual void SetCourseCentralityBinning(Bool_t flag) { fCourseCentralityBinning = flag; } |
04af8d15 | 73 | virtual void SetSkipTrigger(Bool_t flag) { fSkipTrigger = flag; } |
1ccd8a0a | 74 | virtual void SetInjectedSignals(Bool_t flag) { fInjectedSignals = flag; } |
e0331fd9 | 75 | |
76 | // histogram settings | |
408d1ac9 | 77 | void SetEfficiencyCorrection(THnF* hist) { fEfficiencyCorrection = hist; } |
e0331fd9 | 78 | |
79 | // for event QA | |
80 | void SetTracksInVertex( Int_t val ){ fnTracksVertex = val; } | |
81 | void SetZVertex( Double_t val ) { fZVertex = val; } | |
82 | ||
83 | // track cuts | |
84 | void SetTrackEtaCut( Double_t val ) { fTrackEtaCut = val; } | |
9da2f080 | 85 | void SetOnlyOneEtaSide(Int_t flag) { fOnlyOneEtaSide = flag; } |
e0331fd9 | 86 | void SetPtMin(Double_t val) { fPtMin = val; } |
87 | void SetFilterBit( UInt_t val ) { fFilterBit = val; } | |
2a910c25 | 88 | |
e0331fd9 | 89 | void SetEventSelectionBit( UInt_t val ) { fSelectBit = val; } |
90 | void SetUseChargeHadrons( Bool_t val ) { fUseChargeHadrons = val; } | |
2a910c25 | 91 | void SetSelectCharge(Int_t selectCharge) { fSelectCharge = selectCharge; } |
7a77d480 | 92 | void SetSelectTriggerCharge(Int_t selectCharge) { fTriggerSelectCharge = selectCharge; } |
d38fa455 | 93 | void SetTriggerRestrictEta(Float_t eta) { fTriggerRestrictEta = eta; } |
00b6f3c6 | 94 | void SetEtaOrdering(Bool_t flag) { fEtaOrdering = flag; } |
b0d56b29 | 95 | void SetPairCuts(Bool_t conversions, Bool_t resonances) { fCutConversions = conversions; fCutResonances = resonances; } |
2a910c25 | 96 | void SetCentralityMethod(const char* method) { fCentralityMethod = method; } |
c05ff6be | 97 | void SetFillpT(Bool_t flag) { fFillpT = flag; } |
a26093ba | 98 | void SetStepsFillSkip(Bool_t step0, Bool_t step6) { fFillOnlyStep0 = step0; fSkipStep6 = step6; } |
5e053cad | 99 | void SetRejectCentralityOutliers(Bool_t flag = kTRUE) { fRejectCentralityOutliers = flag; } |
f613255f | 100 | void SetRemoveWeakDecays(Bool_t flag = kTRUE) { fRemoveWeakDecays = flag; } |
e0331fd9 | 101 | |
102 | private: | |
103 | AliAnalysisTaskPhiCorrelations(const AliAnalysisTaskPhiCorrelations &det); | |
104 | AliAnalysisTaskPhiCorrelations& operator=(const AliAnalysisTaskPhiCorrelations &det); | |
105 | void AddSettingsTree(); // add list of settings to output list | |
106 | // Analysis methods | |
107 | void AnalyseCorrectionMode(); // main algorithm to get correction maps | |
108 | void AnalyseDataMode(); // main algorithm to get raw distributions | |
109 | void Initialize(); // initialize some common pointer | |
7a028750 | 110 | TObjArray* CloneAndReduceTrackList(TObjArray* tracks); |
e0331fd9 | 111 | |
e0331fd9 | 112 | // General configuration |
113 | Int_t fDebug; // Debug flag | |
114 | Int_t fMode; // fMode = 0: data-like analysis | |
115 | // fMode = 1: corrections analysis | |
116 | Bool_t fReduceMemoryFootprint; // reduce memory consumption by writing less debug histograms | |
eed401dc | 117 | Bool_t fFillMixed; // enable event mixing (default: ON) |
ac647b0f | 118 | Int_t fMixingTracks; // size of track buffer for event mixing |
c3294f09 | 119 | Bool_t fCompareCentralities; // use the z vtx axis for a centrality comparison |
1bba939a | 120 | Bool_t fTwoTrackEfficiencyStudy; // two-track efficiency study on |
d4b3dbfc | 121 | Float_t fTwoTrackEfficiencyCut; // enable two-track efficiency cut |
3bbad7c1 | 122 | Int_t fUseVtxAxis; // use z vtx as axis (needs 7-10 times more memory!) |
9894bedd | 123 | Bool_t fCourseCentralityBinning; // less centrality bins |
04af8d15 | 124 | Bool_t fSkipTrigger; // skip trigger selection |
1ccd8a0a | 125 | Bool_t fInjectedSignals; // check header to skip injected signals in MC |
e0331fd9 | 126 | |
127 | // Pointers to external UE classes | |
128 | AliAnalyseLeadingTrackUE* fAnalyseUE; //! points to class containing common analysis algorithms | |
129 | AliUEHistograms* fHistos; //! points to class to handle histograms/containers | |
130 | AliUEHistograms* fHistosMixed; //! points to class to handle mixed histograms/containers | |
131 | ||
408d1ac9 | 132 | THnF* fEfficiencyCorrection; // if non-0 this efficiency correction is applied on the fly to the filling for associated particles. The factor is multiplicative, i.e. should contain 1/efficiency. Axes: eta, pT, centrality, z-vtx |
e0331fd9 | 133 | |
134 | // Handlers and events | |
135 | AliAODEvent* fAOD; //! AOD Event | |
85bfac17 | 136 | AliESDEvent* fESD; //! ESD Event |
e0331fd9 | 137 | TClonesArray* fArrayMC; //! Array of MC particles |
138 | AliInputEventHandler* fInputHandler; //! Generic InputEventHandler | |
139 | AliMCEvent* fMcEvent; //! MC event | |
140 | AliMCEventHandler* fMcHandler; //! MCEventHandler | |
2a910c25 | 141 | AliEventPoolManager* fPoolMgr; //! event pool manager |
e0331fd9 | 142 | |
143 | // Histogram settings | |
144 | TList* fListOfHistos; // Output list of containers | |
145 | ||
146 | // Event QA cuts | |
c05ff6be | 147 | Int_t fnTracksVertex; // QA tracks pointing to principal vertex |
e0331fd9 | 148 | Double_t fZVertex; // Position of Vertex in Z direction |
2a910c25 | 149 | TString fCentralityMethod; // Method to determine centrality |
e0331fd9 | 150 | |
151 | // Track cuts | |
152 | Double_t fTrackEtaCut; // Eta cut on particles | |
9da2f080 | 153 | Int_t fOnlyOneEtaSide; // decides that only trigger particle from one eta side are considered (0 = all; -1 = negative, 1 = positive) |
e0331fd9 | 154 | Double_t fPtMin; // Min pT to start correlations |
155 | UInt_t fFilterBit; // Select tracks from an specific track cut (default 0xFF all track selected) | |
156 | UInt_t fSelectBit; // Select events according to AliAnalysisTaskJetServices bit maps | |
157 | Bool_t fUseChargeHadrons; // Only use charge hadrons | |
158 | ||
2a910c25 | 159 | Int_t fSelectCharge; // (un)like sign selection when building correlations: 0: no selection; 1: unlike sign; 2: like sign |
7a77d480 | 160 | Int_t fTriggerSelectCharge; // select charge of trigger particle |
d38fa455 | 161 | Float_t fTriggerRestrictEta; // restrict eta range for trigger particle (default: -1 [off]) |
00b6f3c6 | 162 | Bool_t fEtaOrdering; // eta ordering, see AliUEHistograms.h for documentation |
b0d56b29 | 163 | Bool_t fCutConversions; // cut on conversions (inv mass) |
164 | Bool_t fCutResonances; // cut on resonances (inv mass) | |
a26093ba | 165 | Bool_t fFillOnlyStep0; // fill only step 0 |
166 | Bool_t fSkipStep6; // skip step 6 when filling | |
5e053cad | 167 | Bool_t fRejectCentralityOutliers; // enable rejection of outliers in centrality vs no track correlation |
f613255f | 168 | Bool_t fRemoveWeakDecays; // remove secondaries from weak decays from tracks and particles |
d38fa455 | 169 | |
c05ff6be | 170 | Bool_t fFillpT; // fill sum pT instead of number density |
2a910c25 | 171 | |
f613255f | 172 | ClassDef( AliAnalysisTaskPhiCorrelations, 18); // Analysis task for delta phi correlations |
e0331fd9 | 173 | }; |
174 | ||
a1c31636 | 175 | class AliDPhiBasicParticle : public AliVParticle |
176 | { | |
177 | public: | |
178 | AliDPhiBasicParticle(Float_t eta, Float_t phi, Float_t pt, Short_t charge) | |
179 | : fEta(eta), fPhi(phi), fpT(pt), fCharge(charge) | |
180 | { | |
181 | } | |
182 | ~AliDPhiBasicParticle() {} | |
183 | ||
184 | // kinematics | |
185 | virtual Double_t Px() const { AliFatal("Not implemented"); return 0; } | |
186 | virtual Double_t Py() const { AliFatal("Not implemented"); return 0; } | |
187 | virtual Double_t Pz() const { AliFatal("Not implemented"); return 0; } | |
188 | virtual Double_t Pt() const { return fpT; } | |
189 | virtual Double_t P() const { AliFatal("Not implemented"); return 0; } | |
190 | virtual Bool_t PxPyPz(Double_t[3]) const { AliFatal("Not implemented"); return 0; } | |
191 | ||
192 | virtual Double_t Xv() const { AliFatal("Not implemented"); return 0; } | |
193 | virtual Double_t Yv() const { AliFatal("Not implemented"); return 0; } | |
194 | virtual Double_t Zv() const { AliFatal("Not implemented"); return 0; } | |
195 | virtual Bool_t XvYvZv(Double_t[3]) const { AliFatal("Not implemented"); return 0; } | |
196 | ||
197 | virtual Double_t OneOverPt() const { AliFatal("Not implemented"); return 0; } | |
198 | virtual Double_t Phi() const { return fPhi; } | |
199 | virtual Double_t Theta() const { AliFatal("Not implemented"); return 0; } | |
200 | ||
201 | ||
202 | virtual Double_t E() const { AliFatal("Not implemented"); return 0; } | |
203 | virtual Double_t M() const { AliFatal("Not implemented"); return 0; } | |
204 | ||
205 | virtual Double_t Eta() const { return fEta; } | |
206 | virtual Double_t Y() const { AliFatal("Not implemented"); return 0; } | |
207 | ||
208 | virtual Short_t Charge() const { return fCharge; } | |
209 | virtual Int_t GetLabel() const { AliFatal("Not implemented"); return 0; } | |
210 | // PID | |
211 | virtual Int_t PdgCode() const { AliFatal("Not implemented"); return 0; } | |
212 | virtual const Double_t *PID() const { AliFatal("Not implemented"); return 0; } | |
213 | ||
214 | private: | |
215 | Float_t fEta; // eta | |
216 | Float_t fPhi; // phi | |
217 | Float_t fpT; // pT | |
218 | Short_t fCharge; // charge | |
219 | ||
220 | ClassDef( AliDPhiBasicParticle, 1); // class which contains only quantities requires for this analysis to reduce memory consumption for event mixing | |
221 | }; | |
222 | ||
e0331fd9 | 223 | #endif |
224 | ||
225 |