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