AliPhysics  vAN-20150425 (2dcf1b0)
 All Classes Namespaces Files Functions Variables Enumerations Enumerator Macros
AliAnaInsideClusterInvariantMass.h
Go to the documentation of this file.
1 #ifndef ALIANAINSIDECLUSTERINVARIANTMASS_H
2 #define ALIANAINSIDECLUSTERINVARIANTMASS_H
3 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
4  * See cxx source for full Copyright notice */
5 
6 //_________________________________________________________________________
12 //
20 
21 // --- ROOT system ---
22 class TList ;
23 class TObjString;
24 class TLorentzVector;
25 
26 // --- ANALYSIS system ---
27 class AliAODCaloCluster;
28 
30 
32 
33  public:
34 
36 
39 
40  void CheckLocalMaximaMCOrigin(AliVCluster* cluster, Int_t mcindex, Int_t noverlaps,
41  Float_t e1, Float_t e2, Float_t mass);
42  //, Float_t m02, TLorentzVector l1, TLorentzVector l2);
43 
44  TObjString * GetAnalysisCuts();
45 
46  TList * GetCreateOutputObjects();
47 
48  void GetMCIndex(AliVCluster * cluster, Int_t & mcindex, Int_t & mcTag);
49 
50  void GetMCPrimaryKine(AliVCluster* cluster, Int_t mcindex, Int_t mcTag, Bool_t matched,
51  Float_t & eprim, Float_t & asymGen, Float_t & angleGen, Int_t & noverlaps );
52 
53  void FillAngleHistograms(Int_t nMax, Bool_t matched, Int_t mcindex,
54  Float_t en, Float_t e1 , Float_t e2,
55  Float_t angle, Float_t mass, Float_t anglePrim,
56  Float_t m02, Float_t asym, Int_t pid, Int_t noverlaps);
57 
58  void FillArmenterosHistograms(Int_t nMax, Int_t ebin, Int_t mcindex,
59  Float_t pi0E, Float_t m02, Int_t pid);
60 
61  void FillThetaStarHistograms(Int_t nMax, Bool_t matched, Int_t mcindex,
62  Float_t pi0E, Float_t m02, Int_t pid);
63 
64  void FillEBinHistograms(Int_t ebin, Int_t nMax, Int_t mcindex, Float_t splitFrac,
65  Float_t mass, Float_t asym, Float_t l0);
66 
67  void FillMCHistograms(Float_t en, Float_t e1 , Float_t e2,
68  Int_t ebin, Int_t mcindex,Int_t noverlaps,
69  Float_t l0, Float_t mass,
70  Int_t nMax, Bool_t matched,
71  Float_t splitFrac, Float_t asym,
72  Float_t eprim, Float_t asymGen);
73 
74  void FillMCOverlapHistograms(Float_t en, Float_t enprim,
75  Int_t nc, Float_t mass, Float_t l0,
76  Float_t asym, Float_t splitFrac,
77  Int_t nlm, Int_t ebin, Bool_t matched,
78  Int_t mcindex, Int_t noverlaps);
79 
80  void FillSSWeightHistograms(AliVCluster *cluster, Int_t nlm, Int_t absId1, Int_t absId2);
81 
82  void FillSSExtraHistograms(AliVCluster *cluster, Int_t nMax,
83  Bool_t matched, Int_t mcindex,
84  Float_t mass , Int_t ebin);
85 
86  void FillNLMDiffCutHistograms(AliVCluster *cluster, AliVCaloCells *cells, Bool_t matched);
87 
88  void FillNCellHistograms(Int_t ncells, Float_t energy, Int_t nMax,
89  Bool_t matched, Int_t mcindex,
90  Float_t mass , Float_t l0);
91 
92  void FillTrackMatchingHistograms(AliVCluster * cluster,Int_t nMax, Int_t mcindex);
93 
94  void FillHistograms1(Float_t en, Float_t e1, Float_t e2,
95  Int_t nMax, Float_t mass, Float_t l0,
96  Float_t eta, Float_t phi,
97  Bool_t matched, Int_t mcindex);
98 
99 
100  void FillHistograms2(Float_t en, Float_t eprim,
101  Float_t e1, Float_t e2, Int_t nMax,
102  Float_t mass, Float_t l0,
103  Bool_t matched, Int_t mcindex);
104 
105  void FillIdPi0Histograms(Float_t en, Float_t e1, Float_t e2,
106  Int_t nc, Int_t nMax, Float_t t12diff,
107  Float_t mass, Float_t l0,
108  Float_t eta, Float_t phi,
109  Bool_t matched, Int_t mcindex);
110 
111  void FillIdEtaHistograms(Float_t en, Float_t e1, Float_t e2,
112  Int_t nc, Int_t nMax, Float_t t12diff,
113  Float_t mass, Float_t l0,
114  Float_t eta, Float_t phi,
115  Bool_t matched, Int_t mcindex);
116 
117  void FillIdConvHistograms(Float_t en, Int_t nMax, Float_t asym,
118  Float_t mass, Float_t l0,
119  Bool_t matched, Int_t mcindex);
120 
121  void Init();
122 
123  void InitParameters();
124 
126 
127  void Print(const Option_t * opt) const;
128 
129  void SetMinNCells(Int_t cut) { fMinNCells = cut ; }
130 
131  void SetMinBadChannelDistance(Float_t cut) { fMinBadDist = cut ; }
132 
133  void SetWCorrectionParameter(Int_t i, Float_t p = 0.07) { if( i<2 ) fWSimu[i] = p; }
134 
137 
140 
143 
146 
149 
152 
155 
158 
161 
164 
167 
170 
173 
176 
179 
182 
183  void SetNWeightForShowerShape(Int_t n) { fSSWeightN = n ; }
184  void SetWeightForShowerShape(Int_t i, Float_t v)
185  { if (i < 20) fSSWeight[i] = v ; }
186 
187  void SetNumberOfNLocMaxSettings(Int_t n) { fNLMSettingN = n ; }
188  void SetNLocMaxMinE (Int_t i, Float_t v) { if (i < 5) fNLMMinE [i] = v ; }
189  void SetNLocMaxMinDiff(Int_t i, Float_t v) { if (i < 5) fNLMMinDiff[i] = v ; }
190 
191 
192  void SetNECellCutForShowerShape(Int_t n) { fSSECellCutN = n ; }
193  void SetECellCutForShowerShape(Int_t i, Float_t v)
194  { if (i < 20) fSSECellCut[i] = v ; }
195 
196 
197  void RecalculateClusterShowerShapeParametersWithCellCut(const AliEMCALGeometry * geom, AliVCaloCells* cells, AliVCluster * cluster,
198  Float_t & l0, Float_t & l1,
199  Float_t & disp, Float_t & dEta, Float_t & dPhi,
200  Float_t & sEta, Float_t & sPhi, Float_t & sEtaPhi,Float_t eCellMin = 0.);
201 
202 
204  enum mcTypes { kmcPhoton = 1, kmcConversion = 2, kmcPi0 = 3, kmcPi0Conv = 4,
205  kmcEta = 5, kmcHadron = 6 };
206 
207  private:
208 
209  Int_t fMinNCells ;
210  Float_t fMinBadDist ;
211  Float_t fHistoECut ;
213 
215  Bool_t fFillTMHisto ;
218  Bool_t fFillMCHisto ;
221  Bool_t fFillEbinHisto ;
223  Bool_t fFillNCellHisto ;
225  Bool_t fFillIdEtaHisto ;
229 
230  Float_t fSSWeight[20];
231  Int_t fSSWeightN;
232 
233  Float_t fSSECellCut[20];
234  Int_t fSSECellCutN;
235 
236  Float_t fNLMMinE [5];
237  Float_t fNLMMinDiff[5];
238  Int_t fNLMSettingN;
239 
240  Float_t fWSimu[2];
241 
243  TLorentzVector fClusterMomentum;
244  TLorentzVector fSubClusterMom1;
245  TLorentzVector fSubClusterMom2;
246  TLorentzVector fSubClusterMomSum;
247  TLorentzVector fSubClusterMomBoost;
248 
249  TLorentzVector fPrimaryMom;
250  TLorentzVector fGrandMotherMom;
251  TLorentzVector fMCDaughMom1;
252  TLorentzVector fMCDaughMom2;
253  TVector3 fProdVertex;
254 
255  //Histograms
256 
257  TH2F * fhMassNLocMax1[7][2] ;
258  TH2F * fhMassNLocMax2[7][2] ;
259  TH2F * fhMassNLocMaxN[7][2] ;
260 
261  TH2F * fhMassSplitENLocMax1[7][2] ;
262  TH2F * fhMassSplitENLocMax2[7][2] ;
263  TH2F * fhMassSplitENLocMaxN[7][2] ;
264 
265  TH2F * fhAsymNLocMax1[7][2] ;
266  TH2F * fhAsymNLocMax2[7][2] ;
267  TH2F * fhAsymNLocMaxN[7][2] ;
268 
272 
276 
280 
284 
288 
292 
296 
300 
304 
308 
309  TH2F * fhMassM02NLocMax1[7][2] ;
310  TH2F * fhMassM02NLocMax2[7][2] ;
311  TH2F * fhMassM02NLocMaxN[7][2] ;
312 
316 
320 
324 
327 
328  TH2F * fhMassDispEtaNLocMax1[7][2] ;
329  TH2F * fhMassDispEtaNLocMax2[7][2] ;
330  TH2F * fhMassDispEtaNLocMaxN[7][2] ;
331 
335 
336  TH2F * fhMassDispPhiNLocMax1[7][2] ;
337  TH2F * fhMassDispPhiNLocMax2[7][2] ;
338  TH2F * fhMassDispPhiNLocMaxN[7][2] ;
339 
343 
344  TH2F * fhMassDispAsyNLocMax1[7][2] ;
345  TH2F * fhMassDispAsyNLocMax2[7][2] ;
346  TH2F * fhMassDispAsyNLocMaxN[7][2] ;
347 
351 
352  TH2F * fhNLocMax [7][2] ;
353  TH2F * fhNLocMaxM02Cut[7][2] ;
354  TH2F * fhNLocMaxIdPi0 [7][2] ;
355 
356  TH2F * fhSplitClusterENLocMax[7][2] ;
358 
359  TH2F * fhLM1NLocMax [7][2] ;
360  TH2F * fhLM1NLocMaxM02Cut[7][2] ;
361  TH2F * fhLM1NLocMaxIdPi0 [7][2] ;
362 
363  TH2F * fhLM2NLocMax [7][2] ;
364  TH2F * fhLM2NLocMaxM02Cut[7][2] ;
365  TH2F * fhLM2NLocMaxIdPi0 [7][2] ;
366 
367  TH2F * fhM02NLocMax1 [7][2] ;
368  TH2F * fhM02NLocMax2 [7][2] ;
369  TH2F * fhM02NLocMaxN [7][2] ;
370 
374 
375  TH2F * fhMCGenFracNLocMax1[7][2] ;
376  TH2F * fhMCGenFracNLocMax2[7][2] ;
377  TH2F * fhMCGenFracNLocMaxN[7][2] ;
378 
382 
386 
390 
394 
398 
402 
403  TH2F * fhMCGenEvsSplitENLocMax1[7][2] ;
404  TH2F * fhMCGenEvsSplitENLocMax2[7][2] ;
405  TH2F * fhMCGenEvsSplitENLocMaxN[7][2] ;
406 
407  TH2F * fhMCGenFracNLocMaxEbin[7][4] ;
409 
413 
417 
418  TH2F * fhNCellNLocMax1[7][2] ;
419  TH2F * fhNCellNLocMax2[7][2] ;
420  TH2F * fhNCellNLocMaxN[7][2] ;
421 
426 
431 
436 
437  TH2F * fhM02Pi0NLocMax1[7][2] ;
438  TH2F * fhM02EtaNLocMax1[7][2] ;
439  TH2F * fhM02ConNLocMax1[7][2] ;
440 
441  TH2F * fhM02Pi0NLocMax2[7][2] ;
442  TH2F * fhM02EtaNLocMax2[7][2] ;
443  TH2F * fhM02ConNLocMax2[7][2] ;
444 
445  TH2F * fhM02Pi0NLocMaxN[7][2] ;
446  TH2F * fhM02EtaNLocMaxN[7][2] ;
447  TH2F * fhM02ConNLocMaxN[7][2] ;
448 
449  TH2F * fhMassPi0NLocMax1[7][2] ;
450  TH2F * fhMassEtaNLocMax1[7][2] ;
451  TH2F * fhMassConNLocMax1[7][2] ;
452 
453  TH2F * fhMassPi0NLocMax2[7][2] ;
454  TH2F * fhMassEtaNLocMax2[7][2] ;
455  TH2F * fhMassConNLocMax2[7][2] ;
456 
457  TH2F * fhMassPi0NLocMaxN[7][2] ;
458  TH2F * fhMassEtaNLocMaxN[7][2] ;
459  TH2F * fhMassConNLocMaxN[7][2] ;
460 
461  TH2F * fhNCellPi0NLocMax1[7][2] ;
462  TH2F * fhNCellEtaNLocMax1[7][2] ;
463  TH2F * fhNCellPi0NLocMax2[7][2] ;
464  TH2F * fhNCellEtaNLocMax2[7][2] ;
465  TH2F * fhNCellPi0NLocMaxN[7][2] ;
466  TH2F * fhNCellEtaNLocMaxN[7][2] ;
467 
468  TH2F * fhMassAfterCutsNLocMax1[7][2] ;
469  TH2F * fhMassAfterCutsNLocMax2[7][2] ;
470  TH2F * fhMassAfterCutsNLocMaxN[7][2] ;
471 
475 
476  TH2F * fhMassSplitEPi0NLocMax1[7][2] ;
477  TH2F * fhMassSplitEPi0NLocMax2[7][2] ;
478  TH2F * fhMassSplitEPi0NLocMaxN[7][2] ;
479 
480  TH2F * fhAsyPi0NLocMax1[7][2] ;
481  TH2F * fhAsyEtaNLocMax1[7][2] ;
482  TH2F * fhAsyConNLocMax1[7][2] ;
483 
484  TH2F * fhAsyPi0NLocMax2[7][2] ;
485  TH2F * fhAsyEtaNLocMax2[7][2] ;
486  TH2F * fhAsyConNLocMax2[7][2] ;
487 
488  TH2F * fhAsyPi0NLocMaxN[7][2] ;
489  TH2F * fhAsyEtaNLocMaxN[7][2] ;
490  TH2F * fhAsyConNLocMaxN[7][2] ;
491 
492  TH2F * fhSplitEFractionNLocMax1[7][2] ;
493  TH2F * fhSplitEFractionNLocMax2[7][2] ;
494  TH2F * fhSplitEFractionNLocMaxN[7][2] ;
495 
499 
503 
504  TH2F * fhAnglePairNLocMax1[7][2] ;
505  TH2F * fhAnglePairNLocMax2[7][2] ;
506  TH2F * fhAnglePairNLocMaxN[7][2] ;
507 
511 
512  TH2F * fhAnglePairPi0NLocMax1[7][2] ;
513  TH2F * fhAnglePairPi0NLocMax2[7][2] ;
514  TH2F * fhAnglePairPi0NLocMaxN[7][2] ;
515 
516  TH2F * fhAnglePairMassNLocMax1[7][2] ;
517  TH2F * fhAnglePairMassNLocMax2[7][2] ;
518  TH2F * fhAnglePairMassNLocMaxN[7][2] ;
519 
520  TH2F * fhAnglePairM02NLocMax1[7][2] ;
521  TH2F * fhAnglePairM02NLocMax2[7][2] ;
522  TH2F * fhAnglePairM02NLocMaxN[7][2] ;
523 
527 
531 
535 
539 
543 
544  TH2F * fhArmNLocMax1[7][4] ;
545  TH2F * fhArmNLocMax2[7][4] ;
546  TH2F * fhArmNLocMaxN[7][4] ;
547 
548  TH2F * fhArmAfterCutsNLocMax1[7][4] ;
549  TH2F * fhArmAfterCutsNLocMax2[7][4] ;
550  TH2F * fhArmAfterCutsNLocMaxN[7][4] ;
551 
552  TH2F * fhArmPi0NLocMax1[7][4] ;
553  TH2F * fhArmPi0NLocMax2[7][4] ;
554  TH2F * fhArmPi0NLocMaxN[7][4] ;
555 
556  TH2F * fhCosThStarNLocMax1[7][2] ;
557  TH2F * fhCosThStarNLocMax2[7][2] ;
558  TH2F * fhCosThStarNLocMaxN[7][2] ;
559 
563 
564  TH2F * fhCosThStarPi0NLocMax1[7][2] ;
565  TH2F * fhCosThStarPi0NLocMax2[7][2] ;
566  TH2F * fhCosThStarPi0NLocMaxN[7][2] ;
567 
574 
581 
588 
595 
602 
612 
613  TH2F * fhPi0CellE[3] ;
614  TH2F * fhPi0CellEFrac[3] ;
615  TH2F * fhPi0CellLogEFrac[3] ;
619  TH2F * fhPi0CellEMaxFrac [3];
620  TH2F * fhPi0CellEMax2Frac [3];
621 
622  TH2F * fhM02WeightPi0[3][20] ;
623  TH2F * fhM02ECellCutPi0[3][20] ;
624 
631 
632  TH2F * fhMCEM02Overlap0[3][7];
633  TH2F * fhMCEM02Overlap1[3][7];
634  TH2F * fhMCEM02OverlapN[3][7];
635  TH2F * fhMCEM02Overlap0Match[3][7];
636  TH2F * fhMCEM02Overlap1Match[3][7];
637  TH2F * fhMCEM02OverlapNMatch[3][7];
638 
639  TH2F * fhMCEMassOverlap0[3][7];
640  TH2F * fhMCEMassOverlap1[3][7];
641  TH2F * fhMCEMassOverlapN[3][7];
642  TH2F * fhMCEMassOverlap0Match[3][7];
643  TH2F * fhMCEMassOverlap1Match[3][7];
644  TH2F * fhMCEMassOverlapNMatch[3][7];
645 
646  TH2F * fhMCESplitEFracOverlap0[3][7];
647  TH2F * fhMCESplitEFracOverlap1[3][7];
648  TH2F * fhMCESplitEFracOverlapN[3][7];
652 
653  TH2F * fhMCEAsymOverlap0[3][7];
654  TH2F * fhMCEAsymOverlap1[3][7];
655  TH2F * fhMCEAsymOverlapN[3][7];
656  TH2F * fhMCEAsymOverlap0Match[3][7];
657  TH2F * fhMCEAsymOverlap1Match[3][7];
658  TH2F * fhMCEAsymOverlapNMatch[3][7];
659 
660  TH2F * fhMCENCellOverlap0[3][7];
661  TH2F * fhMCENCellOverlap1[3][7];
662  TH2F * fhMCENCellOverlapN[3][7];
663  TH2F * fhMCENCellOverlap0Match[3][7];
664  TH2F * fhMCENCellOverlap1Match[3][7];
665  TH2F * fhMCENCellOverlapNMatch[3][7];
666 
667  TH2F * fhMCEEpriOverlap0[3][7];
668  TH2F * fhMCEEpriOverlap1[3][7];
669  TH2F * fhMCEEpriOverlapN[3][7];
670  TH2F * fhMCEEpriOverlap0Match[3][7];
671  TH2F * fhMCEEpriOverlap1Match[3][7];
672  TH2F * fhMCEEpriOverlapNMatch[3][7];
673 
674  TH2F * fhMCEEpriOverlap0IdPi0[3][7];
675  TH2F * fhMCEEpriOverlap1IdPi0[3][7];
676  TH2F * fhMCEEpriOverlapNIdPi0[3][7];
677 
678  TH2F * fhMCPi0MassM02Overlap0[3][4];
679  TH2F * fhMCPi0MassM02Overlap1[3][4];
680  TH2F * fhMCPi0MassM02OverlapN[3][4];
684 
685  TH2F * fhMCENOverlaps[3][7];
686  TH2F * fhMCENOverlapsMatch[3][7];
687 
692  // No match between highest energy local maxima and highest energy MC particle
697 
702  // No match between highest energy local maxima and highest energy MC particle
707 
714 
721 
726 
731 
736 
741 
742 
747 
752 
757 
762 
763 
770 
777 
778 
781 
782  TH2F * fhMassBadDistClose[3];
783  TH2F * fhM02BadDistClose[3];
784  TH2F * fhMassOnBorder[3];
785  TH2F * fhM02OnBorder[3];
786 
787 
788  TH2F * fhNLocMaxDiffCut [5][5] [2] ;
789  TH2F * fhM02NLocMaxDiffCut[5][5][3][2] ;
790  TH2F * fhMassNLocMaxDiffCut[5][5][3][2] ;
791 
792  TH2F * fhNLocMaxDiffCutPi0 [5][5] [2] ;
793  TH2F * fhM02NLocMaxDiffCutPi0[5][5][3][2] ;
794  TH2F * fhMassNLocMaxDiffCutPi0[5][5][3][2] ;
795 
798 
801 
803  ClassDef(AliAnaInsideClusterInvariantMass,30) ;
805 
806 } ;
807 
808 #endif //ALIANAINSIDECLUSTERINVARIANTMASS_H
809 
810 
811 
TH2F * fhMCENCellOverlap0[3][7]
! E vs sum of splitted cluster energy asymmetry for different MC origin, no other MC particles contri...
TH2F * fhMassAsyCutNLocMaxN
! Mass(E) asym selection, not matched, Mass of split clusters, NLM > 2
Bool_t fFillMCHisto
Fill MC energy fraction histos.
Bool_t fFillArmenterosHisto
Fill armenteros type histo.
TH2F * fhMCEM02OverlapNMatch[3][7]
! E vs M02 for different MC origin, N other MC particles contributes, charged cluster ...
Float_t fSSWeight[20]
List of weights to test.
TH2F * fhMCGenEFracvsSplitEFracNLocMaxN[7][2]
! E generated particle / E reconstructed vs E1+E2 reconstructed / E reconstructed for N max in cluste...
TH2F * fhMCPi0DecayPhotonAdjHighLMOverlapMass[3]
! E vs Mass when cluster originated in pi0 merging and MC photon decay hit the adjacent cell local ma...
Bool_t fFillAngleHisto
Fill splitted clusters angle histograms.
TH2F * fhMassAsyCutNLocMax1
! Mass(E) asym selection, not matched, Mass of split clusters, NLM = 1
TH2F * fhMCPi0DecayPhotonAdjHighLM
! E vs NLM when cluster originated in pi0 merging and MC photon decay hit the adjacent cell local max...
TH2F * fhMCGenFracNLocMax2[7][2]
! E generated particle / E reconstructed vs E reconstructed for N max in cluster = 2...
TH2F * fhCosThStarNLocMax1[7][2]
! cos(theta^star) vs E, NLM=1
TH2F * fhMCGenSplitEFracNLocMax2NoOverlap[7][2]
! E generated particle / E1+E2 reconstructed vs E reconstructed for N max in cluster = 2...
TH2F * fhMCGenSplitEFracNLocMaxNNoOverlap[7][2]
! E generated particle / E1+E2 reconstructed vs E reconstructed for N max in cluster > 2...
TH2F * fhMCPi0DecayPhotonHitHighLMDiffELM2[3]
! E vs Ephoton-Esplit when cluster originated in pi0 merging and MC photon decay hit the cell local m...
TH2F * fhMCPi0DecayPhotonHitOtherLM
! E vs NLM when cluster originated in pi0 merging and MC photon decay hit the cell local maximas...
TH2F * fhEtaEPairDiffTimeNLM2
! E vs Pair of clusters time difference vs E, for selected eta, NLM=2
TH2F * fhAnglePairPrimPi0vsRecoNLocMax1
! Pair opening angle pi0 generated vs reconstructed
TH2F * fhMCPi0MassM02Overlap0[3][4]
! MC Pi0 M02 vs Mass for different MC origin, no other MC particles contributes, neutral cluster...
TH2F * fhMassSplitEFractionNLocMax1Ebin[7][4]
! Mass vs sum of splitted cluster energy / cluster energy for N max in cluster = 1, 1-6 for different MC particle types, not track matched
Bool_t fFillTMHisto
Fill track matching histos,.
TH2F * fhMassSplitEAfterCutsNLocMax1[7][2]
! Split Inv Mass vs E1+E2, NLM=1, after M02, asymmetry cuts, different MC particle types...
TH2F * fhSplitEFractionAfterCutsNLocMax2[7][2]
! sum of splitted cluster energy / cluster energy for N Local Maxima = 2, after M02 and asymmetry cut...
TH2F * fhCosThStarPi0NLocMax2[7][2]
! cos(theta^star) vs E, after M02, asymmetry and pi0 mass cuts, NLM=2
TH2F * fhAnglePairMassNLocMax2[7][2]
! Pair opening angle vs Mass for E > 7 GeV
void GetMCPrimaryKine(AliVCluster *cluster, Int_t mcindex, Int_t mcTag, Bool_t matched, Float_t &eprim, Float_t &asymGen, Float_t &angleGen, Int_t &noverlaps)
Check origin of the candidates, get primary kinematics if overlapped meson decay. ...
TH2F * fhSplitEFractionNLocMaxN[7][2]
! sum of splitted cluster energy / cluster energy for N Local Maxima > 2
TH2F * fhAnglePairOverM02NLocMaxN[7][2]
! Pair opening angle / m02 vs E, NLM=N
TH2F * fhSplitEFractionvsAsyNLocMax2[2]
! sum of splitted cluster energy / cluster energy for N Local Maxima = 2 vs |A|
TH2F * fhMassM02NLocMax2[7][2]
! Mass of splitted clusters when 2 local max vs M02, for E > 8 GeV, 1-6 for different MC particle typ...
TH2F * fhAnglePairOverM02NLocMax1[7][2]
! Pair opening angle / m02 vs E, NLM=1
TH2F * fhM02WeightPi0[3][20]
! M02 for selected pi0 with different weight, with NLM = 1, = 2, > 2
TH2F * fhMCGenEvsSplitENLocMaxN[7][2]
! E generated particle vs E1+E2 for N max in cluster > 2, 1-6 for different MC particle types ...
TH2F * fhMCESplitEFracOverlap1Match[3][7]
! E vs sum of splitted cluster energy / cluster energy for different MC origin, 1 other MC particles ...
TH2F * fhMCPi0DecayPhotonAdjOtherLM
! E vs NLM when cluster originated in pi0 merging and MC photon decay do not hit the adjacent cell lo...
TH2F * fhMassEnCutNLocMax2
! Mass(E) E sub-cluster cut selection, not matched, Mass of split clusters, NLM = 1 ...
TH2F * fhMCPi0DecayPhotonAdjacentMass[3]
! E vs Mass when cluster originated in pi0 merging and MC photon decay hit adjacen cells...
Int_t fNLMSettingN
Total number of NLM settings to test.
void FillIdPi0Histograms(Float_t en, Float_t e1, Float_t e2, Int_t nc, Int_t nMax, Float_t t12diff, Float_t mass, Float_t l0, Float_t eta, Float_t phi, Bool_t matched, Int_t mcindex)
Fill histograms for clusters passing the pi0 selection.
TH2F * fhArmPi0NLocMax2[7][4]
! Armenteros for selected pi0, N Local Maxima = 2
Bool_t fFillIdConvHisto
Fill histograms for clusters identified as conversion.
TH2F * fhMCPi0DecayPhotonAdjHighLMOverlapDiffELM2[3]
! E vs Ephoton-Esplit when cluster originated in pi0 merging and MC photon decay hit the adjacent cel...
TH2F * fhM02NLocMax2[7][2]
! M02 vs E for N max in cluster = 2, 1-6 for different MC particle types
TH2F * fhArmNLocMaxN[7][4]
! Armenteros of >2 cells local maxima vs E, 1-6 for different MC particle types
TH2F * fhMassNLocMax2[7][2]
! Split Inv Mass vs cluster E, NLM=2, different MC particle types, track matching on/off ...
Int_t fSSWeightN
Total number of weights to test.
TH2F * fhMCPi0DecayPhotonHitHighLMOverlapMass[3]
! E vs Mass when cluster originated in pi0 merging and MC photon decay hit the cell local maxima...
void FillAngleHistograms(Int_t nMax, Bool_t matched, Int_t mcindex, Float_t en, Float_t e1, Float_t e2, Float_t angle, Float_t mass, Float_t anglePrim, Float_t m02, Float_t asym, Int_t pid, Int_t noverlaps)
Fill histograms related to opening angle of the split clusters.
TH2F * fhMCGenEvsSplitENLocMax1[7][2]
! E generated particle vs E1+E2 for N max in cluster = 1, 1-6 for different MC particle types ...
TH2F * fhMCPi0MassM02Overlap1Match[3][4]
! MC Pi0 M02 vs Mass for different MC origin, 1 other MC particles contributes, charged cluster...
TH2F * fhMCEMassOverlap0Match[3][7]
! E vs Mass for different MC origin, no other MC particles contributes, charged cluster ...
TH2F * fhMCPi0DecayPhotonAdjHighLMOverlapDiffELM2vsELM2[3]
! E vs Ephoton-Esplit when cluster originated in pi0 merging and MC photon decay hit the adjacent cel...
TH2F * fhMassSplitECutNLocMaxN
! 85% of split energy, not matched, Mass of split clusters, NLM > 2
TH2F * fhMassM02NLocMax1[7][2]
! Mass of splitted clusters when 1 local max vs M02, for E > 8 GeV, 1-6 for different MC particle typ...
TH2F * fhMCPi0DecayPhotonAdjOtherLMOverlap
! E vs NLM when cluster originated in pi0 merging and MC photon decay do not hit the adjacent cell lo...
TH2F * fhMCGenSplitEFracNLocMax1NoOverlap[7][2]
! E generated particle / E1+E2 reconstructed vs E reconstructed for N max in cluster = 1...
TH2F * fhAnglePairM02NLocMaxN[7][2]
! Pair opening angle vs M02 for E > 7 GeV
TH2F * fhMassDispEtaNLocMaxN[7][2]
! Mass of >2 cells local maxima, vs M02, for E > 8 GeV, 1-6 for different MC particle types ...
TH2F * fhCosThStarNLocMax2[7][2]
! cos(theta^star) vs E, NLM=2
TH2F * fhEtaEPairDiffTimeNLMN
! E vs Pair of clusters time difference vs E, for selected eta, NLM>2
TH2F * fhNCellM02ELowNLocMax1MCPi0
! n cells in cluster vs m02 for low energy clusters, for N max in cluster = 1, for MC pi0 ...
TH2F * fhMCEOverlapTypeMatch
! What particles overlap with pi0, charged clusters
TH2F * fhAsyEtaNLocMaxN[7][2]
! Asy for Mass around eta, N Local Maxima > 2
TH2F * fhMCPi0DecayPhotonHitOtherLMDiffELM2[3]
! E vs Ephoton-Esplit when cluster originated in pi0 merging and MC photon decay hit the cell local m...
TH2F * fhArmNLocMax2[7][4]
! Armenteros of 2 cells local maxima vs E, 1-6 for different MC particle types
TH2F * fhPi0EPairDiffTimeNLM2
! E vs Pair of clusters time difference vs E, for selected pi0, NLM=2
TH2F * fhMCENCellOverlap1Match[3][7]
! E vs sum of splitted cluster energy asymmetry for different MC origin, 1 other MC particles contrib...
TH2F * fhMCGenEFracvsSplitEFracNLocMax1[7][2]
! E generated particle / E reconstructed vs E1+E2 reconstructed / E reconstructed for N max in cluste...
TH2F * fhAnglePairPi0NLocMax1[7][2]
! Pair opening angle vs E
TLorentzVector fGrandMotherMom
! Primary momentum
TH2F * fhMCESplitEFracOverlap0Match[3][7]
! E vs sum of splitted cluster energy / cluster energy for different MC origin, no other MC particles...
void FillMCHistograms(Float_t en, Float_t e1, Float_t e2, Int_t ebin, Int_t mcindex, Int_t noverlaps, Float_t l0, Float_t mass, Int_t nMax, Bool_t matched, Float_t splitFrac, Float_t asym, Float_t eprim, Float_t asymGen)
Fill histograms depending on MC input.
TH2F * fhTrackMatchedDPhiNLocMax2[7]
! Phi distance between track and cluster vs cluster E, 2 local maximum
TH2F * fhM02ConNLocMaxN[7][2]
! M02 for Mass around close to 0, N Local Maxima > 2
TH2F * fhNLocMaxDiffCutPi0[5][5][2]
! Number of maxima for different values of min Loc Max value and min difference between cells...
TH2F * fhMCPi0DecayPhotonAdjOtherLMOverlapDiffELM1vsELM1[3]
! E vs Ephoton-Esplit when cluster originated in pi0 merging and MC photon decay do not hit the adjac...
TH2F * fhM02NLocMaxDiffCut[5][5][3][2]
! M02 for 3 kinds of number of maxima for different values of min Loc Max value and min difference be...
TH2F * fhMassSplitEPi0NLocMaxN[7][2]
! Split Inv Mass vs E1+E2, NLM>2, after pi0 selection, different MC particle types, track matching on/off
TH2F * fhTrackMatchedDPhiNLocMax2Neg[7]
! Phi distance between track and cluster vs cluster E, 2 local maximum
Float_t fMinBadDist
Minimal distance to bad channel to accept cluster.
TH2F * fhMassDispPhiNLocMax2[7][2]
! Mass of 2 cells local maxima, vs M02, for E > 8 GeV, 1-6 for different MC particle types ...
TH2F * fhAnglePairNLocMax1[7][2]
! Pair opening angle vs E
TH2F * fhPi0CellEMaxClusterFrac[3]
! pi0's energy vs energy fraction of main LM and cluster energy with NLM = 1, = 2, > 2
TH2F * fhNCellMassELowNLocMax2MCPi0
! n cells in cluster vs mass for low energy clusters, for N max in cluster = 2, for MC pi0 ...
TH2F * fhMassNLocMaxDiffCut[5][5][3][2]
! Mass for 3 kinds of number of maxima for different values of min Loc Max value and min difference b...
TH2F * fhPi0CellEMax2ClusterFrac[3]
! pi0's energy vs energy fraction of second LM and cluster energy with NLM = 1, = 2...
TH2F * fhCosThStarPi0NLocMax1[7][2]
! cos(theta^star) vs E, after M02, asymmetry and pi0 mass cuts, NLM=1
TH2F * fhMCPi0DecayPhotonHitHighLMOverlapDiffELM1[3]
! E vs Ephoton-Esplit cluster when cluster originated in pi0 merging and MC photon decay hit the cell...
TH2F * fhMCEM02OverlapN[3][7]
! E vs M02 for different MC origin, N other MC particles contributes, neutral cluster ...
TH2F * fhTrackMatchedDPhiNLocMaxNNeg[7]
! Phi distance between track and cluster vs cluster E, more than 2 local maximum
TH2F * fhMCEAsymOverlap0Match[3][7]
! E vs sum of splitted cluster energy asymmetry for different MC origin, no other MC particles contri...
TH2F * fhSplitEFracEnCutNLocMaxN
! Split E fraction (E) E sub-cluster cut selection, not matched, Mass of split clusters, NLM > 2
TH2F * fhNCellMassEHighNLocMax1MCPi0
! n cells in cluster vs mass for high energy clusters, for N max in cluster = 1, for MC pi0 ...
TH2F * fhMCENOverlaps[3][7]
! E vs number of Overlaps in MC, neutral cluster
TH2F * fhM02Pi0NLocMax1[7][2]
! M02 for Mass around pi0, N Local Maxima = 1
TH2F * fhPi0CellEFrac[3]
! pi0's energy vs cluster cell energy fraction with NLM = 1, = 2, > 2
TH2F * fhMassSplitEPi0NLocMax2[7][2]
! Split Inv Mass vs E1+E2, NLM=2, after pi0 selection, different MC particle types, track matching on/off
TH2F * fhMCEEpriOverlap1[3][7]
! E reco vs primary for different MC origin, 1 other MC particles contributes, neutral cluster ...
TH2F * fhClusterEtaPhiNLocMax2
! Eta vs Phi of clusters with N Local Maxima = 2, E > 8 GeV
TH2F * fhMCPi0DecayPhotonAdjHighLMOverlap
! E vs NLM when cluster originated in pi0 merging and MC photon decay hit the adjacent cell local max...
TH2F * fhM02ConNLocMax1[7][2]
! M02 for Mass around close to 0, N Local Maxima = 1
TH2F * fhMCPi0NoneNLMPairNoMCMatchOverlap
! E vs NLM when cluster originated in pi0 merging and a both no NLM corresponds to the photons...
TLorentzVector fClusterMomentum
! Cluster momentum
TH2F * fhMCENCellOverlapN[3][7]
! E vs sum of splitted cluster energy asymmetry for different MC origin, N other MC particles contrib...
TH2F * fhMCPi0DecayPhotonAdjHighLMDiffELM2[3]
! E vs Ephoton-Esplit when cluster originated in pi0 merging and MC photon decay hit the adjacent cel...
TH2F * fhSplitEFractionvsAsyNLocMaxN[2]
! sum of splitted cluster energy / cluster energy for N Local Maxima > 2 vs |A|
TH2F * fhMassM02NLocMaxN[7][2]
! Mass of splitted clusters when >2 local max vs M02, for E > 8 GeV, 1-6 for different MC particle ty...
TH2F * fhMCPi0DecayPhotonHitOtherLMDiffELM1[3]
! E vs Ephoton-Esplit when cluster originated in pi0 merging and MC photon decay hit the cell local m...
TH2F * fhMassConNLocMaxN[7][2]
! Mass for selected around close to 0, N Local Maxima > 2
TH2F * fhAsymNLocMax1[7][2]
! Asymmetry of 2 highest energy cells when 1 local max vs E, 1-6 for different MC particle types ...
TH2F * fhEtaEtaPhiNLocMax2
! Eta vs Phi of eta's with N Local Maxima = 2, E > 8 GeV
TH2F * fhMCGenFracAfterCutsNLocMax2MCPi0
! E generated particle / E reconstructed vs E reconstructed for N max in cluster = 2...
TH2F * fhAsyConNLocMax1[7][2]
! Asy for Mass around close to 0, N Local Maxima = 1
TH2F * fhArmNLocMax1[7][4]
! Armenteros of 2 highest energy cells when 1 local max vs E, 1-6 for different MC particle types ...
TH2F * fhMassDispEtaNLocMax1Ebin[4]
! Mass of 2 highest energy cells when 1 local max, vs M02, 4 E bins, neutral clusters ...
TH2F * fhCentralityEtaNLocMax1
! Centrality for selected eta, N Local Maxima = 1
TH2F * fhSplitEFracEnCutNLocMax2
! Split E fraction (E) E sub-cluster cut selection, not matched, Mass of split clusters, NLM = 1
TH2F * fhNCellNLocMaxN[7][2]
! n cells in cluster vs E for N max in cluster > 2, 1-6 for different MC particle types ...
TH2F * fhAnglePairMassNLocMaxN[7][2]
! Pair opening angle vs Mass for E > 7 GeV
TH2F * fhMCEAsymOverlapN[3][7]
! E vs sum of splitted cluster energy asymmetry for different MC origin, N other MC particles contrib...
TH2F * fhMCPi0DecayPhotonHitHighLM
! E vs NLM when cluster originated in pi0 merging and MC photon decay hit the cell local maxima ...
TH2F * fhCentralityPi0NLocMaxN
! Centrality for selected pi0, N Local Maxima > 2
void FillArmenterosHistograms(Int_t nMax, Int_t ebin, Int_t mcindex, Float_t pi0E, Float_t m02, Int_t pid)
Fill Armeteros type histograms of split clusters.
TH2F * fhMassM02CutNLocMax1
! Mass(E) M02 selection, not matched, Mass of split clusters, NLM = 1
TH2F * fhMassM02NLocMax2Ebin[4]
! Mass of splitted clusters when 2 local max vs M02, 4 E bins, neutral clusters
void SetWCorrectionParameter(Int_t i, Float_t p=0.07)
TH2F * fhSplitClusterEPi0NLocMax[7][2]
! Number of maxima in cluster vs E of splitted clusters when cluster id as pi0, 1-6 for different MC ...
TVector3 fProdVertex
! primary production vertex
TH2F * fhMCPi0DecayPhotonHitOtherLMOverlapDiffELM2[3]
! E vs Ephoton-Esplit when cluster originated in pi0 merging and MC photon decay hit the cell local m...
TH2F * fhMCPi0DecayPhotonHitHighLMOverlap
! E vs NLM when cluster originated in pi0 merging and MC photon decay hit the cell local maxima...
TH2F * fhMCPi0DecayPhotonAdjOtherLMDiffELM1vsELM1[3]
! E vs Ephoton-Esplit when cluster originated in pi0 merging and MC photon decay do not hit the adjac...
TH2F * fhNCellEtaNLocMax2[7][2]
! n cells for selected around eta, N Local Maxima = 2
TH2F * fhMCEEpriOverlapNMatch[3][7]
! E reco vs primary for different MC origin, N other MC particles contributes, charged cluster ...
void FillEBinHistograms(Int_t ebin, Int_t nMax, Int_t mcindex, Float_t splitFrac, Float_t mass, Float_t asym, Float_t l0)
Fill some histograms integrating in few energy bins.
TH2F * fhMCPi0DecayPhotonAdjOtherLMDiffELM2vsELM2[3]
! E vs Ephoton-Esplit when cluster originated in pi0 merging and MC photon decay do not hit the adjac...
TH2F * fhSplitClusterENLocMax[7][2]
! Number of maxima in cluster vs E of splitted clusters, 1-6 for different MC particle types ...
void FillHistograms2(Float_t en, Float_t eprim, Float_t e1, Float_t e2, Int_t nMax, Float_t mass, Float_t l0, Bool_t matched, Int_t mcindex)
Fill histograms for clusters passing the first M02 selection.
TH2F * fhAsyPi0NLocMax2[7][2]
! Asy for Mass around pi0, N Local Maxima = 2
TH2F * fhMassAsyCutNLocMax2
! Mass(E) asym selection, not matched, Mass of split clusters, NLM = 1
TH2F * fhMCPi0DecayPhotonHitHighLMDiffELM1[3]
! E vs Ephoton-Esplit cluster when cluster originated in pi0 merging and MC photon decay hit the cell...
TH2F * fhNLocMaxIdPi0[7][2]
! Number of maxima in cluster vs E, 1-6 for different MC particle types, after pi0 selection ...
TH2F * fhMCPi0MassM02OverlapN[3][4]
! MC Pi0 M02 vs Mass for different MC origin, N other MC particles contributes, neutral cluster...
TH2F * fhLM2NLocMaxIdPi0[7][2]
! Split cluster 2 E distribution vs Number of maxima in cluster vs E, 1-6 for different MC particle t...
TH2F * fhLM2NLocMaxM02Cut[7][2]
! Split cluster 2 E distribution vs Number of maxima in cluster vs E, 1-6 for different MC particle t...
TH2F * fhTrackMatchedDEtaNLocMax1Neg[7]
! Eta distance between track and cluster vs cluster E, 1 local maximum
TLorentzVector fMCDaughMom2
! Primary momentum
TH2F * fhAnglePairPrimPi0vsRecoNLocMaxN
! Pair opening angle pi0 generated vs reconstructed
TLorentzVector fMCDaughMom1
! Primary momentum
Float_t fSSECellCut[20]
List of cell min energy cuts to test.
TH2F * fhMassSplitEPi0NLocMax1[7][2]
! Split Inv Mass vs E1+E2, NLM=1, after pi0 selection, different MC particle types, track matching on/off
TH2F * fhAnglePairMassNLocMax1[7][2]
! Pair opening angle vs Mass for E > 7 GeV
TH2F * fhMCGenFracNLocMaxEbinMatched[7][4]
! NLM vs E generated particle / E reconstructed vs E reconstructed 1-6 for different MC particle type...
TH2F * fhMCENOverlapsMatch[3][7]
! E vs number of Overlaps in MC, charged cluster
TH2F * fhLM1NLocMax[7][2]
! Split cluster 1 E distribution vs Number of maxima in cluster vs E, 1-6 for different MC particle t...
void FillIdConvHistograms(Float_t en, Int_t nMax, Float_t asym, Float_t mass, Float_t l0, Bool_t matched, Int_t mcindex)
Fill histograms for clusters passing the photon selection.
TH2F * fhMCPi0DecayPhotonAdjHighLMDiffELM2vsELM2[3]
! E vs Ephoton-Esplit when cluster originated in pi0 merging and MC photon decay hit the adjacent cel...
TH2F * fhCentralityEtaNLocMax2
! Centrality for selected eta, N Local Maxima = 2
TH2F * fhMassDispPhiNLocMaxN[7][2]
! Mass of >2 cells local maxima, vs M02, for E > 8 GeV, 1-6 for different MC particle types ...
TH2F * fhPi0EPairDiffTimeNLMN
! E vs Pair of clusters time difference vs E, for selected pi0, NLM>2
TH2F * fhMCPi0DecayPhotonHitOtherLMOverlapDiffELM2vsELM2[3]
! E vs Ephoton-Esplit when cluster originated in pi0 merging and MC photon decay hit the cell local m...
TH2F * fhTrackMatchedDEtaNLocMaxNPos[7]
! Eta distance between track and cluster vs cluster E, more than 2 local maximum
TH2F * fhMCAsymM02NLocMaxNMCPi0Ebin[4]
! M02 vs decay asymmetry for N max in cluster > 2, for 4 energy bins
TH2F * fhAsymEnCutNLocMaxN
! Asym(E) E sub-cluster cut selection, not matched, Mass of split clusters, NLM > 2 ...
TH2F * fhMCPi0NoneNLMPairNoMCMatch
! E vs NLM when cluster originated in pi0 merging and a both no NLM corresponds to the photons ...
TH2F * fhMassM02NLocMaxNEbin[4]
! Mass of splitted clusters when >2 local max vs M02, 4 E bins, neutral clusters
TH2F * fhMassDispAsyNLocMax1[7][2]
! Mass of 2 highest energy cells when 1 local max, vs M02, for E > 8 GeV, 1-6 for different MC partic...
TH2F * fhCosThStarPi0NLocMaxN[7][2]
! cos(theta^star) vs E, after M02, asymmetry and pi0 mass cuts, NLM>2
TH2F * fhMCGenEvsSplitENLocMax2[7][2]
! E generated particle vs E1+E2 for N max in cluster = 2, 1-6 for different MC particle types ...
TH2F * fhM02Pi0NLocMax2[7][2]
! M02 for Mass around pi0, N Local Maxima = 2
TH2F * fhMassEtaNLocMaxN[7][2]
! Mass for selected around eta, N Local Maxima > 2
TH2F * fhMCGenFracNLocMax1[7][2]
! E generated particle / E reconstructed vs E reconstructed for N max in cluster = 1...
void FillMCOverlapHistograms(Float_t en, Float_t enprim, Int_t nc, Float_t mass, Float_t l0, Float_t asym, Float_t splitFrac, Int_t nlm, Int_t ebin, Bool_t matched, Int_t mcindex, Int_t noverlaps)
Fill histograms depending on number of overlaps.
TH2F * fhAsyPi0NLocMaxN[7][2]
! Asy for Mass around pi0, N Local Maxima > 2
TH2F * fhMCPi0DecayPhotonAdjHighLMDiffELM1[3]
! E vs Ephoton-Esplit cluster when cluster originated in pi0 merging and MC photon decay hit the adja...
TH2F * fhMCEEpriOverlap0IdPi0[3][7]
! E reco vs primary for different MC origin, no other MC particles contributes, neutral cluster...
TH2F * fhMCPi0DecayPhotonHitHighLMOverlapDiffELM2vsELM2[3]
! E vs Ephoton-Esplit cluster when cluster originated in pi0 merging and MC photon decay hit the cell...
TH2F * fhNCellMassEHighNLocMaxNMCPi0
! n cells in cluster vs mass for high energy clusters, for N max in cluster > 2, for MC pi0 ...
TH2F * fhPi0EtaPhiNLocMaxN
! Eta vs Phi of pi0's with N Local Maxima > N, E > 8 GeV
TH2F * fhMCPi0DecayPhotonHitOtherLMOverlapDiffELM1vsELM1[3]
! E vs Ephoton-Esplit when cluster originated in pi0 merging and MC photon decay hit the cell local m...
TH2F * fhEventPlaneEtaNLocMax1
! Event plane for selected eta, N Local Maxima = 1
TH2F * fhMCPi0LowNLMPairOverlap
! E vs NLM when cluster originated in pi0 merging and a pair of local maxima except highest energy co...
TH2F * fhMCPi0DecayPhotonAdjOtherLMDiffELM2[3]
! E vs Ephoton-Esplit when cluster originated in pi0 merging and MC photon decay do not hit the adjac...
TH2F * fhM02NLocMaxDiffCutPi0[5][5][3][2]
! M02 for 3 kinds of number of maxima for different values of min Loc Max value and min difference be...
TH2F * fhMassNLocMaxDiffCutPi0[5][5][3][2]
! M02 for 3 kinds of number of maxima for different values of min Loc Max value and min difference be...
TH2F * fhMCESplitEFracOverlap0[3][7]
! E vs sum of splitted cluster energy / cluster energy for different MC origin, no other MC particles...
TH2F * fhAsyMCGenRecoNLocMax1EbinPi0[4]
! Generated vs reconstructed asymmetry of splitted clusters from pi0 when 1 local max...
Float_t fNLMMinE[5]
List of local maxima min energy.
TH2F * fhM02MCGenFracNLocMax1Ebin[7][4]
! M02 vs E generated particle / E reconstructed vs E reconstructed for N max in cluster = 1...
Base class for CaloTrackCorr analysis algorithms.
TH2F * fhAnglePairPrimPi0RecoNLocMaxN
! Pair opening angle pi0 generated/reconstructed vs E
TH2F * fhAsyPi0NLocMax1[7][2]
! Asy for Mass around pi0, N Local Maxima = 1
void FillTrackMatchingHistograms(AliVCluster *cluster, Int_t nMax, Int_t mcindex)
Fill histograms related to cluster-track matching.
TH2F * fhMCEEpriOverlapNIdPi0[3][7]
! E reco vs primary for different MC origin, 1 other MC particles contributes, neutral cluster...
TH2F * fhPi0CellE[3]
! pi0's energy vs cluster cell energy with NLM = 1, = 2, > 2
TLorentzVector fSubClusterMom2
! Sub-Cluster momentum
TH2F * fhAsyMCGenRecoNLocMax2EbinPi0[4]
! Generated vs reconstructed asymmetry of splitted clusters from pi0 when 2 local max...
TH2F * fhArmPi0NLocMax1[7][4]
! Armenteros for selected pi0, N Local Maxima = 1
TH2F * fhEventPlaneEtaNLocMax2
! Event plane for selected eta, N Local Maxima = 2
TH2F * fhTrackMatchedDEtaNLocMax1Pos[7]
! Eta distance between track and cluster vs cluster E, 1 local maximum
TH2F * fhAnglePairOverM02NLocMax2[7][2]
! Pair opening angle / m02 vs E, NLM=2
TH2F * fhLM1NLocMaxM02Cut[7][2]
! Split cluster 1 E distribution vs Number of maxima in cluster vs E, 1-6 for different MC particle t...
TH2F * fhAnglePairPi0NLocMaxN[7][2]
! Pair opening angle vs E
TH2F * fhNLocMaxDiffCut[5][5][2]
! Number of maxima for different values of min Loc Max value and min difference between cells...
void GetMCIndex(AliVCluster *cluster, Int_t &mcindex, Int_t &mcTag)
Assign MC index depending on MC bit set, to be used in MC histograms arrays.
void RecalculateClusterShowerShapeParametersWithCellCut(const AliEMCALGeometry *geom, AliVCaloCells *cells, AliVCluster *cluster, Float_t &l0, Float_t &l1, Float_t &disp, Float_t &dEta, Float_t &dPhi, Float_t &sEta, Float_t &sPhi, Float_t &sEtaPhi, Float_t eCellMin=0.)
TH2F * fhM02AsyCutNLocMax2
! M02(E) asym selection, not matched, M02, NLM = 2
TH2F * fhMCPi0DecayPhotonAdjHighLMOverlapDiffELM1[3]
! E vs Ephoton-Esplit when cluster originated in pi0 merging and MC photon decay hit the adjacent cel...
TH2F * fhMassEtaNLocMax2[7][2]
! Mass for selected around eta, N Local Maxima = 2
Bool_t fFillHighMultHisto
Fill centrality/event plane histograms.
TH2F * fhAnglePairM02NLocMax1[7][2]
! Pair opening angle vs M02 for E > 7 GeV
TH2F * fhMassMCGenFracNLocMax1Ebin[7][4]
! Mass vs E generated particle / E reconstructed vs E reconstructed for N max in cluster = 1...
TH2F * fhMassSplitENLocMax1[7][2]
! Split Inv Mass vs E1+E2, NLM=1, different MC particle types, track matching on/off ...
TH2F * fhAnglePairPrimPi0RecoNLocMax2
! Pair opening angle pi0 generated/reconstructed vs E
TH2F * fhTrackMatchedDPhiNLocMax1[7]
! Phi distance between track and cluster vs cluster E, 1 local maximum
TH2F * fhArmAfterCutsNLocMax2[7][4]
! Armenteros after M02, asymmetry cuts for pi0, N Local Maxima = 2
TH2F * fhPi0EtaPhiNLocMax2
! Eta vs Phi of pi0's with N Local Maxima = 2, E > 8 GeV
TH2F * fhTrackMatchedDPhiNLocMaxNPos[7]
! Phi distance between track and cluster vs cluster E, more than 2 local maximum
TH2F * fhTrackMatchedDEtaNLocMaxN[7]
! Eta distance between track and cluster vs cluster E, more than 2 local maximum
TH2F * fhMCPi0DecayPhotonHitOtherLMOverlapMass[3]
! E vs Mass when cluster originated in pi0 merging and MC photon decay hit the cell local maximas...
TH2F * fhMCPi0MassM02Overlap0Match[3][4]
! MC Pi0 M02 vs Mass for different MC origin, no other MC particles contributes, charged cluster...
TH2F * fhMCPi0DecayPhotonHitHighLMDiffELM2vsELM2[3]
! E vs Ephoton-Esplit when cluster originated in pi0 merging and MC photon decay hit the cell local m...
void FillHistograms1(Float_t en, Float_t e1, Float_t e2, Int_t nMax, Float_t mass, Float_t l0, Float_t eta, Float_t phi, Bool_t matched, Int_t mcindex)
Fill histograms for clusters before any selection after spliting.
TH2F * fhNCellMassELowNLocMaxNMCPi0
! n cells in cluster vs mass for low energy clusters, for N max in cluster > 2, for MC pi0 ...
TH2F * fhMCPi0DecayPhotonAdjOtherLMOverlapDiffELM2vsELM2[3]
! E vs Ephoton-Esplit when cluster originated in pi0 merging and MC photon decay do not hit the adjac...
TH2F * fhPi0CellLogEFrac[3]
! pi0's energy vs cluster log cell energy fraction with NLM = 1, = 2, > 2
TH2F * fhAsyConNLocMaxN[7][2]
! Asy for Mass around close to 0, N Local Maxima > 2
TH2F * fhMCPi0AnyNLMPairNoMCMatch
! E vs NLM when cluster originated in pi0 merging and a both highest energy pairs and other pairs cor...
TH2F * fhAsymEnCutNLocMax2
! Asym(E) E sub-cluster cut selection, not matched, Mass of split clusters, NLM = 1 ...
TH2F * fhTrackMatchedDEtaNLocMaxNNeg[7]
! Eta distance between track and cluster vs cluster E, more than 2 local maximum
TH2F * fhAnglePairAfterCutsNLocMax2[7][2]
! Pair opening angle vs E
TH2F * fhTrackMatchedDEtaNLocMax2Pos[7]
! Eta distance between track and cluster vs cluster E, 2 local maximum
TH2F * fhMCPi0NoneNLMPairOverlap
! E vs NLM when cluster originated in pi0 merging and a both no NLM corresponds to the photons...
TH2F * fhMCPi0MassM02OverlapNMatch[3][4]
! MC Pi0 M02 vs Mass for different MC origin, N other MC particles contributes, charged cluster...
TH2F * fhNLocMaxM02Cut[7][2]
! Number of maxima in cluster vs E, 1-6 for different MC particle types, after SS cut ...
TH2F * fhMCPi0HighNLMPairOverlap
! E vs NLM when cluster originated in pi0 merging and highest energy local maxima correspond to 2 pho...
TH2F * fhCentralityPi0NLocMax2
! Centrality for selected pi0, N Local Maxima = 2
TH2F * fhMassAsyNLocMax2Ebin[4]
! Mass of Mass of splitted clusters when 2 local max vs asymmetry, 4 E bins, neutral clusters ...
TLorentzVector fPrimaryMom
! Primary momentum
Bool_t fCheckSplitDistToBad
Check the distance to bad channel and to EMCal borders of split clusters.
TH2F * fhMassSplitEAfterCutsNLocMaxN[7][2]
! Split Inv Mass vs E1+E2, NLM>2, after M02, asymmetry cuts, different MC particle types...
TH2F * fhMassSplitEFractionNLocMaxNEbin[7][4]
! Mass vs sum of splitted cluster energy / cluster energy for N max in cluster > 2, 1-6 for different MC particle types, not track matched
TH2F * fhMCENCellOverlap1[3][7]
! E vs sum of splitted cluster energy asymmetry for different MC origin, 1 other MC particles contrib...
TH2F * fhMCPi0LowNLMPairNoMCMatch
! E vs NLM when cluster originated in pi0 merging and a pair of local maxima except highest energy co...
TH2F * fhPi0CellEMax2Frac[3]
! pi0's energy vs energy fraction of second LM and cluster cell energy with NLM = 1...
TH2F * fhSplitEFractionNLocMax1[7][2]
! sum of splitted cluster energy / cluster energy for N Local Maxima = 1
TH2F * fhSplitEFractionAfterCutsNLocMaxN[7][2]
! sum of splitted cluster energy / cluster energy for N Local Maxima > 2, after M02 and asymmetry cut...
TH2F * fhMCPi0DecayPhotonHitHighLMOverlapDiffELM1vsELM1[3]
! E vs Ephoton-Esplit cluster when cluster originated in pi0 merging and MC photon decay hit the cell...
void CheckLocalMaximaMCOrigin(AliVCluster *cluster, Int_t mcindex, Int_t noverlaps, Float_t e1, Float_t e2, Float_t mass)
TH2F * fhTrackMatchedDPhiNLocMax2Pos[7]
! Phi distance between track and cluster vs cluster E, 2 local maximum
TH2F * fhAnglePairNLocMaxN[7][2]
! Pair opening angle vs E
TH2F * fhMCPi0DecayPhotonAdjacentOverlap
! E vs NLM when cluster originated in pi0 merging and MC photon decay hit adjacen cells...
TH2F * fhMassDispPhiNLocMax1Ebin[4]
! Mass of 2 highest energy cells when 1 local max, vs M02, 4 E bins, neutral clusters ...
TH2F * fhM02EnCutNLocMaxN
! M02(E) E sub-cluster cut selection, not matched, Mass of split clusters, NLM > 2 ...
TH2F * fhTrackMatchedDPhiNLocMaxN[7]
! Phi distance between track and cluster vs cluster E, more than 2 local maximum
TH2F * fhMCGenFracAfterCutsNLocMax1MCPi0
! E generated particle / E reconstructed vs E reconstructed for N max in cluster = 1...
TH2F * fhMassBadDistClose[3]
! Split mass of clusters with second LM close to bad channel
TH2F * fhMCEMassOverlap0[3][7]
! E vs Mass for different MC origin, no other MC particles contributes, neutral cluster ...
TH2F * fhNCellPi0NLocMax2[7][2]
! n cells for selected around pi0, N Local Maxima = 2
TH2F * fhAsymNLocMaxN[7][2]
! Asymmetry of >2 cells local maxima vs E, 1-6 for different MC particle types
TH2F * fhMCGenFracNLocMaxNNoOverlap[7][2]
! E generated particle / E reconstructed vs E reconstructed for N max in cluster > 2...
TH2F * fhNCellM02ELowNLocMaxNMCPi0
! n cells in cluster vs m02 for low energy clusters, for N max in cluster > 2, for MC pi0 ...
TH2F * fhMassDispPhiNLocMax2Ebin[4]
! Mass of 2 cells local maxima, vs M02, 4 E bins, neutral clusters
TH2F * fhSplitEFractionNLocMax2[7][2]
! sum of splitted cluster energy / cluster energy for N Local Maxima = 2
Float_t fHistoECut
Fixed E cut for some histograms.
mcTypes
Enumerate indeces for MC histograms depending on the particle ID that generated the cluster...
TH2F * fhMCPi0HighNLMPair
! E vs NLM when cluster originated in pi0 merging and highest energy local maxima correspond to 2 pho...
TH2F * fhArmAfterCutsNLocMax1[7][4]
! Armenteros after M02, asymmetry cuts for pi0, N Local Maxima = 1
TH2F * fhCosThStarAfterCutsNLocMaxN[7][2]
! cos(theta^star) vs E, after M02, asymmetry cuts, NLM>2
TH2F * fhMassNLocMax1[7][2]
! Split Inv Mass vs cluster E, NLM=1, different MC particle types, track matching on/off ...
TH2F * fhMCEEpriOverlap1IdPi0[3][7]
! E reco vs primary for different MC origin, 1 other MC particles contributes, neutral cluster...
TH2F * fhLM2NLocMax[7][2]
! Split cluster 2 E distribution vs Number of maxima in cluster vs E, 1-6 for different MC particle t...
TH2F * fhPi0EtaPhiNLocMax1
! Eta vs Phi of pi0's with N Local Maxima = 1, E > 8 GeV
void FillThetaStarHistograms(Int_t nMax, Bool_t matched, Int_t mcindex, Float_t pi0E, Float_t m02, Int_t pid)
TH2F * fhMCPi0DecayPhotonHitNoLMOverlapMass[3]
! E vs Mass when cluster originated in pi0 merging and MC photon decay do not hit the cell local maxi...
TH2F * fhM02EtaNLocMax2[7][2]
! M02 for Mass around eta, N Local Maxima = 2
TH2F * fhNCellM02EHighNLocMax1MCPi0
! n cells in cluster vs m02 for high energy clusters, for N max in cluster = 1, for MC pi0 ...
TH2F * fhMassConNLocMax1[7][2]
! Mass for selected around close to 0, N Local Maxima = 1
TH2F * fhMCEAsymOverlapNMatch[3][7]
! E vs sum of splitted cluster energy asymmetry for different MC origin, N other MC particles contrib...
TH2F * fhMCESplitEFracOverlap1[3][7]
! E vs sum of splitted cluster energy / cluster energy for different MC origin, 1 other MC particles ...
TH2F * fhMCESplitEFracOverlapN[3][7]
! E vs sum of splitted cluster energy / cluster energy for different MC origin, N other MC particles ...
TH2F * fhMCEEpriOverlap0Match[3][7]
! E reco vs primary for different MC origin, no other MC particles contributes, charged cluster ...
TH2F * fhSplitEFractionAfterCutsNLocMax1[7][2]
! sum of splitted cluster energy / cluster energy for N Local Maxima = 1, after M02 and asymmetry cut...
TH2F * fhNCellNLocMax2[7][2]
! n cells in cluster vs E for N max in cluster = 2, 1-6 for different MC particle types ...
Bool_t fFillNLMDiffCutHisto
Fill NLM histograms for different settings.
TH2F * fhM02BadDistClose[3]
! M02 of clusters with second LM close to bad channel
TLorentzVector fSubClusterMom1
! Sub-Cluster momentum
TH2F * fhAnglePairOverM02NLocMaxNOverlap0[7][2]
! Pair opening angle / m02 vs E, NLM=N
AliAnaInsideClusterInvariantMass & operator=(const AliAnaInsideClusterInvariantMass &split)
Assignment operator not implemented.
TH2F * fhMassEnCutNLocMax1
! Mass(E) E sub-cluster cut selection, not matched, Mass of split clusters, NLM = 1 ...
TH2F * fhEtaEtaPhiNLocMaxN
! Eta vs Phi of eta's with N Local Maxima > N, E > 8 GeV
TH2F * fhEventPlanePi0NLocMax2
! Event plane for selected pi0, N Local Maxima = 2
TH2F * fhLM1NLocMaxIdPi0[7][2]
! Split cluster 1 E distribution vs Number of maxima in cluster vs E, 1-6 for different MC particle t...
TH2F * fhM02OnBorder[3]
! M02 of clusters with second LM close to EMCAL border
TH2F * fhMassOnBorder[3]
! Split mass of clusters with second LM on EMCAL border
TH2F * fhMassSplitECutNLocMax1
! 85% of split energy, not matched, Mass of split clusters, NLM = 1
TH2F * fhMCEMassOverlap1Match[3][7]
! E vs Mass for different MC origin, 1 other MC particles contributes, charged cluster ...
TH2F * fhMCPi0DecayPhotonHitNoLMMass[3]
! E vs Mass when cluster originated in pi0 merging and MC photon decay do not hit the cell local maxi...
TH2F * fhMCPi0MassM02Overlap1[3][4]
! MC Pi0 M02 vs Mass for different MC origin, 1 other MC particles contributes, neutral cluster...
TH2F * fhTrackMatchedDEtaNLocMax2[7]
! Eta distance between track and cluster vs cluster E, 2 local maximum
TH2F * fhAnglePairPi0NLocMax2[7][2]
! Pair opening angle vs E
TH2F * fhMCAsymM02NLocMax1MCPi0Ebin[4]
! M02 vs decay asymmetry for N max in cluster = 1, for 4 energy bins
TH2F * fhAnglePairAfterCutsNLocMax1[7][2]
! Pair opening angle vs E
TH2F * fhMCPi0DecayPhotonAdjHighLMMass[3]
! E vs Mass when cluster originated in pi0 merging and MC photon decay hit the adjacent cell local ma...
TH2F * fhMCPi0DecayPhotonHitNoLMOverlap
! E vs NLM when cluster originated in pi0 merging and MC photon decay do not hit the cell local maxim...
Bool_t fFillIdEtaHisto
Fill histograms for clusters identified as Eta.
TH2F * fhMCEMassOverlap1[3][7]
! E vs Mass for different MC origin, 1 other MC particles contributes, neutral cluster ...
TH2F * fhMCPi0AnyNLMPairOverlap
! E vs NLM when cluster originated in pi0 merging and a both highest energy pairs and other pairs cor...
TH2F * fhAsymM02CutNLocMax2
! Asym(E) M02 selection, not matched, energy asymmetry of split clusters, NLM = 2 ...
TH2F * fhMCGenSplitEFracNLocMax1[7][2]
! E generated particle / E1+E2 reconstructed vs E reconstructed for N max in cluster = 1...
TH2F * fhMCEM02Overlap0Match[3][7]
! E vs M02 for different MC origin, no other MC particles contributes, charged cluster ...
TH2F * fhMCPi0DecayPhotonHitOtherLMOverlapDiffELM1[3]
! E vs Ephoton-Esplit when cluster originated in pi0 merging and MC photon decay hit the cell local m...
TH2F * fhCosThStarAfterCutsNLocMax2[7][2]
! cos(theta^star) vs E, after M02, asymmetry cuts, NLM=2
TLorentzVector fSubClusterMomSum
! Sub-Cluster momentum sum, armenteros
TH2F * fhMCPi0DecayPhotonAdjacentOverlapMass[3]
! E vs Mass when cluster originated in pi0 merging and MC photon decay hit adjacen cells...
TH2F * fhMassPi0NLocMax1[7][2]
! Mass for selected pi0, N Local Maxima = 1
Bool_t fFillSSExtraHisto
Fill shower shape extra histos.
Bool_t fFillTMResidualHisto
Fill track matching histos, residuals.
TH2F * fhMassDispAsyNLocMax2[7][2]
! Mass of 2 cells local maxima, vs M02, for E > 8 GeV, 1-6 for different MC particle types ...
TH2F * fhAsyMCGenRecoDiffMCPi0Conv[3]
! reconstructed-generated asymmetry of splitted clusters vs E from converted pi0, for 3 NLM cases ...
TH2F * fhMassDispEtaNLocMax2[7][2]
! Mass of 2 cells local maxima, vs M02, for E > 8 GeV, 1-6 for different MC particle types ...
TH2F * fhM02EnCutNLocMax1
! M02(E) E sub-cluster cut selection, not matched, Mass of split clusters, NLM = 1 ...
TH2F * fhPi0EPairDiffTimeNLM1
! E vs Pair of clusters time difference vs E, for selected pi0, NLM=1
TH2F * fhMCPi0HighNLMPairNoMCMatchOverlap
! E vs NLM when cluster originated in pi0 merging and highest energy local maxima correspond to 2 pho...
TH2F * fhAnglePairPrimPi0OverM02NLocMaxN
! Pair opening angle / m02 vs E, NLM=N, prim pi0
TH2F * fhM02NLocMaxN[7][2]
! M02 vs E for N max in cluster > 2, 1-6 for different MC particle types
TH2F * fhMCAsymM02NLocMax2MCPi0Ebin[4]
! M02 vs decay asymmetry for N max in cluster = 2, for 4 energy bins
TH2F * fhM02ECellCutPi0[3][20]
! M02 for selected pi0 with different cut on cell energy, with NLM = 1, = 2, > 2
TH2F * fhMCGenFracAfterCutsNLocMaxNMCPi0
! E generated particle / E reconstructed vs E reconstructed for N max in cluster > 2...
TH2F * fhMassAsyNLocMax1Ebin[4]
! Mass of Mass of splitted clusters when 1 local max vs asymmetry, 4 E bins, neutral clusters ...
TH2F * fhNCellMassELowNLocMax1MCPi0
! n cells in cluster vs mass for low energy clusters, for N max in cluster = 1, for MC pi0 ...
void FillNCellHistograms(Int_t ncells, Float_t energy, Int_t nMax, Bool_t matched, Int_t mcindex, Float_t mass, Float_t l0)
Fill optional histograms depending on the number of cells of the cluster.
Bool_t fFillNCellHisto
Fill n cells in cluster dependent histograms.
TH2F * fhNCellPi0NLocMax1[7][2]
! n cells for selected around pi0, N Local Maxima = 1
TH2F * fhMCEAsymOverlap1Match[3][7]
! E vs sum of splitted cluster energy asymmetry for different MC origin, 1 other MC particles contrib...
Bool_t fFillSSWeightHisto
Fill weigth histograms.
TH2F * fhMassAfterCutsNLocMax2[7][2]
! Mass after M02, asymmetry cuts for MC part, N Local Maxima = 2
TH2F * fhNCellNLocMax1[7][2]
! n cells in cluster vs E for N max in cluster = 1, 1-6 for different MC particle types ...
TH2F * fhMCPi0DecayPhotonHitOtherLMDiffELM1vsELM1[3]
! E vs Ephoton-Esplit when cluster originated in pi0 merging and MC photon decay hit the cell local m...
TH2F * fhCentralityPi0NLocMax1
! Centrality for selected pi0, N Local Maxima = 1
TH2F * fhTrackMatchedDEtaNLocMax2Neg[7]
! Eta distance between track and cluster vs cluster E, 2 local maximum
TH2F * fhAnglePairPrimPi0vsRecoNLocMax2
! Pair opening angle pi0 generated vs reconstructed
TH2F * fhM02EtaNLocMaxN[7][2]
! M02 for Mass around eta, N Local Maxima > 2
TH2F * fhM02AsyCutNLocMax1
! M02(E) asym selection, not matched, M02, NLM = 1
TH2F * fhMassDispEtaNLocMax2Ebin[4]
! Mass of 2 cells local maxima, vs M02, 4 E bins, neutral clusters
TH2F * fhPi0CellEMaxFrac[3]
! pi0's energy vs energy fraction of main LM and cluster cell energy with NLM = 1, = 2, > 2
TH2F * fhM02Pi0NLocMaxN[7][2]
! M02 for Mass around pi0, N Local Maxima > 2
TH2F * fhMCGenFracNLocMaxN[7][2]
! E generated particle / E reconstructed vs E reconstructed for N max in cluster > 2...
TH2F * fhMCPi0DecayPhotonAdjHighLMOverlapDiffELM1vsELM1[3]
! E vs Ephoton-Esplit when cluster originated in pi0 merging and MC photon decay hit the adjacent cel...
TH2F * fhMCEOverlapType
! What particles overlap with pi0, neutral clusters
TH2F * fhNCellMassEHighNLocMax2MCPi0
! n cells in cluster vs mass for high energy clusters, for N max in cluster = 2, for MC pi0 ...
TH2F * fhMassSplitEFractionNLocMax2Ebin[7][4]
! Mass vs sum of splitted cluster energy / cluster energy for N max in cluster = 2, 1-6 for different MC particle types, not track matched
TH2F * fhMCEAsymOverlap1[3][7]
! E vs sum of splitted cluster energy asymmetry for different MC origin, 1 other MC particles contrib...
TH2F * fhMCPi0DecayPhotonHitOtherLMMass[3]
! E vs Mass when cluster originated in pi0 merging and MC photon decay hit the cell local maximas...
TH2F * fhMCPi0DecayPhotonAdjOtherLMOverlapMass[3]
! E vs Mass when cluster originated in pi0 merging and MC photon decay do not hit the adjacent cell l...
TH2F * fhMCPi0NoneNLMPair
! E vs NLM when cluster originated in pi0 merging and a both no NLM corresponds to the photons ...
TH2F * fhAsyMCGenRecoNLocMaxNEbinPi0[4]
! Generated vs reconstructed asymmetry of splitted clusters from pi0 when >2 local max...
TH2F * fhMCEM02Overlap1Match[3][7]
! E vs M02 for different MC origin, 1 other MC particles contributes, charged cluster ...
TH2F * fhMCEM02Overlap0[3][7]
! E vs M02 for different MC origin, no other MC particles contributes, neutral cluster ...
TH2F * fhArmAfterCutsNLocMaxN[7][4]
! Armenteros after M02, asymmetry cuts for pi0, N Local Maxima > 2
TH2F * fhM02NLocMax1[7][2]
! M02 vs E for N max in cluster = 1, 1-6 for different MC particle types
TH2F * fhMCEMassOverlapN[3][7]
! E vs Mass for different MC origin, N other MC particles contributes, neutral cluster ...
TH2F * fhAsymM02CutNLocMaxN
! Asym(E) M02 selection, not matched, energy asymmetry of split clusters, NLM > 2 ...
TH2F * fhMassDispAsyNLocMax1Ebin[4]
! Mass of 2 highest energy cells when 1 local max, vs M02, 4 E bins, neutral clusters ...
TH2F * fhCosThStarNLocMaxN[7][2]
! cos(theta^star) vs E, NLM>2
TH2F * fhEventPlanePi0NLocMax1
! Event plane for selected pi0, N Local Maxima = 1
TH2F * fhMassM02NLocMax1Ebin[4]
! Mass of splitted clusters when 1 local max vs M02, 4 E bins, neutral clusters
TH2F * fhMCGenSplitEFracNLocMax2[7][2]
! E generated particle / E1+E2 reconstructed vs E reconstructed for N max in cluster = 2...
TH2F * fhMassPi0NLocMaxN[7][2]
! Mass for selected around pi0, N Local Maxima > 2
TH2F * fhMCPi0DecayPhotonHitOtherLMOverlap
! E vs NLM when cluster originated in pi0 merging and MC photon decay hit the cell local maximas...
void FillSSExtraHistograms(AliVCluster *cluster, Int_t nMax, Bool_t matched, Int_t mcindex, Float_t mass, Int_t ebin)
Fill optional histograms with more SS parameters.
TH2F * fhMassSplitENLocMax2[7][2]
! Split Inv Mass vs E1+E2, NLM=2, different MC particle types, track matching on/off ...
Bool_t fFillMCOverlapHisto
Fill MC particles overlap histograms.
TH2F * fhNCellM02EHighNLocMaxNMCPi0
! n cells in cluster vs m02 for high energy clusters, for N max in cluster > 2, for MC pi0 ...
TH2F * fhAsyEtaNLocMax2[7][2]
! Asy for Mass around eta, N Local Maxima = 2
TH2F * fhAsymNLocMax2[7][2]
! Asymmetry of 2 cells local maxima vs E, 1-6 for different MC particle types
void InitParameters()
Initialize the parameters of the analysis with default values.
TH2F * fhMCEEpriOverlap0[3][7]
! E reco vs primary for different MC origin, no other MC particles contributes, neutral cluster ...
TH2F * fhMassDispPhiNLocMaxNEbin[4]
! Mass of >2 cells local maxima, vs M02, 4 E bins, neutral clusters
TH2F * fhMCENCellOverlap0Match[3][7]
! E vs sum of splitted cluster energy asymmetry for different MC origin, no other MC particles contri...
TH2F * fhMCPi0DecayPhotonHitHighLMOverlapDiffELM2[3]
! E vs Ephoton-Esplit cluster when cluster originated in pi0 merging and MC photon decay hit the cell...
TH2F * fhMassDispAsyNLocMaxN[7][2]
! Mass of >2 cells local maxima, vs M02, for E > 8 GeV, 1-6 for different MC particle types ...
Int_t fMinNCells
Study clusters with ncells larger than cut.
void FillNLMDiffCutHistograms(AliVCluster *cluster, AliVCaloCells *cells, Bool_t matched)
Calculate NLM for different settings and recalculate splitting.
TH2F * fhNCellM02ELowNLocMax2MCPi0
! n cells in cluster vs m02 for low energy clusters, for N max in cluster = 2, for MC pi0 ...
TH2F * fhAsyConNLocMax2[7][2]
! Asy for Mass around close to 0, N Local Maxima = 2
TH2F * fhMassSplitENLocMaxN[7][2]
! Split Inv Mass vs E1+E2, NLM>2, different MC particle types, track matching on/off ...
TH2F * fhEtaEPairDiffTimeNLM1
! E vs Pair of clusters time difference vs E, for selected eta, NLM=1
TH2F * fhAnglePairM02NLocMax2[7][2]
! Pair opening angle vs M02 for E > 7 GeV
TH2F * fhMCGenFracNLocMax1NoOverlap[7][2]
! E generated particle / E reconstructed vs E reconstructed for N max in cluster = 1...
TH2F * fhMassAfterCutsNLocMax1[7][2]
! Mass after M02, asymmetry cuts for MC part, N Local Maxima = 1
TH2F * fhEtaEtaPhiNLocMax1
! Eta vs Phi of eta's with N Local Maxima = 1, E > 8 GeV
TH2F * fhMCPi0DecayPhotonAdjOtherLMMass[3]
! E vs Mass when cluster originated in pi0 merging and MC photon decay do not hit the adjacent cell l...
TH2F * fhClusterEtaPhiNLocMax1
! Eta vs Phi of clusters with N Local Maxima = 1, E > 8 GeV
TH2F * fhM02EtaNLocMax1[7][2]
! M02 for Mass around eta, N Local Maxima = 1
TH2F * fhAnglePairAfterCutsNLocMaxN[7][2]
! Pair opening angle vs E
TH2F * fhAnglePairOverM02NLocMax2Overlap0[7][2]
! Pair opening angle / m02 vs E, NLM=2
TH2F * fhMassPi0NLocMax2[7][2]
! Mass for selected around pi0, N Local Maxima = 2
TH2F * fhMCPi0DecayPhotonHitHighLMDiffELM1vsELM1[3]
! E vs Ephoton-Esplit cluster when cluster originated in pi0 merging and MC photon decay hit the cell...
void FillSSWeightHistograms(AliVCluster *cluster, Int_t nlm, Int_t absId1, Int_t absId2)
Recalculate shower shape for different weights predifined in an array and fill histograms.
TH2F * fhMCPi0DecayPhotonAdjOtherLMOverlapDiffELM1[3]
! E vs Ephoton-Esplit when cluster originated in pi0 merging and MC photon decay do not hit the adjac...
TH2F * fhMCESplitEFracOverlapNMatch[3][7]
! E vs sum of splitted cluster energy / cluster energy for different MC origin, N other MC particles ...
TH2F * fhMCPi0DecayPhotonAdjOtherLMDiffELM1[3]
! E vs Ephoton-Esplit when cluster originated in pi0 merging and MC photon decay do not hit the adjac...
TH2F * fhMCGenSplitEFracAfterCutsNLocMax1MCPi0
! E generated particle / E1+E2 reconstructed vs E reconstructed for N max in cluster = 1...
TH2F * fhTrackMatchedDEtaNLocMax1[7]
! Eta distance between track and cluster vs cluster E, 1 local maximum
TH2F * fhMCEM02Overlap1[3][7]
! E vs M02 for different MC origin, 1 other MC particles contributes, neutral cluster ...
TH2F * fhM02MCGenFracNLocMaxNEbin[7][4]
! M02 vs E generated particle / E reconstructed vs E reconstructed for N max in cluster > 2...
TH2F * fhMassM02CutNLocMaxN
! Mass(E) M02 selection, not matched, Mass of split clusters, NLM > 2
void MakeAnalysisFillHistograms()
Main method. Search for pi0/eta in GetCalorimeter() with shower shape + split analysis.
TH2F * fhMCPi0DecayPhotonAdjHighLMDiffELM1vsELM1[3]
! E vs Ephoton-Esplit cluster when cluster originated in pi0 merging and MC photon decay hit the adja...
TH2F * fhMCEAsymOverlap0[3][7]
! E vs sum of splitted cluster energy asymmetry for different MC origin, no other MC particles contri...
TH2F * fhNCellPi0NLocMaxN[7][2]
! n cells for selected around pi0, N Local Maxima > 2
TH2F * fhMassDispPhiNLocMax1[7][2]
! Mass of 2 highest energy cells when 1 local max, vs M02, for E > 8 GeV, 1-6 for different MC partic...
TH2F * fhAsyEtaNLocMax1[7][2]
! Asy for Mass around eta, N Local Maxima = 1
TH2F * fhMCGenFracNLocMax2NoOverlap[7][2]
! E generated particle / E reconstructed vs E reconstructed for N max in cluster = 2...
TH2F * fhMCPi0DecayPhotonAdjOtherLMOverlapDiffELM2[3]
! E vs Ephoton-Esplit when cluster originated in pi0 merging and MC photon decay do not hit the adjac...
TH2F * fhMassDispEtaNLocMax1[7][2]
! Mass of 2 highest energy cells when 1 local max, vs M02, for E > 8 GeV, 1-6 for different MC partic...
TH2F * fhSplitEFractionvsAsyNLocMax1[2]
! sum of splitted cluster energy / cluster energy for N Local Maxima = 1 vs |A|
TH2F * fhTrackMatchedDPhiNLocMax1Pos[7]
! Phi distance between track and cluster vs cluster E, 1 local maximum
TH2F * fhMCPi0AnyNLMPairNoMCMatchOverlap
! E vs NLM when cluster originated in pi0 merging and a both highest energy pairs and other pairs cor...
TH2F * fhMCPi0HighNLMPairNoMCMatch
! E vs NLM when cluster originated in pi0 merging and highest energy local maxima correspond to 2 pho...
TH2F * fhMCEEpriOverlap1Match[3][7]
! E reco vs primary for different MC origin, 1 other MC particles contributes, charged cluster ...
Select clusters with large shower shape, split them and tag them as pi0/eta via invariant mass...
TH2F * fhNCellEtaNLocMaxN[7][2]
! n cells for selected around eta, N Local Maxima > 2
TH2F * fhMassNLocMaxN[7][2]
! Split Inv Mass vs cluster E, NLM>2, different MC particle types, track matching on/off ...
TH2F * fhMCGenEFracvsSplitEFracNLocMax2[7][2]
! E generated particle / E reconstructed vs E1+E2 reconstructed / E reconstructed for N max in cluste...
TH2F * fhAnglePairPrimPi0OverM02NLocMax2
! Pair opening angle / m02 vs E, NLM=2, prim pi0
TH2F * fhArmPi0NLocMaxN[7][4]
! Armenteros for selected pi0, N Local Maxima > 2
void FillIdEtaHistograms(Float_t en, Float_t e1, Float_t e2, Int_t nc, Int_t nMax, Float_t t12diff, Float_t mass, Float_t l0, Float_t eta, Float_t phi, Bool_t matched, Int_t mcindex)
Fill histograms for clusters passing the eta selection.
TH2F * fhCentralityEtaNLocMaxN
! Centrality for selected eta, N Local Maxima > 2
TH2F * fhMassDispEtaNLocMaxNEbin[4]
! Mass of >2 cells local maxima, vs M02, 4 E bins, neutral clusters
TH2F * fhMassM02CutNLocMax2
! Mass(E) M02 selection, not matched, Mass of split clusters, NLM = 1
TH2F * fhNCellEtaNLocMax1[7][2]
! n cells for selected around eta, N Local Maxima = 1
TH2F * fhAnglePairOverM02NLocMax1Overlap0[7][2]
! Pair opening angle / m02 vs E, NLM=1
TH2F * fhEventPlaneEtaNLocMaxN
! Event plane for selected eta, N Local Maxima > 2
TH2F * fhMCPi0DecayPhotonHitNoLM
! E vs NLM when cluster originated in pi0 merging and MC photon decay do not hit the cell local maxim...
TH2F * fhMCPi0LowNLMPairNoMCMatchOverlap
! E vs NLM when cluster originated in pi0 merging and a pair of local maxima except highest energy co...
TH2F * fhAsymEnCutNLocMax1
! Asym(E) E sub-cluster cut selection, not matched, Mass of split clusters, NLM = 1 ...
TH2F * fhM02MCGenFracNLocMax2Ebin[7][4]
! M02 vs E generated particle / E reconstructed vs E reconstructed for N max in cluster = 2...
Bool_t fFillThetaStarHisto
Fill cosThetaStar histos.
TH2F * fhM02EnCutNLocMax2
! M02(E) E sub-cluster cut selection, not matched, Mass of split clusters, NLM = 1 ...
TH2F * fhNLocMax[7][2]
! Number of maxima in cluster vs E, 1-6 for different MC particle types
TH2F * fhMassAsyNLocMaxNEbin[4]
! Mass of Mass of splitted clusters when >2 local max vs asymmetry, 4 E bins, neutral clusters ...
TH2F * fhMassSplitECutNLocMax2
! 85% of split energy, not matched, Mass of split clusters, NLM = 1
TH2F * fhMCGenSplitEFracNLocMaxN[7][2]
! E generated particle / E1+E2 reconstructed vs E reconstructed for N max in cluster > 2...
TH2F * fhMassMCGenFracNLocMaxNEbin[7][4]
! Mass vs E generated particle / E reconstructed vs E reconstructed for N max in cluster > 2...
TH2F * fhMassAfterCutsNLocMaxN[7][2]
! Mass after M02, asymmetry cuts for MC part, N Local Maxima > 2
TH2F * fhTrackMatchedDPhiNLocMax1Neg[7]
! Phi distance between track and cluster vs cluster E, 1 local maximum
void Print(const Option_t *opt) const
Print some relevant parameters set for the analysis.
TH2F * fhMassEtaNLocMax1[7][2]
! Mass for selected around eta, N Local Maxima = 1
TLorentzVector fSubClusterMomBoost
! Sub-Cluster momentum boosted, armenteros
TH2F * fhMassSplitEAfterCutsNLocMax2[7][2]
! Split Inv Mass vs E1+E2, NLM=2, after M02, asymmetry cuts, different MC particle types...
TH2F * fhCosThStarAfterCutsNLocMax1[7][2]
! cos(theta^star) vs E, after M02, asymmetry cuts, NLM=1
TH2F * fhMassDispAsyNLocMax2Ebin[4]
! Mass of 2 cells local maxima, vs M02, 4 E bins, neutral clusters
TH2F * fhMCPi0DecayPhotonHitHighLMMass[3]
! E vs Mass when cluster originated in pi0 merging and MC photon decay hit the cell local maxima ...
TH2F * fhClusterEtaPhiNLocMaxN
! Eta vs Phi of clusters with N Local Maxima > 2, E > 8 GeV
TH2F * fhMCPi0DecayPhotonHitOtherLMDiffELM2vsELM2[3]
! E vs Ephoton-Esplit when cluster originated in pi0 merging and MC photon decay hit the cell local m...
Int_t fSSECellCutN
Total number of cell min energy cuts to test.
TH2F * fhMassDispAsyNLocMaxNEbin[4]
! Mass of >2 cells local maxima, vs M02, 4 E bins, neutral clusters
Float_t fNLMMinDiff[5]
List of local maxima min difference cell energy.
TH2F * fhMCEMassOverlapNMatch[3][7]
! E vs Mass for different MC origin, N other MC particles contributes, charged cluster ...
TH2F * fhM02ConNLocMax2[7][2]
! M02 for Mass around close to 0, N Local Maxima = 2
TH2F * fhMassConNLocMax2[7][2]
! Mass for selected around close to 0, N Local Maxima = 2
TH2F * fhAnglePairPrimPi0OverM02NLocMax1
! Pair opening angle / m02 vs E, NLM=1, prim pi0
TH2F * fhAsyMCGenRecoDiffMCPi0[3]
! reconstructed-generated asymmetry of splitted clusters vs E from pi0, for 3 NLM cases ...
TH2F * fhEventPlanePi0NLocMaxN
! Event plane for selected pi0, N Local Maxima > 2
TH2F * fhMCPi0AnyNLMPair
! E vs NLM when cluster originated in pi0 merging and a both highest energy pairs and other pairs cor...
TH2F * fhNCellM02EHighNLocMax2MCPi0
! n cells in cluster vs m02 for high energy clusters, for N max in cluster = 2, for MC pi0 ...
TH2F * fhAnglePairPrimPi0RecoNLocMax1
! Pair opening angle pi0 generated/reconstructed vs E
TH2F * fhMCPi0LowNLMPair
! E vs NLM when cluster originated in pi0 merging and a pair of local maxima except highest energy co...
TH2F * fhMassMCGenFracNLocMax2Ebin[7][4]
! Mass vs E generated particle / E reconstructed vs E reconstructed for N max in cluster = 2...
TH2F * fhMCGenSplitEFracAfterCutsNLocMax2MCPi0
! E generated particle / E1+E2 reconstructed vs E reconstructed for N max in cluster = 2...
TH2F * fhAnglePairNLocMax2[7][2]
! Pair opening angle vs E
TH2F * fhSplitEFracEnCutNLocMax1
! Split E fraction (E) E sub-cluster cut selection, not matched, Mass of split clusters, NLM = 1
TH2F * fhM02AsyCutNLocMaxN
! M02(E) asym selection, not matched, M02, NLM > 2
TH2F * fhMCGenFracNLocMaxEbin[7][4]
! NLM vs E generated particle / E reconstructed vs E reconstructed 1-6 for different MC particle type...
TH2F * fhAsymM02CutNLocMax1
! Asym(E) M02 selection, not matched, energy asymmetry of split clusters, NLM = 1 ...
TH2F * fhMCENCellOverlapNMatch[3][7]
! E vs sum of splitted cluster energy asymmetry for different MC origin, N other MC particles contrib...
TH2F * fhMCEEpriOverlapN[3][7]
! E reco vs primary for different MC origin, N other MC particles contributes, neutral cluster ...
virtual ~AliAnaInsideClusterInvariantMass()
Virtual destructor. Not implemeted.
TH2F * fhPi0CellEMaxEMax2Frac[3]
! pi0's energy vs fraction of 2 main maxima energy with NLM = 1, = 2, > 2
TObjString * GetAnalysisCuts()
Save parameters used for analysis.
TH2F * fhMCGenSplitEFracAfterCutsNLocMaxNMCPi0
! E generated particle / E1+E2 reconstructed vs E reconstructed for N max in cluster > 2...
TH2F * fhMassEnCutNLocMaxN
! Mass(E) E sub-cluster cut selection, not matched, Mass of split clusters, NLM > 2 ...
TH2F * fhMCPi0DecayPhotonAdjacent
! E vs NLM when cluster originated in pi0 merging and MC photon decay hit adjacen cells...