changed binning, vertex range, leading, subleading
[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),
ea693273 72fRadioFrac(0.2),
73fMinDist(0.1),
75bf77e3 74fCentMin(0.),
75fCentMax(100.),
76fNInputTracksMin(0),
77fNInputTracksMax(-1),
ea693273 78fAngStructCloseTracks(0),
75bf77e3 79fJetEtaMin(-.5),
80fJetEtaMax(.5),
81fJetPtMin(20.),
82fJetTriggerExcludeMask(AliAODJet::kHighTrackPtTriggered),
83fJetPtFractionMin(0.5),
84fNMatchJets(4),
85fMatchMaxDist(0.8),
86fKeepJets(kFALSE),
75bf77e3 87fkNbranches(2),
88fkEvtClasses(12),
89fOutputList(0x0),
90fbEvent(kTRUE),
91fHistEvtSelection(0x0),
92fHistJetSelection(0x0),
93fh2JetSelection(0x0),
94fh2JetCoreMethod1C10(0x0),
95fh2JetCoreMethod2C10(0x0),
96fh2JetCoreMethod3C10(0x0),
ea693273 97fh2JetCoreMethod1C20(0x0),
98fh2JetCoreMethod2C20(0x0),
99fh2JetCoreMethod3C20(0x0),
75bf77e3 100fh2JetCoreMethod1C30(0x0),
101fh2JetCoreMethod2C30(0x0),
102fh2JetCoreMethod3C30(0x0),
103fh2JetCoreMethod1C60(0x0),
104fh2JetCoreMethod2C60(0x0),
ea693273 105fh2JetCoreMethod3C60(0x0),
a9e585a7 106fh2JetCoreMethod3C10lead(0x0),
107fh2JetCoreMethod3C20lead(0x0),
108fh2JetCoreMethod3C30lead(0x0),
109fh2JetCoreMethod3C60lead(0x0),
110fh2JetCoreMethod3C10sublead(0x0),
111fh2JetCoreMethod3C20sublead(0x0),
112fh2JetCoreMethod3C30sublead(0x0),
113fh2JetCoreMethod3C60sublead(0x0),
ea693273 114fh2SumPtInC10(0x0),
115fh2SumPtInC20(0x0),
116fh2SumPtInC30(0x0),
117fh2SumPtInC60(0x0),
a9e585a7 118fh2SumPtInC10lead(0x0),
119fh2SumPtInC20lead(0x0),
120fh2SumPtInC30lead(0x0),
121fh2SumPtInC60lead(0x0),
122fh2SumPtInC10sublead(0x0),
123fh2SumPtInC20sublead(0x0),
124fh2SumPtInC30sublead(0x0),
125fh2SumPtInC60sublead(0x0),
ea693273 126fh2SumPtOutC10(0x0),
ea693273 127fh2SumPtOutC20(0x0),
128fh2SumPtOutC30(0x0),
129fh2SumPtOutC60(0x0),
a9e585a7 130fh2SumPtOutC10lead(0x0),
131fh2SumPtOutC20lead(0x0),
132fh2SumPtOutC30lead(0x0),
133fh2SumPtOutC60lead(0x0),
134fh2SumPtOutC10sublead(0x0),
135fh2SumPtOutC20sublead(0x0),
136fh2SumPtOutC30sublead(0x0),
137fh2SumPtOutC60sublead(0x0),
ea693273 138fh2SumPtInC10bkg(0x0),
139fh2SumPtInC20bkg(0x0),
140fh2SumPtInC30bkg(0x0),
141fh2SumPtInC60bkg(0x0),
a9e585a7 142fh2SumPtInC10bkglead(0x0),
143fh2SumPtInC20bkglead(0x0),
144fh2SumPtInC30bkglead(0x0),
145fh2SumPtInC60bkglead(0x0),
146fh2SumPtInC10bkgsublead(0x0),
147fh2SumPtInC20bkgsublead(0x0),
148fh2SumPtInC30bkgsublead(0x0),
149fh2SumPtInC60bkgsublead(0x0),
ea693273 150fh2SumPtOutC10bkg(0x0),
151fh2SumPtOutC20bkg(0x0),
152fh2SumPtOutC30bkg(0x0),
153fh2SumPtOutC60bkg(0x0),
a9e585a7 154fh2SumPtOutC10bkglead(0x0),
155fh2SumPtOutC20bkglead(0x0),
156fh2SumPtOutC30bkglead(0x0),
157fh2SumPtOutC60bkglead(0x0),
158fh2SumPtOutC10bkgsublead(0x0),
159fh2SumPtOutC20bkgsublead(0x0),
160fh2SumPtOutC30bkgsublead(0x0),
161fh2SumPtOutC60bkgsublead(0x0),
ea693273 162fh2DeltaRC10pt1(0x0),
163fh2DeltaRC20pt1(0x0),
164fh2DeltaRC30pt1(0x0),
165fh2DeltaRC60pt1(0x0),
166fh2DeltaRC10pt2(0x0),
167fh2DeltaRC20pt2(0x0),
168fh2DeltaRC30pt2(0x0),
169fh2DeltaRC60pt2(0x0),
170fh2DeltaRC10pt3(0x0),
171fh2DeltaRC20pt3(0x0),
172fh2DeltaRC30pt3(0x0),
173fh2DeltaRC60pt3(0x0),
174fh2DeltaRC10pt4(0x0),
175fh2DeltaRC20pt4(0x0),
176fh2DeltaRC30pt4(0x0),
177fh2DeltaRC60pt4(0x0),
178fh2DeltaEtaC10pt1(0x0),
179fh2DeltaEtaC20pt1(0x0),
180fh2DeltaEtaC30pt1(0x0),
181fh2DeltaEtaC60pt1(0x0),
182fh2DeltaEtaC10pt2(0x0),
183fh2DeltaEtaC20pt2(0x0),
184fh2DeltaEtaC30pt2(0x0),
185fh2DeltaEtaC60pt2(0x0),
186fh2DeltaEtaC10pt3(0x0),
187fh2DeltaEtaC20pt3(0x0),
188fh2DeltaEtaC30pt3(0x0),
189fh2DeltaEtaC60pt3(0x0),
190fh2DeltaEtaC10pt4(0x0),
191fh2DeltaEtaC20pt4(0x0),
192fh2DeltaEtaC30pt4(0x0),
193fh2DeltaEtaC60pt4(0x0),
194fh2DeltaPhiC10pt1(0x0),
195fh2DeltaPhiC20pt1(0x0),
196fh2DeltaPhiC30pt1(0x0),
197fh2DeltaPhiC60pt1(0x0),
198fh2DeltaPhiC10pt2(0x0),
199fh2DeltaPhiC20pt2(0x0),
200fh2DeltaPhiC30pt2(0x0),
201fh2DeltaPhiC60pt2(0x0),
202fh2DeltaPhiC10pt3(0x0),
203fh2DeltaPhiC20pt3(0x0),
204fh2DeltaPhiC30pt3(0x0),
205fh2DeltaPhiC60pt3(0x0),
206fh2DeltaPhiC10pt4(0x0),
207fh2DeltaPhiC20pt4(0x0),
208fh2DeltaPhiC30pt4(0x0),
209fh2DeltaPhiC60pt4(0x0),
a9e585a7 210fh2DeltaRC10pt1lead(0x0),
211fh2DeltaRC20pt1lead(0x0),
212fh2DeltaRC30pt1lead(0x0),
213fh2DeltaRC60pt1lead(0x0),
214fh2DeltaRC10pt2lead(0x0),
215fh2DeltaRC20pt2lead(0x0),
216fh2DeltaRC30pt2lead(0x0),
217fh2DeltaRC60pt2lead(0x0),
218fh2DeltaRC10pt3lead(0x0),
219fh2DeltaRC20pt3lead(0x0),
220fh2DeltaRC30pt3lead(0x0),
221fh2DeltaRC60pt3lead(0x0),
222fh2DeltaRC10pt4lead(0x0),
223fh2DeltaRC20pt4lead(0x0),
224fh2DeltaRC30pt4lead(0x0),
225fh2DeltaRC60pt4lead(0x0),
226fh2DeltaEtaC10pt1lead(0x0),
227fh2DeltaEtaC20pt1lead(0x0),
228fh2DeltaEtaC30pt1lead(0x0),
229fh2DeltaEtaC60pt1lead(0x0),
230fh2DeltaEtaC10pt2lead(0x0),
231fh2DeltaEtaC20pt2lead(0x0),
232fh2DeltaEtaC30pt2lead(0x0),
233fh2DeltaEtaC60pt2lead(0x0),
234fh2DeltaEtaC10pt3lead(0x0),
235fh2DeltaEtaC20pt3lead(0x0),
236fh2DeltaEtaC30pt3lead(0x0),
237fh2DeltaEtaC60pt3lead(0x0),
238fh2DeltaEtaC10pt4lead(0x0),
239fh2DeltaEtaC20pt4lead(0x0),
240fh2DeltaEtaC30pt4lead(0x0),
241fh2DeltaEtaC60pt4lead(0x0),
242fh2DeltaPhiC10pt1lead(0x0),
243fh2DeltaPhiC20pt1lead(0x0),
244fh2DeltaPhiC30pt1lead(0x0),
245fh2DeltaPhiC60pt1lead(0x0),
246fh2DeltaPhiC10pt2lead(0x0),
247fh2DeltaPhiC20pt2lead(0x0),
248fh2DeltaPhiC30pt2lead(0x0),
249fh2DeltaPhiC60pt2lead(0x0),
250fh2DeltaPhiC10pt3lead(0x0),
251fh2DeltaPhiC20pt3lead(0x0),
252fh2DeltaPhiC30pt3lead(0x0),
253fh2DeltaPhiC60pt3lead(0x0),
254fh2DeltaPhiC10pt4lead(0x0),
255fh2DeltaPhiC20pt4lead(0x0),
256fh2DeltaPhiC30pt4lead(0x0),
257fh2DeltaPhiC60pt4lead(0x0),
258fh2DeltaRC10pt1sublead(0x0),
259fh2DeltaRC20pt1sublead(0x0),
260fh2DeltaRC30pt1sublead(0x0),
261fh2DeltaRC60pt1sublead(0x0),
262fh2DeltaRC10pt2sublead(0x0),
263fh2DeltaRC20pt2sublead(0x0),
264fh2DeltaRC30pt2sublead(0x0),
265fh2DeltaRC60pt2sublead(0x0),
266fh2DeltaRC10pt3sublead(0x0),
267fh2DeltaRC20pt3sublead(0x0),
268fh2DeltaRC30pt3sublead(0x0),
269fh2DeltaRC60pt3sublead(0x0),
270fh2DeltaRC10pt4sublead(0x0),
271fh2DeltaRC20pt4sublead(0x0),
272fh2DeltaRC30pt4sublead(0x0),
273fh2DeltaRC60pt4sublead(0x0),
274fh2DeltaEtaC10pt1sublead(0x0),
275fh2DeltaEtaC20pt1sublead(0x0),
276fh2DeltaEtaC30pt1sublead(0x0),
277fh2DeltaEtaC60pt1sublead(0x0),
278fh2DeltaEtaC10pt2sublead(0x0),
279fh2DeltaEtaC20pt2sublead(0x0),
280fh2DeltaEtaC30pt2sublead(0x0),
281fh2DeltaEtaC60pt2sublead(0x0),
282fh2DeltaEtaC10pt3sublead(0x0),
283fh2DeltaEtaC20pt3sublead(0x0),
284fh2DeltaEtaC30pt3sublead(0x0),
285fh2DeltaEtaC60pt3sublead(0x0),
286fh2DeltaEtaC10pt4sublead(0x0),
287fh2DeltaEtaC20pt4sublead(0x0),
288fh2DeltaEtaC30pt4sublead(0x0),
289fh2DeltaEtaC60pt4sublead(0x0),
290fh2DeltaPhiC10pt1sublead(0x0),
291fh2DeltaPhiC20pt1sublead(0x0),
292fh2DeltaPhiC30pt1sublead(0x0),
293fh2DeltaPhiC60pt1sublead(0x0),
294fh2DeltaPhiC10pt2sublead(0x0),
295fh2DeltaPhiC20pt2sublead(0x0),
296fh2DeltaPhiC30pt2sublead(0x0),
297fh2DeltaPhiC60pt2sublead(0x0),
298fh2DeltaPhiC10pt3sublead(0x0),
299fh2DeltaPhiC20pt3sublead(0x0),
300fh2DeltaPhiC30pt3sublead(0x0),
301fh2DeltaPhiC60pt3sublead(0x0),
302fh2DeltaPhiC10pt4sublead(0x0),
303fh2DeltaPhiC20pt4sublead(0x0),
304fh2DeltaPhiC30pt4sublead(0x0),
305fh2DeltaPhiC60pt4sublead(0x0),
ea693273 306fh2AngStructpt1C10(0x0),
307fh2AngStructpt2C10(0x0),
308fh2AngStructpt3C10(0x0),
309fh2AngStructpt4C10(0x0),
310fh2AngStructpt1C20(0x0),
311fh2AngStructpt2C20(0x0),
312fh2AngStructpt3C20(0x0),
313fh2AngStructpt4C20(0x0),
314fh2AngStructpt1C30(0x0),
315fh2AngStructpt2C30(0x0),
316fh2AngStructpt3C30(0x0),
317fh2AngStructpt4C30(0x0),
318fh2AngStructpt1C60(0x0),
319fh2AngStructpt2C60(0x0),
320fh2AngStructpt3C60(0x0),
321fh2AngStructpt4C60(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
332AliAnalysisTaskJetCore::AliAnalysisTaskJetCore(const char *name) :
333AliAnalysisTaskSE(name),
334fESD(0x0),
335fAOD(0x0),
ea693273 336fAODExtension(0x0),
75bf77e3 337fBackgroundBranch(""),
ea693273 338fNonStdFile(""),
75bf77e3 339fIsPbPb(kTRUE),
340fOfflineTrgMask(AliVEvent::kAny),
341fMinContribVtx(1),
a9e585a7 342fVtxZMin(-10.),
343fVtxZMax(10.),
75bf77e3 344fEvtClassMin(0),
345fEvtClassMax(4),
ea693273 346fRadioFrac(0.2),
347fMinDist(0.1),
75bf77e3 348fCentMin(0.),
349fCentMax(100.),
350fNInputTracksMin(0),
351fNInputTracksMax(-1),
ea693273 352fAngStructCloseTracks(0),
75bf77e3 353fJetEtaMin(-.5),
354fJetEtaMax(.5),
355fJetPtMin(20.),
356fJetTriggerExcludeMask(AliAODJet::kHighTrackPtTriggered),
357fJetPtFractionMin(0.5),
358fNMatchJets(4),
359fMatchMaxDist(0.8),
360fKeepJets(kFALSE),
75bf77e3 361fkNbranches(2),
362fkEvtClasses(12),
363fOutputList(0x0),
364fbEvent(kTRUE),
365fHistEvtSelection(0x0),
366fHistJetSelection(0x0),
367fh2JetSelection(0x0),
368fh2JetCoreMethod1C10(0x0),
369fh2JetCoreMethod2C10(0x0),
370fh2JetCoreMethod3C10(0x0),
ea693273 371fh2JetCoreMethod1C20(0x0),
372fh2JetCoreMethod2C20(0x0),
373fh2JetCoreMethod3C20(0x0),
75bf77e3 374fh2JetCoreMethod1C30(0x0),
375fh2JetCoreMethod2C30(0x0),
376fh2JetCoreMethod3C30(0x0),
377fh2JetCoreMethod1C60(0x0),
378fh2JetCoreMethod2C60(0x0),
ea693273 379fh2JetCoreMethod3C60(0x0),
a9e585a7 380fh2JetCoreMethod3C10lead(0x0),
381fh2JetCoreMethod3C20lead(0x0),
382fh2JetCoreMethod3C30lead(0x0),
383fh2JetCoreMethod3C60lead(0x0),
384fh2JetCoreMethod3C10sublead(0x0),
385fh2JetCoreMethod3C20sublead(0x0),
386fh2JetCoreMethod3C30sublead(0x0),
387fh2JetCoreMethod3C60sublead(0x0),
ea693273 388fh2SumPtInC10(0x0),
389fh2SumPtInC20(0x0),
390fh2SumPtInC30(0x0),
391fh2SumPtInC60(0x0),
a9e585a7 392fh2SumPtInC10lead(0x0),
393fh2SumPtInC20lead(0x0),
394fh2SumPtInC30lead(0x0),
395fh2SumPtInC60lead(0x0),
396fh2SumPtInC10sublead(0x0),
397fh2SumPtInC20sublead(0x0),
398fh2SumPtInC30sublead(0x0),
399fh2SumPtInC60sublead(0x0),
ea693273 400fh2SumPtOutC10(0x0),
ea693273 401fh2SumPtOutC20(0x0),
402fh2SumPtOutC30(0x0),
403fh2SumPtOutC60(0x0),
a9e585a7 404fh2SumPtOutC10lead(0x0),
405fh2SumPtOutC20lead(0x0),
406fh2SumPtOutC30lead(0x0),
407fh2SumPtOutC60lead(0x0),
408fh2SumPtOutC10sublead(0x0),
409fh2SumPtOutC20sublead(0x0),
410fh2SumPtOutC30sublead(0x0),
411fh2SumPtOutC60sublead(0x0),
ea693273 412fh2SumPtInC10bkg(0x0),
413fh2SumPtInC20bkg(0x0),
414fh2SumPtInC30bkg(0x0),
415fh2SumPtInC60bkg(0x0),
a9e585a7 416fh2SumPtInC10bkglead(0x0),
417fh2SumPtInC20bkglead(0x0),
418fh2SumPtInC30bkglead(0x0),
419fh2SumPtInC60bkglead(0x0),
420fh2SumPtInC10bkgsublead(0x0),
421fh2SumPtInC20bkgsublead(0x0),
422fh2SumPtInC30bkgsublead(0x0),
423fh2SumPtInC60bkgsublead(0x0),
ea693273 424fh2SumPtOutC10bkg(0x0),
425fh2SumPtOutC20bkg(0x0),
426fh2SumPtOutC30bkg(0x0),
427fh2SumPtOutC60bkg(0x0),
a9e585a7 428fh2SumPtOutC10bkglead(0x0),
429fh2SumPtOutC20bkglead(0x0),
430fh2SumPtOutC30bkglead(0x0),
431fh2SumPtOutC60bkglead(0x0),
432fh2SumPtOutC10bkgsublead(0x0),
433fh2SumPtOutC20bkgsublead(0x0),
434fh2SumPtOutC30bkgsublead(0x0),
435fh2SumPtOutC60bkgsublead(0x0),
ea693273 436fh2DeltaRC10pt1(0x0),
437fh2DeltaRC20pt1(0x0),
438fh2DeltaRC30pt1(0x0),
439fh2DeltaRC60pt1(0x0),
440fh2DeltaRC10pt2(0x0),
441fh2DeltaRC20pt2(0x0),
442fh2DeltaRC30pt2(0x0),
443fh2DeltaRC60pt2(0x0),
444fh2DeltaRC10pt3(0x0),
445fh2DeltaRC20pt3(0x0),
446fh2DeltaRC30pt3(0x0),
447fh2DeltaRC60pt3(0x0),
448fh2DeltaRC10pt4(0x0),
449fh2DeltaRC20pt4(0x0),
450fh2DeltaRC30pt4(0x0),
451fh2DeltaRC60pt4(0x0),
452fh2DeltaEtaC10pt1(0x0),
453fh2DeltaEtaC20pt1(0x0),
454fh2DeltaEtaC30pt1(0x0),
455fh2DeltaEtaC60pt1(0x0),
456fh2DeltaEtaC10pt2(0x0),
457fh2DeltaEtaC20pt2(0x0),
458fh2DeltaEtaC30pt2(0x0),
459fh2DeltaEtaC60pt2(0x0),
460fh2DeltaEtaC10pt3(0x0),
461fh2DeltaEtaC20pt3(0x0),
462fh2DeltaEtaC30pt3(0x0),
463fh2DeltaEtaC60pt3(0x0),
464fh2DeltaEtaC10pt4(0x0),
465fh2DeltaEtaC20pt4(0x0),
466fh2DeltaEtaC30pt4(0x0),
467fh2DeltaEtaC60pt4(0x0),
468fh2DeltaPhiC10pt1(0x0),
469fh2DeltaPhiC20pt1(0x0),
470fh2DeltaPhiC30pt1(0x0),
471fh2DeltaPhiC60pt1(0x0),
472fh2DeltaPhiC10pt2(0x0),
473fh2DeltaPhiC20pt2(0x0),
474fh2DeltaPhiC30pt2(0x0),
475fh2DeltaPhiC60pt2(0x0),
476fh2DeltaPhiC10pt3(0x0),
477fh2DeltaPhiC20pt3(0x0),
478fh2DeltaPhiC30pt3(0x0),
479fh2DeltaPhiC60pt3(0x0),
480fh2DeltaPhiC10pt4(0x0),
481fh2DeltaPhiC20pt4(0x0),
482fh2DeltaPhiC30pt4(0x0),
483fh2DeltaPhiC60pt4(0x0),
a9e585a7 484fh2DeltaRC10pt1lead(0x0),
485fh2DeltaRC20pt1lead(0x0),
486fh2DeltaRC30pt1lead(0x0),
487fh2DeltaRC60pt1lead(0x0),
488fh2DeltaRC10pt2lead(0x0),
489fh2DeltaRC20pt2lead(0x0),
490fh2DeltaRC30pt2lead(0x0),
491fh2DeltaRC60pt2lead(0x0),
492fh2DeltaRC10pt3lead(0x0),
493fh2DeltaRC20pt3lead(0x0),
494fh2DeltaRC30pt3lead(0x0),
495fh2DeltaRC60pt3lead(0x0),
496fh2DeltaRC10pt4lead(0x0),
497fh2DeltaRC20pt4lead(0x0),
498fh2DeltaRC30pt4lead(0x0),
499fh2DeltaRC60pt4lead(0x0),
500fh2DeltaEtaC10pt1lead(0x0),
501fh2DeltaEtaC20pt1lead(0x0),
502fh2DeltaEtaC30pt1lead(0x0),
503fh2DeltaEtaC60pt1lead(0x0),
504fh2DeltaEtaC10pt2lead(0x0),
505fh2DeltaEtaC20pt2lead(0x0),
506fh2DeltaEtaC30pt2lead(0x0),
507fh2DeltaEtaC60pt2lead(0x0),
508fh2DeltaEtaC10pt3lead(0x0),
509fh2DeltaEtaC20pt3lead(0x0),
510fh2DeltaEtaC30pt3lead(0x0),
511fh2DeltaEtaC60pt3lead(0x0),
512fh2DeltaEtaC10pt4lead(0x0),
513fh2DeltaEtaC20pt4lead(0x0),
514fh2DeltaEtaC30pt4lead(0x0),
515fh2DeltaEtaC60pt4lead(0x0),
516fh2DeltaPhiC10pt1lead(0x0),
517fh2DeltaPhiC20pt1lead(0x0),
518fh2DeltaPhiC30pt1lead(0x0),
519fh2DeltaPhiC60pt1lead(0x0),
520fh2DeltaPhiC10pt2lead(0x0),
521fh2DeltaPhiC20pt2lead(0x0),
522fh2DeltaPhiC30pt2lead(0x0),
523fh2DeltaPhiC60pt2lead(0x0),
524fh2DeltaPhiC10pt3lead(0x0),
525fh2DeltaPhiC20pt3lead(0x0),
526fh2DeltaPhiC30pt3lead(0x0),
527fh2DeltaPhiC60pt3lead(0x0),
528fh2DeltaPhiC10pt4lead(0x0),
529fh2DeltaPhiC20pt4lead(0x0),
530fh2DeltaPhiC30pt4lead(0x0),
531fh2DeltaPhiC60pt4lead(0x0),
532fh2DeltaRC10pt1sublead(0x0),
533fh2DeltaRC20pt1sublead(0x0),
534fh2DeltaRC30pt1sublead(0x0),
535fh2DeltaRC60pt1sublead(0x0),
536fh2DeltaRC10pt2sublead(0x0),
537fh2DeltaRC20pt2sublead(0x0),
538fh2DeltaRC30pt2sublead(0x0),
539fh2DeltaRC60pt2sublead(0x0),
540fh2DeltaRC10pt3sublead(0x0),
541fh2DeltaRC20pt3sublead(0x0),
542fh2DeltaRC30pt3sublead(0x0),
543fh2DeltaRC60pt3sublead(0x0),
544fh2DeltaRC10pt4sublead(0x0),
545fh2DeltaRC20pt4sublead(0x0),
546fh2DeltaRC30pt4sublead(0x0),
547fh2DeltaRC60pt4sublead(0x0),
548fh2DeltaEtaC10pt1sublead(0x0),
549fh2DeltaEtaC20pt1sublead(0x0),
550fh2DeltaEtaC30pt1sublead(0x0),
551fh2DeltaEtaC60pt1sublead(0x0),
552fh2DeltaEtaC10pt2sublead(0x0),
553fh2DeltaEtaC20pt2sublead(0x0),
554fh2DeltaEtaC30pt2sublead(0x0),
555fh2DeltaEtaC60pt2sublead(0x0),
556fh2DeltaEtaC10pt3sublead(0x0),
557fh2DeltaEtaC20pt3sublead(0x0),
558fh2DeltaEtaC30pt3sublead(0x0),
559fh2DeltaEtaC60pt3sublead(0x0),
560fh2DeltaEtaC10pt4sublead(0x0),
561fh2DeltaEtaC20pt4sublead(0x0),
562fh2DeltaEtaC30pt4sublead(0x0),
563fh2DeltaEtaC60pt4sublead(0x0),
564fh2DeltaPhiC10pt1sublead(0x0),
565fh2DeltaPhiC20pt1sublead(0x0),
566fh2DeltaPhiC30pt1sublead(0x0),
567fh2DeltaPhiC60pt1sublead(0x0),
568fh2DeltaPhiC10pt2sublead(0x0),
569fh2DeltaPhiC20pt2sublead(0x0),
570fh2DeltaPhiC30pt2sublead(0x0),
571fh2DeltaPhiC60pt2sublead(0x0),
572fh2DeltaPhiC10pt3sublead(0x0),
573fh2DeltaPhiC20pt3sublead(0x0),
574fh2DeltaPhiC30pt3sublead(0x0),
575fh2DeltaPhiC60pt3sublead(0x0),
576fh2DeltaPhiC10pt4sublead(0x0),
577fh2DeltaPhiC20pt4sublead(0x0),
578fh2DeltaPhiC30pt4sublead(0x0),
579fh2DeltaPhiC60pt4sublead(0x0),
ea693273 580fh2AngStructpt1C10(0x0),
581fh2AngStructpt2C10(0x0),
582fh2AngStructpt3C10(0x0),
583fh2AngStructpt4C10(0x0),
584fh2AngStructpt1C20(0x0),
585fh2AngStructpt2C20(0x0),
586fh2AngStructpt3C20(0x0),
587fh2AngStructpt4C20(0x0),
588fh2AngStructpt1C30(0x0),
589fh2AngStructpt2C30(0x0),
590fh2AngStructpt3C30(0x0),
591fh2AngStructpt4C30(0x0),
592fh2AngStructpt1C60(0x0),
593fh2AngStructpt2C60(0x0),
594fh2AngStructpt3C60(0x0),
595fh2AngStructpt4C60(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
609AliAnalysisTaskJetCore::~AliAnalysisTaskJetCore()
610{
611 delete fListJets[0];
612 delete fListJets[1];
613}
614
615void AliAnalysisTaskJetCore::SetBranchNames(const TString &branch1, const TString &branch2)
616{
617 fJetBranchName[0] = branch1;
618 fJetBranchName[1] = branch2;
619}
620
621void 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
631void 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
1198void 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
1824void 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
1843Int_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 1895Int_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 1932Double_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