]>
Commit | Line | Data |
---|---|---|
a9e585a7 | 1 | |
568f8fa2 | 2 | // ****************************************** |
3 | // This task computes several jet observables like | |
4 | // the fraction of energy in inner and outer coronnas, | |
5 | // the distance from track to jet axis and a | |
6 | // correlation strength distribution of particles inside jets. | |
7 | // Author: lcunquei@cern.ch | |
8 | // ******************************************* | |
9 | ||
10 | ||
11 | /************************************************************************** | |
12 | * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. * | |
13 | * * | |
14 | * Author: The ALICE Off-line Project. * | |
15 | * Contributors are mentioned in the code where appropriate. * | |
16 | * * | |
17 | * Permission to use, copy, modify and distribute this software and its * | |
18 | * documentation strictly for non-commercial purposes is hereby granted * | |
19 | * without fee, provided that the above copyright notice appears in all * | |
20 | * copies and that both the copyright notice and this permission notice * | |
21 | * appear in the supporting documentation. The authors make no claims * | |
22 | * about the suitability of this software for any purpose. It is * | |
23 | * provided "as is" without express or implied warranty. * | |
24 | **************************************************************************/ | |
25 | ||
26 | ||
75bf77e3 | 27 | #include "TChain.h" |
28 | #include "TTree.h" | |
29 | #include "TMath.h" | |
30 | #include "TH1F.h" | |
31 | #include "TH2F.h" | |
32 | #include "TH3F.h" | |
33 | #include "THnSparse.h" | |
34 | #include "TCanvas.h" | |
35 | ||
36 | #include "AliLog.h" | |
37 | ||
38 | #include "AliAnalysisTask.h" | |
39 | #include "AliAnalysisManager.h" | |
40 | ||
41 | #include "AliVEvent.h" | |
42 | #include "AliESDEvent.h" | |
43 | #include "AliESDInputHandler.h" | |
44 | #include "AliCentrality.h" | |
45 | #include "AliAnalysisHelperJetTasks.h" | |
46 | #include "AliInputEventHandler.h" | |
47 | #include "AliAODJetEventBackground.h" | |
48 | #include "AliAnalysisTaskFastEmbedding.h" | |
49 | ||
50 | #include "AliAODEvent.h" | |
ea693273 | 51 | #include "AliAODHandler.h" |
75bf77e3 | 52 | #include "AliAODJet.h" |
53 | ||
54 | #include "AliAnalysisTaskJetCore.h" | |
55 | ||
56 | ClassImp(AliAnalysisTaskJetCore) | |
57 | ||
58 | AliAnalysisTaskJetCore::AliAnalysisTaskJetCore() : | |
59 | AliAnalysisTaskSE(), | |
60 | fESD(0x0), | |
61 | fAOD(0x0), | |
ea693273 | 62 | fAODExtension(0x0), |
75bf77e3 | 63 | fBackgroundBranch(""), |
ea693273 | 64 | fNonStdFile(""), |
75bf77e3 | 65 | fIsPbPb(kTRUE), |
66 | fOfflineTrgMask(AliVEvent::kAny), | |
67 | fMinContribVtx(1), | |
a9e585a7 | 68 | fVtxZMin(-10.), |
69 | fVtxZMax(10.), | |
75bf77e3 | 70 | fEvtClassMin(0), |
71 | fEvtClassMax(4), | |
ea693273 | 72 | fRadioFrac(0.2), |
73 | fMinDist(0.1), | |
75bf77e3 | 74 | fCentMin(0.), |
75 | fCentMax(100.), | |
76 | fNInputTracksMin(0), | |
77 | fNInputTracksMax(-1), | |
ea693273 | 78 | fAngStructCloseTracks(0), |
75bf77e3 | 79 | fJetEtaMin(-.5), |
80 | fJetEtaMax(.5), | |
81 | fJetPtMin(20.), | |
82 | fJetTriggerExcludeMask(AliAODJet::kHighTrackPtTriggered), | |
83 | fJetPtFractionMin(0.5), | |
84 | fNMatchJets(4), | |
85 | fMatchMaxDist(0.8), | |
86 | fKeepJets(kFALSE), | |
75bf77e3 | 87 | fkNbranches(2), |
88 | fkEvtClasses(12), | |
89 | fOutputList(0x0), | |
90 | fbEvent(kTRUE), | |
91 | fHistEvtSelection(0x0), | |
92 | fHistJetSelection(0x0), | |
93 | fh2JetSelection(0x0), | |
94 | fh2JetCoreMethod1C10(0x0), | |
95 | fh2JetCoreMethod2C10(0x0), | |
96 | fh2JetCoreMethod3C10(0x0), | |
ea693273 | 97 | fh2JetCoreMethod1C20(0x0), |
98 | fh2JetCoreMethod2C20(0x0), | |
99 | fh2JetCoreMethod3C20(0x0), | |
75bf77e3 | 100 | fh2JetCoreMethod1C30(0x0), |
101 | fh2JetCoreMethod2C30(0x0), | |
102 | fh2JetCoreMethod3C30(0x0), | |
103 | fh2JetCoreMethod1C60(0x0), | |
104 | fh2JetCoreMethod2C60(0x0), | |
ea693273 | 105 | fh2JetCoreMethod3C60(0x0), |
a9e585a7 | 106 | fh2JetCoreMethod3C10lead(0x0), |
107 | fh2JetCoreMethod3C20lead(0x0), | |
108 | fh2JetCoreMethod3C30lead(0x0), | |
109 | fh2JetCoreMethod3C60lead(0x0), | |
110 | fh2JetCoreMethod3C10sublead(0x0), | |
111 | fh2JetCoreMethod3C20sublead(0x0), | |
112 | fh2JetCoreMethod3C30sublead(0x0), | |
113 | fh2JetCoreMethod3C60sublead(0x0), | |
ea693273 | 114 | fh2SumPtInC10(0x0), |
115 | fh2SumPtInC20(0x0), | |
116 | fh2SumPtInC30(0x0), | |
117 | fh2SumPtInC60(0x0), | |
a9e585a7 | 118 | fh2SumPtInC10lead(0x0), |
119 | fh2SumPtInC20lead(0x0), | |
120 | fh2SumPtInC30lead(0x0), | |
121 | fh2SumPtInC60lead(0x0), | |
122 | fh2SumPtInC10sublead(0x0), | |
123 | fh2SumPtInC20sublead(0x0), | |
124 | fh2SumPtInC30sublead(0x0), | |
125 | fh2SumPtInC60sublead(0x0), | |
ea693273 | 126 | fh2SumPtOutC10(0x0), |
ea693273 | 127 | fh2SumPtOutC20(0x0), |
128 | fh2SumPtOutC30(0x0), | |
129 | fh2SumPtOutC60(0x0), | |
a9e585a7 | 130 | fh2SumPtOutC10lead(0x0), |
131 | fh2SumPtOutC20lead(0x0), | |
132 | fh2SumPtOutC30lead(0x0), | |
133 | fh2SumPtOutC60lead(0x0), | |
134 | fh2SumPtOutC10sublead(0x0), | |
135 | fh2SumPtOutC20sublead(0x0), | |
136 | fh2SumPtOutC30sublead(0x0), | |
137 | fh2SumPtOutC60sublead(0x0), | |
ea693273 | 138 | fh2SumPtInC10bkg(0x0), |
139 | fh2SumPtInC20bkg(0x0), | |
140 | fh2SumPtInC30bkg(0x0), | |
141 | fh2SumPtInC60bkg(0x0), | |
a9e585a7 | 142 | fh2SumPtInC10bkglead(0x0), |
143 | fh2SumPtInC20bkglead(0x0), | |
144 | fh2SumPtInC30bkglead(0x0), | |
145 | fh2SumPtInC60bkglead(0x0), | |
146 | fh2SumPtInC10bkgsublead(0x0), | |
147 | fh2SumPtInC20bkgsublead(0x0), | |
148 | fh2SumPtInC30bkgsublead(0x0), | |
149 | fh2SumPtInC60bkgsublead(0x0), | |
ea693273 | 150 | fh2SumPtOutC10bkg(0x0), |
151 | fh2SumPtOutC20bkg(0x0), | |
152 | fh2SumPtOutC30bkg(0x0), | |
153 | fh2SumPtOutC60bkg(0x0), | |
a9e585a7 | 154 | fh2SumPtOutC10bkglead(0x0), |
155 | fh2SumPtOutC20bkglead(0x0), | |
156 | fh2SumPtOutC30bkglead(0x0), | |
157 | fh2SumPtOutC60bkglead(0x0), | |
158 | fh2SumPtOutC10bkgsublead(0x0), | |
159 | fh2SumPtOutC20bkgsublead(0x0), | |
160 | fh2SumPtOutC30bkgsublead(0x0), | |
161 | fh2SumPtOutC60bkgsublead(0x0), | |
ea693273 | 162 | fh2DeltaRC10pt1(0x0), |
163 | fh2DeltaRC20pt1(0x0), | |
164 | fh2DeltaRC30pt1(0x0), | |
165 | fh2DeltaRC60pt1(0x0), | |
166 | fh2DeltaRC10pt2(0x0), | |
167 | fh2DeltaRC20pt2(0x0), | |
168 | fh2DeltaRC30pt2(0x0), | |
169 | fh2DeltaRC60pt2(0x0), | |
170 | fh2DeltaRC10pt3(0x0), | |
171 | fh2DeltaRC20pt3(0x0), | |
172 | fh2DeltaRC30pt3(0x0), | |
173 | fh2DeltaRC60pt3(0x0), | |
174 | fh2DeltaRC10pt4(0x0), | |
175 | fh2DeltaRC20pt4(0x0), | |
176 | fh2DeltaRC30pt4(0x0), | |
177 | fh2DeltaRC60pt4(0x0), | |
178 | fh2DeltaEtaC10pt1(0x0), | |
179 | fh2DeltaEtaC20pt1(0x0), | |
180 | fh2DeltaEtaC30pt1(0x0), | |
181 | fh2DeltaEtaC60pt1(0x0), | |
182 | fh2DeltaEtaC10pt2(0x0), | |
183 | fh2DeltaEtaC20pt2(0x0), | |
184 | fh2DeltaEtaC30pt2(0x0), | |
185 | fh2DeltaEtaC60pt2(0x0), | |
186 | fh2DeltaEtaC10pt3(0x0), | |
187 | fh2DeltaEtaC20pt3(0x0), | |
188 | fh2DeltaEtaC30pt3(0x0), | |
189 | fh2DeltaEtaC60pt3(0x0), | |
190 | fh2DeltaEtaC10pt4(0x0), | |
191 | fh2DeltaEtaC20pt4(0x0), | |
192 | fh2DeltaEtaC30pt4(0x0), | |
193 | fh2DeltaEtaC60pt4(0x0), | |
194 | fh2DeltaPhiC10pt1(0x0), | |
195 | fh2DeltaPhiC20pt1(0x0), | |
196 | fh2DeltaPhiC30pt1(0x0), | |
197 | fh2DeltaPhiC60pt1(0x0), | |
198 | fh2DeltaPhiC10pt2(0x0), | |
199 | fh2DeltaPhiC20pt2(0x0), | |
200 | fh2DeltaPhiC30pt2(0x0), | |
201 | fh2DeltaPhiC60pt2(0x0), | |
202 | fh2DeltaPhiC10pt3(0x0), | |
203 | fh2DeltaPhiC20pt3(0x0), | |
204 | fh2DeltaPhiC30pt3(0x0), | |
205 | fh2DeltaPhiC60pt3(0x0), | |
206 | fh2DeltaPhiC10pt4(0x0), | |
207 | fh2DeltaPhiC20pt4(0x0), | |
208 | fh2DeltaPhiC30pt4(0x0), | |
209 | fh2DeltaPhiC60pt4(0x0), | |
a9e585a7 | 210 | fh2DeltaRC10pt1lead(0x0), |
211 | fh2DeltaRC20pt1lead(0x0), | |
212 | fh2DeltaRC30pt1lead(0x0), | |
213 | fh2DeltaRC60pt1lead(0x0), | |
214 | fh2DeltaRC10pt2lead(0x0), | |
215 | fh2DeltaRC20pt2lead(0x0), | |
216 | fh2DeltaRC30pt2lead(0x0), | |
217 | fh2DeltaRC60pt2lead(0x0), | |
218 | fh2DeltaRC10pt3lead(0x0), | |
219 | fh2DeltaRC20pt3lead(0x0), | |
220 | fh2DeltaRC30pt3lead(0x0), | |
221 | fh2DeltaRC60pt3lead(0x0), | |
222 | fh2DeltaRC10pt4lead(0x0), | |
223 | fh2DeltaRC20pt4lead(0x0), | |
224 | fh2DeltaRC30pt4lead(0x0), | |
225 | fh2DeltaRC60pt4lead(0x0), | |
226 | fh2DeltaEtaC10pt1lead(0x0), | |
227 | fh2DeltaEtaC20pt1lead(0x0), | |
228 | fh2DeltaEtaC30pt1lead(0x0), | |
229 | fh2DeltaEtaC60pt1lead(0x0), | |
230 | fh2DeltaEtaC10pt2lead(0x0), | |
231 | fh2DeltaEtaC20pt2lead(0x0), | |
232 | fh2DeltaEtaC30pt2lead(0x0), | |
233 | fh2DeltaEtaC60pt2lead(0x0), | |
234 | fh2DeltaEtaC10pt3lead(0x0), | |
235 | fh2DeltaEtaC20pt3lead(0x0), | |
236 | fh2DeltaEtaC30pt3lead(0x0), | |
237 | fh2DeltaEtaC60pt3lead(0x0), | |
238 | fh2DeltaEtaC10pt4lead(0x0), | |
239 | fh2DeltaEtaC20pt4lead(0x0), | |
240 | fh2DeltaEtaC30pt4lead(0x0), | |
241 | fh2DeltaEtaC60pt4lead(0x0), | |
242 | fh2DeltaPhiC10pt1lead(0x0), | |
243 | fh2DeltaPhiC20pt1lead(0x0), | |
244 | fh2DeltaPhiC30pt1lead(0x0), | |
245 | fh2DeltaPhiC60pt1lead(0x0), | |
246 | fh2DeltaPhiC10pt2lead(0x0), | |
247 | fh2DeltaPhiC20pt2lead(0x0), | |
248 | fh2DeltaPhiC30pt2lead(0x0), | |
249 | fh2DeltaPhiC60pt2lead(0x0), | |
250 | fh2DeltaPhiC10pt3lead(0x0), | |
251 | fh2DeltaPhiC20pt3lead(0x0), | |
252 | fh2DeltaPhiC30pt3lead(0x0), | |
253 | fh2DeltaPhiC60pt3lead(0x0), | |
254 | fh2DeltaPhiC10pt4lead(0x0), | |
255 | fh2DeltaPhiC20pt4lead(0x0), | |
256 | fh2DeltaPhiC30pt4lead(0x0), | |
257 | fh2DeltaPhiC60pt4lead(0x0), | |
258 | fh2DeltaRC10pt1sublead(0x0), | |
259 | fh2DeltaRC20pt1sublead(0x0), | |
260 | fh2DeltaRC30pt1sublead(0x0), | |
261 | fh2DeltaRC60pt1sublead(0x0), | |
262 | fh2DeltaRC10pt2sublead(0x0), | |
263 | fh2DeltaRC20pt2sublead(0x0), | |
264 | fh2DeltaRC30pt2sublead(0x0), | |
265 | fh2DeltaRC60pt2sublead(0x0), | |
266 | fh2DeltaRC10pt3sublead(0x0), | |
267 | fh2DeltaRC20pt3sublead(0x0), | |
268 | fh2DeltaRC30pt3sublead(0x0), | |
269 | fh2DeltaRC60pt3sublead(0x0), | |
270 | fh2DeltaRC10pt4sublead(0x0), | |
271 | fh2DeltaRC20pt4sublead(0x0), | |
272 | fh2DeltaRC30pt4sublead(0x0), | |
273 | fh2DeltaRC60pt4sublead(0x0), | |
274 | fh2DeltaEtaC10pt1sublead(0x0), | |
275 | fh2DeltaEtaC20pt1sublead(0x0), | |
276 | fh2DeltaEtaC30pt1sublead(0x0), | |
277 | fh2DeltaEtaC60pt1sublead(0x0), | |
278 | fh2DeltaEtaC10pt2sublead(0x0), | |
279 | fh2DeltaEtaC20pt2sublead(0x0), | |
280 | fh2DeltaEtaC30pt2sublead(0x0), | |
281 | fh2DeltaEtaC60pt2sublead(0x0), | |
282 | fh2DeltaEtaC10pt3sublead(0x0), | |
283 | fh2DeltaEtaC20pt3sublead(0x0), | |
284 | fh2DeltaEtaC30pt3sublead(0x0), | |
285 | fh2DeltaEtaC60pt3sublead(0x0), | |
286 | fh2DeltaEtaC10pt4sublead(0x0), | |
287 | fh2DeltaEtaC20pt4sublead(0x0), | |
288 | fh2DeltaEtaC30pt4sublead(0x0), | |
289 | fh2DeltaEtaC60pt4sublead(0x0), | |
290 | fh2DeltaPhiC10pt1sublead(0x0), | |
291 | fh2DeltaPhiC20pt1sublead(0x0), | |
292 | fh2DeltaPhiC30pt1sublead(0x0), | |
293 | fh2DeltaPhiC60pt1sublead(0x0), | |
294 | fh2DeltaPhiC10pt2sublead(0x0), | |
295 | fh2DeltaPhiC20pt2sublead(0x0), | |
296 | fh2DeltaPhiC30pt2sublead(0x0), | |
297 | fh2DeltaPhiC60pt2sublead(0x0), | |
298 | fh2DeltaPhiC10pt3sublead(0x0), | |
299 | fh2DeltaPhiC20pt3sublead(0x0), | |
300 | fh2DeltaPhiC30pt3sublead(0x0), | |
301 | fh2DeltaPhiC60pt3sublead(0x0), | |
302 | fh2DeltaPhiC10pt4sublead(0x0), | |
303 | fh2DeltaPhiC20pt4sublead(0x0), | |
304 | fh2DeltaPhiC30pt4sublead(0x0), | |
305 | fh2DeltaPhiC60pt4sublead(0x0), | |
ea693273 | 306 | fh2AngStructpt1C10(0x0), |
307 | fh2AngStructpt2C10(0x0), | |
308 | fh2AngStructpt3C10(0x0), | |
309 | fh2AngStructpt4C10(0x0), | |
310 | fh2AngStructpt1C20(0x0), | |
311 | fh2AngStructpt2C20(0x0), | |
312 | fh2AngStructpt3C20(0x0), | |
313 | fh2AngStructpt4C20(0x0), | |
314 | fh2AngStructpt1C30(0x0), | |
315 | fh2AngStructpt2C30(0x0), | |
316 | fh2AngStructpt3C30(0x0), | |
317 | fh2AngStructpt4C30(0x0), | |
318 | fh2AngStructpt1C60(0x0), | |
319 | fh2AngStructpt2C60(0x0), | |
320 | fh2AngStructpt3C60(0x0), | |
321 | fh2AngStructpt4C60(0x0) | |
75bf77e3 | 322 | { |
323 | // default Constructor | |
324 | ||
325 | fJetBranchName[0] = ""; | |
326 | fJetBranchName[1] = ""; | |
327 | ||
328 | fListJets[0] = new TList; | |
329 | fListJets[1] = new TList; | |
330 | } | |
331 | ||
332 | AliAnalysisTaskJetCore::AliAnalysisTaskJetCore(const char *name) : | |
333 | AliAnalysisTaskSE(name), | |
334 | fESD(0x0), | |
335 | fAOD(0x0), | |
ea693273 | 336 | fAODExtension(0x0), |
75bf77e3 | 337 | fBackgroundBranch(""), |
ea693273 | 338 | fNonStdFile(""), |
75bf77e3 | 339 | fIsPbPb(kTRUE), |
340 | fOfflineTrgMask(AliVEvent::kAny), | |
341 | fMinContribVtx(1), | |
a9e585a7 | 342 | fVtxZMin(-10.), |
343 | fVtxZMax(10.), | |
75bf77e3 | 344 | fEvtClassMin(0), |
345 | fEvtClassMax(4), | |
ea693273 | 346 | fRadioFrac(0.2), |
347 | fMinDist(0.1), | |
75bf77e3 | 348 | fCentMin(0.), |
349 | fCentMax(100.), | |
350 | fNInputTracksMin(0), | |
351 | fNInputTracksMax(-1), | |
ea693273 | 352 | fAngStructCloseTracks(0), |
75bf77e3 | 353 | fJetEtaMin(-.5), |
354 | fJetEtaMax(.5), | |
355 | fJetPtMin(20.), | |
356 | fJetTriggerExcludeMask(AliAODJet::kHighTrackPtTriggered), | |
357 | fJetPtFractionMin(0.5), | |
358 | fNMatchJets(4), | |
359 | fMatchMaxDist(0.8), | |
360 | fKeepJets(kFALSE), | |
75bf77e3 | 361 | fkNbranches(2), |
362 | fkEvtClasses(12), | |
363 | fOutputList(0x0), | |
364 | fbEvent(kTRUE), | |
365 | fHistEvtSelection(0x0), | |
366 | fHistJetSelection(0x0), | |
367 | fh2JetSelection(0x0), | |
368 | fh2JetCoreMethod1C10(0x0), | |
369 | fh2JetCoreMethod2C10(0x0), | |
370 | fh2JetCoreMethod3C10(0x0), | |
ea693273 | 371 | fh2JetCoreMethod1C20(0x0), |
372 | fh2JetCoreMethod2C20(0x0), | |
373 | fh2JetCoreMethod3C20(0x0), | |
75bf77e3 | 374 | fh2JetCoreMethod1C30(0x0), |
375 | fh2JetCoreMethod2C30(0x0), | |
376 | fh2JetCoreMethod3C30(0x0), | |
377 | fh2JetCoreMethod1C60(0x0), | |
378 | fh2JetCoreMethod2C60(0x0), | |
ea693273 | 379 | fh2JetCoreMethod3C60(0x0), |
a9e585a7 | 380 | fh2JetCoreMethod3C10lead(0x0), |
381 | fh2JetCoreMethod3C20lead(0x0), | |
382 | fh2JetCoreMethod3C30lead(0x0), | |
383 | fh2JetCoreMethod3C60lead(0x0), | |
384 | fh2JetCoreMethod3C10sublead(0x0), | |
385 | fh2JetCoreMethod3C20sublead(0x0), | |
386 | fh2JetCoreMethod3C30sublead(0x0), | |
387 | fh2JetCoreMethod3C60sublead(0x0), | |
ea693273 | 388 | fh2SumPtInC10(0x0), |
389 | fh2SumPtInC20(0x0), | |
390 | fh2SumPtInC30(0x0), | |
391 | fh2SumPtInC60(0x0), | |
a9e585a7 | 392 | fh2SumPtInC10lead(0x0), |
393 | fh2SumPtInC20lead(0x0), | |
394 | fh2SumPtInC30lead(0x0), | |
395 | fh2SumPtInC60lead(0x0), | |
396 | fh2SumPtInC10sublead(0x0), | |
397 | fh2SumPtInC20sublead(0x0), | |
398 | fh2SumPtInC30sublead(0x0), | |
399 | fh2SumPtInC60sublead(0x0), | |
ea693273 | 400 | fh2SumPtOutC10(0x0), |
ea693273 | 401 | fh2SumPtOutC20(0x0), |
402 | fh2SumPtOutC30(0x0), | |
403 | fh2SumPtOutC60(0x0), | |
a9e585a7 | 404 | fh2SumPtOutC10lead(0x0), |
405 | fh2SumPtOutC20lead(0x0), | |
406 | fh2SumPtOutC30lead(0x0), | |
407 | fh2SumPtOutC60lead(0x0), | |
408 | fh2SumPtOutC10sublead(0x0), | |
409 | fh2SumPtOutC20sublead(0x0), | |
410 | fh2SumPtOutC30sublead(0x0), | |
411 | fh2SumPtOutC60sublead(0x0), | |
ea693273 | 412 | fh2SumPtInC10bkg(0x0), |
413 | fh2SumPtInC20bkg(0x0), | |
414 | fh2SumPtInC30bkg(0x0), | |
415 | fh2SumPtInC60bkg(0x0), | |
a9e585a7 | 416 | fh2SumPtInC10bkglead(0x0), |
417 | fh2SumPtInC20bkglead(0x0), | |
418 | fh2SumPtInC30bkglead(0x0), | |
419 | fh2SumPtInC60bkglead(0x0), | |
420 | fh2SumPtInC10bkgsublead(0x0), | |
421 | fh2SumPtInC20bkgsublead(0x0), | |
422 | fh2SumPtInC30bkgsublead(0x0), | |
423 | fh2SumPtInC60bkgsublead(0x0), | |
ea693273 | 424 | fh2SumPtOutC10bkg(0x0), |
425 | fh2SumPtOutC20bkg(0x0), | |
426 | fh2SumPtOutC30bkg(0x0), | |
427 | fh2SumPtOutC60bkg(0x0), | |
a9e585a7 | 428 | fh2SumPtOutC10bkglead(0x0), |
429 | fh2SumPtOutC20bkglead(0x0), | |
430 | fh2SumPtOutC30bkglead(0x0), | |
431 | fh2SumPtOutC60bkglead(0x0), | |
432 | fh2SumPtOutC10bkgsublead(0x0), | |
433 | fh2SumPtOutC20bkgsublead(0x0), | |
434 | fh2SumPtOutC30bkgsublead(0x0), | |
435 | fh2SumPtOutC60bkgsublead(0x0), | |
ea693273 | 436 | fh2DeltaRC10pt1(0x0), |
437 | fh2DeltaRC20pt1(0x0), | |
438 | fh2DeltaRC30pt1(0x0), | |
439 | fh2DeltaRC60pt1(0x0), | |
440 | fh2DeltaRC10pt2(0x0), | |
441 | fh2DeltaRC20pt2(0x0), | |
442 | fh2DeltaRC30pt2(0x0), | |
443 | fh2DeltaRC60pt2(0x0), | |
444 | fh2DeltaRC10pt3(0x0), | |
445 | fh2DeltaRC20pt3(0x0), | |
446 | fh2DeltaRC30pt3(0x0), | |
447 | fh2DeltaRC60pt3(0x0), | |
448 | fh2DeltaRC10pt4(0x0), | |
449 | fh2DeltaRC20pt4(0x0), | |
450 | fh2DeltaRC30pt4(0x0), | |
451 | fh2DeltaRC60pt4(0x0), | |
452 | fh2DeltaEtaC10pt1(0x0), | |
453 | fh2DeltaEtaC20pt1(0x0), | |
454 | fh2DeltaEtaC30pt1(0x0), | |
455 | fh2DeltaEtaC60pt1(0x0), | |
456 | fh2DeltaEtaC10pt2(0x0), | |
457 | fh2DeltaEtaC20pt2(0x0), | |
458 | fh2DeltaEtaC30pt2(0x0), | |
459 | fh2DeltaEtaC60pt2(0x0), | |
460 | fh2DeltaEtaC10pt3(0x0), | |
461 | fh2DeltaEtaC20pt3(0x0), | |
462 | fh2DeltaEtaC30pt3(0x0), | |
463 | fh2DeltaEtaC60pt3(0x0), | |
464 | fh2DeltaEtaC10pt4(0x0), | |
465 | fh2DeltaEtaC20pt4(0x0), | |
466 | fh2DeltaEtaC30pt4(0x0), | |
467 | fh2DeltaEtaC60pt4(0x0), | |
468 | fh2DeltaPhiC10pt1(0x0), | |
469 | fh2DeltaPhiC20pt1(0x0), | |
470 | fh2DeltaPhiC30pt1(0x0), | |
471 | fh2DeltaPhiC60pt1(0x0), | |
472 | fh2DeltaPhiC10pt2(0x0), | |
473 | fh2DeltaPhiC20pt2(0x0), | |
474 | fh2DeltaPhiC30pt2(0x0), | |
475 | fh2DeltaPhiC60pt2(0x0), | |
476 | fh2DeltaPhiC10pt3(0x0), | |
477 | fh2DeltaPhiC20pt3(0x0), | |
478 | fh2DeltaPhiC30pt3(0x0), | |
479 | fh2DeltaPhiC60pt3(0x0), | |
480 | fh2DeltaPhiC10pt4(0x0), | |
481 | fh2DeltaPhiC20pt4(0x0), | |
482 | fh2DeltaPhiC30pt4(0x0), | |
483 | fh2DeltaPhiC60pt4(0x0), | |
a9e585a7 | 484 | fh2DeltaRC10pt1lead(0x0), |
485 | fh2DeltaRC20pt1lead(0x0), | |
486 | fh2DeltaRC30pt1lead(0x0), | |
487 | fh2DeltaRC60pt1lead(0x0), | |
488 | fh2DeltaRC10pt2lead(0x0), | |
489 | fh2DeltaRC20pt2lead(0x0), | |
490 | fh2DeltaRC30pt2lead(0x0), | |
491 | fh2DeltaRC60pt2lead(0x0), | |
492 | fh2DeltaRC10pt3lead(0x0), | |
493 | fh2DeltaRC20pt3lead(0x0), | |
494 | fh2DeltaRC30pt3lead(0x0), | |
495 | fh2DeltaRC60pt3lead(0x0), | |
496 | fh2DeltaRC10pt4lead(0x0), | |
497 | fh2DeltaRC20pt4lead(0x0), | |
498 | fh2DeltaRC30pt4lead(0x0), | |
499 | fh2DeltaRC60pt4lead(0x0), | |
500 | fh2DeltaEtaC10pt1lead(0x0), | |
501 | fh2DeltaEtaC20pt1lead(0x0), | |
502 | fh2DeltaEtaC30pt1lead(0x0), | |
503 | fh2DeltaEtaC60pt1lead(0x0), | |
504 | fh2DeltaEtaC10pt2lead(0x0), | |
505 | fh2DeltaEtaC20pt2lead(0x0), | |
506 | fh2DeltaEtaC30pt2lead(0x0), | |
507 | fh2DeltaEtaC60pt2lead(0x0), | |
508 | fh2DeltaEtaC10pt3lead(0x0), | |
509 | fh2DeltaEtaC20pt3lead(0x0), | |
510 | fh2DeltaEtaC30pt3lead(0x0), | |
511 | fh2DeltaEtaC60pt3lead(0x0), | |
512 | fh2DeltaEtaC10pt4lead(0x0), | |
513 | fh2DeltaEtaC20pt4lead(0x0), | |
514 | fh2DeltaEtaC30pt4lead(0x0), | |
515 | fh2DeltaEtaC60pt4lead(0x0), | |
516 | fh2DeltaPhiC10pt1lead(0x0), | |
517 | fh2DeltaPhiC20pt1lead(0x0), | |
518 | fh2DeltaPhiC30pt1lead(0x0), | |
519 | fh2DeltaPhiC60pt1lead(0x0), | |
520 | fh2DeltaPhiC10pt2lead(0x0), | |
521 | fh2DeltaPhiC20pt2lead(0x0), | |
522 | fh2DeltaPhiC30pt2lead(0x0), | |
523 | fh2DeltaPhiC60pt2lead(0x0), | |
524 | fh2DeltaPhiC10pt3lead(0x0), | |
525 | fh2DeltaPhiC20pt3lead(0x0), | |
526 | fh2DeltaPhiC30pt3lead(0x0), | |
527 | fh2DeltaPhiC60pt3lead(0x0), | |
528 | fh2DeltaPhiC10pt4lead(0x0), | |
529 | fh2DeltaPhiC20pt4lead(0x0), | |
530 | fh2DeltaPhiC30pt4lead(0x0), | |
531 | fh2DeltaPhiC60pt4lead(0x0), | |
532 | fh2DeltaRC10pt1sublead(0x0), | |
533 | fh2DeltaRC20pt1sublead(0x0), | |
534 | fh2DeltaRC30pt1sublead(0x0), | |
535 | fh2DeltaRC60pt1sublead(0x0), | |
536 | fh2DeltaRC10pt2sublead(0x0), | |
537 | fh2DeltaRC20pt2sublead(0x0), | |
538 | fh2DeltaRC30pt2sublead(0x0), | |
539 | fh2DeltaRC60pt2sublead(0x0), | |
540 | fh2DeltaRC10pt3sublead(0x0), | |
541 | fh2DeltaRC20pt3sublead(0x0), | |
542 | fh2DeltaRC30pt3sublead(0x0), | |
543 | fh2DeltaRC60pt3sublead(0x0), | |
544 | fh2DeltaRC10pt4sublead(0x0), | |
545 | fh2DeltaRC20pt4sublead(0x0), | |
546 | fh2DeltaRC30pt4sublead(0x0), | |
547 | fh2DeltaRC60pt4sublead(0x0), | |
548 | fh2DeltaEtaC10pt1sublead(0x0), | |
549 | fh2DeltaEtaC20pt1sublead(0x0), | |
550 | fh2DeltaEtaC30pt1sublead(0x0), | |
551 | fh2DeltaEtaC60pt1sublead(0x0), | |
552 | fh2DeltaEtaC10pt2sublead(0x0), | |
553 | fh2DeltaEtaC20pt2sublead(0x0), | |
554 | fh2DeltaEtaC30pt2sublead(0x0), | |
555 | fh2DeltaEtaC60pt2sublead(0x0), | |
556 | fh2DeltaEtaC10pt3sublead(0x0), | |
557 | fh2DeltaEtaC20pt3sublead(0x0), | |
558 | fh2DeltaEtaC30pt3sublead(0x0), | |
559 | fh2DeltaEtaC60pt3sublead(0x0), | |
560 | fh2DeltaEtaC10pt4sublead(0x0), | |
561 | fh2DeltaEtaC20pt4sublead(0x0), | |
562 | fh2DeltaEtaC30pt4sublead(0x0), | |
563 | fh2DeltaEtaC60pt4sublead(0x0), | |
564 | fh2DeltaPhiC10pt1sublead(0x0), | |
565 | fh2DeltaPhiC20pt1sublead(0x0), | |
566 | fh2DeltaPhiC30pt1sublead(0x0), | |
567 | fh2DeltaPhiC60pt1sublead(0x0), | |
568 | fh2DeltaPhiC10pt2sublead(0x0), | |
569 | fh2DeltaPhiC20pt2sublead(0x0), | |
570 | fh2DeltaPhiC30pt2sublead(0x0), | |
571 | fh2DeltaPhiC60pt2sublead(0x0), | |
572 | fh2DeltaPhiC10pt3sublead(0x0), | |
573 | fh2DeltaPhiC20pt3sublead(0x0), | |
574 | fh2DeltaPhiC30pt3sublead(0x0), | |
575 | fh2DeltaPhiC60pt3sublead(0x0), | |
576 | fh2DeltaPhiC10pt4sublead(0x0), | |
577 | fh2DeltaPhiC20pt4sublead(0x0), | |
578 | fh2DeltaPhiC30pt4sublead(0x0), | |
579 | fh2DeltaPhiC60pt4sublead(0x0), | |
ea693273 | 580 | fh2AngStructpt1C10(0x0), |
581 | fh2AngStructpt2C10(0x0), | |
582 | fh2AngStructpt3C10(0x0), | |
583 | fh2AngStructpt4C10(0x0), | |
584 | fh2AngStructpt1C20(0x0), | |
585 | fh2AngStructpt2C20(0x0), | |
586 | fh2AngStructpt3C20(0x0), | |
587 | fh2AngStructpt4C20(0x0), | |
588 | fh2AngStructpt1C30(0x0), | |
589 | fh2AngStructpt2C30(0x0), | |
590 | fh2AngStructpt3C30(0x0), | |
591 | fh2AngStructpt4C30(0x0), | |
592 | fh2AngStructpt1C60(0x0), | |
593 | fh2AngStructpt2C60(0x0), | |
594 | fh2AngStructpt3C60(0x0), | |
595 | fh2AngStructpt4C60(0x0) | |
75bf77e3 | 596 | |
597 | { | |
598 | // Constructor | |
599 | ||
600 | fJetBranchName[0] = ""; | |
601 | fJetBranchName[1] = ""; | |
602 | ||
603 | fListJets[0] = new TList; | |
604 | fListJets[1] = new TList; | |
605 | ||
606 | DefineOutput(1, TList::Class()); | |
607 | } | |
608 | ||
609 | AliAnalysisTaskJetCore::~AliAnalysisTaskJetCore() | |
610 | { | |
611 | delete fListJets[0]; | |
612 | delete fListJets[1]; | |
613 | } | |
614 | ||
615 | void AliAnalysisTaskJetCore::SetBranchNames(const TString &branch1, const TString &branch2) | |
616 | { | |
617 | fJetBranchName[0] = branch1; | |
618 | fJetBranchName[1] = branch2; | |
619 | } | |
620 | ||
621 | void AliAnalysisTaskJetCore::Init() | |
622 | { | |
623 | ||
624 | // check for jet branches | |
625 | if(!strlen(fJetBranchName[0].Data()) || !strlen(fJetBranchName[1].Data())){ | |
626 | AliError("Jet branch name not set."); | |
627 | } | |
628 | ||
629 | } | |
630 | ||
631 | void AliAnalysisTaskJetCore::UserCreateOutputObjects() | |
632 | { | |
633 | // Create histograms | |
634 | // Called once | |
635 | OpenFile(1); | |
636 | if(!fOutputList) fOutputList = new TList; | |
637 | fOutputList->SetOwner(kTRUE); | |
638 | ||
639 | Bool_t oldStatus = TH1::AddDirectoryStatus(); | |
640 | TH1::AddDirectory(kFALSE); | |
641 | ||
642 | ||
643 | fHistEvtSelection = new TH1I("fHistEvtSelection", "event selection", 6, -0.5, 5.5); | |
644 | fHistEvtSelection->GetXaxis()->SetBinLabel(1,"ACCEPTED"); | |
645 | fHistEvtSelection->GetXaxis()->SetBinLabel(2,"events IN"); | |
646 | fHistEvtSelection->GetXaxis()->SetBinLabel(3,"event selection (rejected)"); | |
647 | fHistEvtSelection->GetXaxis()->SetBinLabel(4,"vertex cut (rejected)"); | |
648 | fHistEvtSelection->GetXaxis()->SetBinLabel(5,"centrality (rejected)"); | |
649 | fHistEvtSelection->GetXaxis()->SetBinLabel(6,"multiplicity (rejected)"); | |
650 | ||
651 | fHistJetSelection = new TH1I("fHistJetSelection", "jet selection", 8, -0.5, 7.5); | |
652 | fHistJetSelection->GetXaxis()->SetBinLabel(1,"ACCEPTED"); | |
653 | fHistJetSelection->GetXaxis()->SetBinLabel(2,"probes IN"); | |
654 | fHistJetSelection->GetXaxis()->SetBinLabel(3,"no matching jet"); | |
655 | fHistJetSelection->GetXaxis()->SetBinLabel(4,"not in list"); | |
656 | fHistJetSelection->GetXaxis()->SetBinLabel(5,"fraction cut"); | |
657 | fHistJetSelection->GetXaxis()->SetBinLabel(6,"acceptance cut"); | |
658 | fHistJetSelection->GetXaxis()->SetBinLabel(7,"p_{T} cut"); | |
659 | fHistJetSelection->GetXaxis()->SetBinLabel(8,"trigger exclude mask"); | |
660 | ||
661 | fh2JetSelection = new TH2F("fh2JetSelection", "jet selection", 8, -0.5, 7.5,100,0.,200.); | |
662 | fh2JetSelection->GetXaxis()->SetBinLabel(1,"ACCEPTED"); | |
663 | fh2JetSelection->GetXaxis()->SetBinLabel(2,"probes IN"); | |
664 | fh2JetSelection->GetXaxis()->SetBinLabel(3,"no matching jet"); | |
665 | fh2JetSelection->GetXaxis()->SetBinLabel(4,"not in list"); | |
666 | fh2JetSelection->GetXaxis()->SetBinLabel(5,"fraction cut"); | |
667 | fh2JetSelection->GetXaxis()->SetBinLabel(6,"acceptance cut"); | |
668 | fh2JetSelection->GetXaxis()->SetBinLabel(7,"p_{T} cut"); | |
669 | fh2JetSelection->GetXaxis()->SetBinLabel(8,"trigger exclude mask"); | |
670 | ||
671 | ||
ea693273 | 672 | //UInt_t entries = 0; // bit coded, see GetDimParams() below |
673 | //UInt_t opt = 0; // bit coded, default (0) or high resolution (1) | |
674 | ||
675 | // Int_t bins[5]={10,10,150,150,15}; | |
676 | //Double_t xmin[5]={0.,0.,0.,0.,0.}; | |
677 | //Double_t xmax[5]={100.,6.5,150.,1.5,1.5}; | |
678 | //fhnDeltaRjets = new THnSparseF("fhnDeltaRjets", "fhnDeltaRjets",5,bins,xmin,xmax); | |
679 | ||
680 | ||
75bf77e3 | 681 | |
682 | fh2JetCoreMethod1C10 = new TH2F("JetCoreMethod1C10","",150, 0., 150.,100, 0., 1.5); | |
683 | fh2JetCoreMethod2C10 = new TH2F("JetCoreMethod2C10","",150, 0., 150.,100, 0., 1.5); | |
684 | fh2JetCoreMethod3C10 = new TH2F("JetCoreMethod3C10","",150, 0., 150.,100, 0., 1.5); | |
ea693273 | 685 | fh2JetCoreMethod1C20 = new TH2F("JetCoreMethod1C20","",150, 0., 150.,100, 0., 1.5); |
686 | fh2JetCoreMethod2C20 = new TH2F("JetCoreMethod2C20","",150, 0., 150.,100, 0., 1.5); | |
687 | fh2JetCoreMethod3C20 = new TH2F("JetCoreMethod3C20","",150, 0., 150.,100, 0., 1.5); | |
75bf77e3 | 688 | fh2JetCoreMethod1C30 = new TH2F("JetCoreMethod1C30","",150, 0., 150.,100, 0., 1.5); |
689 | fh2JetCoreMethod2C30 = new TH2F("JetCoreMethod2C30","",150, 0., 150.,100, 0., 1.5); | |
690 | fh2JetCoreMethod3C30 = new TH2F("JetCoreMethod3C30","",150, 0., 150.,100, 0., 1.5); | |
691 | fh2JetCoreMethod1C60 = new TH2F("JetCoreMethod1C60","",150, 0., 150.,100, 0., 1.5); | |
692 | fh2JetCoreMethod2C60 = new TH2F("JetCoreMethod2C60","",150, 0., 150.,100, 0., 1.5); | |
693 | fh2JetCoreMethod3C60 = new TH2F("JetCoreMethod3C60","",150, 0., 150.,100, 0., 1.5); | |
a9e585a7 | 694 | fh2JetCoreMethod3C10lead = new TH2F("JetCoreMethod3C10lead","",150, 0., 150.,100, 0., 1.5); |
695 | fh2JetCoreMethod3C20lead = new TH2F("JetCoreMethod3C20lead","",150, 0., 150.,100, 0., 1.5); | |
696 | fh2JetCoreMethod3C30lead = new TH2F("JetCoreMethod3C30lead","",150, 0., 150.,100, 0., 1.5); | |
697 | fh2JetCoreMethod3C60lead = new TH2F("JetCoreMethod3C60lead","",150, 0., 150.,100, 0., 1.5); | |
698 | fh2JetCoreMethod3C10sublead = new TH2F("JetCoreMethod3C10sublead","",150, 0., 150.,100, 0., 1.5); | |
699 | fh2JetCoreMethod3C20sublead = new TH2F("JetCoreMethod3C20sublead","",150, 0., 150.,100, 0., 1.5); | |
700 | fh2JetCoreMethod3C30sublead = new TH2F("JetCoreMethod3C30sublead","",150, 0., 150.,100, 0., 1.5); | |
701 | fh2JetCoreMethod3C60sublead = new TH2F("JetCoreMethod3C60sublead","",150, 0., 150.,100, 0., 1.5); | |
ea693273 | 702 | |
a9e585a7 | 703 | fh2SumPtInC10 = new TH2F("PtSumInC10","",150, 0., 150.,5000, 0., 25.); |
704 | fh2SumPtInC20 = new TH2F("PtSumInC20","",150, 0., 150.,5000, 0., 25.); | |
705 | fh2SumPtInC30 = new TH2F("PtSumInC30","",150, 0., 150.,5000, 0., 25.); | |
706 | fh2SumPtInC60 = new TH2F("PtSumInC60","",150, 0., 150.,5000, 0., 25.); | |
707 | fh2SumPtInC10lead = new TH2F("PtSumInC10lead","",150, 0., 150.,5000, 0., 25.); | |
708 | fh2SumPtInC20lead = new TH2F("PtSumInC20lead","",150, 0., 150.,5000, 0., 25.); | |
709 | fh2SumPtInC30lead = new TH2F("PtSumInC30lead","",150, 0., 150.,5000, 0., 25.); | |
710 | fh2SumPtInC60lead = new TH2F("PtSumInC60lead","",150, 0., 150.,5000, 0., 25.); | |
711 | fh2SumPtInC10sublead = new TH2F("PtSumInC10sublead","",150, 0., 150.,5000, 0., 25.); | |
712 | fh2SumPtInC20sublead = new TH2F("PtSumInC20sublead","",150, 0., 150.,5000, 0., 25.); | |
713 | fh2SumPtInC30sublead = new TH2F("PtSumInC30sublead","",150, 0., 150.,5000, 0., 25.); | |
714 | fh2SumPtInC60sublead = new TH2F("PtSumInC60sublead","",150, 0., 150.,5000, 0., 25.); | |
715 | ||
716 | fh2SumPtOutC10 = new TH2F("PtSumOutC10","",150, 0., 150.,5000, 0., 25.); | |
717 | fh2SumPtOutC20 = new TH2F("PtSumOutC20","",150, 0., 150.,5000, 0., 25.); | |
718 | fh2SumPtOutC30 = new TH2F("PtSumOutC30","",150, 0., 150.,5000, 0., 25.); | |
719 | fh2SumPtOutC60 = new TH2F("PtSumOutC60","",150, 0., 150.,5000, 0., 25.); | |
720 | fh2SumPtOutC10lead = new TH2F("PtSumOutC10lead","",150, 0., 150.,5000, 0., 25.); | |
721 | fh2SumPtOutC20lead = new TH2F("PtSumOutC20lead","",150, 0., 150.,5000, 0., 25.); | |
722 | fh2SumPtOutC30lead = new TH2F("PtSumOutC30lead","",150, 0., 150.,5000, 0., 25.); | |
723 | fh2SumPtOutC60lead = new TH2F("PtSumOutC60lead","",150, 0., 150.,5000, 0., 25.); | |
724 | fh2SumPtOutC10sublead = new TH2F("PtSumOutC10sublead","",150, 0., 150.,5000, 0., 25.); | |
725 | fh2SumPtOutC20sublead = new TH2F("PtSumOutC20sublead","",150, 0., 150.,5000, 0., 25.); | |
726 | fh2SumPtOutC30sublead = new TH2F("PtSumOutC30sublead","",150, 0., 150.,5000, 0., 25.); | |
727 | fh2SumPtOutC60sublead = new TH2F("PtSumOutC60sublead","",150, 0., 150.,5000, 0., 25.); | |
728 | ||
729 | ||
730 | fh2SumPtInC10bkg = new TH2F("PtSumInC10bkg","",150, 0., 150.,5000, 0., 25.); | |
731 | fh2SumPtInC20bkg = new TH2F("PtSumInC20bkg","",150, 0., 150.,5000, 0., 25.); | |
732 | fh2SumPtInC30bkg = new TH2F("PtSumInC30bkg","",150, 0., 150.,5000, 0., 25.); | |
733 | fh2SumPtInC60bkg = new TH2F("PtSumInC60bkg","",150, 0., 150.,5000, 0., 25.); | |
734 | fh2SumPtInC10bkglead = new TH2F("PtSumInC10bkglead","",150, 0., 150.,5000, 0., 25.); | |
735 | fh2SumPtInC20bkglead = new TH2F("PtSumInC20bkglead","",150, 0., 150.,5000, 0., 25.); | |
736 | fh2SumPtInC30bkglead = new TH2F("PtSumInC30bkglead","",150, 0., 150.,5000, 0., 25.); | |
737 | fh2SumPtInC60bkglead = new TH2F("PtSumInC60bkglead","",150, 0., 150.,5000, 0., 25.); | |
738 | fh2SumPtInC10bkgsublead = new TH2F("PtSumInC10bkgsublead","",150, 0., 150.,5000, 0., 25.); | |
739 | fh2SumPtInC20bkgsublead = new TH2F("PtSumInC20bkgsublead","",150, 0., 150.,5000, 0., 25.); | |
740 | fh2SumPtInC30bkgsublead = new TH2F("PtSumInC30bkgsublead","",150, 0., 150.,5000, 0., 25.); | |
741 | fh2SumPtInC60bkgsublead = new TH2F("PtSumInC60bkgsublead","",150, 0., 150.,5000, 0., 25.); | |
742 | ||
743 | ||
744 | fh2SumPtOutC10bkg = new TH2F("PtSumOutC10bkg","",150, 0., 150.,5000, 0., 25.); | |
745 | fh2SumPtOutC20bkg = new TH2F("PtSumOutC20bkg","",150, 0., 150.,5000, 0., 25.); | |
746 | fh2SumPtOutC30bkg = new TH2F("PtSumOutC30bkg","",150, 0., 150.,5000, 0., 25.); | |
747 | fh2SumPtOutC60bkg = new TH2F("PtSumOutC60bkg","",150, 0., 150.,5000, 0., 25.); | |
748 | fh2SumPtOutC10bkglead = new TH2F("PtSumOutC10bkglead","",150, 0., 150.,5000, 0., 25.); | |
749 | fh2SumPtOutC20bkglead = new TH2F("PtSumOutC20bkglead","",150, 0., 150.,5000, 0., 25.); | |
750 | fh2SumPtOutC30bkglead = new TH2F("PtSumOutC30bkglead","",150, 0., 150.,5000, 0., 25.); | |
751 | fh2SumPtOutC60bkglead = new TH2F("PtSumOutC60bkglead","",150, 0., 150.,5000, 0., 25.); | |
752 | fh2SumPtOutC10bkgsublead = new TH2F("PtSumOutC10bkgsublead","",150, 0., 150.,5000, 0., 25.); | |
753 | fh2SumPtOutC20bkgsublead = new TH2F("PtSumOutC20bkgsublead","",150, 0., 150.,5000, 0., 25.); | |
754 | fh2SumPtOutC30bkgsublead = new TH2F("PtSumOutC30bkgsublead","",150, 0., 150.,5000, 0., 25.); | |
755 | fh2SumPtOutC60bkgsublead = new TH2F("PtSumOutC60bkgsublead","",150, 0., 150.,5000, 0., 25.); | |
ea693273 | 756 | |
757 | ||
758 | ||
759 | fh2DeltaRC10pt1 = new TH2F("DeltaRC10pt1","",150, 0., 150.,100,0.,1.5); | |
760 | fh2DeltaRC20pt1 = new TH2F("DeltaRC20pt1","",150, 0., 150.,100,0.,1.5); | |
761 | fh2DeltaRC30pt1 = new TH2F("DeltaRC30pt1","",150, 0., 150.,100,0.,1.5); | |
762 | fh2DeltaRC60pt1 = new TH2F("DeltaRC60pt1","",150, 0., 150.,100,0.,1.5); | |
763 | fh2DeltaRC10pt2 = new TH2F("DeltaRC10pt2","",150, 0., 150.,100,0.,1.5); | |
764 | fh2DeltaRC20pt2 = new TH2F("DeltaRC20pt2","",150, 0., 150.,100,0.,1.5); | |
765 | fh2DeltaRC30pt2 = new TH2F("DeltaRC30pt2","",150, 0., 150.,100,0.,1.5); | |
766 | fh2DeltaRC60pt2 = new TH2F("DeltaRC60pt2","",150, 0., 150.,100,0.,1.5); | |
767 | fh2DeltaRC10pt3 = new TH2F("DeltaRC10pt3","",150, 0., 150.,100,0.,1.5); | |
768 | fh2DeltaRC20pt3 = new TH2F("DeltaRC20pt3","",150, 0., 150.,100,0.,1.5); | |
769 | fh2DeltaRC30pt3 = new TH2F("DeltaRC30pt3","",150, 0., 150.,100,0.,1.5); | |
770 | fh2DeltaRC60pt3 = new TH2F("DeltaRC60pt3","",150, 0., 150.,100,0.,1.5); | |
771 | fh2DeltaRC10pt4 = new TH2F("DeltaRC10pt4","",150, 0., 150.,100,0.,1.5); | |
772 | fh2DeltaRC20pt4 = new TH2F("DeltaRC20pt4","",150, 0., 150.,100,0.,1.5); | |
773 | fh2DeltaRC30pt4 = new TH2F("DeltaRC30pt4","",150, 0., 150.,100,0.,1.5); | |
774 | fh2DeltaRC60pt4 = new TH2F("DeltaRC60pt4","",150, 0., 150.,100,0.,1.5); | |
ea693273 | 775 | fh2DeltaEtaC10pt1 = new TH2F("DeltaEtaC10pt1","",150, 0., 150.,100,-1.5,1.5); |
776 | fh2DeltaEtaC20pt1 = new TH2F("DeltaEtaC20pt1","",150, 0., 150.,100,-1.5,1.5); | |
777 | fh2DeltaEtaC30pt1 = new TH2F("DeltaEtaC30pt1","",150, 0., 150.,100,-1.5,1.5); | |
778 | fh2DeltaEtaC60pt1 = new TH2F("DeltaEtaC60pt1","",150, 0., 150.,100,-1.5,1.5); | |
779 | fh2DeltaEtaC10pt2 = new TH2F("DeltaEtaC10pt2","",150, 0., 150.,100,-1.5,1.5); | |
780 | fh2DeltaEtaC20pt2 = new TH2F("DeltaEtaC20pt2","",150, 0., 150.,100,-1.5,1.5); | |
781 | fh2DeltaEtaC30pt2 = new TH2F("DeltaEtaC30pt2","",150, 0., 150.,100,-1.5,1.5); | |
782 | fh2DeltaEtaC60pt2 = new TH2F("DeltaEtaC60pt2","",150, 0., 150.,100,-1.5,1.5); | |
783 | fh2DeltaEtaC10pt3 = new TH2F("DeltaEtaC10pt3","",150, 0., 150.,100,-1.5,1.5); | |
784 | fh2DeltaEtaC20pt3 = new TH2F("DeltaEtaC20pt3","",150, 0., 150.,100,-1.5,1.5); | |
785 | fh2DeltaEtaC30pt3 = new TH2F("DeltaEtaC30pt3","",150, 0., 150.,100,-1.5,1.5); | |
786 | fh2DeltaEtaC60pt3 = new TH2F("DeltaEtaC60pt3","",150, 0., 150.,100,-1.5,1.5); | |
787 | fh2DeltaEtaC10pt4 = new TH2F("DeltaEtaC10pt4","",150, 0., 150.,100,-1.5,1.5); | |
788 | fh2DeltaEtaC20pt4 = new TH2F("DeltaEtaC20pt4","",150, 0., 150.,100,-1.5,1.5); | |
789 | fh2DeltaEtaC30pt4 = new TH2F("DeltaEtaC30pt4","",150, 0., 150.,100,-1.5,1.5); | |
790 | fh2DeltaEtaC60pt4 = new TH2F("DeltaEtaC60pt4","",150, 0., 150.,100,-1.5,1.5); | |
791 | fh2DeltaPhiC10pt1 = new TH2F("DeltaPhiC10pt1","",150, 0., 150.,100,-6.5,6.5); | |
792 | fh2DeltaPhiC20pt1 = new TH2F("DeltaPhiC20pt1","",150, 0., 150.,100,-6.5,6.5); | |
793 | fh2DeltaPhiC30pt1 = new TH2F("DeltaPhiC30pt1","",150, 0., 150.,100,-6.5,6.5); | |
794 | fh2DeltaPhiC60pt1 = new TH2F("DeltaPhiC60pt1","",150, 0., 150.,100,-6.5,6.5); | |
795 | fh2DeltaPhiC10pt2 = new TH2F("DeltaPhiC10pt2","",150, 0., 150.,100,-6.5,6.5); | |
796 | fh2DeltaPhiC20pt2 = new TH2F("DeltaPhiC20pt2","",150, 0., 150.,100,-6.5,6.5); | |
797 | fh2DeltaPhiC30pt2 = new TH2F("DeltaPhiC30pt2","",150, 0., 150.,100,-6.5,6.5); | |
798 | fh2DeltaPhiC60pt2 = new TH2F("DeltaPhiC60pt2","",150, 0., 150.,100,-6.5,6.5); | |
799 | fh2DeltaPhiC10pt3 = new TH2F("DeltaPhiC10pt3","",150, 0., 150.,100,-6.5,6.5); | |
800 | fh2DeltaPhiC20pt3 = new TH2F("DeltaPhiC20pt3","",150, 0., 150.,100,-6.5,6.5); | |
801 | fh2DeltaPhiC30pt3 = new TH2F("DeltaPhiC30pt3","",150, 0., 150.,100,-6.5,6.5); | |
802 | fh2DeltaPhiC60pt3 = new TH2F("DeltaPhiC60pt3","",150, 0., 150.,100,-6.5,6.5); | |
803 | fh2DeltaPhiC10pt4 = new TH2F("DeltaPhiC10pt4","",150, 0., 150.,100,-6.5,6.5); | |
804 | fh2DeltaPhiC20pt4 = new TH2F("DeltaPhiC20pt4","",150, 0., 150.,100,-6.5,6.5); | |
805 | fh2DeltaPhiC30pt4 = new TH2F("DeltaPhiC30pt4","",150, 0., 150.,100,-6.5,6.5); | |
806 | fh2DeltaPhiC60pt4 = new TH2F("DeltaPhiC60pt4","",150, 0., 150.,100,-6.5,6.5); | |
807 | ||
a9e585a7 | 808 | fh2DeltaRC10pt1lead = new TH2F("DeltaRC10pt1lead","",150, 0., 150.,100,0.,1.5); |
809 | fh2DeltaRC20pt1lead = new TH2F("DeltaRC20pt1lead","",150, 0., 150.,100,0.,1.5); | |
810 | fh2DeltaRC30pt1lead = new TH2F("DeltaRC30pt1lead","",150, 0., 150.,100,0.,1.5); | |
811 | fh2DeltaRC60pt1lead = new TH2F("DeltaRC60pt1lead","",150, 0., 150.,100,0.,1.5); | |
812 | fh2DeltaRC10pt2lead = new TH2F("DeltaRC10pt2lead","",150, 0., 150.,100,0.,1.5); | |
813 | fh2DeltaRC20pt2lead = new TH2F("DeltaRC20pt2lead","",150, 0., 150.,100,0.,1.5); | |
814 | fh2DeltaRC30pt2lead = new TH2F("DeltaRC30pt2lead","",150, 0., 150.,100,0.,1.5); | |
815 | fh2DeltaRC60pt2lead = new TH2F("DeltaRC60pt2lead","",150, 0., 150.,100,0.,1.5); | |
816 | fh2DeltaRC10pt3lead = new TH2F("DeltaRC10pt3lead","",150, 0., 150.,100,0.,1.5); | |
817 | fh2DeltaRC20pt3lead = new TH2F("DeltaRC20pt3lead","",150, 0., 150.,100,0.,1.5); | |
818 | fh2DeltaRC30pt3lead = new TH2F("DeltaRC30pt3lead","",150, 0., 150.,100,0.,1.5); | |
819 | fh2DeltaRC60pt3lead = new TH2F("DeltaRC60pt3lead","",150, 0., 150.,100,0.,1.5); | |
820 | fh2DeltaRC10pt4lead = new TH2F("DeltaRC10pt4lead","",150, 0., 150.,100,0.,1.5); | |
821 | fh2DeltaRC20pt4lead = new TH2F("DeltaRC20pt4lead","",150, 0., 150.,100,0.,1.5); | |
822 | fh2DeltaRC30pt4lead = new TH2F("DeltaRC30pt4lead","",150, 0., 150.,100,0.,1.5); | |
823 | fh2DeltaRC60pt4lead = new TH2F("DeltaRC60pt4lead","",150, 0., 150.,100,0.,1.5); | |
824 | fh2DeltaEtaC10pt1lead = new TH2F("DeltaEtaC10pt1lead","",150, 0., 150.,100,-1.5,1.5); | |
825 | fh2DeltaEtaC20pt1lead = new TH2F("DeltaEtaC20pt1lead","",150, 0., 150.,100,-1.5,1.5); | |
826 | fh2DeltaEtaC30pt1lead = new TH2F("DeltaEtaC30pt1lead","",150, 0., 150.,100,-1.5,1.5); | |
827 | fh2DeltaEtaC60pt1lead = new TH2F("DeltaEtaC60pt1lead","",150, 0., 150.,100,-1.5,1.5); | |
828 | fh2DeltaEtaC10pt2lead = new TH2F("DeltaEtaC10pt2lead","",150, 0., 150.,100,-1.5,1.5); | |
829 | fh2DeltaEtaC20pt2lead = new TH2F("DeltaEtaC20pt2lead","",150, 0., 150.,100,-1.5,1.5); | |
830 | fh2DeltaEtaC30pt2lead = new TH2F("DeltaEtaC30pt2lead","",150, 0., 150.,100,-1.5,1.5); | |
831 | fh2DeltaEtaC60pt2lead = new TH2F("DeltaEtaC60pt2lead","",150, 0., 150.,100,-1.5,1.5); | |
832 | fh2DeltaEtaC10pt3lead = new TH2F("DeltaEtaC10pt3lead","",150, 0., 150.,100,-1.5,1.5); | |
833 | fh2DeltaEtaC20pt3lead = new TH2F("DeltaEtaC20pt3lead","",150, 0., 150.,100,-1.5,1.5); | |
834 | fh2DeltaEtaC30pt3lead = new TH2F("DeltaEtaC30pt3lead","",150, 0., 150.,100,-1.5,1.5); | |
835 | fh2DeltaEtaC60pt3lead = new TH2F("DeltaEtaC60pt3lead","",150, 0., 150.,100,-1.5,1.5); | |
836 | fh2DeltaEtaC10pt4lead = new TH2F("DeltaEtaC10pt4lead","",150, 0., 150.,100,-1.5,1.5); | |
837 | fh2DeltaEtaC20pt4lead = new TH2F("DeltaEtaC20pt4lead","",150, 0., 150.,100,-1.5,1.5); | |
838 | fh2DeltaEtaC30pt4lead = new TH2F("DeltaEtaC30pt4lead","",150, 0., 150.,100,-1.5,1.5); | |
839 | fh2DeltaEtaC60pt4lead = new TH2F("DeltaEtaC60pt4lead","",150, 0., 150.,100,-1.5,1.5); | |
840 | fh2DeltaPhiC10pt1lead = new TH2F("DeltaPhiC10pt1lead","",150, 0., 150.,100,-6.5,6.5); | |
841 | fh2DeltaPhiC20pt1lead = new TH2F("DeltaPhiC20pt1lead","",150, 0., 150.,100,-6.5,6.5); | |
842 | fh2DeltaPhiC30pt1lead = new TH2F("DeltaPhiC30pt1lead","",150, 0., 150.,100,-6.5,6.5); | |
843 | fh2DeltaPhiC60pt1lead = new TH2F("DeltaPhiC60pt1lead","",150, 0., 150.,100,-6.5,6.5); | |
844 | fh2DeltaPhiC10pt2lead = new TH2F("DeltaPhiC10pt2lead","",150, 0., 150.,100,-6.5,6.5); | |
845 | fh2DeltaPhiC20pt2lead = new TH2F("DeltaPhiC20pt2lead","",150, 0., 150.,100,-6.5,6.5); | |
846 | fh2DeltaPhiC30pt2lead = new TH2F("DeltaPhiC30pt2lead","",150, 0., 150.,100,-6.5,6.5); | |
847 | fh2DeltaPhiC60pt2lead = new TH2F("DeltaPhiC60pt2lead","",150, 0., 150.,100,-6.5,6.5); | |
848 | fh2DeltaPhiC10pt3lead = new TH2F("DeltaPhiC10pt3lead","",150, 0., 150.,100,-6.5,6.5); | |
849 | fh2DeltaPhiC20pt3lead = new TH2F("DeltaPhiC20pt3lead","",150, 0., 150.,100,-6.5,6.5); | |
850 | fh2DeltaPhiC30pt3lead = new TH2F("DeltaPhiC30pt3lead","",150, 0., 150.,100,-6.5,6.5); | |
851 | fh2DeltaPhiC60pt3lead = new TH2F("DeltaPhiC60pt3lead","",150, 0., 150.,100,-6.5,6.5); | |
852 | fh2DeltaPhiC10pt4lead = new TH2F("DeltaPhiC10pt4lead","",150, 0., 150.,100,-6.5,6.5); | |
853 | fh2DeltaPhiC20pt4lead = new TH2F("DeltaPhiC20pt4lead","",150, 0., 150.,100,-6.5,6.5); | |
854 | fh2DeltaPhiC30pt4lead = new TH2F("DeltaPhiC30pt4lead","",150, 0., 150.,100,-6.5,6.5); | |
855 | fh2DeltaPhiC60pt4lead = new TH2F("DeltaPhiC60pt4lead","",150, 0., 150.,100,-6.5,6.5); | |
856 | ||
857 | ||
858 | fh2DeltaRC10pt1sublead = new TH2F("DeltaRC10pt1sublead","",150, 0., 150.,100,0.,1.5); | |
859 | fh2DeltaRC20pt1sublead = new TH2F("DeltaRC20pt1sublead","",150, 0., 150.,100,0.,1.5); | |
860 | fh2DeltaRC30pt1sublead = new TH2F("DeltaRC30pt1sublead","",150, 0., 150.,100,0.,1.5); | |
861 | fh2DeltaRC60pt1sublead = new TH2F("DeltaRC60pt1sublead","",150, 0., 150.,100,0.,1.5); | |
862 | fh2DeltaRC10pt2sublead = new TH2F("DeltaRC10pt2sublead","",150, 0., 150.,100,0.,1.5); | |
863 | fh2DeltaRC20pt2sublead = new TH2F("DeltaRC20pt2sublead","",150, 0., 150.,100,0.,1.5); | |
864 | fh2DeltaRC30pt2sublead = new TH2F("DeltaRC30pt2sublead","",150, 0., 150.,100,0.,1.5); | |
865 | fh2DeltaRC60pt2sublead = new TH2F("DeltaRC60pt2sublead","",150, 0., 150.,100,0.,1.5); | |
866 | fh2DeltaRC10pt3sublead = new TH2F("DeltaRC10pt3sublead","",150, 0., 150.,100,0.,1.5); | |
867 | fh2DeltaRC20pt3sublead = new TH2F("DeltaRC20pt3sublead","",150, 0., 150.,100,0.,1.5); | |
868 | fh2DeltaRC30pt3sublead = new TH2F("DeltaRC30pt3sublead","",150, 0., 150.,100,0.,1.5); | |
869 | fh2DeltaRC60pt3sublead = new TH2F("DeltaRC60pt3sublead","",150, 0., 150.,100,0.,1.5); | |
870 | fh2DeltaRC10pt4sublead = new TH2F("DeltaRC10pt4sublead","",150, 0., 150.,100,0.,1.5); | |
871 | fh2DeltaRC20pt4sublead = new TH2F("DeltaRC20pt4sublead","",150, 0., 150.,100,0.,1.5); | |
872 | fh2DeltaRC30pt4sublead = new TH2F("DeltaRC30pt4sublead","",150, 0., 150.,100,0.,1.5); | |
873 | fh2DeltaRC60pt4sublead = new TH2F("DeltaRC60pt4sublead","",150, 0., 150.,100,0.,1.5); | |
874 | fh2DeltaEtaC10pt1sublead = new TH2F("DeltaEtaC10pt1sublead","",150, 0., 150.,100,-1.5,1.5); | |
875 | fh2DeltaEtaC20pt1sublead = new TH2F("DeltaEtaC20pt1sublead","",150, 0., 150.,100,-1.5,1.5); | |
876 | fh2DeltaEtaC30pt1sublead = new TH2F("DeltaEtaC30pt1sublead","",150, 0., 150.,100,-1.5,1.5); | |
877 | fh2DeltaEtaC60pt1sublead = new TH2F("DeltaEtaC60pt1sublead","",150, 0., 150.,100,-1.5,1.5); | |
878 | fh2DeltaEtaC10pt2sublead = new TH2F("DeltaEtaC10pt2sublead","",150, 0., 150.,100,-1.5,1.5); | |
879 | fh2DeltaEtaC20pt2sublead = new TH2F("DeltaEtaC20pt2sublead","",150, 0., 150.,100,-1.5,1.5); | |
880 | fh2DeltaEtaC30pt2sublead = new TH2F("DeltaEtaC30pt2sublead","",150, 0., 150.,100,-1.5,1.5); | |
881 | fh2DeltaEtaC60pt2sublead = new TH2F("DeltaEtaC60pt2sublead","",150, 0., 150.,100,-1.5,1.5); | |
882 | fh2DeltaEtaC10pt3sublead = new TH2F("DeltaEtaC10pt3sublead","",150, 0., 150.,100,-1.5,1.5); | |
883 | fh2DeltaEtaC20pt3sublead = new TH2F("DeltaEtaC20pt3sublead","",150, 0., 150.,100,-1.5,1.5); | |
884 | fh2DeltaEtaC30pt3sublead = new TH2F("DeltaEtaC30pt3sublead","",150, 0., 150.,100,-1.5,1.5); | |
885 | fh2DeltaEtaC60pt3sublead = new TH2F("DeltaEtaC60pt3sublead","",150, 0., 150.,100,-1.5,1.5); | |
886 | fh2DeltaEtaC10pt4sublead = new TH2F("DeltaEtaC10pt4sublead","",150, 0., 150.,100,-1.5,1.5); | |
887 | fh2DeltaEtaC20pt4sublead = new TH2F("DeltaEtaC20pt4sublead","",150, 0., 150.,100,-1.5,1.5); | |
888 | fh2DeltaEtaC30pt4sublead = new TH2F("DeltaEtaC30pt4sublead","",150, 0., 150.,100,-1.5,1.5); | |
889 | fh2DeltaEtaC60pt4sublead = new TH2F("DeltaEtaC60pt4sublead","",150, 0., 150.,100,-1.5,1.5); | |
890 | fh2DeltaPhiC10pt1sublead = new TH2F("DeltaPhiC10pt1sublead","",150, 0., 150.,100,-6.5,6.5); | |
891 | fh2DeltaPhiC20pt1sublead = new TH2F("DeltaPhiC20pt1sublead","",150, 0., 150.,100,-6.5,6.5); | |
892 | fh2DeltaPhiC30pt1sublead = new TH2F("DeltaPhiC30pt1sublead","",150, 0., 150.,100,-6.5,6.5); | |
893 | fh2DeltaPhiC60pt1sublead = new TH2F("DeltaPhiC60pt1sublead","",150, 0., 150.,100,-6.5,6.5); | |
894 | fh2DeltaPhiC10pt2sublead = new TH2F("DeltaPhiC10pt2sublead","",150, 0., 150.,100,-6.5,6.5); | |
895 | fh2DeltaPhiC20pt2sublead = new TH2F("DeltaPhiC20pt2sublead","",150, 0., 150.,100,-6.5,6.5); | |
896 | fh2DeltaPhiC30pt2sublead = new TH2F("DeltaPhiC30pt2sublead","",150, 0., 150.,100,-6.5,6.5); | |
897 | fh2DeltaPhiC60pt2sublead = new TH2F("DeltaPhiC60pt2sublead","",150, 0., 150.,100,-6.5,6.5); | |
898 | fh2DeltaPhiC10pt3sublead = new TH2F("DeltaPhiC10pt3sublead","",150, 0., 150.,100,-6.5,6.5); | |
899 | fh2DeltaPhiC20pt3sublead = new TH2F("DeltaPhiC20pt3sublead","",150, 0., 150.,100,-6.5,6.5); | |
900 | fh2DeltaPhiC30pt3sublead = new TH2F("DeltaPhiC30pt3sublead","",150, 0., 150.,100,-6.5,6.5); | |
901 | fh2DeltaPhiC60pt3sublead = new TH2F("DeltaPhiC60pt3sublead","",150, 0., 150.,100,-6.5,6.5); | |
902 | fh2DeltaPhiC10pt4sublead = new TH2F("DeltaPhiC10pt4sublead","",150, 0., 150.,100,-6.5,6.5); | |
903 | fh2DeltaPhiC20pt4sublead = new TH2F("DeltaPhiC20pt4sublead","",150, 0., 150.,100,-6.5,6.5); | |
904 | fh2DeltaPhiC30pt4sublead = new TH2F("DeltaPhiC30pt4sublead","",150, 0., 150.,100,-6.5,6.5); | |
905 | fh2DeltaPhiC60pt4sublead = new TH2F("DeltaPhiC60pt4sublead","",150, 0., 150.,100,-6.5,6.5); | |
906 | ||
907 | ||
908 | ||
909 | ||
910 | ||
911 | if(fAngStructCloseTracks>0){ | |
912 | fh2AngStructpt1C10 = new TH2F("Ang struct pt1 C10","",15,0.,1.5,150,0.,10.); | |
913 | fh2AngStructpt2C10 = new TH2F("Ang struct pt2 C10","",15,0.,1.5,150,0.,10.); | |
914 | fh2AngStructpt3C10 = new TH2F("Ang struct pt3 C10","",15,0.,1.5,150,0.,10.); | |
915 | fh2AngStructpt4C10 = new TH2F("Ang struct pt4 C10","",15,0.,1.5,150,0.,10.); | |
916 | fh2AngStructpt1C20 = new TH2F("Ang struct pt1 C20","",15,0.,1.5,150,0.,10.); | |
917 | fh2AngStructpt2C20 = new TH2F("Ang struct pt2 C20","",15,0.,1.5,150,0.,10.); | |
918 | fh2AngStructpt3C20 = new TH2F("Ang struct pt3 C20","",15,0.,1.5,150,0.,10.); | |
919 | fh2AngStructpt4C20 = new TH2F("Ang struct pt4 C20","",15,0.,1.5,150,0.,10.); | |
920 | fh2AngStructpt1C30 = new TH2F("Ang struct pt1 C30","",15,0.,1.5,150,0.,10.); | |
921 | fh2AngStructpt2C30 = new TH2F("Ang struct pt2 C30","",15,0.,1.5,150,0.,10.); | |
922 | fh2AngStructpt3C30 = new TH2F("Ang struct pt3 C30","",15,0.,1.5,150,0.,10.); | |
923 | fh2AngStructpt4C30 = new TH2F("Ang struct pt4 C30","",15,0.,1.5,150,0.,10.); | |
924 | fh2AngStructpt1C60 = new TH2F("Ang struct pt1 C60","",15,0.,1.5,150,0.,10.); | |
925 | fh2AngStructpt2C60 = new TH2F("Ang struct pt2 C60","",15,0.,1.5,150,0.,10.); | |
926 | fh2AngStructpt3C60 = new TH2F("Ang struct pt3 C60","",15,0.,1.5,150,0.,10.); | |
927 | fh2AngStructpt4C60 = new TH2F("Ang struct pt4 C60","",15,0.,1.5,150,0.,10.); } | |
75bf77e3 | 928 | |
929 | fOutputList->Add(fHistEvtSelection); | |
930 | fOutputList->Add(fHistJetSelection); | |
931 | fOutputList->Add(fh2JetSelection); | |
a9e585a7 | 932 | |
75bf77e3 | 933 | fOutputList->Add(fh2JetCoreMethod1C10); |
934 | fOutputList->Add(fh2JetCoreMethod2C10); | |
935 | fOutputList->Add(fh2JetCoreMethod3C10); | |
ea693273 | 936 | fOutputList->Add(fh2JetCoreMethod1C20); |
937 | fOutputList->Add(fh2JetCoreMethod2C20); | |
938 | fOutputList->Add(fh2JetCoreMethod3C20); | |
75bf77e3 | 939 | fOutputList->Add(fh2JetCoreMethod1C30); |
940 | fOutputList->Add(fh2JetCoreMethod2C30); | |
941 | fOutputList->Add(fh2JetCoreMethod3C30); | |
942 | fOutputList->Add(fh2JetCoreMethod1C60); | |
943 | fOutputList->Add(fh2JetCoreMethod2C60); | |
944 | fOutputList->Add(fh2JetCoreMethod3C60); | |
a9e585a7 | 945 | |
946 | fOutputList->Add(fh2JetCoreMethod3C10lead); | |
947 | fOutputList->Add(fh2JetCoreMethod3C20lead); | |
948 | fOutputList->Add(fh2JetCoreMethod3C30lead); | |
949 | fOutputList->Add(fh2JetCoreMethod3C60lead); | |
950 | fOutputList->Add(fh2JetCoreMethod3C10sublead); | |
951 | fOutputList->Add(fh2JetCoreMethod3C20sublead); | |
952 | fOutputList->Add(fh2JetCoreMethod3C30sublead); | |
953 | fOutputList->Add(fh2JetCoreMethod3C60sublead); | |
954 | ||
ea693273 | 955 | fOutputList->Add(fh2SumPtInC10); |
956 | fOutputList->Add(fh2SumPtInC20); | |
957 | fOutputList->Add(fh2SumPtInC30); | |
958 | fOutputList->Add(fh2SumPtInC60); | |
a9e585a7 | 959 | fOutputList->Add(fh2SumPtInC10lead); |
960 | fOutputList->Add(fh2SumPtInC20lead); | |
961 | fOutputList->Add(fh2SumPtInC30lead); | |
962 | fOutputList->Add(fh2SumPtInC60lead); | |
963 | fOutputList->Add(fh2SumPtInC10sublead); | |
964 | fOutputList->Add(fh2SumPtInC20sublead); | |
965 | fOutputList->Add(fh2SumPtInC30sublead); | |
966 | fOutputList->Add(fh2SumPtInC60sublead); | |
967 | ||
ea693273 | 968 | fOutputList->Add(fh2SumPtOutC10); |
ea693273 | 969 | fOutputList->Add(fh2SumPtOutC20); |
970 | fOutputList->Add(fh2SumPtOutC30); | |
971 | fOutputList->Add(fh2SumPtOutC60); | |
a9e585a7 | 972 | fOutputList->Add(fh2SumPtOutC10lead); |
973 | fOutputList->Add(fh2SumPtOutC20lead); | |
974 | fOutputList->Add(fh2SumPtOutC30lead); | |
975 | fOutputList->Add(fh2SumPtOutC60lead); | |
976 | fOutputList->Add(fh2SumPtOutC10sublead); | |
977 | fOutputList->Add(fh2SumPtOutC20sublead); | |
978 | fOutputList->Add(fh2SumPtOutC30sublead); | |
979 | fOutputList->Add(fh2SumPtOutC60sublead); | |
ea693273 | 980 | |
981 | fOutputList->Add(fh2SumPtInC10bkg); | |
982 | fOutputList->Add(fh2SumPtInC20bkg); | |
983 | fOutputList->Add(fh2SumPtInC30bkg); | |
984 | fOutputList->Add(fh2SumPtInC60bkg); | |
a9e585a7 | 985 | fOutputList->Add(fh2SumPtInC10bkglead); |
986 | fOutputList->Add(fh2SumPtInC20bkglead); | |
987 | fOutputList->Add(fh2SumPtInC30bkglead); | |
988 | fOutputList->Add(fh2SumPtInC60bkglead); | |
989 | fOutputList->Add(fh2SumPtInC10bkgsublead); | |
990 | fOutputList->Add(fh2SumPtInC20bkgsublead); | |
991 | fOutputList->Add(fh2SumPtInC30bkgsublead); | |
992 | fOutputList->Add(fh2SumPtInC60bkgsublead); | |
993 | ||
ea693273 | 994 | fOutputList->Add(fh2SumPtOutC10bkg); |
995 | fOutputList->Add(fh2SumPtOutC20bkg); | |
996 | fOutputList->Add(fh2SumPtOutC30bkg); | |
997 | fOutputList->Add(fh2SumPtOutC60bkg); | |
a9e585a7 | 998 | fOutputList->Add(fh2SumPtOutC10bkglead); |
999 | fOutputList->Add(fh2SumPtOutC20bkglead); | |
1000 | fOutputList->Add(fh2SumPtOutC30bkglead); | |
1001 | fOutputList->Add(fh2SumPtOutC60bkglead); | |
1002 | fOutputList->Add(fh2SumPtOutC10bkgsublead); | |
1003 | fOutputList->Add(fh2SumPtOutC20bkgsublead); | |
1004 | fOutputList->Add(fh2SumPtOutC30bkgsublead); | |
1005 | fOutputList->Add(fh2SumPtOutC60bkgsublead); | |
1006 | ||
ea693273 | 1007 | fOutputList->Add(fh2DeltaRC10pt1); |
1008 | fOutputList->Add(fh2DeltaRC20pt1); | |
1009 | fOutputList->Add(fh2DeltaRC30pt1); | |
1010 | fOutputList->Add(fh2DeltaRC60pt1); | |
1011 | fOutputList->Add(fh2DeltaRC10pt2); | |
1012 | fOutputList->Add(fh2DeltaRC20pt2); | |
1013 | fOutputList->Add(fh2DeltaRC30pt2); | |
1014 | fOutputList->Add(fh2DeltaRC60pt2); | |
1015 | fOutputList->Add(fh2DeltaRC10pt3); | |
1016 | fOutputList->Add(fh2DeltaRC20pt3); | |
1017 | fOutputList->Add(fh2DeltaRC30pt3); | |
1018 | fOutputList->Add(fh2DeltaRC60pt3); | |
1019 | fOutputList->Add(fh2DeltaRC10pt4); | |
1020 | fOutputList->Add(fh2DeltaRC20pt4); | |
1021 | fOutputList->Add(fh2DeltaRC30pt4); | |
1022 | fOutputList->Add(fh2DeltaRC60pt4); | |
ea693273 | 1023 | fOutputList->Add(fh2DeltaEtaC10pt1); |
1024 | fOutputList->Add(fh2DeltaEtaC20pt1); | |
1025 | fOutputList->Add(fh2DeltaEtaC30pt1); | |
1026 | fOutputList->Add(fh2DeltaEtaC60pt1); | |
1027 | fOutputList->Add(fh2DeltaEtaC10pt2); | |
1028 | fOutputList->Add(fh2DeltaEtaC20pt2); | |
1029 | fOutputList->Add(fh2DeltaEtaC30pt2); | |
1030 | fOutputList->Add(fh2DeltaEtaC60pt2); | |
1031 | fOutputList->Add(fh2DeltaEtaC10pt3); | |
1032 | fOutputList->Add(fh2DeltaEtaC20pt3); | |
1033 | fOutputList->Add(fh2DeltaEtaC30pt3); | |
1034 | fOutputList->Add(fh2DeltaEtaC60pt3); | |
1035 | fOutputList->Add(fh2DeltaEtaC10pt4); | |
1036 | fOutputList->Add(fh2DeltaEtaC20pt4); | |
1037 | fOutputList->Add(fh2DeltaEtaC30pt4); | |
1038 | fOutputList->Add(fh2DeltaEtaC60pt4); | |
1039 | fOutputList->Add(fh2DeltaPhiC10pt1); | |
1040 | fOutputList->Add(fh2DeltaPhiC20pt1); | |
1041 | fOutputList->Add(fh2DeltaPhiC30pt1); | |
1042 | fOutputList->Add(fh2DeltaPhiC60pt1); | |
1043 | fOutputList->Add(fh2DeltaPhiC10pt2); | |
1044 | fOutputList->Add(fh2DeltaPhiC20pt2); | |
1045 | fOutputList->Add(fh2DeltaPhiC30pt2); | |
1046 | fOutputList->Add(fh2DeltaPhiC60pt2); | |
1047 | fOutputList->Add(fh2DeltaPhiC10pt3); | |
1048 | fOutputList->Add(fh2DeltaPhiC20pt3); | |
1049 | fOutputList->Add(fh2DeltaPhiC30pt3); | |
1050 | fOutputList->Add(fh2DeltaPhiC60pt3); | |
1051 | fOutputList->Add(fh2DeltaPhiC10pt4); | |
1052 | fOutputList->Add(fh2DeltaPhiC20pt4); | |
1053 | fOutputList->Add(fh2DeltaPhiC30pt4); | |
1054 | fOutputList->Add(fh2DeltaPhiC60pt4); | |
a9e585a7 | 1055 | |
1056 | fOutputList->Add(fh2DeltaRC10pt1lead); | |
1057 | fOutputList->Add(fh2DeltaRC20pt1lead); | |
1058 | fOutputList->Add(fh2DeltaRC30pt1lead); | |
1059 | fOutputList->Add(fh2DeltaRC60pt1lead); | |
1060 | fOutputList->Add(fh2DeltaRC10pt2lead); | |
1061 | fOutputList->Add(fh2DeltaRC20pt2lead); | |
1062 | fOutputList->Add(fh2DeltaRC30pt2lead); | |
1063 | fOutputList->Add(fh2DeltaRC60pt2lead); | |
1064 | fOutputList->Add(fh2DeltaRC10pt3lead); | |
1065 | fOutputList->Add(fh2DeltaRC20pt3lead); | |
1066 | fOutputList->Add(fh2DeltaRC30pt3lead); | |
1067 | fOutputList->Add(fh2DeltaRC60pt3lead); | |
1068 | fOutputList->Add(fh2DeltaRC10pt4lead); | |
1069 | fOutputList->Add(fh2DeltaRC20pt4lead); | |
1070 | fOutputList->Add(fh2DeltaRC30pt4lead); | |
1071 | fOutputList->Add(fh2DeltaRC60pt4lead); | |
1072 | fOutputList->Add(fh2DeltaEtaC10pt1lead); | |
1073 | fOutputList->Add(fh2DeltaEtaC20pt1lead); | |
1074 | fOutputList->Add(fh2DeltaEtaC30pt1lead); | |
1075 | fOutputList->Add(fh2DeltaEtaC60pt1lead); | |
1076 | fOutputList->Add(fh2DeltaEtaC10pt2lead); | |
1077 | fOutputList->Add(fh2DeltaEtaC20pt2lead); | |
1078 | fOutputList->Add(fh2DeltaEtaC30pt2lead); | |
1079 | fOutputList->Add(fh2DeltaEtaC60pt2lead); | |
1080 | fOutputList->Add(fh2DeltaEtaC10pt3lead); | |
1081 | fOutputList->Add(fh2DeltaEtaC20pt3lead); | |
1082 | fOutputList->Add(fh2DeltaEtaC30pt3lead); | |
1083 | fOutputList->Add(fh2DeltaEtaC60pt3lead); | |
1084 | fOutputList->Add(fh2DeltaEtaC10pt4lead); | |
1085 | fOutputList->Add(fh2DeltaEtaC20pt4lead); | |
1086 | fOutputList->Add(fh2DeltaEtaC30pt4lead); | |
1087 | fOutputList->Add(fh2DeltaEtaC60pt4lead); | |
1088 | fOutputList->Add(fh2DeltaPhiC10pt1lead); | |
1089 | fOutputList->Add(fh2DeltaPhiC20pt1lead); | |
1090 | fOutputList->Add(fh2DeltaPhiC30pt1lead); | |
1091 | fOutputList->Add(fh2DeltaPhiC60pt1lead); | |
1092 | fOutputList->Add(fh2DeltaPhiC10pt2lead); | |
1093 | fOutputList->Add(fh2DeltaPhiC20pt2lead); | |
1094 | fOutputList->Add(fh2DeltaPhiC30pt2lead); | |
1095 | fOutputList->Add(fh2DeltaPhiC60pt2lead); | |
1096 | fOutputList->Add(fh2DeltaPhiC10pt3lead); | |
1097 | fOutputList->Add(fh2DeltaPhiC20pt3lead); | |
1098 | fOutputList->Add(fh2DeltaPhiC30pt3lead); | |
1099 | fOutputList->Add(fh2DeltaPhiC60pt3lead); | |
1100 | fOutputList->Add(fh2DeltaPhiC10pt4lead); | |
1101 | fOutputList->Add(fh2DeltaPhiC20pt4lead); | |
1102 | fOutputList->Add(fh2DeltaPhiC30pt4lead); | |
1103 | fOutputList->Add(fh2DeltaPhiC60pt4lead); | |
1104 | ||
1105 | ||
1106 | ||
1107 | fOutputList->Add(fh2DeltaRC10pt1sublead); | |
1108 | fOutputList->Add(fh2DeltaRC20pt1sublead); | |
1109 | fOutputList->Add(fh2DeltaRC30pt1sublead); | |
1110 | fOutputList->Add(fh2DeltaRC60pt1sublead); | |
1111 | fOutputList->Add(fh2DeltaRC10pt2sublead); | |
1112 | fOutputList->Add(fh2DeltaRC20pt2sublead); | |
1113 | fOutputList->Add(fh2DeltaRC30pt2sublead); | |
1114 | fOutputList->Add(fh2DeltaRC60pt2sublead); | |
1115 | fOutputList->Add(fh2DeltaRC10pt3sublead); | |
1116 | fOutputList->Add(fh2DeltaRC20pt3sublead); | |
1117 | fOutputList->Add(fh2DeltaRC30pt3sublead); | |
1118 | fOutputList->Add(fh2DeltaRC60pt3sublead); | |
1119 | fOutputList->Add(fh2DeltaRC10pt4sublead); | |
1120 | fOutputList->Add(fh2DeltaRC20pt4sublead); | |
1121 | fOutputList->Add(fh2DeltaRC30pt4sublead); | |
1122 | fOutputList->Add(fh2DeltaRC60pt4sublead); | |
1123 | fOutputList->Add(fh2DeltaEtaC10pt1sublead); | |
1124 | fOutputList->Add(fh2DeltaEtaC20pt1sublead); | |
1125 | fOutputList->Add(fh2DeltaEtaC30pt1sublead); | |
1126 | fOutputList->Add(fh2DeltaEtaC60pt1sublead); | |
1127 | fOutputList->Add(fh2DeltaEtaC10pt2sublead); | |
1128 | fOutputList->Add(fh2DeltaEtaC20pt2sublead); | |
1129 | fOutputList->Add(fh2DeltaEtaC30pt2sublead); | |
1130 | fOutputList->Add(fh2DeltaEtaC60pt2sublead); | |
1131 | fOutputList->Add(fh2DeltaEtaC10pt3sublead); | |
1132 | fOutputList->Add(fh2DeltaEtaC20pt3sublead); | |
1133 | fOutputList->Add(fh2DeltaEtaC30pt3sublead); | |
1134 | fOutputList->Add(fh2DeltaEtaC60pt3sublead); | |
1135 | fOutputList->Add(fh2DeltaEtaC10pt4sublead); | |
1136 | fOutputList->Add(fh2DeltaEtaC20pt4sublead); | |
1137 | fOutputList->Add(fh2DeltaEtaC30pt4sublead); | |
1138 | fOutputList->Add(fh2DeltaEtaC60pt4sublead); | |
1139 | fOutputList->Add(fh2DeltaPhiC10pt1sublead); | |
1140 | fOutputList->Add(fh2DeltaPhiC20pt1sublead); | |
1141 | fOutputList->Add(fh2DeltaPhiC30pt1sublead); | |
1142 | fOutputList->Add(fh2DeltaPhiC60pt1sublead); | |
1143 | fOutputList->Add(fh2DeltaPhiC10pt2sublead); | |
1144 | fOutputList->Add(fh2DeltaPhiC20pt2sublead); | |
1145 | fOutputList->Add(fh2DeltaPhiC30pt2sublead); | |
1146 | fOutputList->Add(fh2DeltaPhiC60pt2sublead); | |
1147 | fOutputList->Add(fh2DeltaPhiC10pt3sublead); | |
1148 | fOutputList->Add(fh2DeltaPhiC20pt3sublead); | |
1149 | fOutputList->Add(fh2DeltaPhiC30pt3sublead); | |
1150 | fOutputList->Add(fh2DeltaPhiC60pt3sublead); | |
1151 | fOutputList->Add(fh2DeltaPhiC10pt4sublead); | |
1152 | fOutputList->Add(fh2DeltaPhiC20pt4sublead); | |
1153 | fOutputList->Add(fh2DeltaPhiC30pt4sublead); | |
1154 | fOutputList->Add(fh2DeltaPhiC60pt4sublead); | |
75bf77e3 | 1155 | |
a9e585a7 | 1156 | |
1157 | ||
1158 | ||
1159 | ||
1160 | if(fAngStructCloseTracks>0){ | |
ea693273 | 1161 | fOutputList->Add(fh2AngStructpt1C10); |
1162 | fOutputList->Add(fh2AngStructpt2C10); | |
1163 | fOutputList->Add(fh2AngStructpt3C10); | |
1164 | fOutputList->Add(fh2AngStructpt4C10); | |
1165 | fOutputList->Add(fh2AngStructpt1C20); | |
1166 | fOutputList->Add(fh2AngStructpt2C20); | |
1167 | fOutputList->Add(fh2AngStructpt3C20); | |
1168 | fOutputList->Add(fh2AngStructpt4C20); | |
1169 | fOutputList->Add(fh2AngStructpt1C30); | |
1170 | fOutputList->Add(fh2AngStructpt2C30); | |
1171 | fOutputList->Add(fh2AngStructpt3C30); | |
1172 | fOutputList->Add(fh2AngStructpt4C30); | |
1173 | fOutputList->Add(fh2AngStructpt1C60); | |
1174 | fOutputList->Add(fh2AngStructpt2C60); | |
1175 | fOutputList->Add(fh2AngStructpt3C60); | |
a9e585a7 | 1176 | fOutputList->Add(fh2AngStructpt4C60);} |
ea693273 | 1177 | |
1178 | ||
1179 | ||
75bf77e3 | 1180 | |
1181 | // =========== Switch on Sumw2 for all histos =========== | |
1182 | for (Int_t i=0; i<fOutputList->GetEntries(); ++i) { | |
1183 | TH1 *h1 = dynamic_cast<TH1*>(fOutputList->At(i)); | |
1184 | if (h1){ | |
1185 | h1->Sumw2(); | |
1186 | continue; | |
1187 | } | |
ea693273 | 1188 | THnSparse *hn = dynamic_cast<THnSparse*>(fOutputList->At(i)); |
1189 | if (hn){ | |
1190 | hn->Sumw2(); | |
1191 | } | |
75bf77e3 | 1192 | } |
1193 | TH1::AddDirectory(oldStatus); | |
1194 | ||
1195 | PostData(1, fOutputList); | |
1196 | } | |
1197 | ||
1198 | void AliAnalysisTaskJetCore::UserExec(Option_t *) | |
1199 | { | |
1200 | ||
1201 | ||
1202 | if(!strlen(fJetBranchName[0].Data()) || !strlen(fJetBranchName[1].Data())){ | |
1203 | AliError("Jet branch name not set."); | |
1204 | return; | |
1205 | } | |
1206 | ||
1207 | fESD=dynamic_cast<AliESDEvent*>(InputEvent()); | |
1208 | if (!fESD) { | |
1209 | AliError("ESD not available"); | |
1210 | fAOD = dynamic_cast<AliAODEvent*>(InputEvent()); | |
1211 | } else { | |
1212 | fAOD = dynamic_cast<AliAODEvent*>(AODEvent()); | |
1213 | } | |
ea693273 | 1214 | |
1215 | if(fNonStdFile.Length()!=0){ | |
1216 | // case that we have an AOD extension we need can fetch the jets from the extended output | |
1217 | AliAODHandler *aodH = dynamic_cast<AliAODHandler*>(AliAnalysisManager::GetAnalysisManager()->GetOutputEventHandler()); | |
1218 | fAODExtension = (aodH?aodH->GetExtension(fNonStdFile.Data()):0); | |
1219 | if(!fAODExtension){ | |
1220 | if(fDebug>1)Printf("AODExtension found for %s",fNonStdFile.Data()); | |
1221 | }} | |
1222 | ||
1223 | ||
1224 | ||
1225 | ||
75bf77e3 | 1226 | |
1227 | // -- event selection -- | |
1228 | fHistEvtSelection->Fill(1); // number of events before event selection | |
1229 | ||
1230 | // physics selection | |
1231 | AliInputEventHandler* inputHandler = (AliInputEventHandler*) | |
1232 | ((AliAnalysisManager::GetAnalysisManager())->GetInputEventHandler()); | |
1233 | if(!(inputHandler->IsEventSelected() & fOfflineTrgMask)){ | |
1234 | if(fDebug) Printf(" Trigger Selection: event REJECTED ... "); | |
1235 | fHistEvtSelection->Fill(2); | |
1236 | PostData(1, fOutputList); | |
1237 | return; | |
1238 | } | |
1239 | ||
1240 | // vertex selection | |
1241 | AliAODVertex* primVtx = fAOD->GetPrimaryVertex(); | |
1242 | Int_t nTracksPrim = primVtx->GetNContributors(); | |
1243 | if ((nTracksPrim < fMinContribVtx) || | |
1244 | (primVtx->GetZ() < fVtxZMin) || | |
1245 | (primVtx->GetZ() > fVtxZMax) ){ | |
1246 | if(fDebug) Printf("%s:%d primary vertex z = %f: event REJECTED...",(char*)__FILE__,__LINE__,primVtx->GetZ()); | |
1247 | fHistEvtSelection->Fill(3); | |
1248 | PostData(1, fOutputList); | |
1249 | return; | |
1250 | } | |
1251 | ||
1252 | // event class selection (from jet helper task) | |
1253 | Int_t eventClass = AliAnalysisHelperJetTasks::EventClass(); | |
1254 | if(fDebug) Printf("Event class %d", eventClass); | |
1255 | if (eventClass < fEvtClassMin || eventClass > fEvtClassMax){ | |
1256 | fHistEvtSelection->Fill(4); | |
1257 | PostData(1, fOutputList); | |
1258 | return; | |
1259 | } | |
1260 | ||
1261 | // centrality selection | |
1262 | AliCentrality *cent = 0x0; | |
1263 | Float_t centValue = 0.; | |
1264 | if(fESD) cent = fESD->GetCentrality(); | |
1265 | if(cent) centValue = cent->GetCentralityPercentile("V0M"); | |
1266 | if(fDebug) printf("centrality: %f\n", centValue); | |
1267 | if (centValue < fCentMin || centValue > fCentMax){ | |
1268 | fHistEvtSelection->Fill(4); | |
1269 | PostData(1, fOutputList); | |
1270 | return; | |
1271 | } | |
1272 | ||
1273 | ||
1274 | // multiplicity due to input tracks | |
ea693273 | 1275 | //Int_t nInputTracks = GetNInputTracks(); |
1276 | //if (nInputTracks < fNInputTracksMin || (fNInputTracksMax > -1 && nInputTracks > fNInputTracksMax)){ | |
1277 | // fHistEvtSelection->Fill(5); | |
1278 | // PostData(1, fOutputList); | |
1279 | // return; | |
1280 | // } | |
75bf77e3 | 1281 | |
1282 | ||
568f8fa2 | 1283 | fHistEvtSelection->Fill(0); |
1284 | // accepted events | |
75bf77e3 | 1285 | // -- end event selection -- |
ea693273 | 1286 | |
75bf77e3 | 1287 | // get background |
1288 | AliAODJetEventBackground* externalBackground = 0; | |
ea693273 | 1289 | if(fAOD&&!externalBackground&&fBackgroundBranch.Length()){ |
75bf77e3 | 1290 | externalBackground = (AliAODJetEventBackground*)(fAOD->FindListObject(fBackgroundBranch.Data())); |
1291 | if(!externalBackground)Printf("%s:%d Background branch not found %s",(char*)__FILE__,__LINE__,fBackgroundBranch.Data());; | |
1292 | } | |
ea693273 | 1293 | if(fAODExtension&&!externalBackground&&fBackgroundBranch.Length()){ |
1294 | externalBackground = (AliAODJetEventBackground*)(fAODExtension->GetAOD()->FindListObject(fBackgroundBranch.Data())); | |
1295 | if(!externalBackground)Printf("%s:%d Background branch not found %s",(char*)__FILE__,__LINE__,fBackgroundBranch.Data());; | |
1296 | } | |
1297 | ||
75bf77e3 | 1298 | Float_t rho = 0; |
1299 | if(externalBackground)rho = externalBackground->GetBackground(0); | |
1300 | ||
1301 | ||
1302 | // fetch jets | |
1303 | TClonesArray *aodJets[2]; | |
ea693273 | 1304 | aodJets[0]=0; |
1305 | if(fAOD&&!aodJets[0]){ | |
a9e585a7 | 1306 | aodJets[0] = dynamic_cast<TClonesArray*>(fAOD->FindListObject(fJetBranchName[0].Data())); |
ea693273 | 1307 | aodJets[1] = dynamic_cast<TClonesArray*>(fAOD->FindListObject(fJetBranchName[1].Data())); } |
1308 | if(fAODExtension && !aodJets[0]){ | |
a9e585a7 | 1309 | aodJets[0] = dynamic_cast<TClonesArray*>(fAODExtension->GetAOD()->FindListObject(fJetBranchName[0].Data())); |
ea693273 | 1310 | aodJets[1] = dynamic_cast<TClonesArray*>(fAODExtension->GetAOD()->FindListObject(fJetBranchName[1].Data())); } |
1311 | ||
a9e585a7 | 1312 | Double_t ptsub[aodJets[0]->GetEntriesFast()]; |
1313 | Int_t inord[aodJets[0]->GetEntriesFast()]; | |
1314 | for(Int_t n=0;n<aodJets[0]->GetEntriesFast();n++){ | |
1315 | ptsub[n]=0; | |
1316 | inord[n]=0;} | |
1317 | ||
ea693273 | 1318 | TList ParticleList; |
1319 | Int_t nT = GetListOfTracks(&ParticleList); | |
1320 | for (Int_t iJetType = 0; iJetType < 2; iJetType++) { | |
75bf77e3 | 1321 | fListJets[iJetType]->Clear(); |
1322 | if (!aodJets[iJetType]) continue; | |
1323 | ||
1324 | if(fDebug) Printf("%s: %d jets",fJetBranchName[iJetType].Data(),aodJets[iJetType]->GetEntriesFast()); | |
1325 | ||
a9e585a7 | 1326 | |
75bf77e3 | 1327 | for (Int_t iJet = 0; iJet < aodJets[iJetType]->GetEntriesFast(); iJet++) { |
1328 | AliAODJet *jet = dynamic_cast<AliAODJet*>((*aodJets[iJetType])[iJet]); | |
1329 | if (jet) fListJets[iJetType]->Add(jet); | |
a9e585a7 | 1330 | if(iJetType==0){ |
1331 | ptsub[iJet]=jet->Pt()-rho*jet->EffectiveAreaCharged(); | |
1332 | }}} | |
75bf77e3 | 1333 | |
1334 | Double_t etabig=0; | |
1335 | Double_t ptbig=0; | |
1336 | Double_t areabig=0; | |
1337 | Double_t phibig=0.; | |
1338 | Double_t etasmall=0; | |
1339 | Double_t ptsmall=0; | |
1340 | Double_t areasmall=0; | |
1341 | Double_t distr=0.; | |
1342 | Double_t phismall=0.; | |
a9e585a7 | 1343 | Int_t indexlead=-1; |
1344 | Int_t indexsublead=-1; | |
1345 | Int_t indexstop=-1; | |
1346 | ||
1347 | TMath::Sort(fListJets[0]->GetEntries(),ptsub,inord); | |
1348 | ||
1349 | for(Int_t jj=0;jj<fListJets[0]->GetEntries();jj++){ | |
1350 | AliAODJet* jetlead = (AliAODJet*)(fListJets[0]->At(inord[jj])); | |
1351 | if(jetlead->Pt()-rho*jetlead->EffectiveAreaCharged()<=0) continue; | |
1352 | if((jetlead->Eta()<fJetEtaMin)||(jetlead->Eta()>fJetEtaMax)) continue; | |
1353 | indexlead=inord[jj]; | |
1354 | indexstop=jj; | |
1355 | break;} | |
1356 | if((indexstop>-1)&&(indexstop+1<fListJets[0]->GetEntries()-1)){ | |
1357 | for(Int_t k=indexstop+1;k<fListJets[0]->GetEntries();k++){ | |
1358 | AliAODJet* jetsublead = (AliAODJet*)(fListJets[0]->At(inord[k])); | |
1359 | if(jetsublead->Pt()-rho*jetsublead->EffectiveAreaCharged()<=0) continue; | |
1360 | if((jetsublead->Eta()<fJetEtaMin)||(jetsublead->Eta()>fJetEtaMax)) continue; | |
1361 | indexsublead=inord[k]; | |
1362 | break;}} | |
1363 | ||
ea693273 | 1364 | |
1365 | Double_t up1[15]={0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.}; | |
1366 | Double_t up2[15]={0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.}; | |
1367 | Double_t up3[15]={0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.}; | |
1368 | Double_t up4[15]={0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.}; | |
1369 | Double_t down1[15]={0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.}; | |
1370 | Double_t down2[15]={0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.}; | |
1371 | Double_t down3[15]={0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.}; | |
1372 | Double_t down4[15]={0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.}; | |
1373 | ||
1374 | ||
1375 | ||
75bf77e3 | 1376 | for(Int_t i=0; i<fListJets[0]->GetEntries(); ++i){ |
1377 | AliAODJet* jetbig = (AliAODJet*)(fListJets[0]->At(i)); | |
1378 | etabig = jetbig->Eta(); | |
1379 | phibig = jetbig->Phi(); | |
1380 | ptbig = jetbig->Pt(); | |
1381 | if(ptbig==0) continue; | |
1382 | areabig = jetbig->EffectiveAreaCharged(); | |
ea693273 | 1383 | Double_t ptcorr=ptbig-rho*areabig; |
1384 | if(ptcorr<=0) continue; | |
75bf77e3 | 1385 | if((etabig<fJetEtaMin)||(etabig>fJetEtaMax)) continue; |
1386 | Double_t dismin=100.; | |
1387 | Double_t ptmax=-10.; | |
1388 | Int_t index1=-1; | |
1389 | Int_t index2=-1; | |
ea693273 | 1390 | Double_t fracin=0.; |
1391 | Double_t sumPtIn=0.; | |
1392 | Double_t sumPtOut=0.; | |
75bf77e3 | 1393 | //compute sum of the pt of the tracks in a concentric cone |
1394 | TRefArray *genTrackList = jetbig->GetRefTracks(); | |
1395 | Int_t nTracksGenJet = genTrackList->GetEntriesFast(); | |
1396 | AliAODTrack* genTrack; | |
1397 | for(Int_t ir=0; ir<nTracksGenJet; ++ir){ | |
1398 | genTrack = (AliAODTrack*)(genTrackList->At(ir)); | |
1399 | Float_t etr=genTrack->Eta(); | |
1400 | Float_t phir=genTrack->Phi(); | |
1401 | distr=(etr-etabig)*(etr-etabig)+(phir-phibig)*(phir-phibig); | |
1402 | distr=TMath::Sqrt(distr); | |
ea693273 | 1403 | if(distr<=fRadioFrac){ fracin=fracin+genTrack->Pt();}} |
1404 | ||
a9e585a7 | 1405 | if(centValue<10) {fh2JetCoreMethod3C10->Fill(ptcorr,fracin/ptbig); |
1406 | if(i==indexlead) fh2JetCoreMethod3C10lead->Fill(ptcorr,fracin/ptbig); | |
1407 | if(i==indexsublead) fh2JetCoreMethod3C10sublead->Fill(ptcorr,fracin/ptbig);} | |
1408 | ||
ea693273 | 1409 | |
a9e585a7 | 1410 | if((centValue>20)&&(centValue<40)) {fh2JetCoreMethod3C20->Fill(ptcorr,fracin/ptbig); |
1411 | if(i==indexlead) fh2JetCoreMethod3C20lead->Fill(ptcorr,fracin/ptbig); | |
1412 | if(i==indexsublead) fh2JetCoreMethod3C20sublead->Fill(ptcorr,fracin/ptbig);} | |
1413 | if((centValue>30)&&(centValue<60)){ fh2JetCoreMethod3C30->Fill(ptcorr,fracin/ptbig); | |
1414 | if(i==indexlead) fh2JetCoreMethod3C30lead->Fill(ptcorr,fracin/ptbig); | |
1415 | if(i==indexsublead) fh2JetCoreMethod3C30sublead->Fill(ptcorr,fracin/ptbig);} | |
1416 | if(centValue>60){fh2JetCoreMethod3C60->Fill(ptcorr,fracin/ptbig); | |
1417 | if(i==indexlead) fh2JetCoreMethod3C60lead->Fill(ptcorr,fracin/ptbig); | |
1418 | if(i==indexsublead) fh2JetCoreMethod3C60sublead->Fill(ptcorr,fracin/ptbig);} | |
1419 | ||
ea693273 | 1420 | |
75bf77e3 | 1421 | |
1422 | for(Int_t j=0; j<fListJets[1]->GetEntries(); ++j){ | |
1423 | AliAODJet* jetsmall = (AliAODJet*)(fListJets[1]->At(j)); | |
1424 | etasmall = jetsmall->Eta(); | |
1425 | phismall = jetsmall->Phi(); | |
1426 | ptsmall = jetsmall->Pt(); | |
1427 | areasmall = jetsmall->EffectiveAreaCharged(); | |
ea693273 | 1428 | Double_t tmpDeltaR=(phismall-phibig)*(phismall-phibig)+(etasmall-etabig)*(etasmall-etabig); |
1429 | tmpDeltaR=TMath::Sqrt(tmpDeltaR); | |
1430 | //Fraction in the jet core | |
1431 | if((ptsmall>ptmax)&&(tmpDeltaR<=fRadioFrac)){ptmax=ptsmall; | |
75bf77e3 | 1432 | index2=j;} |
ea693273 | 1433 | if(tmpDeltaR<=dismin){ dismin=tmpDeltaR; |
1434 | index1=j;}} //en of loop over R=0.2 jets | |
75bf77e3 | 1435 | //method1:most concentric jet=core |
1436 | if(dismin<fMinDist){ AliAODJet* jetmethod1 = (AliAODJet*)(fListJets[1]->At(index1)); | |
ea693273 | 1437 | if(centValue<10) fh2JetCoreMethod1C10->Fill(ptcorr,jetmethod1->Pt()/ptbig); |
1438 | if((centValue>20)&&(centValue<40)) fh2JetCoreMethod1C20->Fill(ptcorr,jetmethod1->Pt()/ptbig); | |
1439 | if((centValue>30)&&(centValue<60)) fh2JetCoreMethod1C30->Fill(ptcorr,jetmethod1->Pt()/ptbig); | |
1440 | if(centValue>60) fh2JetCoreMethod1C60->Fill(ptcorr,jetmethod1->Pt()/ptbig); } | |
75bf77e3 | 1441 | //method2:hardest contained jet=core |
1442 | if(index2!=-1){ | |
1443 | AliAODJet* jetmethod2 = (AliAODJet*)(fListJets[1]->At(index2)); | |
ea693273 | 1444 | if(centValue<10) fh2JetCoreMethod2C10->Fill(ptcorr,jetmethod2->Pt()/ptbig); |
1445 | if((centValue>20)&&(centValue<40)) fh2JetCoreMethod2C20->Fill(ptcorr,jetmethod2->Pt()/ptbig); | |
1446 | if((centValue>30)&&(centValue<60)) fh2JetCoreMethod2C30->Fill(ptcorr,jetmethod2->Pt()/ptbig); | |
1447 | if(centValue>60) fh2JetCoreMethod2C60->Fill(ptcorr,jetmethod2->Pt()/ptbig); } | |
1448 | ||
1449 | Double_t R=fRadioFrac*2.; | |
1450 | ||
1451 | for(int it = 0;it<nT;++it){ | |
1452 | ||
1453 | AliVParticle *part = (AliVParticle*)ParticleList.At(it); | |
1454 | Float_t deltaR = jetbig->DeltaR(part); | |
1455 | Float_t deltaEta = part->Eta()-etabig; | |
1456 | Float_t deltaPhi = part->Phi()-phibig; | |
1457 | if((deltaR>=0.4)&&(deltaR<=0.6))sumPtIn=sumPtIn+part->Pt(); | |
1458 | if((deltaR>=0.8)&&(deltaR<=1.))sumPtOut=sumPtOut+part->Pt(); | |
1459 | if(centValue<10.){ | |
a9e585a7 | 1460 | |
ea693273 | 1461 | if((ptcorr>=70.)&&(ptcorr<85.)) {fh2DeltaRC10pt1->Fill(part->Pt(),deltaR); |
1462 | if((part->Phi()<=phibig+R)&&(part->Phi()>=phibig-R))fh2DeltaEtaC10pt1->Fill(part->Pt(),deltaEta); | |
1463 | if((part->Eta()<=etabig+R)&&(part->Eta()>=etabig-R))fh2DeltaPhiC10pt1->Fill(part->Pt(),deltaPhi);} | |
1464 | if((ptcorr>=85.)&&(ptcorr<100.)) {fh2DeltaRC10pt2->Fill(part->Pt(),deltaR); | |
1465 | if((part->Phi()<=phibig+R)&&(part->Phi()>=phibig-R)) fh2DeltaEtaC10pt2->Fill(part->Pt(),deltaEta); | |
1466 | if((part->Eta()<=etabig+R)&&(part->Eta()>=etabig-R)) fh2DeltaPhiC10pt2->Fill(part->Pt(),deltaPhi); } | |
1467 | if((ptcorr>=100.)&&(ptcorr<120.)) {fh2DeltaRC10pt3->Fill(part->Pt(),deltaR); | |
1468 | if((part->Phi()<=phibig+R)&&(part->Phi()>=phibig-R)) fh2DeltaEtaC10pt3->Fill(part->Pt(),deltaEta); | |
1469 | if((part->Eta()<=etabig+R)&&(part->Eta()>=etabig-R)) fh2DeltaPhiC10pt3->Fill(part->Pt(),deltaPhi);} | |
1470 | if((ptcorr>=120.)&&(ptcorr<140.)) {fh2DeltaRC10pt4->Fill(part->Pt(),deltaR); | |
1471 | if((part->Phi()<=phibig+R)&&(part->Phi()>=phibig-R)) fh2DeltaEtaC10pt4->Fill(part->Pt(),deltaEta); | |
a9e585a7 | 1472 | if((part->Eta()<=etabig+R)&&(part->Eta()>=etabig-R))fh2DeltaPhiC10pt4->Fill(part->Pt(),deltaPhi); } |
1473 | ||
1474 | if(i==indexlead){ | |
1475 | if((ptcorr>=70.)&&(ptcorr<85.)) {fh2DeltaRC10pt1lead->Fill(part->Pt(),deltaR); | |
1476 | if((part->Phi()<=phibig+R)&&(part->Phi()>=phibig-R))fh2DeltaEtaC10pt1lead->Fill(part->Pt(),deltaEta); | |
1477 | if((part->Eta()<=etabig+R)&&(part->Eta()>=etabig-R))fh2DeltaPhiC10pt1lead->Fill(part->Pt(),deltaPhi);} | |
1478 | if((ptcorr>=85.)&&(ptcorr<100.)) {fh2DeltaRC10pt2lead->Fill(part->Pt(),deltaR); | |
1479 | if((part->Phi()<=phibig+R)&&(part->Phi()>=phibig-R)) fh2DeltaEtaC10pt2lead->Fill(part->Pt(),deltaEta); | |
1480 | if((part->Eta()<=etabig+R)&&(part->Eta()>=etabig-R)) fh2DeltaPhiC10pt2lead->Fill(part->Pt(),deltaPhi); } | |
1481 | if((ptcorr>=100.)&&(ptcorr<120.)) {fh2DeltaRC10pt3lead->Fill(part->Pt(),deltaR); | |
1482 | if((part->Phi()<=phibig+R)&&(part->Phi()>=phibig-R)) fh2DeltaEtaC10pt3lead->Fill(part->Pt(),deltaEta); | |
1483 | if((part->Eta()<=etabig+R)&&(part->Eta()>=etabig-R)) fh2DeltaPhiC10pt3lead->Fill(part->Pt(),deltaPhi);} | |
1484 | if((ptcorr>=120.)&&(ptcorr<140.)) {fh2DeltaRC10pt4lead->Fill(part->Pt(),deltaR); | |
1485 | if((part->Phi()<=phibig+R)&&(part->Phi()>=phibig-R)) fh2DeltaEtaC10pt4lead->Fill(part->Pt(),deltaEta); | |
1486 | if((part->Eta()<=etabig+R)&&(part->Eta()>=etabig-R))fh2DeltaPhiC10pt4lead->Fill(part->Pt(),deltaPhi); }} | |
1487 | ||
1488 | ||
1489 | if(i==indexsublead){ | |
1490 | if((ptcorr>=70.)&&(ptcorr<85.)) {fh2DeltaRC10pt1sublead->Fill(part->Pt(),deltaR); | |
1491 | if((part->Phi()<=phibig+R)&&(part->Phi()>=phibig-R))fh2DeltaEtaC10pt1sublead->Fill(part->Pt(),deltaEta); | |
1492 | if((part->Eta()<=etabig+R)&&(part->Eta()>=etabig-R))fh2DeltaPhiC10pt1sublead->Fill(part->Pt(),deltaPhi);} | |
1493 | if((ptcorr>=85.)&&(ptcorr<100.)) {fh2DeltaRC10pt2sublead->Fill(part->Pt(),deltaR); | |
1494 | if((part->Phi()<=phibig+R)&&(part->Phi()>=phibig-R)) fh2DeltaEtaC10pt2sublead->Fill(part->Pt(),deltaEta); | |
1495 | if((part->Eta()<=etabig+R)&&(part->Eta()>=etabig-R)) fh2DeltaPhiC10pt2sublead->Fill(part->Pt(),deltaPhi); } | |
1496 | if((ptcorr>=100.)&&(ptcorr<120.)) {fh2DeltaRC10pt3sublead->Fill(part->Pt(),deltaR); | |
1497 | if((part->Phi()<=phibig+R)&&(part->Phi()>=phibig-R)) fh2DeltaEtaC10pt3sublead->Fill(part->Pt(),deltaEta); | |
1498 | if((part->Eta()<=etabig+R)&&(part->Eta()>=etabig-R)) fh2DeltaPhiC10pt3sublead->Fill(part->Pt(),deltaPhi);} | |
1499 | if((ptcorr>=120.)&&(ptcorr<140.)) {fh2DeltaRC10pt4sublead->Fill(part->Pt(),deltaR); | |
1500 | if((part->Phi()<=phibig+R)&&(part->Phi()>=phibig-R)) fh2DeltaEtaC10pt4sublead->Fill(part->Pt(),deltaEta); | |
1501 | if((part->Eta()<=etabig+R)&&(part->Eta()>=etabig-R))fh2DeltaPhiC10pt4sublead->Fill(part->Pt(),deltaPhi); }} | |
1502 | ||
1503 | ||
1504 | ||
1505 | ||
1506 | } | |
ea693273 | 1507 | |
1508 | if((centValue>20.)&&(centValue<40.)){ | |
1509 | if((ptcorr>=70.)&&(ptcorr<85.)) {fh2DeltaRC20pt1->Fill(part->Pt(),deltaR); | |
1510 | if((part->Phi()<=phibig+R)&&(part->Phi()>=phibig-R))fh2DeltaEtaC20pt1->Fill(part->Pt(),deltaEta); | |
1511 | if((part->Eta()<=etabig+R)&&(part->Eta()>=etabig-R))fh2DeltaPhiC20pt1->Fill(part->Pt(),deltaPhi);} | |
1512 | if((ptcorr>=85.)&&(ptcorr<100.)) {fh2DeltaRC20pt2->Fill(part->Pt(),deltaR); | |
1513 | if((part->Phi()<=phibig+R)&&(part->Phi()>=phibig-R)) fh2DeltaEtaC20pt2->Fill(part->Pt(),deltaEta); | |
1514 | if((part->Eta()<=etabig+R)&&(part->Eta()>=etabig-R)) fh2DeltaPhiC20pt2->Fill(part->Pt(),deltaPhi); } | |
1515 | if((ptcorr>=100.)&&(ptcorr<120.)) {fh2DeltaRC20pt3->Fill(part->Pt(),deltaR); | |
1516 | if((part->Phi()<=phibig+R)&&(part->Phi()>=phibig-R)) fh2DeltaEtaC20pt3->Fill(part->Pt(),deltaEta); | |
1517 | if((part->Eta()<=etabig+R)&&(part->Eta()>=etabig-R)) fh2DeltaPhiC20pt3->Fill(part->Pt(),deltaPhi);} | |
1518 | if((ptcorr>=120.)&&(ptcorr<140.)) {fh2DeltaRC20pt4->Fill(part->Pt(),deltaR); | |
1519 | if((part->Phi()<=phibig+R)&&(part->Phi()>=phibig-R)) fh2DeltaEtaC20pt4->Fill(part->Pt(),deltaEta); | |
a9e585a7 | 1520 | if((part->Eta()<=etabig+R)&&(part->Eta()>=etabig-R))fh2DeltaPhiC20pt4->Fill(part->Pt(),deltaPhi); } |
1521 | ||
1522 | if(i==indexlead){ | |
1523 | if((ptcorr>=70.)&&(ptcorr<85.)) {fh2DeltaRC20pt1lead->Fill(part->Pt(),deltaR); | |
1524 | if((part->Phi()<=phibig+R)&&(part->Phi()>=phibig-R))fh2DeltaEtaC20pt1lead->Fill(part->Pt(),deltaEta); | |
1525 | if((part->Eta()<=etabig+R)&&(part->Eta()>=etabig-R))fh2DeltaPhiC20pt1lead->Fill(part->Pt(),deltaPhi);} | |
1526 | if((ptcorr>=85.)&&(ptcorr<100.)) {fh2DeltaRC20pt2lead->Fill(part->Pt(),deltaR); | |
1527 | if((part->Phi()<=phibig+R)&&(part->Phi()>=phibig-R)) fh2DeltaEtaC20pt2lead->Fill(part->Pt(),deltaEta); | |
1528 | if((part->Eta()<=etabig+R)&&(part->Eta()>=etabig-R)) fh2DeltaPhiC20pt2lead->Fill(part->Pt(),deltaPhi); } | |
1529 | if((ptcorr>=100.)&&(ptcorr<120.)) {fh2DeltaRC20pt3lead->Fill(part->Pt(),deltaR); | |
1530 | if((part->Phi()<=phibig+R)&&(part->Phi()>=phibig-R)) fh2DeltaEtaC20pt3lead->Fill(part->Pt(),deltaEta); | |
1531 | if((part->Eta()<=etabig+R)&&(part->Eta()>=etabig-R)) fh2DeltaPhiC20pt3lead->Fill(part->Pt(),deltaPhi);} | |
1532 | if((ptcorr>=120.)&&(ptcorr<140.)) {fh2DeltaRC20pt4lead->Fill(part->Pt(),deltaR); | |
1533 | if((part->Phi()<=phibig+R)&&(part->Phi()>=phibig-R)) fh2DeltaEtaC20pt4lead->Fill(part->Pt(),deltaEta); | |
1534 | if((part->Eta()<=etabig+R)&&(part->Eta()>=etabig-R))fh2DeltaPhiC20pt4lead->Fill(part->Pt(),deltaPhi); }} | |
1535 | ||
1536 | if(i==indexsublead){ | |
1537 | if((ptcorr>=70.)&&(ptcorr<85.)) {fh2DeltaRC20pt1sublead->Fill(part->Pt(),deltaR); | |
1538 | if((part->Phi()<=phibig+R)&&(part->Phi()>=phibig-R))fh2DeltaEtaC20pt1sublead->Fill(part->Pt(),deltaEta); | |
1539 | if((part->Eta()<=etabig+R)&&(part->Eta()>=etabig-R))fh2DeltaPhiC20pt1sublead->Fill(part->Pt(),deltaPhi);} | |
1540 | if((ptcorr>=85.)&&(ptcorr<100.)) {fh2DeltaRC20pt2sublead->Fill(part->Pt(),deltaR); | |
1541 | if((part->Phi()<=phibig+R)&&(part->Phi()>=phibig-R)) fh2DeltaEtaC20pt2sublead->Fill(part->Pt(),deltaEta); | |
1542 | if((part->Eta()<=etabig+R)&&(part->Eta()>=etabig-R)) fh2DeltaPhiC20pt2sublead->Fill(part->Pt(),deltaPhi); } | |
1543 | if((ptcorr>=100.)&&(ptcorr<120.)) {fh2DeltaRC20pt3sublead->Fill(part->Pt(),deltaR); | |
1544 | if((part->Phi()<=phibig+R)&&(part->Phi()>=phibig-R)) fh2DeltaEtaC20pt3sublead->Fill(part->Pt(),deltaEta); | |
1545 | if((part->Eta()<=etabig+R)&&(part->Eta()>=etabig-R)) fh2DeltaPhiC20pt3sublead->Fill(part->Pt(),deltaPhi);} | |
1546 | if((ptcorr>=120.)&&(ptcorr<140.)) {fh2DeltaRC20pt4sublead->Fill(part->Pt(),deltaR); | |
1547 | if((part->Phi()<=phibig+R)&&(part->Phi()>=phibig-R)) fh2DeltaEtaC20pt4sublead->Fill(part->Pt(),deltaEta); | |
1548 | if((part->Eta()<=etabig+R)&&(part->Eta()>=etabig-R))fh2DeltaPhiC20pt4sublead->Fill(part->Pt(),deltaPhi); }} | |
1549 | ||
1550 | } | |
ea693273 | 1551 | |
1552 | if((centValue>30.)&&(centValue<60.)){ | |
1553 | ||
1554 | if((ptcorr>=70.)&&(ptcorr<85.)) {fh2DeltaRC30pt1->Fill(part->Pt(),deltaR); | |
1555 | if((part->Phi()<=phibig+R)&&(part->Phi()>=phibig-R))fh2DeltaEtaC30pt1->Fill(part->Pt(),deltaEta); | |
1556 | if((part->Eta()<=etabig+R)&&(part->Eta()>=etabig-R))fh2DeltaPhiC30pt1->Fill(part->Pt(),deltaPhi);} | |
1557 | if((ptcorr>=85.)&&(ptcorr<100.)) {fh2DeltaRC30pt2->Fill(part->Pt(),deltaR); | |
1558 | if((part->Phi()<=phibig+R)&&(part->Phi()>=phibig-R)) fh2DeltaEtaC30pt2->Fill(part->Pt(),deltaEta); | |
1559 | if((part->Eta()<=etabig+R)&&(part->Eta()>=etabig-R)) fh2DeltaPhiC30pt2->Fill(part->Pt(),deltaPhi); } | |
1560 | if((ptcorr>=100.)&&(ptcorr<120.)) {fh2DeltaRC30pt3->Fill(part->Pt(),deltaR); | |
1561 | if((part->Phi()<=phibig+R)&&(part->Phi()>=phibig-R)) fh2DeltaEtaC30pt3->Fill(part->Pt(),deltaEta); | |
1562 | if((part->Eta()<=etabig+R)&&(part->Eta()>=etabig-R)) fh2DeltaPhiC30pt3->Fill(part->Pt(),deltaPhi);} | |
1563 | if((ptcorr>=120.)&&(ptcorr<140.)) {fh2DeltaRC30pt4->Fill(part->Pt(),deltaR); | |
1564 | if((part->Phi()<=phibig+R)&&(part->Phi()>=phibig-R)) fh2DeltaEtaC30pt4->Fill(part->Pt(),deltaEta); | |
a9e585a7 | 1565 | if((part->Eta()<=etabig+R)&&(part->Eta()>=etabig-R))fh2DeltaPhiC30pt4->Fill(part->Pt(),deltaPhi); } |
1566 | ||
1567 | if(i==indexlead){ | |
1568 | if((ptcorr>=70.)&&(ptcorr<85.)) {fh2DeltaRC30pt1lead->Fill(part->Pt(),deltaR); | |
1569 | if((part->Phi()<=phibig+R)&&(part->Phi()>=phibig-R))fh2DeltaEtaC30pt1lead->Fill(part->Pt(),deltaEta); | |
1570 | if((part->Eta()<=etabig+R)&&(part->Eta()>=etabig-R))fh2DeltaPhiC30pt1lead->Fill(part->Pt(),deltaPhi);} | |
1571 | if((ptcorr>=85.)&&(ptcorr<100.)) {fh2DeltaRC30pt2lead->Fill(part->Pt(),deltaR); | |
1572 | if((part->Phi()<=phibig+R)&&(part->Phi()>=phibig-R)) fh2DeltaEtaC30pt2lead->Fill(part->Pt(),deltaEta); | |
1573 | if((part->Eta()<=etabig+R)&&(part->Eta()>=etabig-R)) fh2DeltaPhiC30pt2lead->Fill(part->Pt(),deltaPhi); } | |
1574 | if((ptcorr>=100.)&&(ptcorr<120.)) {fh2DeltaRC30pt3lead->Fill(part->Pt(),deltaR); | |
1575 | if((part->Phi()<=phibig+R)&&(part->Phi()>=phibig-R)) fh2DeltaEtaC30pt3lead->Fill(part->Pt(),deltaEta); | |
1576 | if((part->Eta()<=etabig+R)&&(part->Eta()>=etabig-R)) fh2DeltaPhiC30pt3lead->Fill(part->Pt(),deltaPhi);} | |
1577 | if((ptcorr>=120.)&&(ptcorr<140.)) {fh2DeltaRC30pt4lead->Fill(part->Pt(),deltaR); | |
1578 | if((part->Phi()<=phibig+R)&&(part->Phi()>=phibig-R)) fh2DeltaEtaC30pt4lead->Fill(part->Pt(),deltaEta); | |
1579 | if((part->Eta()<=etabig+R)&&(part->Eta()>=etabig-R))fh2DeltaPhiC30pt4lead->Fill(part->Pt(),deltaPhi); }} | |
1580 | ||
1581 | ||
1582 | if(i==indexsublead){ | |
1583 | if((ptcorr>=70.)&&(ptcorr<85.)) {fh2DeltaRC30pt1sublead->Fill(part->Pt(),deltaR); | |
1584 | if((part->Phi()<=phibig+R)&&(part->Phi()>=phibig-R))fh2DeltaEtaC30pt1sublead->Fill(part->Pt(),deltaEta); | |
1585 | if((part->Eta()<=etabig+R)&&(part->Eta()>=etabig-R))fh2DeltaPhiC30pt1sublead->Fill(part->Pt(),deltaPhi);} | |
1586 | if((ptcorr>=85.)&&(ptcorr<100.)) {fh2DeltaRC30pt2sublead->Fill(part->Pt(),deltaR); | |
1587 | if((part->Phi()<=phibig+R)&&(part->Phi()>=phibig-R)) fh2DeltaEtaC30pt2sublead->Fill(part->Pt(),deltaEta); | |
1588 | if((part->Eta()<=etabig+R)&&(part->Eta()>=etabig-R)) fh2DeltaPhiC30pt2sublead->Fill(part->Pt(),deltaPhi); } | |
1589 | if((ptcorr>=100.)&&(ptcorr<120.)) {fh2DeltaRC30pt3sublead->Fill(part->Pt(),deltaR); | |
1590 | if((part->Phi()<=phibig+R)&&(part->Phi()>=phibig-R)) fh2DeltaEtaC30pt3sublead->Fill(part->Pt(),deltaEta); | |
1591 | if((part->Eta()<=etabig+R)&&(part->Eta()>=etabig-R)) fh2DeltaPhiC30pt3sublead->Fill(part->Pt(),deltaPhi);} | |
1592 | if((ptcorr>=120.)&&(ptcorr<140.)) {fh2DeltaRC30pt4sublead->Fill(part->Pt(),deltaR); | |
1593 | if((part->Phi()<=phibig+R)&&(part->Phi()>=phibig-R)) fh2DeltaEtaC30pt4sublead->Fill(part->Pt(),deltaEta); | |
1594 | if((part->Eta()<=etabig+R)&&(part->Eta()>=etabig-R))fh2DeltaPhiC30pt4sublead->Fill(part->Pt(),deltaPhi); }} | |
1595 | ||
1596 | ||
1597 | ||
1598 | ||
1599 | ||
1600 | } | |
ea693273 | 1601 | |
1602 | ||
1603 | if(centValue>60.){ | |
1604 | if((ptcorr>=70.)&&(ptcorr<85.)) {fh2DeltaRC60pt1->Fill(part->Pt(),deltaR); | |
1605 | if((part->Phi()<=phibig+R)&&(part->Phi()>=phibig-R))fh2DeltaEtaC60pt1->Fill(part->Pt(),deltaEta); | |
1606 | if((part->Eta()<=etabig+R)&&(part->Eta()>=etabig-R))fh2DeltaPhiC60pt1->Fill(part->Pt(),deltaPhi);} | |
1607 | if((ptcorr>=85.)&&(ptcorr<100.)) {fh2DeltaRC60pt2->Fill(part->Pt(),deltaR); | |
1608 | if((part->Phi()<=phibig+R)&&(part->Phi()>=phibig-R)) fh2DeltaEtaC60pt2->Fill(part->Pt(),deltaEta); | |
1609 | if((part->Eta()<=etabig+R)&&(part->Eta()>=etabig-R)) fh2DeltaPhiC60pt2->Fill(part->Pt(),deltaPhi); } | |
1610 | if((ptcorr>=100.)&&(ptcorr<120.)) {fh2DeltaRC60pt3->Fill(part->Pt(),deltaR); | |
1611 | if((part->Phi()<=phibig+R)&&(part->Phi()>=phibig-R)) fh2DeltaEtaC60pt3->Fill(part->Pt(),deltaEta); | |
1612 | if((part->Eta()<=etabig+R)&&(part->Eta()>=etabig-R)) fh2DeltaPhiC60pt3->Fill(part->Pt(),deltaPhi);} | |
1613 | if((ptcorr>=120.)&&(ptcorr<140.)) {fh2DeltaRC60pt4->Fill(part->Pt(),deltaR); | |
1614 | if((part->Phi()<=phibig+R)&&(part->Phi()>=phibig-R)) fh2DeltaEtaC60pt4->Fill(part->Pt(),deltaEta); | |
a9e585a7 | 1615 | if((part->Eta()<=etabig+R)&&(part->Eta()>=etabig-R))fh2DeltaPhiC60pt4->Fill(part->Pt(),deltaPhi);} |
1616 | ||
1617 | if(i==indexlead){ | |
1618 | if((ptcorr>=70.)&&(ptcorr<85.)) {fh2DeltaRC60pt1lead->Fill(part->Pt(),deltaR); | |
1619 | if((part->Phi()<=phibig+R)&&(part->Phi()>=phibig-R))fh2DeltaEtaC60pt1lead->Fill(part->Pt(),deltaEta); | |
1620 | if((part->Eta()<=etabig+R)&&(part->Eta()>=etabig-R))fh2DeltaPhiC60pt1lead->Fill(part->Pt(),deltaPhi);} | |
1621 | if((ptcorr>=85.)&&(ptcorr<100.)) {fh2DeltaRC60pt2lead->Fill(part->Pt(),deltaR); | |
1622 | if((part->Phi()<=phibig+R)&&(part->Phi()>=phibig-R)) fh2DeltaEtaC60pt2lead->Fill(part->Pt(),deltaEta); | |
1623 | if((part->Eta()<=etabig+R)&&(part->Eta()>=etabig-R)) fh2DeltaPhiC60pt2lead->Fill(part->Pt(),deltaPhi); } | |
1624 | if((ptcorr>=100.)&&(ptcorr<120.)) {fh2DeltaRC60pt3lead->Fill(part->Pt(),deltaR); | |
1625 | if((part->Phi()<=phibig+R)&&(part->Phi()>=phibig-R)) fh2DeltaEtaC60pt3lead->Fill(part->Pt(),deltaEta); | |
1626 | if((part->Eta()<=etabig+R)&&(part->Eta()>=etabig-R)) fh2DeltaPhiC60pt3lead->Fill(part->Pt(),deltaPhi);} | |
1627 | if((ptcorr>=120.)&&(ptcorr<140.)) {fh2DeltaRC60pt4lead->Fill(part->Pt(),deltaR); | |
1628 | if((part->Phi()<=phibig+R)&&(part->Phi()>=phibig-R)) fh2DeltaEtaC60pt4lead->Fill(part->Pt(),deltaEta); | |
1629 | if((part->Eta()<=etabig+R)&&(part->Eta()>=etabig-R))fh2DeltaPhiC60pt4lead->Fill(part->Pt(),deltaPhi); }} | |
1630 | ||
1631 | if(i==indexsublead){ | |
1632 | if((ptcorr>=70.)&&(ptcorr<85.)) {fh2DeltaRC60pt1sublead->Fill(part->Pt(),deltaR); | |
1633 | if((part->Phi()<=phibig+R)&&(part->Phi()>=phibig-R))fh2DeltaEtaC60pt1sublead->Fill(part->Pt(),deltaEta); | |
1634 | if((part->Eta()<=etabig+R)&&(part->Eta()>=etabig-R))fh2DeltaPhiC60pt1sublead->Fill(part->Pt(),deltaPhi);} | |
1635 | if((ptcorr>=85.)&&(ptcorr<100.)) {fh2DeltaRC60pt2sublead->Fill(part->Pt(),deltaR); | |
1636 | if((part->Phi()<=phibig+R)&&(part->Phi()>=phibig-R)) fh2DeltaEtaC60pt2sublead->Fill(part->Pt(),deltaEta); | |
1637 | if((part->Eta()<=etabig+R)&&(part->Eta()>=etabig-R)) fh2DeltaPhiC60pt2sublead->Fill(part->Pt(),deltaPhi); } | |
1638 | if((ptcorr>=100.)&&(ptcorr<120.)) {fh2DeltaRC60pt3sublead->Fill(part->Pt(),deltaR); | |
1639 | if((part->Phi()<=phibig+R)&&(part->Phi()>=phibig-R)) fh2DeltaEtaC60pt3sublead->Fill(part->Pt(),deltaEta); | |
1640 | if((part->Eta()<=etabig+R)&&(part->Eta()>=etabig-R)) fh2DeltaPhiC60pt3sublead->Fill(part->Pt(),deltaPhi);} | |
1641 | if((ptcorr>=120.)&&(ptcorr<140.)) {fh2DeltaRC60pt4sublead->Fill(part->Pt(),deltaR); | |
1642 | if((part->Phi()<=phibig+R)&&(part->Phi()>=phibig-R)) fh2DeltaEtaC60pt4sublead->Fill(part->Pt(),deltaEta); | |
1643 | if((part->Eta()<=etabig+R)&&(part->Eta()>=etabig-R))fh2DeltaPhiC60pt4sublead->Fill(part->Pt(),deltaPhi); }} | |
1644 | ||
1645 | } | |
ea693273 | 1646 | |
568f8fa2 | 1647 | } |
1648 | //end of track loop | |
ea693273 | 1649 | Double_t coronain=rho*TMath::Pi()*(1.-0.8*0.8); |
1650 | Double_t coronaout=rho*TMath::Pi()*(0.6*0.6-0.4*0.4); | |
1651 | if(centValue<10.){ | |
1652 | fh2SumPtInC10bkg->Fill(ptcorr,coronain/ptbig); | |
1653 | fh2SumPtOutC10bkg->Fill(ptcorr,coronaout/ptbig); | |
1654 | fh2SumPtInC10->Fill(ptcorr,sumPtIn/ptbig); | |
1655 | fh2SumPtOutC10->Fill(ptcorr,sumPtOut/ptbig); | |
a9e585a7 | 1656 | if(i==indexlead){ fh2SumPtInC10bkglead->Fill(ptcorr,coronain/ptbig); |
1657 | fh2SumPtOutC10bkglead->Fill(ptcorr,coronaout/ptbig); | |
1658 | fh2SumPtInC10lead->Fill(ptcorr,sumPtIn/ptbig); | |
1659 | fh2SumPtOutC10lead->Fill(ptcorr,sumPtOut/ptbig);} | |
1660 | if(i==indexsublead){ fh2SumPtInC10bkgsublead->Fill(ptcorr,coronain/ptbig); | |
1661 | fh2SumPtOutC10bkgsublead->Fill(ptcorr,coronaout/ptbig); | |
1662 | fh2SumPtInC10sublead->Fill(ptcorr,sumPtIn/ptbig); | |
1663 | fh2SumPtOutC10sublead->Fill(ptcorr,sumPtOut/ptbig);} | |
1664 | } | |
ea693273 | 1665 | if((centValue>20.)&&(centValue<40.)){ |
1666 | fh2SumPtInC20bkg->Fill(ptcorr,coronain/ptbig); | |
1667 | fh2SumPtOutC20bkg->Fill(ptcorr,coronaout/ptbig); | |
1668 | fh2SumPtInC20->Fill(ptcorr,sumPtIn/ptbig); | |
a9e585a7 | 1669 | fh2SumPtOutC20->Fill(ptcorr,sumPtOut/ptbig); |
1670 | if(i==indexlead){ fh2SumPtInC20bkglead->Fill(ptcorr,coronain/ptbig); | |
1671 | fh2SumPtOutC20bkglead->Fill(ptcorr,coronaout/ptbig); | |
1672 | fh2SumPtInC20lead->Fill(ptcorr,sumPtIn/ptbig); | |
1673 | fh2SumPtOutC20lead->Fill(ptcorr,sumPtOut/ptbig);} | |
1674 | if(i==indexsublead){ fh2SumPtInC20bkgsublead->Fill(ptcorr,coronain/ptbig); | |
1675 | fh2SumPtOutC20bkgsublead->Fill(ptcorr,coronaout/ptbig); | |
1676 | fh2SumPtInC20sublead->Fill(ptcorr,sumPtIn/ptbig); | |
1677 | fh2SumPtOutC20sublead->Fill(ptcorr,sumPtOut/ptbig);} | |
1678 | ||
1679 | ||
1680 | ||
1681 | } | |
ea693273 | 1682 | if((centValue>30.)&&(centValue<60.)){ |
1683 | fh2SumPtInC30bkg->Fill(ptcorr,coronain/ptbig); | |
1684 | fh2SumPtOutC30bkg->Fill(ptcorr,coronaout/ptbig); | |
1685 | fh2SumPtInC30->Fill(ptcorr,sumPtIn/ptbig); | |
a9e585a7 | 1686 | fh2SumPtOutC30->Fill(ptcorr,sumPtOut/ptbig); |
1687 | ||
1688 | if(i==indexlead){ fh2SumPtInC30bkglead->Fill(ptcorr,coronain/ptbig); | |
1689 | fh2SumPtOutC30bkglead->Fill(ptcorr,coronaout/ptbig); | |
1690 | fh2SumPtInC30lead->Fill(ptcorr,sumPtIn/ptbig); | |
1691 | fh2SumPtOutC30lead->Fill(ptcorr,sumPtOut/ptbig);} | |
1692 | if(i==indexsublead){ fh2SumPtInC30bkgsublead->Fill(ptcorr,coronain/ptbig); | |
1693 | fh2SumPtOutC30bkgsublead->Fill(ptcorr,coronaout/ptbig); | |
1694 | fh2SumPtInC30sublead->Fill(ptcorr,sumPtIn/ptbig); | |
1695 | fh2SumPtOutC30sublead->Fill(ptcorr,sumPtOut/ptbig);} | |
1696 | ||
1697 | } | |
ea693273 | 1698 | if(centValue>60.){ |
1699 | fh2SumPtInC60bkg->Fill(ptcorr,coronain/ptbig); | |
1700 | fh2SumPtOutC60bkg->Fill(ptcorr,coronaout/ptbig); | |
1701 | fh2SumPtInC60->Fill(ptcorr,sumPtIn/ptbig); | |
a9e585a7 | 1702 | fh2SumPtOutC60->Fill(ptcorr,sumPtOut/ptbig); |
1703 | ||
1704 | if(i==indexlead){ fh2SumPtInC60bkglead->Fill(ptcorr,coronain/ptbig); | |
1705 | fh2SumPtOutC60bkglead->Fill(ptcorr,coronaout/ptbig); | |
1706 | fh2SumPtInC60lead->Fill(ptcorr,sumPtIn/ptbig); | |
1707 | fh2SumPtOutC60lead->Fill(ptcorr,sumPtOut/ptbig);} | |
1708 | if(i==indexsublead){ fh2SumPtInC60bkgsublead->Fill(ptcorr,coronain/ptbig); | |
1709 | fh2SumPtOutC60bkgsublead->Fill(ptcorr,coronaout/ptbig); | |
1710 | fh2SumPtInC60sublead->Fill(ptcorr,sumPtIn/ptbig); | |
1711 | fh2SumPtOutC60sublead->Fill(ptcorr,sumPtOut/ptbig);} | |
1712 | ||
1713 | } | |
ea693273 | 1714 | |
1715 | //////////////////ANGULAR STRUCTURE////////////////////////////////////// | |
1716 | ||
1717 | //tracks up to R=0.8 distant from the jet axis | |
1718 | if(fAngStructCloseTracks==1){ | |
1719 | TList CloseTrackList; | |
1720 | Int_t nn=GetListOfTracksCloseToJet(&CloseTrackList,jetbig); | |
1721 | Double_t difR=0.04; | |
1722 | for(Int_t l=0;l<15;l++){ | |
1723 | Double_t rr=l*0.1+0.1; | |
1724 | for(int it = 0;it<nn;++it){ | |
1725 | AliVParticle *part1 = (AliVParticle*)CloseTrackList.At(it); | |
1726 | for(int itu=it+1;itu<CloseTrackList.GetEntries();itu++){ | |
1727 | AliVParticle *part2 = (AliVParticle*)CloseTrackList.At(itu); | |
1728 | Double_t ptm=part1->Pt(); | |
1729 | Double_t ptn=part2->Pt(); | |
1730 | Double_t Rnm = (part1->Eta()-part2->Eta())*(part1->Eta()-part2->Eta())+(part1->Phi()-part2->Phi())*(part1->Phi()-part2->Phi()); | |
1731 | Rnm=TMath::Sqrt(Rnm); | |
1732 | Double_t deltag=(1./(TMath::Sqrt(2*TMath::Pi())*difR))*TMath::Exp(-1.*(rr-Rnm)*(rr-Rnm)/(2.*difR*difR)); | |
1733 | Double_t stepf=0.5*(1.+TMath::Erf((rr-Rnm)/(TMath::Sqrt(2.)*difR))); | |
1734 | if((ptcorr<85.) && (ptcorr>=70.)){up1[l]=up1[l]+ptm*ptn*Rnm*Rnm*deltag; | |
1735 | down1[l]=down1[l]+ptm*ptn*Rnm*Rnm*stepf;} | |
1736 | if((ptcorr<100.) && (ptcorr>=85.)){up2[l]=up2[l]+ptm*ptn*Rnm*Rnm*deltag; | |
1737 | down2[l]=down2[l]+ptm*ptn*Rnm*Rnm*stepf;} | |
1738 | if((ptcorr<120.) && (ptcorr>=100.)){up3[l]=up3[l]+ptm*ptn*Rnm*Rnm*deltag; | |
1739 | down3[l]=down3[l]+ptm*ptn*Rnm*Rnm*stepf;} | |
1740 | if((ptcorr<140.) && (ptcorr>=120.)){up4[l]=up4[l]+ptm*ptn*Rnm*Rnm*deltag; | |
1741 | down4[l]=down4[l]+ptm*ptn*Rnm*Rnm*stepf;}}}} | |
1742 | } | |
1743 | ||
1744 | //only jet constituents | |
a9e585a7 | 1745 | if(fAngStructCloseTracks==2){ |
ea693273 | 1746 | |
1747 | Double_t difR=0.04; | |
1748 | for(Int_t l=0;l<15;l++){ | |
1749 | Double_t rr=l*0.1+0.1; | |
1750 | ||
1751 | ||
1752 | AliAODTrack* part1; | |
1753 | AliAODTrack* part2; | |
1754 | for(Int_t it=0; it<nTracksGenJet; ++it){ | |
1755 | part1 = (AliAODTrack*)(genTrackList->At(it)); | |
568f8fa2 | 1756 | for(Int_t itu=0; itu<nTracksGenJet; ++itu){ |
ea693273 | 1757 | part2 = (AliAODTrack*)(genTrackList->At(itu)); |
1758 | Double_t ptm=part1->Pt(); | |
1759 | Double_t ptn=part2->Pt(); | |
1760 | Double_t Rnm = (part1->Eta()-part2->Eta())*(part1->Eta()-part2->Eta())+(part1->Phi()-part2->Phi())*(part1->Phi()-part2->Phi()); | |
1761 | Rnm=TMath::Sqrt(Rnm); | |
1762 | Double_t deltag=(1./(TMath::Sqrt(2*TMath::Pi())*difR))*TMath::Exp(-1.*(rr-Rnm)*(rr-Rnm)/(2.*difR*difR)); | |
1763 | Double_t stepf=0.5*(1.+TMath::Erf((rr-Rnm)/(TMath::Sqrt(2.)*difR))); | |
1764 | if((ptcorr<85.) && (ptcorr>=70.)){up1[l]=up1[l]+ptm*ptn*Rnm*Rnm*deltag; | |
1765 | down1[l]=down1[l]+ptm*ptn*Rnm*Rnm*stepf;} | |
1766 | if((ptcorr<100.) && (ptcorr>=85.)){up2[l]=up2[l]+ptm*ptn*Rnm*Rnm*deltag; | |
1767 | down2[l]=down2[l]+ptm*ptn*Rnm*Rnm*stepf;} | |
1768 | if((ptcorr<120.) && (ptcorr>=100.)){up3[l]=up3[l]+ptm*ptn*Rnm*Rnm*deltag; | |
1769 | down3[l]=down3[l]+ptm*ptn*Rnm*Rnm*stepf;} | |
1770 | if((ptcorr<140.) && (ptcorr>=120.)){up4[l]=up4[l]+ptm*ptn*Rnm*Rnm*deltag; | |
1771 | down4[l]=down4[l]+ptm*ptn*Rnm*Rnm*stepf;}}}}} | |
75bf77e3 | 1772 | |
75bf77e3 | 1773 | |
1774 | ||
ea693273 | 1775 | |
1776 | ||
1777 | ||
1778 | ||
1779 | ||
1780 | ||
1781 | ||
1782 | ||
1783 | ||
1784 | ||
1785 | ||
1786 | ||
75bf77e3 | 1787 | } |
ea693273 | 1788 | |
1789 | ||
1790 | //end loop over R=0.4 jets | |
a9e585a7 | 1791 | if(fAngStructCloseTracks>0){ |
ea693273 | 1792 | for(Int_t l=0;l<15;l++){ |
1793 | Double_t rr=l*0.1+0.1; | |
1794 | if(down1[l]!=0){ | |
1795 | if(centValue<10.)fh2AngStructpt1C10->Fill(rr,rr*up1[l]/down1[l]); | |
1796 | if(centValue>20. && centValue<40.) fh2AngStructpt1C20->Fill(rr,rr*up1[l]/down1[l]); | |
1797 | if(centValue>30. && centValue<60.) fh2AngStructpt1C30->Fill(rr,rr*up1[l]/down1[l]); | |
1798 | if(centValue>60.) fh2AngStructpt1C60->Fill(rr,rr*up1[l]/down1[l]);} | |
1799 | if(down2[l]!=0){ | |
1800 | if(centValue<10.) fh2AngStructpt2C10->Fill(rr,rr*up2[l]/down2[l]); | |
1801 | if(centValue>20. && centValue<40.) fh2AngStructpt2C20->Fill(rr,rr*up2[l]/down2[l]); | |
1802 | if(centValue>30. && centValue<60.) fh2AngStructpt2C30->Fill(rr,rr*up2[l]/down2[l]); | |
1803 | if(centValue>60.) fh2AngStructpt2C60->Fill(rr,rr*up2[l]/down2[l]);} | |
1804 | if(down3[l]!=0){ | |
1805 | if(centValue<10.) fh2AngStructpt3C10->Fill(rr,rr*up3[l]/down3[l]); | |
1806 | if(centValue>20. && centValue<40.) fh2AngStructpt3C20->Fill(rr,rr*up3[l]/down3[l]); | |
1807 | if(centValue>30. && centValue<60.) fh2AngStructpt3C30->Fill(rr,rr*up3[l]/down3[l]); | |
1808 | if(centValue>60.) fh2AngStructpt3C60->Fill(rr,rr*up3[l]/down3[l]);} | |
1809 | if(down4[l]!=0){ | |
1810 | if(centValue<10.) fh2AngStructpt4C10->Fill(rr,rr*up4[l]/down4[l]); | |
1811 | if(centValue>20. && centValue<40.) fh2AngStructpt4C20->Fill(rr,rr*up4[l]/down4[l]); | |
1812 | if(centValue>30. && centValue<60.) fh2AngStructpt4C30->Fill(rr,rr*up4[l]/down4[l]); | |
a9e585a7 | 1813 | if(centValue>60.) fh2AngStructpt4C60->Fill(rr,rr*up4[l]/down4[l]);}}} |
ea693273 | 1814 | |
1815 | ||
1816 | ||
1817 | ||
1818 | ||
75bf77e3 | 1819 | |
1820 | ||
1821 | PostData(1, fOutputList); | |
a9e585a7 | 1822 | } |
75bf77e3 | 1823 | |
1824 | void AliAnalysisTaskJetCore::Terminate(const Option_t *) | |
1825 | { | |
1826 | // Draw result to the screen | |
1827 | // Called once at the end of the query | |
1828 | ||
1829 | if (!GetOutputData(1)) | |
1830 | return; | |
1831 | } | |
1832 | ||
ea693273 | 1833 | |
1834 | ||
1835 | ||
1836 | ||
1837 | ||
1838 | ||
1839 | ||
1840 | ||
1841 | ||
1842 | ||
1843 | Int_t AliAnalysisTaskJetCore::GetListOfTracks(TList *list){ | |
1844 | ||
1845 | Int_t iCount = 0; | |
1846 | ||
1847 | Int_t fFilterMask=272; | |
1848 | for(int it = 0;it < fAOD->GetNumberOfTracks();++it){ | |
1849 | AliAODTrack *tr = fAOD->GetTrack(it); | |
1850 | if((fFilterMask>0)&&!(tr->TestFilterBit(fFilterMask)))continue; | |
1851 | if(TMath::Abs(tr->Eta())>0.9)continue; | |
1852 | if(tr->Pt()<0.15)continue; | |
1853 | list->Add(tr); | |
1854 | //cout<<fAOD->GetNumberOfTracks()<<" "<<tr->Pt()<<endl; | |
1855 | iCount++; | |
1856 | } | |
1857 | ||
1858 | list->Sort(); | |
1859 | return iCount; | |
1860 | ||
1861 | } | |
1862 | ||
1863 | Int_t AliAnalysisTaskJetCore::GetListOfTracksCloseToJet(TList *list,AliAODJet *jetbig){ | |
1864 | ||
1865 | Int_t iCount = 0; | |
1866 | ||
1867 | Int_t fFilterMask=272; | |
1868 | for(int it = 0;it < fAOD->GetNumberOfTracks();++it){ | |
1869 | AliAODTrack *tr = fAOD->GetTrack(it); | |
1870 | if((fFilterMask>0)&&!(tr->TestFilterBit(fFilterMask)))continue; | |
1871 | if(TMath::Abs(tr->Eta())>0.9)continue; | |
1872 | if(tr->Pt()<0.15)continue; | |
1873 | Double_t disR=jetbig->DeltaR(tr); | |
1874 | if(disR>0.8) continue; | |
1875 | list->Add(tr); | |
1876 | //cout<<fAOD->GetNumberOfTracks()<<" "<<tr->Pt()<<endl; | |
1877 | iCount++; | |
1878 | } | |
1879 | ||
1880 | list->Sort(); | |
1881 | return iCount; | |
1882 | ||
1883 | } | |
1884 | ||
1885 | ||
1886 | ||
1887 | ||
1888 | ||
1889 | ||
1890 | ||
1891 | ||
1892 | ||
1893 | ||
1894 | ||
75bf77e3 | 1895 | Int_t AliAnalysisTaskJetCore::GetNInputTracks() |
1896 | { | |
1897 | ||
1898 | Int_t nInputTracks = 0; | |
1899 | ||
1900 | TString jbname(fJetBranchName[1]); | |
1901 | //needs complete event, use jets without background subtraction | |
1902 | for(Int_t i=1; i<=3; ++i){ | |
1903 | if(jbname.Contains(Form("B%d",i))) jbname.ReplaceAll(Form("B%d",i),"B0"); | |
1904 | } | |
1905 | // use only HI event | |
1906 | if(jbname.Contains("AODextraonly")) jbname.ReplaceAll("AODextraonly","AOD"); | |
1907 | if(jbname.Contains("AODextra")) jbname.ReplaceAll("AODextra","AOD"); | |
1908 | ||
1909 | if(fDebug) Printf("Multiplicity from jet branch %s", jbname.Data()); | |
1910 | TClonesArray *tmpAODjets = dynamic_cast<TClonesArray*>(fAOD->FindListObject(jbname.Data())); | |
1911 | if(!tmpAODjets){ | |
1912 | Printf("Jet branch %s not found", jbname.Data()); | |
1913 | Printf("AliAnalysisTaskJetCore::GetNInputTracks FAILED"); | |
1914 | return -1; | |
1915 | } | |
1916 | ||
1917 | for (Int_t iJet=0; iJet<tmpAODjets->GetEntriesFast(); iJet++){ | |
1918 | AliAODJet *jet = dynamic_cast<AliAODJet*>((*tmpAODjets)[iJet]); | |
1919 | if(!jet) continue; | |
1920 | TRefArray *trackList = jet->GetRefTracks(); | |
1921 | Int_t nTracks = trackList->GetEntriesFast(); | |
1922 | nInputTracks += nTracks; | |
1923 | if(fDebug) Printf("#jet%d: %d tracks", iJet, nTracks); | |
1924 | } | |
1925 | if(fDebug) Printf("---> input tracks: %d", nInputTracks); | |
1926 | ||
1927 | return nInputTracks; | |
1928 | } | |
1929 | ||
1930 | ||
1931 | ||
ea693273 | 1932 | Double_t AliAnalysisTaskJetCore::RelativePhi(Double_t mphi,Double_t vphi){ |
1933 | ||
1934 | if (vphi < -1*TMath::Pi()) vphi += (2*TMath::Pi()); | |
1935 | else if (vphi > TMath::Pi()) vphi -= (2*TMath::Pi()); | |
1936 | if (mphi < -1*TMath::Pi()) mphi += (2*TMath::Pi()); | |
1937 | else if (mphi > TMath::Pi()) mphi -= (2*TMath::Pi()); | |
1938 | double dphi = mphi-vphi; | |
1939 | if (dphi < -1*TMath::Pi()) dphi += (2*TMath::Pi()); | |
1940 | else if (dphi > TMath::Pi()) dphi -= (2*TMath::Pi()); | |
1941 | return dphi;//dphi in [-Pi, Pi] | |
1942 | } | |
1943 | ||
75bf77e3 | 1944 | |
1945 |