AliPhysics  a9863a5 (a9863a5)
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
AliAnalysisTaskEmcalJetCDF.cxx
Go to the documentation of this file.
1 #include <cmath>
2 
3 #include <Rtypes.h>
4 #include <TMath.h>
5 #include <TMathBase.h>
6 #include <TClonesArray.h>
7 #include <TArray.h>
8 #include <TArrayI.h>
9 
10 #include <TAxis.h>
11 #include <TH1.h>
12 #include <TH1D.h>
13 #include <TProfile.h>
14 #include <TSeqCollection.h>
15 #include <TCollection.h>
16 #include <TList.h>
17 #include <TLorentzVector.h>
18 #include <TVector2.h>
19 #include <TAttMarker.h>
20 
21 #include <AliVCluster.h>
22 #include <AliAODCaloCluster.h>
23 #include <AliESDCaloCluster.h>
24 #include <AliVTrack.h>
25 #include <AliEmcalJet.h>
26 #include <AliVEvent.h>
27 #include <AliRhoParameter.h>
28 #include <AliLog.h>
29 #include <AliJetContainer.h>
30 #include <AliParticleContainer.h>
31 #include <AliClusterContainer.h>
32 #include <AliVEvent.h>
33 #include <AliVParticle.h>
34 #include <AliLog.h>
35 
37 
38 using std::cout;
39 using std::endl;
40 
44 
45 //________________________________________________________________________
47  : AliAnalysisTaskEmcalJet ( "AliAnalysisTaskEmcalJetCDF", kTRUE ),
48  fH1( NULL ),
49  fH2( NULL ),
50  fH3( NULL ),
51  fH4( NULL ),
52  fH4c( NULL ),
53  fH5( NULL ),
54  fH6( NULL ),
55  fH6c( NULL ),
56  fH8( NULL ),
57  fH8xi( NULL ),
58  fH8_all( NULL ),
59  fH8xi_all( NULL ),
60  fH15_bin( NULL ),
61  fH15_bin_n90( NULL ),
62  fH15_bin_n85( NULL ),
63  fH15_bin_n80( NULL ),
64  fH15_bin_n75( NULL ),
65  fH15_bin_n70( NULL ),
66  fH15_bin_pt90( NULL ),
67  fH15_bin_pt85( NULL ),
68  fH15_bin_pt80( NULL ),
69  fH15_bin_pt75( NULL ),
70  fH15_bin_pt70( NULL ),
71  fH15all_bin( NULL ),
72  fH15all_bin_n90( NULL ),
73  fH15all_bin_n85( NULL ),
74  fH15all_bin_n80( NULL ),
75  fH15all_bin_n75( NULL ),
76  fH15all_bin_n70( NULL ),
77  fH15all_bin_pt90( NULL ),
78  fH15all_bin_pt85( NULL ),
79  fH15all_bin_pt80( NULL ),
80  fH15all_bin_pt75( NULL ),
81  fH15all_bin_pt70( NULL ),
82  fH20( NULL ),
83  fH20_n90( NULL ),
84  fH20_n85( NULL ),
85  fH20_n80( NULL ),
86  fH20_n75( NULL ),
87  fH20_n70( NULL ),
88  fH20_pt90( NULL ),
89  fH20_pt85( NULL ),
90  fH20_pt80( NULL ),
91  fH20_pt75( NULL ),
92  fH20_pt70( NULL ),
93  fH20all( NULL ),
94  fH20all_n90( NULL ),
95  fH20all_n85( NULL ),
96  fH20all_n80( NULL ),
97  fH20all_n75( NULL ),
98  fH20all_n70( NULL ),
99  fH20all_pt90( NULL ),
100  fH20all_pt85( NULL ),
101  fH20all_pt80( NULL ),
102  fH20all_pt75( NULL ),
103  fH20all_pt70( NULL ),
104  fHg( NULL ),
105  fHg_n90( NULL ),
106  fHg_n85( NULL ),
107  fHg_n80( NULL ),
108  fHg_n75( NULL ),
109  fHg_n70( NULL ),
110  fHg_pt90( NULL ),
111  fHg_pt85( NULL ),
112  fHg_pt80( NULL ),
113  fHg_pt75( NULL ),
114  fHg_pt70( NULL ),
115  fHptd( NULL ),
116  fHptd_n90( NULL ),
117  fHptd_n85( NULL ),
118  fHptd_n80( NULL ),
119  fHptd_n75( NULL ),
120  fHptd_n70( NULL ),
121  fHptd_pt90( NULL ),
122  fHptd_pt85( NULL ),
123  fHptd_pt80( NULL ),
124  fHptd_pt75( NULL ),
125  fHptd_pt70( NULL ),
126  fJetsCont ( NULL ),
127  fTracksCont ( NULL ),
128  fCaloClustersCont ( NULL ),
129  fTracksContArray ( NULL ),
130  idx_jetcont(0),
131  fNJets_accepted ( 0 ),
132  fNaccPart ( 0 ),
133  fNaccClus ( 0 ),
134  fHistManager("AliAnalysisTaskEmcalJetCDF")
135  {
136  // Default constructor.
137  fDebug = AliLog::GetGlobalDebugLevel();
138  // SetMakeGeneralHistograms ( kTRUE );
139  }
140 
141 //________________________________________________________________________
143  : AliAnalysisTaskEmcalJet ( name, kTRUE ),
144  fH1( NULL ),
145  fH2( NULL ),
146  fH3( NULL ),
147  fH4( NULL ),
148  fH4c( NULL ),
149  fH5( NULL ),
150  fH6( NULL ),
151  fH6c( NULL ),
152  fH8( NULL ),
153  fH8xi( NULL ),
154  fH8_all( NULL ),
155  fH8xi_all( NULL ),
156  fH15_bin( NULL ),
157  fH15_bin_n90( NULL ),
158  fH15_bin_n85( NULL ),
159  fH15_bin_n80( NULL ),
160  fH15_bin_n75( NULL ),
161  fH15_bin_n70( NULL ),
162  fH15_bin_pt90( NULL ),
163  fH15_bin_pt85( NULL ),
164  fH15_bin_pt80( NULL ),
165  fH15_bin_pt75( NULL ),
166  fH15_bin_pt70( NULL ),
167  fH15all_bin( NULL ),
168  fH15all_bin_n90( NULL ),
169  fH15all_bin_n85( NULL ),
170  fH15all_bin_n80( NULL ),
171  fH15all_bin_n75( NULL ),
172  fH15all_bin_n70( NULL ),
173  fH15all_bin_pt90( NULL ),
174  fH15all_bin_pt85( NULL ),
175  fH15all_bin_pt80( NULL ),
176  fH15all_bin_pt75( NULL ),
177  fH15all_bin_pt70( NULL ),
178  fH20( NULL ),
179  fH20_n90( NULL ),
180  fH20_n85( NULL ),
181  fH20_n80( NULL ),
182  fH20_n75( NULL ),
183  fH20_n70( NULL ),
184  fH20_pt90( NULL ),
185  fH20_pt85( NULL ),
186  fH20_pt80( NULL ),
187  fH20_pt75( NULL ),
188  fH20_pt70( NULL ),
189  fH20all( NULL ),
190  fH20all_n90( NULL ),
191  fH20all_n85( NULL ),
192  fH20all_n80( NULL ),
193  fH20all_n75( NULL ),
194  fH20all_n70( NULL ),
195  fH20all_pt90( NULL ),
196  fH20all_pt85( NULL ),
197  fH20all_pt80( NULL ),
198  fH20all_pt75( NULL ),
199  fH20all_pt70( NULL ),
200  fHg( NULL ),
201  fHg_n90( NULL ),
202  fHg_n85( NULL ),
203  fHg_n80( NULL ),
204  fHg_n75( NULL ),
205  fHg_n70( NULL ),
206  fHg_pt90( NULL ),
207  fHg_pt85( NULL ),
208  fHg_pt80( NULL ),
209  fHg_pt75( NULL ),
210  fHg_pt70( NULL ),
211  fHptd( NULL ),
212  fHptd_n90( NULL ),
213  fHptd_n85( NULL ),
214  fHptd_n80( NULL ),
215  fHptd_n75( NULL ),
216  fHptd_n70( NULL ),
217  fHptd_pt90( NULL ),
218  fHptd_pt85( NULL ),
219  fHptd_pt80( NULL ),
220  fHptd_pt75( NULL ),
221  fHptd_pt70( NULL ),
222  fJetsCont ( NULL ),
223  fTracksCont ( NULL ),
224  fCaloClustersCont ( NULL ),
225  fTracksContArray ( NULL ),
226  idx_jetcont(0),
227  fNJets_accepted ( 0 ),
228  fNaccPart ( 0 ),
229  fNaccClus ( 0 ),
230  fHistManager("AliAnalysisTaskEmcalJetCDF")
231  {
232  // Standard constructor.
233  fDebug = AliLog::GetGlobalDebugLevel();
234  // SetMakeGeneralHistograms ( kTRUE );
235  }
236 
237 //________________________________________________________________________
239  {
240  // Destructor.
241  if ( fJetsCont ) { delete fJetsCont ; fJetsCont = NULL; }
242  if ( fTracksCont ) { delete fTracksCont ; fTracksCont = NULL; }
243  if ( fCaloClustersCont ) { delete fCaloClustersCont ; fCaloClustersCont = NULL; }
244 
245  if ( fTracksContArray ) { delete fTracksContArray ; fTracksContArray = NULL; }
247 
248  }
249 
250 //________________________________________________________________________
252  {
253  // Run analysis code here, if needed. It will be executed before FillHistograms().
255  }
256 
257 //________________________________________________________________________
258 Bool_t AliAnalysisTaskEmcalJetCDF::ProcessJetContainer(Int_t idx_jet_container)
259  {
260  fJetsCont = GetJetContainer ( idx_jet_container );
262 
263  if ( !fJetsCont )
264  {
265  std::cout << "ERROR :: Jet Container not found!!!" << std::endl;
266  return kFALSE;
267  }
268 
269  // get particles connected to jets
271  if (fTracksCont)
272  {
274  fTracksCont->SetClassName ( "AliVTrack" );
275 
276  // get array of tracks
278 
279  // Multiplicity in event - accepted tracks in tracks container
281  }
282 
283  // get clusters connected to jets
285  if (fCaloClustersCont)
286  {
287  fCaloClustersCont->SetClassName ( "AliVCluster" );
288 
289  // get array of constituents
291 
292  // accepted clusters in cluster container
294  }
295 
296  // Number of Jets found in event - accepted cuts applied by JetContainer
298 
299  // protection
300  if ( ( fNJets_accepted < 1 ) || ( fNaccPart < 1 ) )
301  {
302  if ( fDebug > 1 )
303  {
304  std::cout << "accepted (fNJets || fNPart) == 0" << std::endl;
305  }
306  return kFALSE;
307  }
308 
309  if ( fDebug > 1 )
310  {
311  std::cout << "fNJets = " << fNJets_accepted << " ; fNPart = " << fNaccPart << std::endl;
312  }
313 
314  // Run analysis code here, if needed. It will be executed before FillHistograms().
315  return kTRUE; // If return kFALSE FillHistogram() will NOT be executed.
316  }
317 
318 //________________________________________________________________________
320  {
321  AliVParticle* track = NULL;
322  AliEmcalJet* jet = NULL;
323 
324  fH5->Fill ( fNJets_accepted ); // Distribution of jets in events;
325 
326  UShort_t counter_part = 0; Double_t counter_pt = 0.; // counter for npart and pt recording
327 
328  UShort_t jet_n90 = -99 ; Double_t jet_pt90 = -99.99 ;
329  UShort_t jet_n85 = -99 ; Double_t jet_pt85 = -99.99 ;
330  UShort_t jet_n80 = -99 ; Double_t jet_pt80 = -99.99 ;
331  UShort_t jet_n75 = -99 ; Double_t jet_pt75 = -99.99 ;
332  UShort_t jet_n70 = -99 ; Double_t jet_pt70 = -99.99 ;
333 
334  // variables used to compute g and ptD
335  Double_t g_tot = 0.; Double_t sum_part_pt_tot = 0.; Double_t sum_part_pt2_tot = 0.;
336 
337  Double_t g_n90 = 0.; Double_t sum_part_pt_n90 = 0.; Double_t sum_part_pt2_n90 = 0.;
338  Double_t g_n85 = 0.; Double_t sum_part_pt_n85 = 0.; Double_t sum_part_pt2_n85 = 0.;
339  Double_t g_n80 = 0.; Double_t sum_part_pt_n80 = 0.; Double_t sum_part_pt2_n80 = 0.;
340  Double_t g_n75 = 0.; Double_t sum_part_pt_n75 = 0.; Double_t sum_part_pt2_n75 = 0.;
341  Double_t g_n70 = 0.; Double_t sum_part_pt_n70 = 0.; Double_t sum_part_pt2_n70 = 0.;
342 
343  Double_t g_pt90 = 0.; Double_t sum_part_pt_pt90 = 0.; Double_t sum_part_pt2_pt90 = 0.;
344  Double_t g_pt85 = 0.; Double_t sum_part_pt_pt85 = 0.; Double_t sum_part_pt2_pt85 = 0.;
345  Double_t g_pt80 = 0.; Double_t sum_part_pt_pt80 = 0.; Double_t sum_part_pt2_pt80 = 0.;
346  Double_t g_pt75 = 0.; Double_t sum_part_pt_pt75 = 0.; Double_t sum_part_pt2_pt75 = 0.;
347  Double_t g_pt70 = 0.; Double_t sum_part_pt_pt70 = 0.; Double_t sum_part_pt2_pt70 = 0.;
348 
349 // **************************************************************
350 // ALL JETS
351 // **************************************************************
352  // jet propreties
353  Double_t jet_pt = 0. ; UShort_t jet_npart = 0; UShort_t jet_nconst = 0;
354 
355  // vector of sorted indexes of particles in jet
356  std::vector< int > jet_sorted_idxvec ;
357 
359  while ( (jet = fJetsCont->GetNextAcceptJet()) )
360  {
361  Int_t track_idx = -999; // index variable for tracks
362 
363  // jet : Sorting by p_T jet constituents
364  jet_sorted_idxvec = jet->SortConstituentsPt ( fTracksContArray );
365 
366  jet_pt = jet->Pt();
367  jet_npart = jet->GetNumberOfTracks();
368  jet_nconst = jet->GetNumberOfConstituents();
369 
370  g_tot = 0.; sum_part_pt_tot = 0.; sum_part_pt2_tot = 0.;
371 
372  g_n90 = 0.; sum_part_pt_n90 = 0.; sum_part_pt2_n90 = 0.;
373  g_n85 = 0.; sum_part_pt_n85 = 0.; sum_part_pt2_n85 = 0.;
374  g_n80 = 0.; sum_part_pt_n80 = 0.; sum_part_pt2_n80 = 0.;
375  g_n75 = 0.; sum_part_pt_n75 = 0.; sum_part_pt2_n75 = 0.;
376  g_n70 = 0.; sum_part_pt_n70 = 0.; sum_part_pt2_n70 = 0.;
377 
378  g_pt90 = 0.; sum_part_pt_pt90 = 0.; sum_part_pt2_pt90 = 0.;
379  g_pt85 = 0.; sum_part_pt_pt85 = 0.; sum_part_pt2_pt85 = 0.;
380  g_pt80 = 0.; sum_part_pt_pt80 = 0.; sum_part_pt2_pt80 = 0.;
381  g_pt75 = 0.; sum_part_pt_pt75 = 0.; sum_part_pt2_pt75 = 0.;
382  g_pt70 = 0.; sum_part_pt_pt70 = 0.; sum_part_pt2_pt70 = 0.;
383 
384  jet_n90 = ( UShort_t ) ( 0.9 * jet_npart );
385  jet_pt90 = 0.9 * jet_pt;
386 
387  jet_n85 = ( UShort_t ) ( 0.85 * jet_npart );
388  jet_pt85 = 0.85 * jet_pt;
389 
390  jet_n80 = ( UShort_t ) ( 0.8 * jet_npart );
391  jet_pt80 = 0.8 * jet_pt;
392 
393  jet_n75 = ( UShort_t ) ( 0.75 * jet_npart );
394  jet_pt75 = 0.75 * jet_pt;
395 
396  jet_n70 = ( UShort_t ) ( 0.7 * jet_npart );
397  jet_pt70 = 0.7 * jet_pt;
398 
399  fH1->Fill ( jet_pt ); // Pt distribution of jets
400  fH2->Fill ( jet->Eta() ); // Eta distribution of jets
401  fH3->Fill ( jet->Phi() ); // Phi distribution of jets
402  fH4->Fill ( jet_npart ); // Multiplicity of jets
403  fH4c->Fill ( jet_nconst ); // Multiplicity of jets - all constituents
404  fH7all->Fill ( jet_pt, jet_nconst ); // N(jet) vs P_{T} - all jets
405 
406  counter_part = 0; counter_pt = 0.; // reset counters
407 
408  for ( Size_t i = 0; i < jet_npart; i++ )
409  {
410  track_idx = jet_sorted_idxvec.at (i);
411  track = jet->TrackAt ( track_idx, fTracksContArray );
412  if (!track) { std::cout << "Parsing tracks of jets :: track not defined but it should!!!" << std::endl; continue; }
413 
414  Double_t dpart = jet->DeltaR ( track );
415  Double_t track_pt = track->Pt();
416 
417  fH8_all->Fill ( jet->GetZ ( track ) ); // Momentum distribution for jets (FF)
418  fH8xi_all->Fill ( jet->GetXi ( track ) ); // Momentum distribution for jets (FF) xi
419 
420  Double_t z_p = Z_ptot(jet,track);
421  fH8_all_p->Fill ( z_p ); // Momentum distribution for jets (FF)
422  fH8xi_all_p->Fill ( Xi (z_p) ); // Momentum distribution for jets (FF) xi
423 
424  Double_t z_pt = Z_pt(jet,track);
425  fH8_all_pt->Fill ( z_pt ); // Momentum distribution for jets (FF)
426  fH8xi_all_pt->Fill ( Xi (z_pt) ); // Momentum distribution for jets (FF) xi
427 
428  fH20all->Fill ( dpart ); // Distribution of R in leading jet
429 
430  fH15all->Fill ( dpart, track_pt ); // p_{T} track vs the Distance R from Jet1
431  fH15all_bin->Fill ( dpart, track_pt ); // p_{T} track vs the Distance R from jet axis
432 
433  fH23all->Fill ( jet_pt, dpart ); // Jet1 Size vs P_{T} - all jets
434 
435  // computing components for g and ptD in the jet tracks loop
436  g_tot += (track_pt * dpart)/jet_pt;
437  sum_part_pt_tot += track_pt;
438  sum_part_pt2_tot += TMath::Power( track_pt, 2 );
439 
440 //#############################################################################################
441  if ( counter_part <= jet_n90 )
442  {
443  fH15all_n90->Fill ( dpart, track_pt ); // p_{T} track vs the Distance R from Jet - 80% of particles
444  fH15all_bin_n90->Fill ( dpart, track_pt ); // p_{T} track vs the Distance R from Jet - 80% of particles
445  fH20all_n90->Fill ( dpart ); // Distribution of R in leading jet
446 
447  // computing components for g and ptD in the jet tracks loop
448  g_n90 += (track_pt * dpart)/jet_pt;
449  sum_part_pt_n90 += track_pt;
450  sum_part_pt2_n90 += TMath::Power( track_pt, 2 );
451  }
452 
453  if ( counter_pt <= jet_pt90 )
454  {
455  fH15all_pt90->Fill ( dpart, track_pt ); // p_{T} track vs the Distance R from Jet - 80% of pt
456  fH15all_bin_pt90->Fill ( dpart, track_pt ); // p_{T} track vs the Distance R from Jet - 80% of pt
457  fH20all_pt90->Fill ( dpart ); // Distribution of R in leading jet
458 
459  // computing components for g and ptD in the jet tracks loop
460  g_pt90 += (track_pt * dpart)/jet_pt;
461  sum_part_pt_pt90 += track_pt;
462  sum_part_pt2_pt90 += TMath::Power( track_pt, 2 );
463  }
464 
465 //#############################################################################################
466  if ( counter_part <= jet_n85 )
467  {
468  fH15all_n85->Fill ( dpart, track_pt ); // p_{T} track vs the Distance R from Jet - 80% of particles
469  fH15all_bin_n85->Fill ( dpart, track_pt ); // p_{T} track vs the Distance R from Jet - 80% of particles
470  fH20all_n85->Fill ( dpart ); // Distribution of R in leading jet
471 
472  // computing components for g and ptD in the jet tracks loop
473  g_n85 += (track_pt * dpart)/jet_pt;
474  sum_part_pt_n85 += track_pt;
475  sum_part_pt2_n85 += TMath::Power( track_pt, 2 );
476  }
477 
478  if ( counter_pt <= jet_pt85 )
479  {
480  fH15all_pt85->Fill ( dpart, track_pt ); // p_{T} track vs the Distance R from Jet - 80% of pt
481  fH15all_bin_pt85->Fill ( dpart, track_pt ); // p_{T} track vs the Distance R from Jet - 80% of pt
482  fH20all_pt85->Fill ( dpart ); // Distribution of R in leading jet
483 
484  // computing components for g and ptD in the jet tracks loop
485  g_pt85 += (track_pt * dpart)/jet_pt;
486  sum_part_pt_pt85 += track_pt;
487  sum_part_pt2_pt85 += TMath::Power( track_pt, 2 );
488  }
489 
490 //#############################################################################################
491  if ( counter_part <= jet_n80 )
492  {
493  fH15all_n80->Fill ( dpart, track_pt ); // p_{T} track vs the Distance R from Jet - 80% of particles
494  fH15all_bin_n80->Fill ( dpart, track_pt ); // p_{T} track vs the Distance R from Jet - 80% of particles
495  fH20all_n80->Fill ( dpart ); // Distribution of R in leading jet
496  fH24all->Fill ( jet_pt, dpart ); // Jet1 Size vs P_{T} - all jets - 80% of particles
497 
498  // computing components for g and ptD in the jet tracks loop
499  g_n80 += (track_pt * dpart)/jet_pt;
500  sum_part_pt_n80 += track_pt;
501  sum_part_pt2_n80 += TMath::Power( track_pt, 2 );
502  }
503 
504  if ( counter_pt <= jet_pt80 )
505  {
506  fH15all_pt80->Fill ( dpart, track_pt ); // p_{T} track vs the Distance R from Jet - 80% of pt
507  fH15all_bin_pt80->Fill ( dpart, track_pt ); // p_{T} track vs the Distance R from Jet - 80% of pt
508  fH20all_pt80->Fill ( dpart ); // Distribution of R in leading jet
509  fH25all->Fill ( jet_pt, dpart ); // Jet1 Size vs P_{T} - all jets - 80% of Pt
510 
511  // computing components for g and ptD in the jet tracks loop
512  g_pt80 += (track_pt * dpart)/jet_pt;
513  sum_part_pt_pt80 += track_pt;
514  sum_part_pt2_pt80 += TMath::Power( track_pt, 2 );
515  }
516 
517 //#############################################################################################
518  if ( counter_part <= jet_n75 )
519  {
520  fH15all_n75->Fill ( dpart, track_pt ); // p_{T} track vs the Distance R from Jet - 80% of particles
521  fH15all_bin_n75->Fill ( dpart, track_pt ); // p_{T} track vs the Distance R from Jet - 80% of particles
522  fH20all_n75->Fill ( dpart ); // Distribution of R in leading jet
523 
524  // computing components for g and ptD in the jet tracks loop
525  g_n75 += (track_pt * dpart)/jet_pt;
526  sum_part_pt_n75 += track_pt;
527  sum_part_pt2_n75 += TMath::Power( track_pt, 2 );
528  }
529 
530  if ( counter_pt <= jet_pt75 )
531  {
532  fH15all_pt75->Fill ( dpart, track_pt ); // p_{T} track vs the Distance R from Jet - 80% of pt
533  fH15all_bin_pt75->Fill ( dpart, track_pt ); // p_{T} track vs the Distance R from Jet - 80% of pt
534  fH20all_pt75->Fill ( dpart ); // Distribution of R in leading jet
535 
536  // computing components for g and ptD in the jet tracks loop
537  g_pt75 += (track_pt * dpart)/jet_pt;
538  sum_part_pt_pt75 += track_pt;
539  sum_part_pt2_pt75 += TMath::Power( track_pt, 2 );
540  }
541 
542 //#############################################################################################
543  if ( counter_part <= jet_n70 )
544  {
545  fH15all_n70->Fill ( dpart, track_pt ); // p_{T} track vs the Distance R from Jet - 80% of particles
546  fH15all_bin_n70->Fill ( dpart, track_pt ); // p_{T} track vs the Distance R from Jet - 80% of particles
547  fH20all_n80->Fill ( dpart ); // Distribution of R in leading jet
548 
549  // computing components for g and ptD in the jet tracks loop
550  g_n70 += (track_pt * dpart)/jet_pt;
551  sum_part_pt_n70 += track_pt;
552  sum_part_pt2_n70 += TMath::Power( track_pt, 2 );
553  }
554 
555  if ( counter_pt <= jet_pt70 )
556  {
557  fH15all_pt70->Fill ( dpart, track_pt ); // p_{T} track vs the Distance R from Jet - 80% of pt
558  fH15all_bin_pt70->Fill ( dpart, track_pt ); // p_{T} track vs the Distance R from Jet - 80% of pt
559  fH20all_pt70->Fill ( dpart ); // Distribution of R in leading jet
560 
561  // computing components for g and ptD in the jet tracks loop
562  g_pt70 += (track_pt * dpart)/jet_pt;
563  sum_part_pt_pt70 += track_pt;
564  sum_part_pt2_pt70 += TMath::Power( track_pt, 2 );
565  }
566  ++counter_part; counter_pt += track_pt;
567  } // end of loop over jet tracks
568 
569  fHg->Fill ( g_tot );
570  fHg_n70->Fill ( g_n70 ); fHg_pt70->Fill ( g_pt70 );
571  fHg_n75->Fill ( g_n75 ); fHg_pt75->Fill ( g_pt75 );
572  fHg_n80->Fill ( g_n80 ); fHg_pt80->Fill ( g_pt80 );
573  fHg_n85->Fill ( g_n85 ); fHg_pt85->Fill ( g_pt85 );
574  fHg_n90->Fill ( g_n90 ); fHg_pt90->Fill ( g_pt90 );
575 
576  fHptd->Fill( TMath::Sqrt(sum_part_pt2_tot)/sum_part_pt_tot );
577 
578  fHptd->Fill( TMath::Sqrt(sum_part_pt2_n70)/sum_part_pt_n70 );
579  fHptd->Fill( TMath::Sqrt(sum_part_pt2_n75)/sum_part_pt_n75 );
580  fHptd->Fill( TMath::Sqrt(sum_part_pt2_n80)/sum_part_pt_n80 );
581  fHptd->Fill( TMath::Sqrt(sum_part_pt2_n85)/sum_part_pt_n85 );
582  fHptd->Fill( TMath::Sqrt(sum_part_pt2_n90)/sum_part_pt_n90 );
583 
584  fHptd->Fill( TMath::Sqrt(sum_part_pt2_pt70)/sum_part_pt_pt70 );
585  fHptd->Fill( TMath::Sqrt(sum_part_pt2_pt75)/sum_part_pt_pt75 );
586  fHptd->Fill( TMath::Sqrt(sum_part_pt2_pt80)/sum_part_pt_pt80 );
587  fHptd->Fill( TMath::Sqrt(sum_part_pt2_pt85)/sum_part_pt_pt85 );
588  fHptd->Fill( TMath::Sqrt(sum_part_pt2_pt90)/sum_part_pt_pt90 );
589 
590  jet = fJetsCont->GetNextAcceptJet();
591  jet_sorted_idxvec.clear();
592  } // end of loop over jets
593 
594  jet = NULL; track = NULL;
595 
596 // **************************************************************
597 // LEADING JETS
598 // **************************************************************
599  AliEmcalJet* jet1 = fJetsCont->GetLeadingJet(); // internaly checked for AcceptedJet
600 
601  if ( fDebug > 1 )
602  {
603  if (jet1)
604  {
605  std::cout << "+++++++++++++++++>>>>>>>>> Leading jet found" << std::endl;
606  jet1->Print();
607  }
608  else
609  { Printf ( "Jet1 not found (did not survive cuts?)\n" ); }
610  }
611 
612  if (jet1)
613  {
614  //vector of sorted indexes of particles in leading jet
615  std::vector< int > jet1_sorted_idxvec;
616 
617  // jet1 : Sorting by p_T jet constituents
618  jet1_sorted_idxvec = jet1->SortConstituentsPt ( fTracksContArray );
619 
620  Int_t track_idx = -999; // index variable for tracks
621  Double_t jet1_pt = jet1->Pt();
622  UInt_t jet1_npart = jet1->GetNumberOfTracks();
623  UInt_t jet1_nconst = jet1->GetNumberOfConstituents();
624 
625  UShort_t jet1_n90 = ( UShort_t ) ( 0.9 * jet1_npart );
626  Double_t jet1_pt90 = 0.9 * jet1_pt;
627 
628  UShort_t jet1_n85 = ( UShort_t ) ( 0.85 * jet1_npart );
629  Double_t jet1_pt85 = 0.85 * jet1_pt;
630 
631  UShort_t jet1_n80 = ( UShort_t ) ( 0.8 * jet1_npart );
632  Double_t jet1_pt80 = 0.8 * jet1_pt;
633 
634  UShort_t jet1_n75 = ( UShort_t ) ( 0.75 * jet1_npart );
635  Double_t jet1_pt75 = 0.75 * jet1_pt;
636 
637  UShort_t jet1_n70 = ( UShort_t ) ( 0.7 * jet1_npart );
638  Double_t jet1_pt70 = 0.7 * jet1_pt;
639 
640  fH6->Fill ( jet1_npart ); // Multiplicity of jet1 - charged tracks
641  fH6c->Fill ( jet1_nconst ); // Multiplicity of jet1 - all constituents
642  fH7->Fill ( jet1_pt, jet1_nconst ); // N(jet) vs P_{T}(jet1)
643 
644  counter_part = 0; counter_pt = 0.; // reset counters
645 
646  //___________________________________________________________________________
647  // parsing tracks of jet1 (leading jet) in decreasing order of Pt
648 
649  for ( Size_t i = 0; i < jet1_npart; i++ )
650  {
651  track_idx = jet1_sorted_idxvec.at (i);
652  track = jet1->TrackAt ( track_idx, fTracksContArray );
653  if (!track) { std::cout << "Parsing tracks of jet1 :: track not defined but it should!!!" << std::endl; continue; }
654 
655  Double_t dpart = jet1->DeltaR ( track );
656  Double_t track_pt = track->Pt();
657 
658  fH8->Fill ( jet1->GetZ ( track ) ); // Momentum distribution for leading jet (FF)
659  fH8xi->Fill ( jet1->GetXi ( track ) ); // Momentum distribution for leading jet (FF) xi
660 
661  Double_t z_p = Z_ptot(jet1, track);
662  fH8_p->Fill ( z_p ); // Momentum distribution for jets (FF)
663  fH8xi_p->Fill ( Xi (z_p) ); // Momentum distribution for jets (FF) xi
664 
665  Double_t z_pt = Z_pt(jet1, track);
666  fH8_pt->Fill ( z_pt ); // Momentum distribution for jets (FF)
667  fH8xi_pt->Fill ( Xi (z_pt) ); // Momentum distribution for jets (FF) xi
668 
669  fH20->Fill ( dpart ); // Distribution of R in leading jet
670  fH15->Fill ( dpart, track_pt ); // <p_{T}> track vs the Distance R from Jet1
671  fH15_bin->Fill ( dpart, track_pt ); // p_{T} sum track vs the Distance R from Jet1
672 
673  fH23->Fill ( jet1_pt, dpart ); // Jet1 Size vs P_{T}(jet1)
674 
675  // fill histograms for 70% of particles with highest pt
676  if ( counter_part <= jet1_n70 )
677  {
678  fH15_n70->Fill ( dpart, track_pt ); // <p_{T}> track vs the Distance R from Jet1 - 80% of particles
679  fH15_bin_n70->Fill ( dpart, track_pt ); // p_{T} sum track vs the Distance R from Jet1 - 80% of particles
680  fH20_n70->Fill ( dpart ); // Distribution of R in leading jet
681  }
682  // fill histograms for 75% of particles with highest pt
683  if ( counter_part <= jet1_n75 )
684  {
685  fH15_n75->Fill ( dpart, track_pt ); // <p_{T}> track vs the Distance R from Jet1 - 80% of particles
686  fH15_bin_n75->Fill ( dpart, track_pt ); // p_{T} sum track vs the Distance R from Jet1 - 80% of particles
687  fH20_n75->Fill ( dpart ); // Distribution of R in leading jet
688  }
689  // fill histograms for 80% of particles with highest pt
690  if ( counter_part <= jet1_n80 )
691  {
692  fH15_n80->Fill ( dpart, track_pt ); // <p_{T}> track vs the Distance R from Jet1 - 80% of particles
693  fH15_bin_n80->Fill ( dpart, track_pt ); // p_{T} sum track vs the Distance R from Jet1 - 80% of particles
694  fH20_n80->Fill ( dpart ); // Distribution of R in leading jet
695  fH24->Fill ( jet1_pt, dpart ); // Jet1 Size vs P_{T}(jet1) - 80% of particles
696  }
697  // fill histograms for 85% of particles with highest pt
698  if ( counter_part <= jet1_n85 )
699  {
700  fH15_n85->Fill ( dpart, track_pt ); // <p_{T}> track vs the Distance R from Jet1 - 80% of particles
701  fH15_bin_n85->Fill ( dpart, track_pt ); // p_{T} sum track vs the Distance R from Jet1 - 80% of particles
702  fH20_n85->Fill ( dpart ); // Distribution of R in leading jet
703  }
704  // fill histograms for 90% of particles with highest pt
705  if ( counter_part <= jet1_n90 )
706  {
707  fH15_n90->Fill ( dpart, track_pt ); // <p_{T}> track vs the Distance R from Jet1 - 80% of particles
708  fH15_bin_n90->Fill ( dpart, track_pt ); // p_{T} sum track vs the Distance R from Jet1 - 80% of particles
709  fH20_n90->Fill ( dpart ); // Distribution of R in leading jet
710  }
711 
712  // fill histograms for particles that have first 70% of pt
713  if ( counter_pt <= jet1_pt70 )
714  {
715  fH15_pt70->Fill ( dpart, track_pt ); // <p_{T}> track vs the Distance R from Jet1 - 80% of pt
716  fH15_bin_pt70->Fill ( dpart, track_pt ); // p_{T} sum track vs the Distance R from Jet1 - 80% of pt
717  fH20_pt70->Fill ( dpart ); // Distribution of R in leading jet
718  }
719  // fill histograms for particles that have first 75% of pt
720  if ( counter_pt <= jet1_pt75 )
721  {
722  fH15_pt75->Fill ( dpart, track_pt ); // <p_{T}> track vs the Distance R from Jet1 - 80% of pt
723  fH15_bin_pt75->Fill ( dpart, track_pt ); // p_{T} sum track vs the Distance R from Jet1 - 80% of pt
724  fH20_pt75->Fill ( dpart ); // Distribution of R in leading jet
725  }
726  // fill histograms for particles that have first 80% of pt
727  if ( counter_pt <= jet1_pt80 )
728  {
729  fH15_pt80->Fill ( dpart, track_pt ); // <p_{T}> track vs the Distance R from Jet1 - 80% of pt
730  fH15_bin_pt80->Fill ( dpart, track_pt ); // p_{T} sum track vs the Distance R from Jet1 - 80% of pt
731  fH20_pt80->Fill ( dpart ); // Distribution of R in leading jet
732  fH25->Fill ( jet1_pt, dpart ); // Jet1 Size vs P_{T}(jet1) - 80% of Pt
733  }
734  // fill histograms for particles that have first 80% of pt
735  if ( counter_pt <= jet1_pt85 )
736  {
737  fH15_pt85->Fill ( dpart, track_pt ); // <p_{T}> track vs the Distance R from Jet1 - 80% of pt
738  fH15_bin_pt85->Fill ( dpart, track_pt ); // p_{T} sum track vs the Distance R from Jet1 - 80% of pt
739  fH20_pt85->Fill ( dpart ); // Distribution of R in leading jet
740  }
741  // fill histograms for particles that have first 80% of pt
742  if ( counter_pt <= jet1_pt90 )
743  {
744  fH15_pt90->Fill ( dpart, track_pt ); // <p_{T}> track vs the Distance R from Jet1 - 80% of pt
745  fH15_bin_pt90->Fill ( dpart, track_pt ); // p_{T} sum track vs the Distance R from Jet1 - 80% of pt
746  fH20_pt90->Fill ( dpart ); // Distribution of R in leading jet
747  }
748 
749  ++counter_part; counter_pt += track_pt;
750  } // end of loop over jet1 tracks
751  jet1_sorted_idxvec.clear();
752  } // end of jet1 (leading jet) processing
753 
754 
755  track = NULL; jet1 = NULL;
756 
757  // post data at every processing
758  PostData ( 1, fOutput ); // Post data for ALL output slots > 0 here.
759  return kTRUE;
760  }
761 
762 //________________________________________________________________________
764  {
765  // Create user output.
766 
768 
769  // Histograms
770  fOutput->SetOwner ( kTRUE );
771 
772  // Create the list of histograms. Only the list is owned.
773 
774  //____________________________________________________________________________________
775  Int_t h1_nbin = 300;
776  Double_t h1_binwidth = 1;
777  Double_t h1_low = 0;
778  Double_t h1_high = h1_low + h1_binwidth * h1_nbin; // 1GeV/bin
779  fH1 = new TH1D ( "histo1", "p_{T} distribution of jets (accepted)", h1_nbin, h1_low, h1_high );
780  fH1->SetStats ( kTRUE );
781  fH1->GetXaxis()->SetTitle ( "p_{T,jet} in GeV/c" );
782  fH1->GetYaxis()->SetTitle ( "Jets" );
783  fH1->GetXaxis()->SetTitleColor ( 1 );
784  fH1->SetMarkerStyle ( kFullCircle );
785  fOutput->Add ( fH1 );
786 
787  Int_t h2_nbin = 200;
788  Double_t h2_binwidth = 0.01;
789  Double_t h2_low = -1;
790  Double_t h2_high = h2_low + h2_binwidth * h2_nbin;
791  fH2 = new TH1D ( "histo2", "#eta distribution of jets (accepted)", h2_nbin, h2_low,
792  h2_high ); // 1 unit of rapidity / 100 bin
793  fH2->SetStats ( kTRUE );
794  fH2->GetXaxis()->SetTitle ( "#eta_{jet}" );
795  fH2->GetYaxis()->SetTitle ( "Jets" );
796  fH2->GetXaxis()->SetTitleColor ( 1 );
797  fH2->SetMarkerStyle ( kFullCircle );
798  fOutput->Add ( fH2 );
799 
800  Int_t h3_nbin = 126;
801  Double_t h3_binwidth = 0.05;
802  Double_t h3_low = 0.;
803  Double_t h3_high = h3_low + h3_binwidth * h3_nbin;
804  fH3 = new TH1D ( "histo3", "#phi distribution of jets (accepted)", h3_nbin, h3_low, h3_high );
805  fH3->SetStats ( kTRUE );
806  fH3->GetXaxis()->SetTitle ( "#phi_{jet}" );
807  fH3->GetYaxis()->SetTitle ( "Jets" );
808  fH3->GetXaxis()->SetTitleColor ( 1 );
809  fH3->SetMarkerStyle ( kFullCircle );
810  fOutput->Add ( fH3 );
811 
812  //____________________________________________________________________________________
813  Int_t h4_nbin = 100;
814  Double_t h4_binwidth = 1;
815  Double_t h4_low = 0;
816  Double_t h4_high = h4_low + h4_binwidth * h4_nbin;
817  fH4 = new TH1D ( "histo4", "Multiplicity of jets (accepted) - charged tracks", h4_nbin, h4_low, h4_high ); // 1 unit of multiplicity /bin
818  fH4->SetStats ( kTRUE );
819  fH4->GetXaxis()->SetTitle ( "N_{tracks}(jet)" );
820  fH4->GetYaxis()->SetTitle ( "Jets" );
821  fH4->GetXaxis()->SetTitleColor ( 1 );
822  fH4->SetMarkerStyle ( kFullCircle );
823  fOutput->Add ( fH4 );
824 
825  fH4c = new TH1D ( "histo4c", "Multiplicity of jets (accepted) - all constituents", h4_nbin, h4_low, h4_high ); // 1 unit of multiplicity /bin
826  fH4c->SetStats ( kTRUE );
827  fH4c->GetXaxis()->SetTitle ( "N_{tracks}(jet)" );
828  fH4c->GetYaxis()->SetTitle ( "Jets" );
829  fH4c->GetXaxis()->SetTitleColor ( 1 );
830  fH4c->SetMarkerStyle ( kFullCircle );
831  fOutput->Add ( fH4c );
832 
833  //____________________________________________________________________________________
834  Int_t h5_nbin = 100;
835  Double_t h5_binwidth = 1;
836  Double_t h5_low = 0;
837  Double_t h5_high = h5_low + h5_binwidth * h5_nbin;
838  fH5 = new TH1D ( "histo5", "Distribution of jets in events", h5_nbin, h5_low, h5_high );
839  fH5->SetStats ( kTRUE );
840  fH5->GetXaxis()->SetTitle ( "N_{jets}" );
841  fH5->GetYaxis()->SetTitle ( "Events" );
842  fH5->GetXaxis()->SetTitleColor ( 1 );
843  fH5->SetMarkerStyle ( kFullCircle );
844  fOutput->Add ( fH5 );
845 
846  //____________________________________________________________________________________
847  Int_t h6_nbin = 100;
848  Double_t h6_binwidth = 1;
849  Double_t h6_low = 0;
850  Double_t h6_high = h6_low + h6_binwidth * h6_nbin;
851  fH6 = new TH1D ( "histo6", "Jet1 Multiplicity Distribution - charged tracks", h6_nbin, h6_low, h6_high );
852  fH6->SetStats ( kTRUE );
853  fH6->GetXaxis()->SetTitle ( "N_{tracks}(jet1)" );
854  fH6->GetYaxis()->SetTitle ( "1/N_{jet1}" );
855  fH6->GetXaxis()->SetTitleColor ( 1 );
856  fH6->SetMarkerStyle ( kFullCircle );
857  fOutput->Add ( fH6 );
858 
859  fH6c = new TH1D ( "histo6c", "Jet1 Multiplicity Distribution - all constituents", h6_nbin, h6_low, h6_high );
860  fH6c->SetStats ( kTRUE );
861  fH6c->GetXaxis()->SetTitle ( "N_{tracks}(jet1)" );
862  fH6c->GetYaxis()->SetTitle ( "1/N_{jet1}" );
863  fH6c->GetXaxis()->SetTitleColor ( 1 );
864  fH6c->SetMarkerStyle ( kFullCircle );
865  fOutput->Add ( fH6c );
866 
867  //____________________________________________________________________________________
868  Int_t h7_nbin = 400;
869  Double_t h7_binwidth = 1;
870  Double_t h7_xlow = 0;
871  Double_t h7_xhigh = h7_xlow + h7_binwidth * h7_nbin;
872  fH7 = new TProfile ( "histo7", "N(jet1) vs P_{T} - jet1", h7_nbin, h7_xlow, h7_xhigh );
873  fH7->SetStats ( kTRUE );
874  fH7->GetXaxis()->SetTitle ( "p_{T}(jet1) (GeV/c)" );
875  fH7->GetYaxis()->SetTitle ( "<N_{tracks}(jet1)> in 1 GeV/c bin" );
876  fH7->GetXaxis()->SetTitleColor ( 1 );
877  fH7->SetMarkerStyle ( kFullCircle );
878  fOutput->Add ( fH7 );
879 
880  fH7all = new TProfile ( "histo7all", "N(jet1) vs P_{T} - all jets", h7_nbin, h7_xlow, h7_xhigh );
881  fH7all->SetStats ( kTRUE );
882  fH7all->GetXaxis()->SetTitle ( "p_{T}(jet) (GeV/c)" );
883  fH7all->GetYaxis()->SetTitle ( "<N_{tracks}(jet1)> in 1 GeV/c bin" );
884  fH7all->GetXaxis()->SetTitleColor ( 1 );
885  fH7all->SetMarkerStyle ( kFullCircle );
886  fOutput->Add ( fH7all );
887 
888 //__________________________________________________________________________________________________
889  Int_t h8_nbin = 101;
890  Double_t h8_binwidth = 0.01;
891  Double_t h8_low = 0;
892  Double_t h8_high = h8_low + h8_binwidth * h8_nbin;
893 //########################################################
894  // Standard implementation of Z
895  fH8 = new TH1D ( "histo8", "Momentum distribution for jet1 (FF)", h8_nbin, h8_low, h8_high );
896  fH8->SetStats ( kTRUE );
897  fH8->GetXaxis()->SetTitle ( "z" );
898  fH8->GetYaxis()->SetTitle ( "F(Z) = 1/N_{jets1} dN/dz" );
899  fH8->GetXaxis()->SetTitleColor ( 1 );
900  fH8->SetMarkerStyle ( kFullCircle );
901  fOutput->Add ( fH8 );
902 
903  fH8_all = new TH1D ( "histo8_all", "Momentum distribution for jets (FF)", h8_nbin, h8_low, h8_high );
904  fH8_all->SetStats ( kTRUE );
905  fH8_all->GetXaxis()->SetTitle ( "z" );
906  fH8_all->GetYaxis()->SetTitle ( "F(Z) = 1/N_{jets} dN/dz" );
907  fH8_all->GetXaxis()->SetTitleColor ( 1 );
908  fH8_all->SetMarkerStyle ( kFullCircle );
909  fOutput->Add ( fH8_all );
910 //########################################################
911  // P_tot implementation of Z
912  fH8_p = new TH1D ( "histo8_p", "Momentum distribution for jet1 (FF) - P_tot", h8_nbin, h8_low, h8_high );
913  fH8_p->SetStats ( kTRUE );
914  fH8_p->GetXaxis()->SetTitle ( "z = p_{track}/p_{jet1}" );
915  fH8_p->GetYaxis()->SetTitle ( "F(Z) = 1/N_{jets1} dN/dz" );
916  fH8_p->GetXaxis()->SetTitleColor ( 1 );
917  fH8_p->SetMarkerStyle ( kFullCircle );
918  fOutput->Add ( fH8_p );
919 
920  fH8_all_p = new TH1D ( "histo8_all_p", "Momentum distribution for jets (FF) - P_tot", h8_nbin, h8_low, h8_high );
921  fH8_all_p->SetStats ( kTRUE );
922  fH8_all_p->GetXaxis()->SetTitle ( "z = p_{track}/p_{jet}" );
923  fH8_all_p->GetYaxis()->SetTitle ( "F(Z) = 1/N_{jets1} dN/dz" );
924  fH8_all_p->GetXaxis()->SetTitleColor ( 1 );
925  fH8_all_p->SetMarkerStyle ( kFullCircle );
926  fOutput->Add ( fH8_all_p );
927 //########################################################
928  // Pt implementation of Z
929  fH8_pt = new TH1D ( "histo8_pt", "Momentum distribution for jet1 (FF) - Pt", h8_nbin, h8_low, h8_high );
930  fH8_pt->SetStats ( kTRUE );
931  fH8_pt->GetXaxis()->SetTitle ( "z = p_{T,track}/p_{T,jet1}" );
932  fH8_pt->GetYaxis()->SetTitle ( "F(Z) = 1/N_{jets1} dN/dz" );
933  fH8_pt->GetXaxis()->SetTitleColor ( 1 );
934  fH8_pt->SetMarkerStyle ( kFullCircle );
935  fOutput->Add ( fH8_pt );
936 
937  fH8_all_pt = new TH1D ( "histo8_all_pt", "Momentum distribution for jets (FF) - Pt", h8_nbin, h8_low, h8_high );
938  fH8_all_pt->SetStats ( kTRUE );
939  fH8_all_pt->GetXaxis()->SetTitle ( "z = p_{T,track}/p_{T,jet1}" );
940  fH8_all_pt->GetYaxis()->SetTitle ( "F(Z) = 1/N_{jets} dN/dz" );
941  fH8_all_pt->GetXaxis()->SetTitleColor ( 1 );
942  fH8_all_pt->SetMarkerStyle ( kFullCircle );
943  fOutput->Add ( fH8_all_pt );
944 //########################################################
945 
946 //__________________________________________________________________________________________________
947  Int_t h8xi_nbin = 300;
948  Double_t h8xi_binwidth = 0.05;
949  Double_t h8xi_low = 0;
950  Double_t h8xi_high = h8xi_low + h8xi_binwidth * h8xi_nbin;
951 
952 //########################################################
953  fH8xi = new TH1D ( "histo8xi", "Momentum distribution for jet1 (FF)", h8xi_nbin, h8xi_low, h8xi_high );
954  fH8xi->SetStats ( kTRUE );
955  fH8xi->GetXaxis()->SetTitle ( "#xi = ln(1/z)" );
956  fH8xi->GetYaxis()->SetTitle ( "D(#xi) = 1/N_{jets1} dN/d#xi" );
957  fH8xi->GetXaxis()->SetTitleColor ( 1 );
958  fH8xi->SetMarkerStyle ( kFullCircle );
959  fOutput->Add ( fH8xi );
960 
961  fH8xi_all = new TH1D ( "histo8xi_all", "Momentum distribution for all jets (FF)", h8xi_nbin, h8xi_low, h8xi_high );
962  fH8xi_all->SetStats ( kTRUE );
963  fH8xi_all->GetXaxis()->SetTitle ( "#xi = ln(1/z)" );
964  fH8xi_all->GetYaxis()->SetTitle ( "D(#xi) = 1/N_{jets1} dN/d#xi" );
965  fH8xi_all->GetXaxis()->SetTitleColor ( 1 );
966  fH8xi_all->SetMarkerStyle ( kFullCircle );
967  fOutput->Add ( fH8xi_all );
968 
969 //########################################################
970  fH8xi_p = new TH1D ( "histo8xi_p", "Momentum distribution for jet1 (FF) - P_tot", h8xi_nbin, h8xi_low, h8xi_high );
971  fH8xi_p->SetStats ( kTRUE );
972  fH8xi_p->GetXaxis()->SetTitle ( "#xi = ln(1/z)" );
973  fH8xi_p->GetYaxis()->SetTitle ( "D(#xi) = 1/N_{jets1} dN/d#xi" );
974  fH8xi_p->GetXaxis()->SetTitleColor ( 1 );
975  fH8xi_p->SetMarkerStyle ( kFullCircle );
976  fOutput->Add ( fH8xi_p );
977 
978  fH8xi_all_p = new TH1D ( "histo8xi_all_p", "Momentum distribution for all jets (FF) - P_tot", h8xi_nbin, h8xi_low, h8xi_high );
979  fH8xi_all_p->SetStats ( kTRUE );
980  fH8xi_all_p->GetXaxis()->SetTitle ( "#xi = ln(1/z)" );
981  fH8xi_all_p->GetYaxis()->SetTitle ( "D(#xi) = 1/N_{jets} dN/d#xi" );
982  fH8xi_all_p->GetXaxis()->SetTitleColor ( 1 );
983  fH8xi_all_p->SetMarkerStyle ( kFullCircle );
984  fOutput->Add ( fH8xi_all_p );
985 
986 //########################################################
987  fH8xi_pt = new TH1D ( "histo8xi_pt", "Momentum distribution for jet1 (FF) - Pt", h8xi_nbin, h8xi_low, h8xi_high );
988  fH8xi_pt->SetStats ( kTRUE );
989  fH8xi_pt->GetXaxis()->SetTitle ( "#xi = ln(1/z)" );
990  fH8xi_pt->GetYaxis()->SetTitle ( "D(#xi) = 1/N_{jets1} dN/d#xi" );
991  fH8xi_pt->GetXaxis()->SetTitleColor ( 1 );
992  fH8xi_pt->SetMarkerStyle ( kFullCircle );
993  fOutput->Add ( fH8xi_pt );
994 
995  fH8xi_all_pt = new TH1D ( "histo8xi_all_pt", "Momentum distribution for all jets (FF) - Pt", h8xi_nbin, h8xi_low, h8xi_high );
996  fH8xi_all_pt->SetStats ( kTRUE );
997  fH8xi_all_pt->GetXaxis()->SetTitle ( "#xi = ln(1/z)" );
998  fH8xi_all_pt->GetYaxis()->SetTitle ( "D(#xi) = 1/N_{jets} dN/d#xi" );
999  fH8xi_all_pt->GetXaxis()->SetTitleColor ( 1 );
1000  fH8xi_all_pt->SetMarkerStyle ( kFullCircle );
1001  fOutput->Add ( fH8xi_all_pt );
1002 
1003 //____________________________________________________________________________________
1004  Int_t h15_nbin = 100;
1005  Double_t h15_binwidth = 0.01;
1006  Double_t h15_xlow = 0.;
1007  Double_t h15_xhigh = h15_xlow + h15_binwidth * h15_nbin;
1008 
1009  fH15 = new TProfile ( "histo15", "<p_{T}> track in dR(jet1)", h15_nbin, h15_xlow, h15_xhigh );
1010  fH15->SetStats ( kTRUE );
1011  fH15->GetXaxis()->SetTitle ( "R" );
1012  fH15->GetYaxis()->SetTitle ( "<p_{T}> track" );
1013  fH15->GetXaxis()->SetTitleColor ( 1 );
1014  fH15->SetMarkerStyle ( kFullCircle );
1015  fOutput->Add ( fH15 );
1016 
1017 //=====================================================================================
1018  fH15_n70 = new TProfile ( "histo15_n70", "<p_{T}> track in dR(jet1) - 70% of particles", h15_nbin, h15_xlow, h15_xhigh );
1019  fH15_n70->SetStats ( kTRUE );
1020  fH15_n70->GetXaxis()->SetTitle ( "R" );
1021  fH15_n70->GetYaxis()->SetTitle ( "<p_{T}> track" );
1022  fH15_n70->GetXaxis()->SetTitleColor ( 1 );
1023  fH15_n70->SetMarkerStyle ( kFullCircle );
1024  fOutput->Add ( fH15_n70 );
1025 
1026  fH15_n75 = new TProfile ( "histo15_n75", "<p_{T}> track in dR(jet1) - 75% of particles", h15_nbin, h15_xlow, h15_xhigh );
1027  fH15_n75->SetStats ( kTRUE );
1028  fH15_n75->GetXaxis()->SetTitle ( "R" );
1029  fH15_n75->GetYaxis()->SetTitle ( "<p_{T}> track" );
1030  fH15_n75->GetXaxis()->SetTitleColor ( 1 );
1031  fH15_n75->SetMarkerStyle ( kFullCircle );
1032  fOutput->Add ( fH15_n75 );
1033 
1034  fH15_n80 = new TProfile ( "histo15_n80", "<p_{T}> track in dR(jet1) - 80% of particles", h15_nbin, h15_xlow, h15_xhigh );
1035  fH15_n80->SetStats ( kTRUE );
1036  fH15_n80->GetXaxis()->SetTitle ( "R" );
1037  fH15_n80->GetYaxis()->SetTitle ( "<p_{T}> track" );
1038  fH15_n80->GetXaxis()->SetTitleColor ( 1 );
1039  fH15_n80->SetMarkerStyle ( kFullCircle );
1040  fOutput->Add ( fH15_n80 );
1041 
1042  fH15_n85 = new TProfile ( "histo15_n85", "<p_{T}> track in dR(jet1) - 85% of particles", h15_nbin, h15_xlow, h15_xhigh );
1043  fH15_n85->SetStats ( kTRUE );
1044  fH15_n85->GetXaxis()->SetTitle ( "R" );
1045  fH15_n85->GetYaxis()->SetTitle ( "<p_{T}> track" );
1046  fH15_n85->GetXaxis()->SetTitleColor ( 1 );
1047  fH15_n85->SetMarkerStyle ( kFullCircle );
1048  fOutput->Add ( fH15_n85 );
1049 
1050  fH15_n90 = new TProfile ( "histo15_n90", "<p_{T}> track in dR(jet1) - 90% of particles", h15_nbin, h15_xlow, h15_xhigh );
1051  fH15_n90->SetStats ( kTRUE );
1052  fH15_n90->GetXaxis()->SetTitle ( "R" );
1053  fH15_n90->GetYaxis()->SetTitle ( "<p_{T}> track" );
1054  fH15_n90->GetXaxis()->SetTitleColor ( 1 );
1055  fH15_n90->SetMarkerStyle ( kFullCircle );
1056  fOutput->Add ( fH15_n90 );
1057 
1058 
1059 //=====================================================================================
1060  fH15_pt70 = new TProfile ( "histo15_pt70", "<p_{T}> track in dR(jet1) - 70% of Pt", h15_nbin, h15_xlow, h15_xhigh );
1061  fH15_pt70->SetStats ( kTRUE );
1062  fH15_pt70->GetXaxis()->SetTitle ( "R" );
1063  fH15_pt70->GetYaxis()->SetTitle ( "<p_{T}> track" );
1064  fH15_pt70->GetXaxis()->SetTitleColor ( 1 );
1065  fH15_pt70->SetMarkerStyle ( kFullCircle );
1066  fOutput->Add ( fH15_pt70 );
1067 
1068  fH15_pt75 = new TProfile ( "histo15_pt75", "<p_{T}> track in dR(jet1) - 75% of Pt", h15_nbin, h15_xlow, h15_xhigh );
1069  fH15_pt75->SetStats ( kTRUE );
1070  fH15_pt75->GetXaxis()->SetTitle ( "R" );
1071  fH15_pt75->GetYaxis()->SetTitle ( "<p_{T}> track" );
1072  fH15_pt75->GetXaxis()->SetTitleColor ( 1 );
1073  fH15_pt75->SetMarkerStyle ( kFullCircle );
1074  fOutput->Add ( fH15_pt75 );
1075 
1076  fH15_pt80 = new TProfile ( "histo15_pt80", "<p_{T}> track in dR(jet1) - 80% of Pt", h15_nbin, h15_xlow, h15_xhigh );
1077  fH15_pt80->SetStats ( kTRUE );
1078  fH15_pt80->GetXaxis()->SetTitle ( "R" );
1079  fH15_pt80->GetYaxis()->SetTitle ( "<p_{T}> track" );
1080  fH15_pt80->GetXaxis()->SetTitleColor ( 1 );
1081  fH15_pt80->SetMarkerStyle ( kFullCircle );
1082  fOutput->Add ( fH15_pt80 );
1083 
1084  fH15_pt85 = new TProfile ( "histo15_pt85", "<p_{T}> track in dR(jet1) - 85% of Pt", h15_nbin, h15_xlow, h15_xhigh );
1085  fH15_pt85->SetStats ( kTRUE );
1086  fH15_pt85->GetXaxis()->SetTitle ( "R" );
1087  fH15_pt85->GetYaxis()->SetTitle ( "<p_{T}> track" );
1088  fH15_pt85->GetXaxis()->SetTitleColor ( 1 );
1089  fH15_pt85->SetMarkerStyle ( kFullCircle );
1090  fOutput->Add ( fH15_pt85 );
1091 
1092  fH15_pt90 = new TProfile ( "histo15_pt90", "<p_{T}> track in dR(jet1) - 90% of Pt", h15_nbin, h15_xlow, h15_xhigh );
1093  fH15_pt90->SetStats ( kTRUE );
1094  fH15_pt90->GetXaxis()->SetTitle ( "R" );
1095  fH15_pt90->GetYaxis()->SetTitle ( "<p_{T}> track" );
1096  fH15_pt90->GetXaxis()->SetTitleColor ( 1 );
1097  fH15_pt90->SetMarkerStyle ( kFullCircle );
1098  fOutput->Add ( fH15_pt90 );
1099 
1100 
1101 
1102 //__________________________________________________________________________________________________
1103  fH15all = new TProfile ( "histo15all", "<p_{T}> track in dR(jet1)", h15_nbin, h15_xlow, h15_xhigh );
1104  fH15all->SetStats ( kTRUE );
1105  fH15all->GetXaxis()->SetTitle ( "R" );
1106  fH15all->GetYaxis()->SetTitle ( "<p_{T}> track" );
1107  fH15all->GetXaxis()->SetTitleColor ( 1 );
1108  fH15all->SetMarkerStyle ( kFullCircle );
1109  fOutput->Add ( fH15all );
1110 
1111 //=====================================================================================
1112  fH15all_n70 = new TProfile ( "histo15all_n70", "<p_{T}> track in dR(jet1) - 70% of particles", h15_nbin, h15_xlow, h15_xhigh );
1113  fH15all_n70->SetStats ( kTRUE );
1114  fH15all_n70->GetXaxis()->SetTitle ( "R" );
1115  fH15all_n70->GetYaxis()->SetTitle ( "<p_{T}> track" );
1116  fH15all_n70->GetXaxis()->SetTitleColor ( 1 );
1117  fH15all_n70->SetMarkerStyle ( kFullCircle );
1118  fOutput->Add ( fH15all_n70 );
1119 
1120  fH15all_n75 = new TProfile ( "histo15all_n75", "<p_{T}> track in dR(jet1) - 75% of particles", h15_nbin, h15_xlow, h15_xhigh );
1121  fH15all_n75->SetStats ( kTRUE );
1122  fH15all_n75->GetXaxis()->SetTitle ( "R" );
1123  fH15all_n75->GetYaxis()->SetTitle ( "<p_{T}> track" );
1124  fH15all_n75->GetXaxis()->SetTitleColor ( 1 );
1125  fH15all_n75->SetMarkerStyle ( kFullCircle );
1126  fOutput->Add ( fH15all_n75 );
1127 
1128  fH15all_n80 = new TProfile ( "histo15all_n80", "<p_{T}> track in dR(jet1) - 80% of particles", h15_nbin, h15_xlow, h15_xhigh );
1129  fH15all_n80->SetStats ( kTRUE );
1130  fH15all_n80->GetXaxis()->SetTitle ( "R" );
1131  fH15all_n80->GetYaxis()->SetTitle ( "<p_{T}> track" );
1132  fH15all_n80->GetXaxis()->SetTitleColor ( 1 );
1133  fH15all_n80->SetMarkerStyle ( kFullCircle );
1134  fOutput->Add ( fH15all_n80 );
1135 
1136  fH15all_n85 = new TProfile ( "histo15all_n85", "<p_{T}> track in dR(jet1) - 85% of particles", h15_nbin, h15_xlow, h15_xhigh );
1137  fH15all_n85->SetStats ( kTRUE );
1138  fH15all_n85->GetXaxis()->SetTitle ( "R" );
1139  fH15all_n85->GetYaxis()->SetTitle ( "<p_{T}> track" );
1140  fH15all_n85->GetXaxis()->SetTitleColor ( 1 );
1141  fH15all_n85->SetMarkerStyle ( kFullCircle );
1142  fOutput->Add ( fH15all_n85 );
1143 
1144  fH15all_n90 = new TProfile ( "histo15all_n90", "<p_{T}> track in dR(jet1) - 90% of particles", h15_nbin, h15_xlow, h15_xhigh );
1145  fH15all_n90->SetStats ( kTRUE );
1146  fH15all_n90->GetXaxis()->SetTitle ( "R" );
1147  fH15all_n90->GetYaxis()->SetTitle ( "<p_{T}> track" );
1148  fH15all_n90->GetXaxis()->SetTitleColor ( 1 );
1149  fH15all_n90->SetMarkerStyle ( kFullCircle );
1150  fOutput->Add ( fH15all_n90 );
1151 
1152 //=====================================================================================
1153  fH15all_pt70 = new TProfile ( "histo15all_pt70", "<p_{T}> track in dR(jet1) - 70% of Pt", h15_nbin, h15_xlow, h15_xhigh );
1154  fH15all_pt70->SetStats ( kTRUE );
1155  fH15all_pt70->GetXaxis()->SetTitle ( "R" );
1156  fH15all_pt70->GetYaxis()->SetTitle ( "<p_{T}> track" );
1157  fH15all_pt70->GetXaxis()->SetTitleColor ( 1 );
1158  fH15all_pt70->SetMarkerStyle ( kFullCircle );
1159  fOutput->Add ( fH15all_pt70 );
1160 
1161  fH15all_pt75 = new TProfile ( "histo15all_pt75", "<p_{T}> track in dR(jet1) - 75% of Pt", h15_nbin, h15_xlow, h15_xhigh );
1162  fH15all_pt75->SetStats ( kTRUE );
1163  fH15all_pt75->GetXaxis()->SetTitle ( "R" );
1164  fH15all_pt75->GetYaxis()->SetTitle ( "<p_{T}> track" );
1165  fH15all_pt75->GetXaxis()->SetTitleColor ( 1 );
1166  fH15all_pt75->SetMarkerStyle ( kFullCircle );
1167  fOutput->Add ( fH15all_pt75 );
1168 
1169  fH15all_pt80 = new TProfile ( "histo15all_pt80", "<p_{T}> track in dR(jet1) - 80% of Pt", h15_nbin, h15_xlow, h15_xhigh );
1170  fH15all_pt80->SetStats ( kTRUE );
1171  fH15all_pt80->GetXaxis()->SetTitle ( "R" );
1172  fH15all_pt80->GetYaxis()->SetTitle ( "<p_{T}> track" );
1173  fH15all_pt80->GetXaxis()->SetTitleColor ( 1 );
1174  fH15all_pt80->SetMarkerStyle ( kFullCircle );
1175  fOutput->Add ( fH15all_pt80 );
1176 
1177  fH15all_pt85 = new TProfile ( "histo15all_pt85", "<p_{T}> track in dR(jet1) - 85% of Pt", h15_nbin, h15_xlow, h15_xhigh );
1178  fH15all_pt85->SetStats ( kTRUE );
1179  fH15all_pt85->GetXaxis()->SetTitle ( "R" );
1180  fH15all_pt85->GetYaxis()->SetTitle ( "<p_{T}> track" );
1181  fH15all_pt85->GetXaxis()->SetTitleColor ( 1 );
1182  fH15all_pt85->SetMarkerStyle ( kFullCircle );
1183  fOutput->Add ( fH15all_pt85 );
1184 
1185  fH15all_pt90 = new TProfile ( "histo15all_pt90", "<p_{T}> track in dR(jet1) - 90% of Pt", h15_nbin, h15_xlow, h15_xhigh );
1186  fH15all_pt90->SetStats ( kTRUE );
1187  fH15all_pt90->GetXaxis()->SetTitle ( "R" );
1188  fH15all_pt90->GetYaxis()->SetTitle ( "<p_{T}> track" );
1189  fH15all_pt90->GetXaxis()->SetTitleColor ( 1 );
1190  fH15all_pt90->SetMarkerStyle ( kFullCircle );
1191  fOutput->Add ( fH15all_pt90 );
1192 
1193 //____________________________________________________________________________________
1194  fH15_bin = new TH1D ( "histo15_bin", "p_{T} SUM (track) in dR(jet1)", h15_nbin, h15_xlow, h15_xhigh );
1195  fH15_bin->SetStats ( kTRUE );
1196  fH15_bin->GetXaxis()->SetTitle ( " R" );
1197  fH15_bin->GetYaxis()->SetTitle ( "p_{T} SUM (track)" );
1198  fH15_bin->GetXaxis()->SetTitleColor ( 1 );
1199  fH15_bin->SetMarkerStyle ( kFullCircle );
1200  fOutput->Add ( fH15_bin );
1201 
1202 
1203 //=====================================================================================
1204  fH15_bin_n70 = new TH1D ( "histo15_bin_n70", "p_{T} SUM (track) in dR(jet1) - 70% of particles", h15_nbin, h15_xlow, h15_xhigh );
1205  fH15_bin_n70->SetStats ( kTRUE );
1206  fH15_bin_n70->GetXaxis()->SetTitle ( "R" );
1207  fH15_bin_n70->GetYaxis()->SetTitle ( "p_{T} SUM (track)" );
1208  fH15_bin_n70->GetXaxis()->SetTitleColor ( 1 );
1209  fH15_bin_n70->SetMarkerStyle ( kFullCircle );
1210  fOutput->Add ( fH15_bin_n70 );
1211 
1212  fH15_bin_n75 = new TH1D ( "histo15_bin_n75", "p_{T} SUM (track) in dR(jet1) - 75% of particles", h15_nbin, h15_xlow, h15_xhigh );
1213  fH15_bin_n75->SetStats ( kTRUE );
1214  fH15_bin_n75->GetXaxis()->SetTitle ( "R" );
1215  fH15_bin_n75->GetYaxis()->SetTitle ( "p_{T} SUM (track)" );
1216  fH15_bin_n75->GetXaxis()->SetTitleColor ( 1 );
1217  fH15_bin_n75->SetMarkerStyle ( kFullCircle );
1218  fOutput->Add ( fH15_bin_n75 );
1219 
1220  fH15_bin_n80 = new TH1D ( "histo15_bin_n80", "p_{T} SUM (track) in dR(jet1) - 80% of particles", h15_nbin, h15_xlow, h15_xhigh );
1221  fH15_bin_n80->SetStats ( kTRUE );
1222  fH15_bin_n80->GetXaxis()->SetTitle ( "R" );
1223  fH15_bin_n80->GetYaxis()->SetTitle ( "p_{T} SUM (track)" );
1224  fH15_bin_n80->GetXaxis()->SetTitleColor ( 1 );
1225  fH15_bin_n80->SetMarkerStyle ( kFullCircle );
1226  fOutput->Add ( fH15_bin_n80 );
1227 
1228  fH15_bin_n85 = new TH1D ( "histo15_bin_n85", "p_{T} SUM (track) in dR(jet1) - 85% of particles", h15_nbin, h15_xlow, h15_xhigh );
1229  fH15_bin_n85->SetStats ( kTRUE );
1230  fH15_bin_n85->GetXaxis()->SetTitle ( "R" );
1231  fH15_bin_n85->GetYaxis()->SetTitle ( "p_{T} SUM (track)" );
1232  fH15_bin_n85->GetXaxis()->SetTitleColor ( 1 );
1233  fH15_bin_n85->SetMarkerStyle ( kFullCircle );
1234  fOutput->Add ( fH15_bin_n85 );
1235 
1236  fH15_bin_n90 = new TH1D ( "histo15_bin_n90", "p_{T} SUM (track) in dR(jet1) - 90% of particles", h15_nbin, h15_xlow, h15_xhigh );
1237  fH15_bin_n90->SetStats ( kTRUE );
1238  fH15_bin_n90->GetXaxis()->SetTitle ( "R" );
1239  fH15_bin_n90->GetYaxis()->SetTitle ( "p_{T} SUM (track)" );
1240  fH15_bin_n90->GetXaxis()->SetTitleColor ( 1 );
1241  fH15_bin_n90->SetMarkerStyle ( kFullCircle );
1242  fOutput->Add ( fH15_bin_n90 );
1243 
1244 //=====================================================================================
1245  fH15_bin_pt70 = new TH1D ( "histo15_bin_pt70", "p_{T} SUM (track) in dR(jet1) - 70% of Pt", h15_nbin, h15_xlow, h15_xhigh );
1246  fH15_bin_pt70->SetStats ( kTRUE );
1247  fH15_bin_pt70->GetXaxis()->SetTitle ( "R" );
1248  fH15_bin_pt70->GetYaxis()->SetTitle ( "p_{T} SUM (track)" );
1249  fH15_bin_pt70->GetXaxis()->SetTitleColor ( 1 );
1250  fH15_bin_pt70->SetMarkerStyle ( kFullCircle );
1251  fOutput->Add ( fH15_bin_pt70 );
1252 
1253  fH15_bin_pt75 = new TH1D ( "histo15_bin_pt75", "p_{T} SUM (track) in dR(jet1) - 75% of Pt", h15_nbin, h15_xlow, h15_xhigh );
1254  fH15_bin_pt75->SetStats ( kTRUE );
1255  fH15_bin_pt75->GetXaxis()->SetTitle ( "R" );
1256  fH15_bin_pt75->GetYaxis()->SetTitle ( "p_{T} SUM (track)" );
1257  fH15_bin_pt75->GetXaxis()->SetTitleColor ( 1 );
1258  fH15_bin_pt75->SetMarkerStyle ( kFullCircle );
1259  fOutput->Add ( fH15_bin_pt75 );
1260 
1261  fH15_bin_pt80 = new TH1D ( "histo15_bin_pt80", "p_{T} SUM (track) in dR(jet1) - 80% of Pt", h15_nbin, h15_xlow, h15_xhigh );
1262  fH15_bin_pt80->SetStats ( kTRUE );
1263  fH15_bin_pt80->GetXaxis()->SetTitle ( "R" );
1264  fH15_bin_pt80->GetYaxis()->SetTitle ( "p_{T} SUM (track)" );
1265  fH15_bin_pt80->GetXaxis()->SetTitleColor ( 1 );
1266  fH15_bin_pt80->SetMarkerStyle ( kFullCircle );
1267  fOutput->Add ( fH15_bin_pt80 );
1268 
1269  fH15_bin_pt85 = new TH1D ( "histo15_bin_pt85", "p_{T} SUM (track) in dR(jet1) - 85% of Pt", h15_nbin, h15_xlow, h15_xhigh );
1270  fH15_bin_pt85->SetStats ( kTRUE );
1271  fH15_bin_pt85->GetXaxis()->SetTitle ( "R" );
1272  fH15_bin_pt85->GetYaxis()->SetTitle ( "p_{T} SUM (track)" );
1273  fH15_bin_pt85->GetXaxis()->SetTitleColor ( 1 );
1274  fH15_bin_pt85->SetMarkerStyle ( kFullCircle );
1275  fOutput->Add ( fH15_bin_pt85 );
1276 
1277  fH15_bin_pt90 = new TH1D ( "histo15_bin_pt90", "p_{T} SUM (track) in dR(jet1) - 90% of Pt", h15_nbin, h15_xlow, h15_xhigh );
1278  fH15_bin_pt90->SetStats ( kTRUE );
1279  fH15_bin_pt90->GetXaxis()->SetTitle ( "R" );
1280  fH15_bin_pt90->GetYaxis()->SetTitle ( "p_{T} SUM (track)" );
1281  fH15_bin_pt90->GetXaxis()->SetTitleColor ( 1 );
1282  fH15_bin_pt90->SetMarkerStyle ( kFullCircle );
1283  fOutput->Add ( fH15_bin_pt90 );
1284 
1285 
1286 //____________________________________________________________________________________
1287  fH15all_bin = new TH1D ( "histo15_bin_all", "p_{T} SUM (track) in dR (jet)", h15_nbin, h15_xlow, h15_xhigh );
1288  fH15all_bin->SetStats ( kTRUE );
1289  fH15all_bin->GetXaxis()->SetTitle ( "R" );
1290  fH15all_bin->GetYaxis()->SetTitle ( "p_{T} SUM (track)" );
1291  fH15all_bin->GetXaxis()->SetTitleColor ( 1 );
1292  fH15all_bin->SetMarkerStyle ( kFullCircle );
1293  fOutput->Add ( fH15all_bin );
1294 
1295 //=====================================================================================
1296  fH15all_bin_n70 = new TH1D ( "histo15_bin_n70_all", "p_{T} SUM (track) in dR (jet) - 70% of particles", h15_nbin, h15_xlow, h15_xhigh );
1297  fH15all_bin_n70->SetStats ( kTRUE );
1298  fH15all_bin_n70->GetXaxis()->SetTitle ( "R" );
1299  fH15all_bin_n70->GetYaxis()->SetTitle ( "p_{T} SUM (track)" );
1300  fH15all_bin_n70->GetXaxis()->SetTitleColor ( 1 );
1301  fH15all_bin_n70->SetMarkerStyle ( kFullCircle );
1302  fOutput->Add ( fH15all_bin_n70 );
1303 
1304  fH15all_bin_n75 = new TH1D ( "histo15_bin_n75_all", "p_{T} SUM (track) in dR (jet) - 75% of particles", h15_nbin, h15_xlow, h15_xhigh );
1305  fH15all_bin_n75->SetStats ( kTRUE );
1306  fH15all_bin_n75->GetXaxis()->SetTitle ( "R" );
1307  fH15all_bin_n75->GetYaxis()->SetTitle ( "p_{T} SUM (track)" );
1308  fH15all_bin_n75->GetXaxis()->SetTitleColor ( 1 );
1309  fH15all_bin_n75->SetMarkerStyle ( kFullCircle );
1310  fOutput->Add ( fH15all_bin_n75 );
1311 
1312  fH15all_bin_n80 = new TH1D ( "histo15_bin_n80_all", "p_{T} SUM (track) in dR (jet) - 80% of particles", h15_nbin, h15_xlow, h15_xhigh );
1313  fH15all_bin_n80->SetStats ( kTRUE );
1314  fH15all_bin_n80->GetXaxis()->SetTitle ( "R" );
1315  fH15all_bin_n80->GetYaxis()->SetTitle ( "p_{T} SUM (track)" );
1316  fH15all_bin_n80->GetXaxis()->SetTitleColor ( 1 );
1317  fH15all_bin_n80->SetMarkerStyle ( kFullCircle );
1318  fOutput->Add ( fH15all_bin_n80 );
1319 
1320  fH15all_bin_n85 = new TH1D ( "histo15_bin_n85_all", "p_{T} SUM (track) in dR (jet) - 85% of particles", h15_nbin, h15_xlow, h15_xhigh );
1321  fH15all_bin_n85->SetStats ( kTRUE );
1322  fH15all_bin_n85->GetXaxis()->SetTitle ( "R" );
1323  fH15all_bin_n85->GetYaxis()->SetTitle ( "p_{T} SUM (track)" );
1324  fH15all_bin_n85->GetXaxis()->SetTitleColor ( 1 );
1325  fH15all_bin_n85->SetMarkerStyle ( kFullCircle );
1326  fOutput->Add ( fH15all_bin_n85 );
1327 
1328  fH15all_bin_n90 = new TH1D ( "histo15_bin_n90_all", "p_{T} SUM (track) in dR (jet) - 90% of particles", h15_nbin, h15_xlow, h15_xhigh );
1329  fH15all_bin_n90->SetStats ( kTRUE );
1330  fH15all_bin_n90->GetXaxis()->SetTitle ( "R" );
1331  fH15all_bin_n90->GetYaxis()->SetTitle ( "p_{T} SUM (track)" );
1332  fH15all_bin_n90->GetXaxis()->SetTitleColor ( 1 );
1333  fH15all_bin_n90->SetMarkerStyle ( kFullCircle );
1334  fOutput->Add ( fH15all_bin_n90 );
1335 
1336 //=====================================================================================
1337  fH15all_bin_pt70 = new TH1D ( "histo15_bin_pt70_all", "p_{T} SUM (track) in dR (jet) - 70% of Pt", h15_nbin, h15_xlow, h15_xhigh );
1338  fH15all_bin_pt70->SetStats ( kTRUE );
1339  fH15all_bin_pt70->GetXaxis()->SetTitle ( "R" );
1340  fH15all_bin_pt70->GetYaxis()->SetTitle ( "p_{T} SUM (track)" );
1341  fH15all_bin_pt70->GetXaxis()->SetTitleColor ( 1 );
1342  fH15all_bin_pt70->SetMarkerStyle ( kFullCircle );
1343  fOutput->Add ( fH15all_bin_pt70 );
1344 
1345  fH15all_bin_pt75 = new TH1D ( "histo15_bin_pt75_all", "p_{T} SUM (track) in dR (jet) - 75% of Pt", h15_nbin, h15_xlow, h15_xhigh );
1346  fH15all_bin_pt75->SetStats ( kTRUE );
1347  fH15all_bin_pt75->GetXaxis()->SetTitle ( "R" );
1348  fH15all_bin_pt75->GetYaxis()->SetTitle ( "p_{T} SUM (track)" );
1349  fH15all_bin_pt75->GetXaxis()->SetTitleColor ( 1 );
1350  fH15all_bin_pt75->SetMarkerStyle ( kFullCircle );
1351  fOutput->Add ( fH15all_bin_pt75 );
1352 
1353  fH15all_bin_pt80 = new TH1D ( "histo15_bin_pt80_all", "p_{T} SUM (track) in dR (jet) - 80% of Pt", h15_nbin, h15_xlow, h15_xhigh );
1354  fH15all_bin_pt80->SetStats ( kTRUE );
1355  fH15all_bin_pt80->GetXaxis()->SetTitle ( "R" );
1356  fH15all_bin_pt80->GetYaxis()->SetTitle ( "p_{T} SUM (track)" );
1357  fH15all_bin_pt80->GetXaxis()->SetTitleColor ( 1 );
1358  fH15all_bin_pt80->SetMarkerStyle ( kFullCircle );
1359  fOutput->Add ( fH15all_bin_pt80 );
1360 
1361  fH15all_bin_pt85 = new TH1D ( "histo15_bin_pt85_all", "p_{T} SUM (track) in dR (jet) - 85% of Pt", h15_nbin, h15_xlow, h15_xhigh );
1362  fH15all_bin_pt85->SetStats ( kTRUE );
1363  fH15all_bin_pt85->GetXaxis()->SetTitle ( "R" );
1364  fH15all_bin_pt85->GetYaxis()->SetTitle ( "p_{T} SUM (track)" );
1365  fH15all_bin_pt85->GetXaxis()->SetTitleColor ( 1 );
1366  fH15all_bin_pt85->SetMarkerStyle ( kFullCircle );
1367  fOutput->Add ( fH15all_bin_pt85 );
1368 
1369  fH15all_bin_pt90 = new TH1D ( "histo15_bin_pt90_all", "p_{T} SUM (track) in dR (jet) - 90% of Pt", h15_nbin, h15_xlow, h15_xhigh );
1370  fH15all_bin_pt90->SetStats ( kTRUE );
1371  fH15all_bin_pt90->GetXaxis()->SetTitle ( "R" );
1372  fH15all_bin_pt90->GetYaxis()->SetTitle ( "p_{T} SUM (track)" );
1373  fH15all_bin_pt90->GetXaxis()->SetTitleColor ( 1 );
1374  fH15all_bin_pt90->SetMarkerStyle ( kFullCircle );
1375  fOutput->Add ( fH15all_bin_pt90 );
1376 
1377 //____________________________________________________________________________________
1378  Int_t h20_nbin = 100;
1379  Double_t h20_binwidth = 0.01;
1380  Double_t h20_low = 0.;
1381  Double_t h20_high = h20_low + h20_binwidth * h20_nbin;
1382  fH20 = new TH1D ( "histo20", "dN/dR (jet1)", h20_nbin, h20_low, h20_high );
1383  fH20->SetStats ( kTRUE );
1384  fH20->GetXaxis()->SetTitle ( "R" );
1385  fH20->GetYaxis()->SetTitle ( "dN/dR" );
1386  fH20->GetXaxis()->SetTitleColor ( 1 );
1387  fH20->SetMarkerStyle ( kFullCircle );
1388  fOutput->Add ( fH20 );
1389 
1390 //=====================================================================================
1391  fH20_n70 = new TH1D ( "histo20_n70", "dN/dR (jet1) - 70% of particles", h20_nbin, h20_low, h20_high );
1392  fH20_n70->SetStats ( kTRUE );
1393  fH20_n70->GetXaxis()->SetTitle ( "R" );
1394  fH20_n70->GetYaxis()->SetTitle ( "dN/dR" );
1395  fH20_n70->GetXaxis()->SetTitleColor ( 1 );
1396  fH20_n70->SetMarkerStyle ( kFullCircle );
1397  fOutput->Add ( fH20_n70 );
1398 
1399  fH20_n75 = new TH1D ( "histo20_n75", "dN/dR (jet1) - 75% of particles", h20_nbin, h20_low, h20_high );
1400  fH20_n75->SetStats ( kTRUE );
1401  fH20_n75->GetXaxis()->SetTitle ( "R" );
1402  fH20_n75->GetYaxis()->SetTitle ( "dN/dR" );
1403  fH20_n75->GetXaxis()->SetTitleColor ( 1 );
1404  fH20_n75->SetMarkerStyle ( kFullCircle );
1405  fOutput->Add ( fH20_n75 );
1406 
1407  fH20_n80 = new TH1D ( "histo20_n80", "dN/dR (jet1) - 80% of particles", h20_nbin, h20_low, h20_high );
1408  fH20_n80->SetStats ( kTRUE );
1409  fH20_n80->GetXaxis()->SetTitle ( "R" );
1410  fH20_n80->GetYaxis()->SetTitle ( "dN/dR" );
1411  fH20_n80->GetXaxis()->SetTitleColor ( 1 );
1412  fH20_n80->SetMarkerStyle ( kFullCircle );
1413  fOutput->Add ( fH20_n80 );
1414 
1415  fH20_n85 = new TH1D ( "histo20_n85", "dN/dR (jet1) - 85% of particles", h20_nbin, h20_low, h20_high );
1416  fH20_n85->SetStats ( kTRUE );
1417  fH20_n85->GetXaxis()->SetTitle ( "R" );
1418  fH20_n85->GetYaxis()->SetTitle ( "dN/dR" );
1419  fH20_n85->GetXaxis()->SetTitleColor ( 1 );
1420  fH20_n85->SetMarkerStyle ( kFullCircle );
1421  fOutput->Add ( fH20_n85 );
1422 
1423  fH20_n90 = new TH1D ( "histo20_n90", "dN/dR (jet1) - 90% of particles", h20_nbin, h20_low, h20_high );
1424  fH20_n90->SetStats ( kTRUE );
1425  fH20_n90->GetXaxis()->SetTitle ( "R" );
1426  fH20_n90->GetYaxis()->SetTitle ( "dN/dR" );
1427  fH20_n90->GetXaxis()->SetTitleColor ( 1 );
1428  fH20_n90->SetMarkerStyle ( kFullCircle );
1429  fOutput->Add ( fH20_n90 );
1430 
1431 //=====================================================================================
1432  fH20_pt70 = new TH1D ( "histo20_pt70", "dN/dR (jet1) - 70% of Pt", h20_nbin, h20_low, h20_high );
1433  fH20_pt70->SetStats ( kTRUE );
1434  fH20_pt70->GetXaxis()->SetTitle ( "R" );
1435  fH20_pt70->GetYaxis()->SetTitle ( "dN/dR" );
1436  fH20_pt70->GetXaxis()->SetTitleColor ( 1 );
1437  fH20_pt70->SetMarkerStyle ( kFullCircle );
1438  fOutput->Add ( fH20_pt70 );
1439 
1440  fH20_pt75 = new TH1D ( "histo20_pt75", "dN/dR (jet1) - 75% of Pt", h20_nbin, h20_low, h20_high );
1441  fH20_pt75->SetStats ( kTRUE );
1442  fH20_pt75->GetXaxis()->SetTitle ( "R" );
1443  fH20_pt75->GetYaxis()->SetTitle ( "dN/dR" );
1444  fH20_pt75->GetXaxis()->SetTitleColor ( 1 );
1445  fH20_pt75->SetMarkerStyle ( kFullCircle );
1446  fOutput->Add ( fH20_pt75 );
1447 
1448  fH20_pt80 = new TH1D ( "histo20_pt80", "dN/dR (jet1) - 80% of Pt", h20_nbin, h20_low, h20_high );
1449  fH20_pt80->SetStats ( kTRUE );
1450  fH20_pt80->GetXaxis()->SetTitle ( "R" );
1451  fH20_pt80->GetYaxis()->SetTitle ( "dN/dR" );
1452  fH20_pt80->GetXaxis()->SetTitleColor ( 1 );
1453  fH20_pt80->SetMarkerStyle ( kFullCircle );
1454  fOutput->Add ( fH20_pt80 );
1455 
1456  fH20_pt85 = new TH1D ( "histo20_pt85", "dN/dR (jet1) - 85% of Pt", h20_nbin, h20_low, h20_high );
1457  fH20_pt85->SetStats ( kTRUE );
1458  fH20_pt85->GetXaxis()->SetTitle ( "R" );
1459  fH20_pt85->GetYaxis()->SetTitle ( "dN/dR" );
1460  fH20_pt85->GetXaxis()->SetTitleColor ( 1 );
1461  fH20_pt85->SetMarkerStyle ( kFullCircle );
1462  fOutput->Add ( fH20_pt85 );
1463 
1464  fH20_pt90 = new TH1D ( "histo20_pt90", "dN/dR (jet1) - 90% of Pt", h20_nbin, h20_low, h20_high );
1465  fH20_pt90->SetStats ( kTRUE );
1466  fH20_pt90->GetXaxis()->SetTitle ( "R" );
1467  fH20_pt90->GetYaxis()->SetTitle ( "dN/dR" );
1468  fH20_pt90->GetXaxis()->SetTitleColor ( 1 );
1469  fH20_pt90->SetMarkerStyle ( kFullCircle );
1470  fOutput->Add ( fH20_pt90 );
1471 
1472 
1473 //____________________________________________________________________________________
1474  fH20all = new TH1D ( "histo20_all", "dN/dR - all jets", h20_nbin, h20_low, h20_high );
1475  fH20all->SetStats ( kTRUE );
1476  fH20all->GetXaxis()->SetTitle ( "R" );
1477  fH20all->GetYaxis()->SetTitle ( "dN/dR" );
1478  fH20all->GetXaxis()->SetTitleColor ( 1 );
1479  fH20all->SetMarkerStyle ( kFullCircle );
1480  fOutput->Add ( fH20all );
1481 
1482 
1483 //=====================================================================================
1484  fH20all_n70 = new TH1D ( "histo20_n70_all", "dN/dR - all jets - 70% of particles", h20_nbin, h20_low, h20_high );
1485  fH20all_n70->SetStats ( kTRUE );
1486  fH20all_n70->GetXaxis()->SetTitle ( "R" );
1487  fH20all_n70->GetYaxis()->SetTitle ( "dN/dR" );
1488  fH20all_n70->GetXaxis()->SetTitleColor ( 1 );
1489  fH20all_n70->SetMarkerStyle ( kFullCircle );
1490  fOutput->Add ( fH20all_n70 );
1491 
1492  fH20all_n75 = new TH1D ( "histo20_n75_all", "dN/dR - all jets - 75% of particles", h20_nbin, h20_low, h20_high );
1493  fH20all_n75->SetStats ( kTRUE );
1494  fH20all_n75->GetXaxis()->SetTitle ( "R" );
1495  fH20all_n75->GetYaxis()->SetTitle ( "dN/dR" );
1496  fH20all_n75->GetXaxis()->SetTitleColor ( 1 );
1497  fH20all_n75->SetMarkerStyle ( kFullCircle );
1498  fOutput->Add ( fH20all_n75 );
1499 
1500  fH20all_n80 = new TH1D ( "histo20_n80_all", "dN/dR - all jets - 80% of particles", h20_nbin, h20_low, h20_high );
1501  fH20all_n80->SetStats ( kTRUE );
1502  fH20all_n80->GetXaxis()->SetTitle ( "R" );
1503  fH20all_n80->GetYaxis()->SetTitle ( "dN/dR" );
1504  fH20all_n80->GetXaxis()->SetTitleColor ( 1 );
1505  fH20all_n80->SetMarkerStyle ( kFullCircle );
1506  fOutput->Add ( fH20all_n80 );
1507 
1508  fH20all_n85 = new TH1D ( "histo20_n85_all", "dN/dR - all jets - 85% of particles", h20_nbin, h20_low, h20_high );
1509  fH20all_n85->SetStats ( kTRUE );
1510  fH20all_n85->GetXaxis()->SetTitle ( "R" );
1511  fH20all_n85->GetYaxis()->SetTitle ( "dN/dR" );
1512  fH20all_n85->GetXaxis()->SetTitleColor ( 1 );
1513  fH20all_n85->SetMarkerStyle ( kFullCircle );
1514  fOutput->Add ( fH20all_n85 );
1515 
1516  fH20all_n90 = new TH1D ( "histo20_n90_all", "dN/dR - all jets - 90% of particles", h20_nbin, h20_low, h20_high );
1517  fH20all_n90->SetStats ( kTRUE );
1518  fH20all_n90->GetXaxis()->SetTitle ( "R" );
1519  fH20all_n90->GetYaxis()->SetTitle ( "dN/dR" );
1520  fH20all_n90->GetXaxis()->SetTitleColor ( 1 );
1521  fH20all_n90->SetMarkerStyle ( kFullCircle );
1522  fOutput->Add ( fH20all_n90 );
1523 
1524 //=====================================================================================
1525  fH20all_pt70 = new TH1D ( "histo20_pt70_all", "dN/dR - all jets - 70% of Pt", h20_nbin, h20_low, h20_high );
1526  fH20all_pt70->SetStats ( kTRUE );
1527  fH20all_pt70->GetXaxis()->SetTitle ( "R" );
1528  fH20all_pt70->GetYaxis()->SetTitle ( "dN/dR" );
1529  fH20all_pt70->GetXaxis()->SetTitleColor ( 1 );
1530  fH20all_pt70->SetMarkerStyle ( kFullCircle );
1531  fOutput->Add ( fH20all_pt70 );
1532 
1533  fH20all_pt75 = new TH1D ( "histo20_pt75_all", "dN/dR - all jets - 75% of Pt", h20_nbin, h20_low, h20_high );
1534  fH20all_pt75->SetStats ( kTRUE );
1535  fH20all_pt75->GetXaxis()->SetTitle ( "R" );
1536  fH20all_pt75->GetYaxis()->SetTitle ( "dN/dR" );
1537  fH20all_pt75->GetXaxis()->SetTitleColor ( 1 );
1538  fH20all_pt75->SetMarkerStyle ( kFullCircle );
1539  fOutput->Add ( fH20all_pt75 );
1540 
1541  fH20all_pt80 = new TH1D ( "histo20_pt80_all", "dN/dR - all jets - 80% of Pt", h20_nbin, h20_low, h20_high );
1542  fH20all_pt80->SetStats ( kTRUE );
1543  fH20all_pt80->GetXaxis()->SetTitle ( "R" );
1544  fH20all_pt80->GetYaxis()->SetTitle ( "dN/dR" );
1545  fH20all_pt80->GetXaxis()->SetTitleColor ( 1 );
1546  fH20all_pt80->SetMarkerStyle ( kFullCircle );
1547  fOutput->Add ( fH20all_pt80 );
1548 
1549  fH20all_pt85 = new TH1D ( "histo20_pt85_all", "dN/dR - all jets - 85% of Pt", h20_nbin, h20_low, h20_high );
1550  fH20all_pt85->SetStats ( kTRUE );
1551  fH20all_pt85->GetXaxis()->SetTitle ( "R" );
1552  fH20all_pt85->GetYaxis()->SetTitle ( "dN/dR" );
1553  fH20all_pt85->GetXaxis()->SetTitleColor ( 1 );
1554  fH20all_pt85->SetMarkerStyle ( kFullCircle );
1555  fOutput->Add ( fH20all_pt85 );
1556 
1557  fH20all_pt90 = new TH1D ( "histo20_pt90_all", "dN/dR - all jets - 90% of Pt", h20_nbin, h20_low, h20_high );
1558  fH20all_pt90->SetStats ( kTRUE );
1559  fH20all_pt90->GetXaxis()->SetTitle ( "R" );
1560  fH20all_pt90->GetYaxis()->SetTitle ( "dN/dR" );
1561  fH20all_pt90->GetXaxis()->SetTitleColor ( 1 );
1562  fH20all_pt90->SetMarkerStyle ( kFullCircle );
1563  fOutput->Add ( fH20all_pt90 );
1564 
1565  //____________________________________________________________________________________
1566  Int_t h23_nbin = 400;
1567  Double_t h23_binwidth = 1.;
1568  Double_t h23_low = 0.;
1569  Double_t h23_high = h23_low + h23_binwidth * h23_nbin;
1570 
1571  fH23 = new TProfile ( "histo23", "Jet1 Size vs P_{T}(jet1)", h23_nbin, h23_low, h23_high );
1572  fH23->SetStats ( kTRUE );
1573  fH23->GetXaxis()->SetTitle ( "p_{T}(jet1) (GeV/c)" );
1574  fH23->GetYaxis()->SetTitle ( "<R(jet1)> in 1 GeV/c bin" );
1575  fH23->GetXaxis()->SetTitleColor ( 1 );
1576  fH23->SetMarkerStyle ( kFullCircle );
1577  fOutput->Add ( fH23 );
1578 
1579  fH24 = new TProfile ( "histo24", "Jet1 Size vs P_{T}(jet1) - 80% of particles", h23_nbin, h23_low, h23_high );
1580  fH24->SetStats ( kTRUE );
1581  fH24->GetXaxis()->SetTitle ( "p_{T}(jet1) (GeV/c)" );
1582  fH24->GetYaxis()->SetTitle ( "<R(jet1)> in 1 GeV/c bin" );
1583  fH24->GetXaxis()->SetTitleColor ( 1 );
1584  fH24->SetMarkerStyle ( kFullCircle );
1585  fOutput->Add ( fH24 );
1586 
1587  fH25 = new TProfile ( "histo25", "Jet1 Size vs P_{T}(jet1) - 80% of Pt", h23_nbin, h23_low, h23_high );
1588  fH25->SetStats ( kTRUE );
1589  fH25->GetXaxis()->SetTitle ( "p_{T}(jet1) (GeV/c)" );
1590  fH25->GetYaxis()->SetTitle ( "<R(jet1)> in 1 GeV/c bin" );
1591  fH25->GetXaxis()->SetTitleColor ( 1 );
1592  fH25->SetMarkerStyle ( kFullCircle );
1593  fOutput->Add ( fH25 );
1594 
1595  fH23all = new TProfile ( "histo23all", "Jet1 Size vs P_{T} - all jets", h23_nbin, h23_low, h23_high );
1596  fH23all->SetStats ( kTRUE );
1597  fH23all->GetXaxis()->SetTitle ( "p_{T}(jet) (GeV/c)" );
1598  fH23all->GetYaxis()->SetTitle ( "<R(jet1)> in 1 GeV/c bin" );
1599  fH23all->GetXaxis()->SetTitleColor ( 1 );
1600  fH23all->SetMarkerStyle ( kFullCircle );
1601  fOutput->Add ( fH23all );
1602 
1603  fH24all = new TProfile ( "histo24all", "Jet1 Size vs P_{T}(jet1) - 80% of particles", h23_nbin, h23_low, h23_high );
1604  fH24all->SetStats ( kTRUE );
1605  fH24all->GetXaxis()->SetTitle ( "p_{T}(jet) (GeV/c)" );
1606  fH24all->GetYaxis()->SetTitle ( "<R(jet1)> in 1 GeV/c bin" );
1607  fH24all->GetXaxis()->SetTitleColor ( 1 );
1608  fH24all->SetMarkerStyle ( kFullCircle );
1609  fOutput->Add ( fH24all );
1610 
1611  fH25all = new TProfile ( "histo25all", "Jet1 Size vs P_{T}(jet1) - 80% of Pt", h23_nbin, h23_low, h23_high );
1612  fH25all->SetStats ( kTRUE );
1613  fH25all->GetXaxis()->SetTitle ( "p_{T}(jet) (GeV/c)" );
1614  fH25all->GetYaxis()->SetTitle ( "<R(jet1)> in 1 GeV/c bin" );
1615  fH25all->GetXaxis()->SetTitleColor ( 1 );
1616  fH25all->SetMarkerStyle ( kFullCircle );
1617  fOutput->Add ( fH25all );
1618 
1619 //=============================================================
1620  Int_t hg_nbin = 100;
1621  Double_t hg_binwidth = 0.01;
1622  Double_t hg_low = 0.;
1623  Double_t hg_high = hg_low + hg_binwidth * hg_nbin;
1624 
1625  fHg = new TH1D ( "histo_g", "dN/dg", hg_nbin, hg_low, hg_high );
1626  fHg->SetStats ( kTRUE );
1627  fHg->GetXaxis()->SetTitle ( "g" );
1628  fHg->GetYaxis()->SetTitle ( "1/N_{jets} dN/dg" );
1629  fHg->GetXaxis()->SetTitleColor ( 1 );
1630  fHg->SetMarkerStyle ( kFullCircle );
1631  fOutput->Add ( fHg );
1632 
1633  fHg_n90 = new TH1D ( "histo_g_n90", "dN/dg", hg_nbin, hg_low, hg_high );
1634  fHg_n90->SetStats ( kTRUE );
1635  fHg_n90->GetXaxis()->SetTitle ( "g" );
1636  fHg_n90->GetYaxis()->SetTitle ( "1/N_{jets} dN/dg" );
1637  fHg_n90->GetXaxis()->SetTitleColor ( 1 );
1638  fHg_n90->SetMarkerStyle ( kFullCircle );
1639  fOutput->Add ( fHg_n90 );
1640 
1641  fHg_n85 = new TH1D ( "histo_g_n85", "dN/dg", hg_nbin, hg_low, hg_high );
1642  fHg_n85->SetStats ( kTRUE );
1643  fHg_n85->GetXaxis()->SetTitle ( "g" );
1644  fHg_n85->GetYaxis()->SetTitle ( "1/N_{jets} dN/dg" );
1645  fHg_n85->GetXaxis()->SetTitleColor ( 1 );
1646  fHg_n85->SetMarkerStyle ( kFullCircle );
1647  fOutput->Add ( fHg_n85 );
1648 
1649  fHg_n80 = new TH1D ( "histo_g_n80", "dN/dg", hg_nbin, hg_low, hg_high );
1650  fHg_n80->SetStats ( kTRUE );
1651  fHg_n80->GetXaxis()->SetTitle ( "g" );
1652  fHg_n80->GetYaxis()->SetTitle ( "1/N_{jets} dN/dg" );
1653  fHg_n80->GetXaxis()->SetTitleColor ( 1 );
1654  fHg_n80->SetMarkerStyle ( kFullCircle );
1655  fOutput->Add ( fHg_n80 );
1656 
1657  fHg_n75 = new TH1D ( "histo_g_n75", "dN/dg", hg_nbin, hg_low, hg_high );
1658  fHg_n75->SetStats ( kTRUE );
1659  fHg_n75->GetXaxis()->SetTitle ( "g" );
1660  fHg_n75->GetYaxis()->SetTitle ( "1/N_{jets} dN/dg" );
1661  fHg_n75->GetXaxis()->SetTitleColor ( 1 );
1662  fHg_n75->SetMarkerStyle ( kFullCircle );
1663  fOutput->Add ( fHg_n75 );
1664 
1665  fHg_n70 = new TH1D ( "histo_g_n70", "dN/dg", hg_nbin, hg_low, hg_high );
1666  fHg_n70->SetStats ( kTRUE );
1667  fHg_n70->GetXaxis()->SetTitle ( "g" );
1668  fHg_n70->GetYaxis()->SetTitle ( "1/N_{jets} dN/dg" );
1669  fHg_n70->GetXaxis()->SetTitleColor ( 1 );
1670  fHg_n70->SetMarkerStyle ( kFullCircle );
1671  fOutput->Add ( fHg_n70 );
1672 
1673  fHg_pt90 = new TH1D ( "histo_g_pt90", "dN/dg", hg_nbin, hg_low, hg_high );
1674  fHg_pt90->SetStats ( kTRUE );
1675  fHg_pt90->GetXaxis()->SetTitle ( "g" );
1676  fHg_pt90->GetYaxis()->SetTitle ( "1/N_{jets} dN/dg" );
1677  fHg_pt90->GetXaxis()->SetTitleColor ( 1 );
1678  fHg_pt90->SetMarkerStyle ( kFullCircle );
1679  fOutput->Add ( fHg_pt90 );
1680 
1681  fHg_pt85 = new TH1D ( "histo_g_pt85", "dN/dg", hg_nbin, hg_low, hg_high );
1682  fHg_pt85->SetStats ( kTRUE );
1683  fHg_pt85->GetXaxis()->SetTitle ( "g" );
1684  fHg_pt85->GetYaxis()->SetTitle ( "1/N_{jets} dN/dg" );
1685  fHg_pt85->GetXaxis()->SetTitleColor ( 1 );
1686  fHg_pt85->SetMarkerStyle ( kFullCircle );
1687  fOutput->Add ( fHg_pt85 );
1688 
1689  fHg_pt80 = new TH1D ( "histo_g_pt80", "dN/dg", hg_nbin, hg_low, hg_high );
1690  fHg_pt80->SetStats ( kTRUE );
1691  fHg_pt80->GetXaxis()->SetTitle ( "g" );
1692  fHg_pt80->GetYaxis()->SetTitle ( "1/N_{jets} dN/dg" );
1693  fHg_pt80->GetXaxis()->SetTitleColor ( 1 );
1694  fHg_pt80->SetMarkerStyle ( kFullCircle );
1695  fOutput->Add ( fHg_pt80 );
1696 
1697  fHg_pt75 = new TH1D ( "histo_g_pt75", "dN/dg", hg_nbin, hg_low, hg_high );
1698  fHg_pt75->SetStats ( kTRUE );
1699  fHg_pt75->GetXaxis()->SetTitle ( "g" );
1700  fHg_pt75->GetYaxis()->SetTitle ( "1/N_{jets} dN/dg" );
1701  fHg_pt75->GetXaxis()->SetTitleColor ( 1 );
1702  fHg_pt75->SetMarkerStyle ( kFullCircle );
1703  fOutput->Add ( fHg_pt75 );
1704 
1705  fHg_pt70 = new TH1D ( "histo_g_pt70", "dN/dg", hg_nbin, hg_low, hg_high );
1706  fHg_pt70->SetStats ( kTRUE );
1707  fHg_pt70->GetXaxis()->SetTitle ( "g" );
1708  fHg_pt70->GetYaxis()->SetTitle ( "1/N_{jets} dN/dg" );
1709  fHg_pt70->GetXaxis()->SetTitleColor ( 1 );
1710  fHg_pt70->SetMarkerStyle ( kFullCircle );
1711  fOutput->Add ( fHg_pt70 );
1712 
1713 
1714 //=============================================================
1715  Int_t hptd_nbin = 100;
1716  Double_t hptd_binwidth = 0.01;
1717  Double_t hptd_low = 0.;
1718  Double_t hptd_high = hptd_low + hptd_binwidth * hptd_nbin;
1719 
1720  fHptd = new TH1D ( "histo_g", "dN/dg", hptd_nbin, hptd_low, hptd_high );
1721  fHptd->SetStats ( kTRUE );
1722  fHptd->GetXaxis()->SetTitle ( "g" );
1723  fHptd->GetYaxis()->SetTitle ( "1/N_{jets} dN/dg" );
1724  fHptd->GetXaxis()->SetTitleColor ( 1 );
1725  fHptd->SetMarkerStyle ( kFullCircle );
1726  fOutput->Add ( fHptd );
1727 
1728  fHptd_n90 = new TH1D ( "histo_ptd_n90", "dN/dg", hptd_nbin, hptd_low, hptd_high );
1729  fHptd_n90->SetStats ( kTRUE );
1730  fHptd_n90->GetXaxis()->SetTitle ( "g" );
1731  fHptd_n90->GetYaxis()->SetTitle ( "1/N_{jets} dN/dg" );
1732  fHptd_n90->GetXaxis()->SetTitleColor ( 1 );
1733  fHptd_n90->SetMarkerStyle ( kFullCircle );
1734  fOutput->Add ( fHptd_n90 );
1735 
1736  fHptd_n85 = new TH1D ( "histo_ptd_n85", "dN/dg", hptd_nbin, hptd_low, hptd_high );
1737  fHptd_n85->SetStats ( kTRUE );
1738  fHptd_n85->GetXaxis()->SetTitle ( "g" );
1739  fHptd_n85->GetYaxis()->SetTitle ( "1/N_{jets} dN/dg" );
1740  fHptd_n85->GetXaxis()->SetTitleColor ( 1 );
1741  fHptd_n85->SetMarkerStyle ( kFullCircle );
1742  fOutput->Add ( fHptd_n85 );
1743 
1744  fHptd_n80 = new TH1D ( "histo_ptd_n80", "dN/dg", hptd_nbin, hptd_low, hptd_high );
1745  fHptd_n80->SetStats ( kTRUE );
1746  fHptd_n80->GetXaxis()->SetTitle ( "g" );
1747  fHptd_n80->GetYaxis()->SetTitle ( "1/N_{jets} dN/dg" );
1748  fHptd_n80->GetXaxis()->SetTitleColor ( 1 );
1749  fHptd_n80->SetMarkerStyle ( kFullCircle );
1750  fOutput->Add ( fHptd_n80 );
1751 
1752  fHptd_n75 = new TH1D ( "histo_ptd_n75", "dN/dg", hptd_nbin, hptd_low, hptd_high );
1753  fHptd_n75->SetStats ( kTRUE );
1754  fHptd_n75->GetXaxis()->SetTitle ( "g" );
1755  fHptd_n75->GetYaxis()->SetTitle ( "1/N_{jets} dN/dg" );
1756  fHptd_n75->GetXaxis()->SetTitleColor ( 1 );
1757  fHptd_n75->SetMarkerStyle ( kFullCircle );
1758  fOutput->Add ( fHptd_n75 );
1759 
1760  fHptd_n70 = new TH1D ( "histo_ptd_n70", "dN/dg", hptd_nbin, hptd_low, hptd_high );
1761  fHptd_n70->SetStats ( kTRUE );
1762  fHptd_n70->GetXaxis()->SetTitle ( "g" );
1763  fHptd_n70->GetYaxis()->SetTitle ( "1/N_{jets} dN/dg" );
1764  fHptd_n70->GetXaxis()->SetTitleColor ( 1 );
1765  fHptd_n70->SetMarkerStyle ( kFullCircle );
1766  fOutput->Add ( fHptd_n70 );
1767 
1768  fHptd_pt90 = new TH1D ( "histo_ptd_pt90", "dN/dg", hptd_nbin, hptd_low, hptd_high );
1769  fHptd_pt90->SetStats ( kTRUE );
1770  fHptd_pt90->GetXaxis()->SetTitle ( "g" );
1771  fHptd_pt90->GetYaxis()->SetTitle ( "1/N_{jets} dN/dg" );
1772  fHptd_pt90->GetXaxis()->SetTitleColor ( 1 );
1773  fHptd_pt90->SetMarkerStyle ( kFullCircle );
1774  fOutput->Add ( fHptd_pt90 );
1775 
1776  fHptd_pt85 = new TH1D ( "histo_ptd_pt85", "dN/dg", hptd_nbin, hptd_low, hptd_high );
1777  fHptd_pt85->SetStats ( kTRUE );
1778  fHptd_pt85->GetXaxis()->SetTitle ( "g" );
1779  fHptd_pt85->GetYaxis()->SetTitle ( "1/N_{jets} dN/dg" );
1780  fHptd_pt85->GetXaxis()->SetTitleColor ( 1 );
1781  fHptd_pt85->SetMarkerStyle ( kFullCircle );
1782  fOutput->Add ( fHptd_pt85 );
1783 
1784  fHptd_pt80 = new TH1D ( "histo_ptd_pt80", "dN/dg", hptd_nbin, hptd_low, hptd_high );
1785  fHptd_pt80->SetStats ( kTRUE );
1786  fHptd_pt80->GetXaxis()->SetTitle ( "g" );
1787  fHptd_pt80->GetYaxis()->SetTitle ( "1/N_{jets} dN/dg" );
1788  fHptd_pt80->GetXaxis()->SetTitleColor ( 1 );
1789  fHptd_pt80->SetMarkerStyle ( kFullCircle );
1790  fOutput->Add ( fHptd_pt80 );
1791 
1792  fHptd_pt75 = new TH1D ( "histo_ptd_pt75", "dN/dg", hptd_nbin, hptd_low, hptd_high );
1793  fHptd_pt75->SetStats ( kTRUE );
1794  fHptd_pt75->GetXaxis()->SetTitle ( "g" );
1795  fHptd_pt75->GetYaxis()->SetTitle ( "1/N_{jets} dN/dg" );
1796  fHptd_pt75->GetXaxis()->SetTitleColor ( 1 );
1797  fHptd_pt75->SetMarkerStyle ( kFullCircle );
1798  fOutput->Add ( fHptd_pt75 );
1799 
1800  fHptd_pt70 = new TH1D ( "histo_ptd_pt70", "dN/dg", hptd_nbin, hptd_low, hptd_high );
1801  fHptd_pt70->SetStats ( kTRUE );
1802  fHptd_pt70->GetXaxis()->SetTitle ( "g" );
1803  fHptd_pt70->GetYaxis()->SetTitle ( "1/N_{jets} dN/dg" );
1804  fHptd_pt70->GetXaxis()->SetTitleColor ( 1 );
1805  fHptd_pt70->SetMarkerStyle ( kFullCircle );
1806  fOutput->Add ( fHptd_pt70 );
1807 
1808 
1809  // =========== Switch on Sumw2 for all histos ===========
1810  for ( Int_t i = 0; i < fOutput->GetEntries(); ++i )
1811  {
1812  TH1 *h1 = dynamic_cast<TH1 *> ( fOutput->At ( i ) );
1813  if ( h1 ) { h1->Sumw2();continue; }
1814 
1815  TProfile *hprof1 = dynamic_cast<TProfile *> ( fOutput->At ( i ) );
1816  if ( hprof1 ) { hprof1->Sumw2(); }
1817  }
1818 
1819  PostData ( 1, fOutput ); // Post data for ALL output slots > 0 here.
1820  }
1821 
1822 //________________________________________________________________________
1823 Double_t AliAnalysisTaskEmcalJetCDF::DeltaR ( const AliVParticle *part1, const AliVParticle *part2 )
1824  {
1825  // Helper function to calculate the distance between two jets or a jet and
1826  // particle
1827  Double_t dPhi = part1->Phi() - part2->Phi();
1828  Double_t dEta = part1->Eta() - part2->Eta();
1829  dPhi = TVector2::Phi_mpi_pi ( dPhi );
1830 
1831  return TMath::Sqrt ( dPhi * dPhi + dEta * dEta );
1832  }
1833 
1834 //__________________________________________________________________________________________________
1835 std::vector<Int_t> AliAnalysisTaskEmcalJetCDF::SortTracksPt ( AliVEvent *event ) const
1836  {
1837  //___________________________________________
1838  // Sorting by p_T (decreasing) event tracks
1839  //^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
1840  Int_t entries = event->GetNumberOfTracks();
1841 
1842  // Create vector for Pt sorting
1843  std::vector<ptidx_pair> pair_list;
1844  pair_list.reserve ( entries );
1845 
1846  for ( Int_t i_entry = 0; i_entry < entries; i_entry++ )
1847  {
1848  AliVParticle *track = event->GetTrack ( i_entry );
1849 
1850  if ( !track )
1851  {
1852  AliError ( Form ( "Unable to find track %d in collection %s", i_entry, event->GetName() ) );
1853  continue;
1854  }
1855 
1856  pair_list.push_back ( std::make_pair ( track->Pt(), i_entry ) );
1857  }
1858 
1859  std::stable_sort ( pair_list.begin(), pair_list.end(), sort_descend() );
1860 
1861  // return an vector of indexes of constituents (sorted descending by pt)
1862  std::vector<Int_t> index_sorted_list;
1863  index_sorted_list.reserve ( entries );
1864 
1865  for ( std::vector< std::pair <Double_t, Int_t> >::iterator it = pair_list.begin(); it != pair_list.end(); ++it )
1866  {
1867  index_sorted_list.push_back ( ( *it ).second );
1868  } // populating the return object with indexes of sorted tracks
1869 
1870  return index_sorted_list;
1871  }
1872 
1873 //__________________________________________________________________________________________________
1874 std::vector<Int_t> AliAnalysisTaskEmcalJetCDF::SortTracksPt ( AliParticleContainer *trackscont ) const
1875  {
1876  //___________________________________________
1877  // Sorting by p_T (decreasing) event tracks
1878  //^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
1879  // Create vector for Pt sorting
1880  std::vector<ptidx_pair> pair_list;
1881 
1882  trackscont->ResetCurrentID();
1883  AliVTrack *track = NULL;
1884  UInt_t i_entry = 0;
1885  while( (track = dynamic_cast<AliVTrack*>(fTracksCont->GetNextAcceptParticle()) ))
1886  {
1887  i_entry++;
1888  pair_list.push_back ( std::make_pair ( track->Pt(), i_entry ) );
1889  }
1890 
1891  std::stable_sort ( pair_list.begin(), pair_list.end(), sort_descend() );
1892 
1893  // return an vector of indexes of constituents (sorted descending by pt)
1894  std::vector<Int_t> index_sorted_list;
1895  index_sorted_list.reserve ( i_entry );
1896 
1897  for ( std::vector< std::pair <Double_t, Int_t> >::iterator it = pair_list.begin(); it != pair_list.end(); ++it )
1898  {
1899  index_sorted_list.push_back ( ( *it ).second );
1900  } // populating the return object with indexes of sorted tracks
1901 
1902  return index_sorted_list;
1903  }
1904 
1905 //________________________________________________________________________
1906 Bool_t AliAnalysisTaskEmcalJetCDF::IdxInArray ( Int_t index, TArrayI &array )
1907  {
1908  for ( Int_t i = 0; i < array.GetSize(); i++ )
1909  {
1910  if ( index == array[i] ) { return kTRUE; }
1911  }
1912 
1913  return kFALSE;
1914  }
1915 
1916 //________________________________________________________________________
1918  {
1920 
1921  if ( fJetsCont ) { delete fJetsCont ; fJetsCont = 0; }
1922  if ( fTracksCont ) { delete fTracksCont ; fTracksCont = 0; }
1923  if ( fCaloClustersCont ) { delete fCaloClustersCont ; fCaloClustersCont = 0; }
1924  }
1925 
1926 //________________________________________________________________________
1928  {
1929  // Called once at the end of the analysis.
1930  // Update pointers reading them from the output slot
1931  fOutput = dynamic_cast<TList *> ( GetOutputData (0) );
1932  }
1933 
1934 //________________________________________________________________________
1935 Double_t AliAnalysisTaskEmcalJetCDF::Z_ptot( const AliEmcalJet* jet, const AliVParticle* trk) const
1936 {
1937  if (trk->P() < 1e-6) return 0.;
1938  return (trk != 0) ? trk->P()/ jet->P() : 0.;
1939 }
1940 
1941 //________________________________________________________________________
1942 Double_t AliAnalysisTaskEmcalJetCDF::Z_pt( const AliEmcalJet* jet, const AliVParticle* trk) const
1943 {
1944  if (trk->P() < 1e-6) return 0.;
1945  return (trk != 0) ? trk->Pt() / jet->Pt() : 0.;
1946 }
1947 
1948 // kate: indent-mode none; indent-width 2; replace-tabs on;
TProfile * fH7
! N(jet1) vs P_{T}(jet1)
TH1D * fH8_all
! Momentum distribution for jets (fragmentation function)
TProfile * fH15_n85
! <p_{T}> track vs the Distance R from Jet1 - 85% of particles
TH1D * fH15_bin
! p_{T} SUM track vs the Distance R from Jet1
TH1D * fHptd_pt85
! Distribution of dispersion pt_D
ClassImp(AliAnalysisTaskTriggerRates) AliAnalysisTaskTriggerRates
TProfile * fH15all_n70
! <p_{T}> track vs the Distance R from Jet1 - 70% of particles
TH1D * fH8xi_all
! Xi distribution for jets (fragmentation function)
virtual AliVParticle * GetNextAcceptParticle()
TH1D * fH8xi_p
! Xi distribution for leading jet (fragmentation function) - Z_ptot
TH1D * fH15_bin_pt90
! p_{T} SUM track vs the Distance R from Jet1 - 80% of Pt
TProfile * fH15_pt80
! <p_{T}> track vs the Distance R from Jet1 - 80% of Pt
TH1D * fH15all_bin_pt70
! p_{T} SUM track vs the Distance R from owner jet - 80% of Pt
TH1D * fH15all_bin_n75
! p_{T} SUM track vs the Distance R from owner jet - 80% of particles
TH1D * fH15_bin_pt80
! p_{T} SUM track vs the Distance R from Jet1 - 80% of Pt
Bool_t IdxInArray(Int_t index, TArrayI &array)
TProfile * fH15all_n90
! <p_{T}> track vs the Distance R from Jet1 - 90% of particles
Double_t GetXi(const AliVParticle *trk) const
TH1D * fH20_pt70
! Distribution of R in leading jet - 70% of Pt
AliJetContainer * fJetsCont
! Jets Container
TH1D * fH20_n85
! Distribution of R in leading jet - 85% of particles
AliJetContainer * GetJetContainer(Int_t i=0) const
Double_t Z_pt(const AliEmcalJet *jet, const AliVParticle *trk) const
Double_t Eta() const
Definition: AliEmcalJet.h:88
TProfile * fH25all
! Jet1 Size vs P_{T}(jet) - 80% of Pt
Double_t Phi() const
Definition: AliEmcalJet.h:84
TH1D * fH15_bin_pt85
! p_{T} SUM track vs the Distance R from Jet1 - 80% of Pt
TProfile * fH15all_n75
! <p_{T}> track vs the Distance R from Jet1 - 75% of particles
TProfile * fH15
! <p_{T}> track vs the Distance R from Jet1
TH1D * fH15all_bin
! p_{T} SUM track vs the Distance R from owner jet
TH1D * fH20all_pt80
! Distribution of R in jets - 80% of Pt
TH1D * fHptd_n70
! Distribution of dispersion pt_D
TH1D * fH8xi_all_pt
! Xi distribution for jets (fragmentation function) - Z_pt
TH1D * fH5
! Distribution of jets in events
TClonesArray * fCaloClustContArray
! the array of clusters from the tracks container
AliClusterContainer * GetClusterContainer() const
TProfile * fH7all
! N(jet1) vs P_{T} - all jets
TH1D * fH15_bin_pt70
! p_{T} SUM track vs the Distance R from Jet1 - 80% of Pt
TProfile * fH24
! Jet1 Size vs P_{T}(jet1) - 80% of particles
TH1D * fHptd_n90
! Distribution of dispersion pt_D
Double_t DeltaR(const AliVParticle *part1, const AliVParticle *part2)
TH1D * fHptd_n80
! Distribution of dispersion pt_D
TH1D * fH20all_n75
! Distribution of R in jets - 80% of particles
TList * fOutput
!output list
TH1D * fH15all_bin_n85
! p_{T} SUM track vs the Distance R from owner jet - 80% of particles
TH1D * fHptd_pt70
! Distribution of dispersion pt_D
TH1D * fH15all_bin_n70
! p_{T} SUM track vs the Distance R from owner jet - 80% of particles
UShort_t GetNumberOfConstituents() const
Definition: AliEmcalJet.h:107
Container for particles within the EMCAL framework.
UShort_t GetNumberOfTracks() const
Definition: AliEmcalJet.h:106
TH1D * fH20all_pt70
! Distribution of R in jets - 80% of Pt
std::ostream & Print(std::ostream &in) const
TH1D * fH20all_n85
! Distribution of R in jets - 80% of particles
Bool_t ProcessJetContainer(Int_t idx_jet_container=0)
TProfile * fH15_n90
! <p_{T}> track vs the Distance R from Jet1 - 90% of particles
AliParticleContainer * GetParticleContainer() const
AliEmcalJet * GetLeadingJet(const char *opt="")
TH1D * fHg_n75
! Distribution of girth (radial girth) g
TProfile * fH24all
! Jet1 Size vs P_{T}(jet) - 80% of particles
TH1D * fHg_n90
! Distribution of girth (radial girth) g
TH1D * fH1
! Pt distribution of all jets
TH1D * fH15all_bin_n90
! p_{T} SUM track vs the Distance R from owner jet - 80% of particles
TH1D * fHg_n70
! Distribution of girth (radial girth) g
Int_t fNaccPart
Multiplicity in event - accepted tracks in tracks container.
TH1D * fH20_n75
! Distribution of R in leading jet - 75% of particles
TH1D * fH15_bin_n90
! p_{T} SUM track vs the Distance R from Jet1 - 90% of particles
TProfile * fH23all
! Jet1 Size vs P_{T}(jet)
TH1D * fH15all_bin_pt90
! p_{T} SUM track vs the Distance R from owner jet - 80% of Pt
Int_t fNaccClus
Multiplicity in event - accepted clusters in cluster container.
TH1D * fH8xi_all_p
! Xi distribution for jets (fragmentation function) - Z_ptot
TH1D * fHg_pt90
! Distribution of girth (radial girth) g
TH1D * fH8xi_pt
! Xi distribution for leading jet (fragmentation function) - Z_pt
std::vector< int > SortConstituentsPt(TClonesArray *tracks) const
TProfile * fH23
! Jet1 Size vs P_{T}(jet1)
Int_t GetNJets() const
TH1D * fH6c
! Jet1 Multiplicity Distribution - all constituents
TH1D * fH15all_bin_pt85
! p_{T} SUM track vs the Distance R from owner jet - 80% of Pt
TH1D * fH15_bin_n80
! p_{T} SUM track vs the Distance R from Jet1 - 80% of particles
AliParticleContainer * fTracksCont
! Tracks Container
TH1D * fHptd_pt90
! Distribution of dispersion pt_D
TH1D * fH20all_pt90
! Distribution of R in jets - 80% of Pt
Int_t GetNAcceptedClusters() const
TProfile * fH25
! Jet1 Size vs P_{T}(jet1) - 80% of Pt
TH1D * fH20_n90
! Distribution of R in leading jet - 90% of particles
TClonesArray * GetArray() const
TH1D * fH15_bin_n70
! p_{T} SUM track vs the Distance R from Jet1 - 70% of particles
TH1D * fH15_bin_pt75
! p_{T} SUM track vs the Distance R from Jet1 - 80% of Pt
TH1D * fH2
! Eta distribution of all jets
TH1D * fH4c
! Multiplicity of all jets - all constituents // 1 unit of multiplicity /bin
TProfile * fH15all_pt75
! <p_{T}> track vs the Distance R from Jet1 - 75% of Pt
Double_t Xi(Double_t z) const
TClonesArray * fTracksContArray
! the array of tracks from the tracks container
TH1D * fH20all
! Distribution of R in jets
TH1D * fH15all_bin_pt75
! p_{T} SUM track vs the Distance R from owner jet - 80% of Pt
TProfile * fH15_pt70
! <p_{T}> track vs the Distance R from Jet1 - 70% of Pt
Double_t GetZ(const Double_t trkPx, const Double_t trkPy, const Double_t trkPz) const
TH1D * fH15all_bin_pt80
! p_{T} SUM track vs the Distance R from owner jet - 80% of Pt
TH1D * fH20_pt80
! Distribution of R in leading jet - 80% of Pt
void SetClassName(const char *clname)
TProfile * fH15all_pt80
! <p_{T}> track vs the Distance R from Jet1 - 80% of Pt
AliEmcalJet * GetNextAcceptJet()
Double_t DeltaR(const AliVParticle *part) const
TH1D * fHptd_pt75
! Distribution of dispersion pt_D
Double_t Pt() const
Definition: AliEmcalJet.h:76
TH1D * fH8
! Momentum distribution for leading jet (fragmentation function)
TH1D * fH20_pt90
! Distribution of R in leading jet - 90% of Pt
TH1D * fH20all_pt75
! Distribution of R in jets - 80% of Pt
TH1D * fHg_n85
! Distribution of girth (radial girth) g
TProfile * fH15_pt75
! <p_{T}> track vs the Distance R from Jet1 - 75% of Pt
Double_t P() const
Definition: AliEmcalJet.h:77
TH1D * fH20all_n80
! Distribution of R in jets - 80% of particles
Short_t TrackAt(Int_t idx) const
Definition: AliEmcalJet.h:127
TH1D * fHptd_n75
! Distribution of dispersion pt_D
TH1D * fH8_all_p
! Momentum distribution for jets (fragmentation function) - Z_ptot
TProfile * fH15_n75
! <p_{T}> track vs the Distance R from Jet1 - 75% of particles
TH1D * fHptd_n85
! Distribution of dispersion pt_D
TH1D * fHg_pt80
! Distribution of girth (radial girth) g
TProfile * fH15all_n80
! <p_{T}> track vs the Distance R from Jet1 - 80% of particles
TH1D * fHg_pt70
! Distribution of girth (radial girth) g
TH1D * fH8_p
! Momentum distribution for leading jet (fragmentation function) - Z_ptot
TProfile * fH15_pt90
! <p_{T}> track vs the Distance R from Jet1 - 80% of Pt
TH1D * fHg_pt85
! Distribution of girth (radial girth) g
TH1D * fHg
! Distribution of girth (radial girth) g = sum_jet_parts ( r_i * ( pt_i/pt_jet ) ) ...
TProfile * fH15_pt85
! <p_{T}> track vs the Distance R from Jet1 - 85% of Pt
Base task in the EMCAL jet framework.
Represent a jet reconstructed using the EMCal jet framework.
Definition: AliEmcalJet.h:44
TH1D * fH15_bin_n85
! p_{T} SUM track vs the Distance R from Jet1 - 85% of particles
TH1D * fHptd_pt80
! Distribution of dispersion pt_D
TH1D * fHg_pt75
! Distribution of girth (radial girth) g
TH1D * fH8xi
! Xi distribution for leading jet (fragmentation function)
TProfile * fH15all_pt70
! <p_{T}> track vs the Distance R from Jet1 - 70% of Pt
TProfile * fH15all_pt85
! <p_{T}> track vs the Distance R from Jet1 - 85% of Pt
TProfile * fH15all_pt90
! <p_{T}> track vs the Distance R from Jet1 - 90% of Pt
TH1D * fH20_n80
! Distribution of R in leading jet - 80% of particles
TH1D * fH20all_n90
! Distribution of R in jets - 80% of particles
Int_t GetNAcceptedParticles() const
TH1D * fHptd
! Distribution of dispersion d pt_D = sqrt ( sum (pt_i^2) )/sum (pt_i)
TH1D * fH15_bin_n75
! p_{T} SUM track vs the Distance R from Jet1 - 75% of particles
TH1D * fH4
! Multiplicity of all jets - tracks // 1 unit of multiplicity /bin
TH1D * fH20_pt75
! Distribution of R in leading jet - 75% of Pt
TH1D * fH3
! Phi distribution of all jets
TProfile * fH15all
! <p_{T}> track vs the Distance R from Jet1
TH1D * fHg_n80
! Distribution of girth (radial girth) g
TProfile * fH15all_n85
! <p_{T}> track vs the Distance R from Jet1 - 85% of particles
AliClusterContainer * fCaloClustersCont
! Clusters Container
Double_t Z_ptot(const AliEmcalJet *jet, const AliVParticle *trk) const
TH1D * fH20_n70
! Distribution of R in leading jet - 70% of particles
functional for sorting pair by first element - descending
Int_t idx_jetcont
index of jet container to be processed
TH1D * fH20all_n70
! Distribution of R in jets - 80% of particles
void ResetCurrentID(Int_t i=-1)
TH1D * fH15all_bin_n80
! p_{T} SUM track vs the Distance R from owner jet - 80% of particles
TProfile * fH15_n70
! <p_{T}> track vs the Distance R from Jet1 - 70% of particles
TProfile * fH15_n80
! <p_{T}> track vs the Distance R from Jet1 - 80% of particles
TH1D * fH8_pt
! Momentum distribution for leading jet (fragmentation function) - Z_pt
TH1D * fH20all_pt85
! Distribution of R in jets - 80% of Pt
TH1D * fH6
! Jet1 Multiplicity Distribution - charged tracks
TH1D * fH8_all_pt
! Momentum distribution for jets (fragmentation function) - Z_pt
Int_t fNJets_accepted
Number of Jets found in event - accepted cuts applied by JetContainer.
TH1D * fH20_pt85
! Distribution of R in leading jet - 85% of Pt
Analysis of leading jets distribution of pt and multiplicity.
std::vector< Int_t > SortTracksPt(AliVEvent *event) const
TH1D * fH20
! Distribution of R in leading jet