PWGGA/CaloTrackCorrelations/AliAnaInsideClusterInvariantMass.h

   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 //_________________________________________________________________________
   7 //
   8 // Split clusters with some criteria and calculate invariant mass
   9 // to identify them as pi0 or conversion
  10 //
  11 //
  12 //-- Author: Gustavo Conesa (LPSC-Grenoble)
  13 //_________________________________________________________________________
  14
  15
  16 // --- ROOT system ---
  17 class TList ;
  18 class TObjString;
  19 class TLorentzVector;
  20
  21 // --- ANALYSIS system ---
  22 class AliAODCaloCluster;
  23
  24 #include "AliAnaCaloTrackCorrBaseClass.h"
  25
  26 class AliAnaInsideClusterInvariantMass : public AliAnaCaloTrackCorrBaseClass {
  27
  28  public:
  29
  30   AliAnaInsideClusterInvariantMass() ; // default ctor
  31   virtual ~AliAnaInsideClusterInvariantMass() { ; } //virtual dtor
  32
  33   void         CheckLocalMaximaMCOrigin(AliVCluster* cluster, Int_t mcindex, Int_t noverlaps,
  34                                         Float_t e1,     Float_t e2,    Float_t mass);
  35                                         //, Float_t m02, TLorentzVector l1, TLorentzVector l2);
  36
  37   TObjString * GetAnalysisCuts();
  38
  39   TList      * GetCreateOutputObjects();
  40
  41   void         GetMCIndex(AliVCluster * cluster, Int_t & mcindex, Int_t & mcTag);
  42
  43   void         GetMCPrimaryKine(AliVCluster* cluster, Int_t mcindex, Int_t mcTag, Bool_t matched,
  44                                 Float_t & eprim, Float_t & asymGen, Float_t & angleGen, Int_t & noverlaps );
  45
  46   void         FillAngleHistograms(Int_t   nMax,  Bool_t  matched, Int_t   mcindex,
  47                                    Float_t en,    Float_t e1  ,    Float_t e2,
  48                                    Float_t angle, Float_t mass,    Float_t anglePrim,
  49                                    Float_t m02,   Float_t asym,    Int_t   pid,    Int_t   noverlaps);
  50
  51   void         FillArmenterosHistograms(Int_t nMax, Int_t ebin, Int_t mcindex,
  52                                         Float_t pi0E, TLorentzVector g1, TLorentzVector g2,
  53                                         Float_t m02, Int_t pid);
  54
  55   void         FillThetaStarHistograms(Int_t nMax, Bool_t matched, Int_t mcindex,
  56                                        Float_t pi0E, TLorentzVector g1, TLorentzVector g2,
  57                                        Float_t m02, Int_t pid);
  58
  59   void         FillEBinHistograms(Int_t ebin, Int_t nMax, Int_t mcindex, Float_t splitFrac,
  60                                   Float_t mass, Float_t asym, Float_t l0);
  61
  62   void         FillMCHistograms(Float_t en,        Float_t e1  , Float_t e2,
  63                                 Int_t ebin,        Int_t mcindex,Int_t noverlaps,
  64                                 Float_t l0,        Float_t mass,
  65                                 Int_t nMax,        Bool_t  matched,
  66                                 Float_t splitFrac, Float_t asym,
  67                                 Float_t eprim,     Float_t asymGen);
  68
  69   void         FillMCOverlapHistograms(Float_t en,      Float_t enprim,
  70                                        Int_t   nc,      Float_t mass,    Float_t l0,
  71                                        Float_t asym,    Float_t splitFrac,
  72                                        Int_t   nlm,     Int_t ebin,   Bool_t matched,
  73                                        Int_t   mcindex, Int_t noverlaps);
  74
  75   void         FillSSWeightHistograms(AliVCluster *cluster, Int_t nlm, Int_t absId1, Int_t absId2);
  76
  77   void         FillSSExtraHistograms(AliVCluster *cluster, Int_t nMax,
  78                                      Bool_t matched, Int_t mcindex,
  79                                      Float_t mass  , Int_t ebin);
  80
  81   void         FillNLMDiffCutHistograms(AliVCluster *cluster, AliVCaloCells *cells, Bool_t matched);
  82
  83   void         FillNCellHistograms(Int_t   ncells,  Float_t energy, Int_t nMax,
  84                                    Bool_t  matched, Int_t mcindex,
  85                                    Float_t mass   , Float_t l0);
  86
  87   void         FillTrackMatchingHistograms(AliVCluster * cluster,Int_t nMax, Int_t mcindex);
  88
  89   void         FillHistograms1(Float_t en,     Float_t e1,     Float_t e2,
  90                                Int_t nMax,     Float_t mass,   Float_t l0,
  91                                Float_t eta,    Float_t phi,
  92                                Bool_t matched, Int_t mcindex);
  93
  94
  95   void         FillHistograms2(Float_t en,     Float_t eprim,
  96                                Float_t e1,     Float_t e2,      Int_t nMax,
  97                                Float_t mass,   Float_t l0,
  98                                Bool_t matched, Int_t mcindex);
  99
 100   void         FillIdPi0Histograms(Float_t en,     Float_t e1,  Float_t e2,
 101                                    Int_t nc,       Int_t nMax,  Float_t t12diff,
 102                                    Float_t mass,   Float_t l0,
 103                                    Float_t eta,    Float_t phi,
 104                                    Bool_t matched, Int_t mcindex);
 105
 106   void         FillIdEtaHistograms(Float_t en,     Float_t e1,  Float_t e2,
 107                                    Int_t nc,       Int_t nMax,  Float_t t12diff,
 108                                    Float_t mass,   Float_t l0,
 109                                    Float_t eta,    Float_t phi,
 110                                    Bool_t matched, Int_t mcindex);
 111
 112   void         FillIdConvHistograms(Float_t en,    Int_t nMax, Float_t asym,
 113                                     Float_t mass,   Float_t l0,
 114                                     Bool_t matched, Int_t mcindex);
 115
 116   void         Init();
 117
 118   void         InitParameters();
 119
 120   void         MakeAnalysisFillHistograms() ;
 121
 122   void         Print(const Option_t * opt) const;
 123
 124   void         SetCalorimeter(TString & det)             { fCalorimeter = det ; }
 125
 126   void         SetMinNCells(Int_t cut)                   { fMinNCells   = cut ; }
 127
 128   void         SetMinBadChannelDistance(Float_t cut)     { fMinBadDist  = cut ; }
 129
 130   void         SetWCorrectionParameter(Int_t i, Float_t p = 0.07) { if( i<2 ) fWSimu[i] = p; }
 131
 132   void         SwitchOnFillAngleHistograms()             { fFillAngleHisto      = kTRUE  ; }
 133   void         SwitchOffFillAngleHistograms()            { fFillAngleHisto      = kFALSE ; }
 134
 135   void         SwitchOnFillArmenterosHistograms()        { fFillArmenterosHisto = kTRUE  ; }
 136   void         SwitchOffFillArmenterosHistograms()       { fFillArmenterosHisto = kFALSE ; }
 137
 138   void         SwitchOnFillThetaStarHistograms()         { fFillThetaStarHisto  = kTRUE  ; }
 139   void         SwitchOffFillThetaStarHistograms()        { fFillThetaStarHisto  = kFALSE ; }
 140
 141   void         SwitchOnFillExtraSSHistograms()           { fFillSSExtraHisto    = kTRUE  ; }
 142   void         SwitchOffFillExtraSSHistograms()          { fFillSSExtraHisto    = kFALSE ; }
 143
 144   void         SwitchOnFillHighMultHistograms()          { fFillHighMultHisto   = kTRUE  ; }
 145   void         SwitchOffFillHighMultHistograms()         { fFillHighMultHisto   = kFALSE ; }
 146
 147   void         SwitchOnFillIdConvHistograms()            { fFillIdConvHisto     = kTRUE  ; }
 148   void         SwitchOffFillIdConvHistograms()           { fFillIdConvHisto     = kFALSE ; }
 149
 150   void         SwitchOnFillIdEtaHistograms()             { fFillIdEtaHisto      = kTRUE  ; }
 151   void         SwitchOffFillIdEtaHistograms()            { fFillIdEtaHisto      = kFALSE ; }
 152
 153   void         SwitchOnFillTMHistograms()                { fFillTMHisto         = kTRUE  ; }
 154   void         SwitchOffFillTMHistograms()               { fFillTMHisto         = kFALSE ; }
 155
 156   void         SwitchOnFillTMResidualHistograms()        { fFillTMResidualHisto = kTRUE  ; }
 157   void         SwitchOffFillTMResidualHistograms()       { fFillTMResidualHisto = kFALSE ; }
 158
 159   void         SwitchOnFillMCPrimaryHistograms()         { fFillMCHisto         = kTRUE  ; }
 160   void         SwitchOffFillMCPrimaryHistograms()        { fFillMCHisto         = kFALSE ; }
 161
 162   void         SwitchOnFillSSWeightHistograms()          { fFillSSWeightHisto   = kTRUE  ; }
 163   void         SwitchOffFillSSWeightHistograms()         { fFillSSWeightHisto   = kFALSE ; }
 164
 165   void         SwitchOnFillNLMDiffCutsHistograms()       { fFillNLMDiffCutHisto = kTRUE  ; }
 166   void         SwitchOffFillNLMDiffCutsHistograms()      { fFillNLMDiffCutHisto = kFALSE ; }
 167
 168   void         SwitchOnFillEbinHistograms()              { fFillEbinHisto       = kTRUE  ; }
 169   void         SwitchOffFillEbinHistograms()             { fFillEbinHisto       = kFALSE ; }
 170
 171   void         SwitchOnFillMCOverlapHistograms()         { fFillMCOverlapHisto  = kTRUE  ; }
 172   void         SwitchOffFillMCOverlapHistograms()        { fFillMCOverlapHisto  = kFALSE ; }
 173
 174   void         SwitchOnFillNCellHistograms()             { fFillNCellHisto      = kTRUE  ; }
 175   void         SwitchOffFillNCellHistograms()            { fFillNCellHisto      = kFALSE ; }
 176
 177   void         SwitchOnSplitClusterDistToBad()           { fCheckSplitDistToBad = kTRUE  ; }
 178   void         SwitchOffSplitClusterDistToBad()          { fCheckSplitDistToBad = kFALSE ; }
 179
 180   void         SetNWeightForShowerShape(Int_t n)         { fSSWeightN = n ; }
 181   void         SetWeightForShowerShape(Int_t i, Float_t v)
 182                                                          { if (i < 20) fSSWeight[i] = v ; }
 183
 184   void         SetNumberOfNLocMaxSettings(Int_t n)       { fNLMSettingN = n ; }
 185   void         SetNLocMaxMinE   (Int_t i, Float_t v)     { if (i < 5) fNLMMinE   [i] = v ; }
 186   void         SetNLocMaxMinDiff(Int_t i, Float_t v)     { if (i < 5) fNLMMinDiff[i] = v ; }
 187
 188
 189   void         SetNECellCutForShowerShape(Int_t n)       { fSSECellCutN = n ; }
 190   void         SetECellCutForShowerShape(Int_t i, Float_t v)
 191                                                          { if (i < 20) fSSECellCut[i] = v ; }
 192
 193
 194   void         RecalculateClusterShowerShapeParametersWithCellCut(const AliEMCALGeometry * geom, AliVCaloCells* cells, AliVCluster * cluster,
 195                                                    Float_t & l0,   Float_t & l1,
 196                                                    Float_t & disp, Float_t & dEta, Float_t & dPhi,
 197                                                    Float_t & sEta, Float_t & sPhi, Float_t & sEtaPhi,Float_t eCellMin = 0.);
 198
 199
 200   //For histograms
 201   enum mcTypes { kmcPhoton = 1, kmcConversion = 2, kmcPi0    = 3,  kmcPi0Conv = 4,
 202                  kmcEta    = 5, kmcHadron = 6  };
 203
 204  private:
 205
 206   TString      fCalorimeter ;          // Calorimeter where the gamma is searched
 207   Int_t        fMinNCells   ;          // Study clusters with ncells larger than cut
 208   Float_t      fMinBadDist  ;          // Minimal distance to bad channel to accept cluster
 209   Float_t      fHistoECut   ;          // Fixed E cut for some histograms
 210   Bool_t       fCheckSplitDistToBad;   // Check the distance to bad channel and to EMCal borders of split clusters
 211
 212   Bool_t       fFillAngleHisto;        // Fill splitted clusters angle histograms
 213   Bool_t       fFillTMHisto ;          // Fill track matching histos,
 214   Bool_t       fFillTMResidualHisto ;  // Fill track matching histos, residuals
 215   Bool_t       fFillSSExtraHisto ;     // Fill shower shape extra histos
 216   Bool_t       fFillMCHisto ;          // Fill MC energy fraction histos
 217   Bool_t       fFillSSWeightHisto ;    // Fill weigth histograms
 218   Bool_t       fFillNLMDiffCutHisto ;  // Fill NLM histograms for different settings
 219   Bool_t       fFillEbinHisto ;        // Fill E bin histograms
 220   Bool_t       fFillMCOverlapHisto ;   // Fill MC particles overlap histograms
 221   Bool_t       fFillNCellHisto ;       // Fill n cells in cluster dependent histograms
 222   Bool_t       fFillIdConvHisto ;      // Fill histograms for clusters identified as conversion
 223   Bool_t       fFillIdEtaHisto ;       // Fill histograms for clusters identified as Eta
 224   Bool_t       fFillHighMultHisto;     // Fill centrality/event plane histograms
 225   Bool_t       fFillArmenterosHisto;   // Fill armenteros type histo
 226   Bool_t       fFillThetaStarHisto;    // Fill cosThetaStar histos
 227
 228   Float_t      fSSWeight[20];          // List of weights to test
 229   Int_t        fSSWeightN;             // Total number of weights to test
 230
 231   Float_t      fSSECellCut[20];        // List of cell min energy cuts to test
 232   Int_t        fSSECellCutN;           // Total number of cell min energy cuts to test
 233
 234   Float_t      fNLMMinE   [5];         // List of local maxima min energy
 235   Float_t      fNLMMinDiff[5];         // List of local maxima min difference cell energy
 236   Int_t        fNLMSettingN;           // Total number of NLM settings to test
 237
 238   Float_t      fWSimu[2];              // Constant and slope of the linear correction factor for the shower
 239                                        // shape weight in simulation, about 1-0.07*w
 240
 241   //Histograms
 242
 243   TH2F       * fhMassNLocMax1[7][2]  ;                  //! Split Inv Mass vs cluster E, NLM=1,  different MC particle types, track matching on/off
 244   TH2F       * fhMassNLocMax2[7][2]  ;                  //! Split Inv Mass vs cluster E, NLM=2,  different MC particle types, track matching on/off
 245   TH2F       * fhMassNLocMaxN[7][2]  ;                  //! Split Inv Mass vs cluster E, NLM>2,  different MC particle types, track matching on/off
 246
 247   TH2F       * fhMassSplitENLocMax1[7][2]  ;            //! Split Inv Mass vs E1+E2, NLM=1,  different MC particle types, track matching  on/off
 248   TH2F       * fhMassSplitENLocMax2[7][2]  ;            //! Split Inv Mass vs E1+E2, NLM=2,  different MC particle types, track matching  on/off
 249   TH2F       * fhMassSplitENLocMaxN[7][2]  ;            //! Split Inv Mass vs E1+E2, NLM>2,  different MC particle types, track matching  on/off
 250
 251   TH2F       * fhAsymNLocMax1[7][2]  ;                  //! Asymmetry of 2 highest energy cells when 1 local max vs E, 1-6 for different MC particle types
 252   TH2F       * fhAsymNLocMax2[7][2]  ;                  //! Asymmetry of 2 cells local maxima vs E,  1-6 for different MC particle types
 253   TH2F       * fhAsymNLocMaxN[7][2]  ;                  //! Asymmetry of >2 cells local maxima vs E, 1-6 for different MC particle types
 254
 255   TH2F       * fhSplitEFractionvsAsyNLocMax1[2] ;       //! sum of splitted cluster energy / cluster energy for N Local Maxima = 1 vs |A|
 256   TH2F       * fhSplitEFractionvsAsyNLocMax2[2] ;       //! sum of splitted cluster energy / cluster energy for N Local Maxima = 2 vs |A|
 257   TH2F       * fhSplitEFractionvsAsyNLocMaxN[2] ;       //! sum of splitted cluster energy / cluster energy for N Local Maxima > 2 vs |A|
 258
 259   TH2F       * fhMassAsyCutNLocMax1  ;                  //! Mass(E) asym selection, not matched, Mass of split clusters, NLM = 1
 260   TH2F       * fhMassAsyCutNLocMax2  ;                  //! Mass(E) asym selection, not matched, Mass of split clusters, NLM = 1
 261   TH2F       * fhMassAsyCutNLocMaxN  ;                  //! Mass(E) asym selection, not matched, Mass of split clusters, NLM > 2
 262
 263   TH2F       * fhM02AsyCutNLocMax1  ;                   //! M02(E) asym selection, not matched, M02, NLM = 1
 264   TH2F       * fhM02AsyCutNLocMax2  ;                   //! M02(E) asym selection, not matched, M02, NLM = 2
 265   TH2F       * fhM02AsyCutNLocMaxN  ;                   //! M02(E) asym selection, not matched, M02, NLM > 2
 266
 267   TH2F       * fhMassM02CutNLocMax1  ;                  //! Mass(E) M02 selection, not matched, Mass of split clusters, NLM = 1
 268   TH2F       * fhMassM02CutNLocMax2  ;                  //! Mass(E) M02 selection, not matched, Mass of split clusters, NLM = 1
 269   TH2F       * fhMassM02CutNLocMaxN  ;                  //! Mass(E) M02 selection, not matched, Mass of split clusters, NLM > 2
 270
 271   TH2F       * fhAsymM02CutNLocMax1  ;                  //! Asym(E) M02 selection, not matched, energy asymmetry of split clusters, NLM = 1
 272   TH2F       * fhAsymM02CutNLocMax2  ;                  //! Asym(E) M02 selection, not matched, energy asymmetry of split clusters, NLM = 2
 273   TH2F       * fhAsymM02CutNLocMaxN  ;                  //! Asym(E) M02 selection, not matched, energy asymmetry of split clusters, NLM > 2
 274
 275   TH2F       * fhMassEnCutNLocMax1  ;                   //! Mass(E) E sub-cluster cut selection, not matched, Mass of split clusters, NLM = 1
 276   TH2F       * fhMassEnCutNLocMax2  ;                   //! Mass(E) E sub-cluster cut selection, not matched, Mass of split clusters, NLM = 1
 277   TH2F       * fhMassEnCutNLocMaxN  ;                   //! Mass(E) E sub-cluster cut selection, not matched, Mass of split clusters, NLM > 2
 278
 279   TH2F       * fhM02EnCutNLocMax1  ;                    //! M02(E) E sub-cluster cut selection, not matched, Mass of split clusters, NLM = 1
 280   TH2F       * fhM02EnCutNLocMax2  ;                    //! M02(E) E sub-cluster cut selection, not matched, Mass of split clusters, NLM = 1
 281   TH2F       * fhM02EnCutNLocMaxN  ;                    //! M02(E) E sub-cluster cut selection, not matched, Mass of split clusters, NLM > 2
 282
 283   TH2F       * fhAsymEnCutNLocMax1  ;                   //! Asym(E) E sub-cluster cut selection, not matched, Mass of split clusters, NLM = 1
 284   TH2F       * fhAsymEnCutNLocMax2  ;                   //! Asym(E) E sub-cluster cut selection, not matched, Mass of split clusters, NLM = 1
 285   TH2F       * fhAsymEnCutNLocMaxN  ;                   //! Asym(E) E sub-cluster cut selection, not matched, Mass of split clusters, NLM > 2
 286
 287   TH2F       * fhSplitEFracEnCutNLocMax1  ;             //! Split E fraction (E) E sub-cluster cut selection, not matched, Mass of split clusters, NLM = 1
 288   TH2F       * fhSplitEFracEnCutNLocMax2  ;             //! Split E fraction (E) E sub-cluster cut selection, not matched, Mass of split clusters, NLM = 1
 289   TH2F       * fhSplitEFracEnCutNLocMaxN  ;             //! Split E fraction (E) E sub-cluster cut selection, not matched, Mass of split clusters, NLM > 2
 290
 291   TH2F       * fhMassSplitECutNLocMax1 ;                //! 85% of split energy, not matched, Mass of split clusters, NLM = 1
 292   TH2F       * fhMassSplitECutNLocMax2 ;                //! 85% of split energy, not matched, Mass of split clusters, NLM = 1
 293   TH2F       * fhMassSplitECutNLocMaxN ;                //! 85% of split energy, not matched, Mass of split clusters, NLM > 2
 294
 295   TH2F       * fhMassM02NLocMax1[7][2]  ;               //! Mass of splitted clusters when 1  local max vs M02, for E > 8 GeV, 1-6 for different MC particle types
 296   TH2F       * fhMassM02NLocMax2[7][2]  ;               //! Mass of splitted clusters when 2  local max vs M02, for E > 8 GeV, 1-6 for different MC particle types
 297   TH2F       * fhMassM02NLocMaxN[7][2]  ;               //! Mass of splitted clusters when >2 local max vs M02, for E > 8 GeV, 1-6 for different MC particle types
 298
 299   TH2F       * fhMassM02NLocMax1Ebin[4] ;               //! Mass of splitted clusters when 1  local max vs M02, 4 E bins, neutral clusters
 300   TH2F       * fhMassM02NLocMax2Ebin[4] ;               //! Mass of splitted clusters when 2  local max vs M02, 4 E bins, neutral clusters
 301   TH2F       * fhMassM02NLocMaxNEbin[4] ;               //! Mass of splitted clusters when >2 local max vs M02, 4 E bins, neutral clusters
 302
 303   TH2F       * fhMassAsyNLocMax1Ebin[4] ;               //! Mass of Mass of splitted clusters when 1  local max vs asymmetry, 4 E bins, neutral clusters
 304   TH2F       * fhMassAsyNLocMax2Ebin[4] ;               //! Mass of Mass of splitted clusters when 2  local max vs asymmetry, 4 E bins, neutral clusters
 305   TH2F       * fhMassAsyNLocMaxNEbin[4] ;               //! Mass of Mass of splitted clusters when >2 local max vs asymmetry, 4 E bins, neutral clusters
 306
 307   TH2F       * fhAsyMCGenRecoNLocMax1EbinPi0[4] ;       //! Generated vs reconstructed asymmetry of splitted clusters from pi0 when 1  local max, 4 E bins, neutral clusters
 308   TH2F       * fhAsyMCGenRecoNLocMax2EbinPi0[4] ;       //! Generated vs reconstructed asymmetry of splitted clusters from pi0 when 2  local max, 4 E bins, neutral clusters
 309   TH2F       * fhAsyMCGenRecoNLocMaxNEbinPi0[4] ;       //! Generated vs reconstructed asymmetry of splitted clusters from pi0 when >2 local max, 4 E bins, neutral clusters
 310
 311   TH2F       * fhAsyMCGenRecoDiffMCPi0[3];              //! reconstructed-generated asymmetry of splitted clusters vs E from pi0, for 3 NLM cases
 312   TH2F       * fhAsyMCGenRecoDiffMCPi0Conv[3];          //! reconstructed-generated asymmetry of splitted clusters vs E from converted pi0, for 3 NLM cases
 313
 314   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 particle types
 315   TH2F       * fhMassDispEtaNLocMax2[7][2]  ;           //! Mass of 2 cells local maxima, vs M02, for E > 8 GeV,  1-6 for different MC particle types
 316   TH2F       * fhMassDispEtaNLocMaxN[7][2]  ;           //! Mass of >2 cells local maxima, vs M02, for E > 8 GeV, 1-6 for different MC particle types
 317
 318   TH2F       * fhMassDispEtaNLocMax1Ebin[4] ;           //! Mass of 2 highest energy cells when 1 local max, vs M02, 4 E bins, neutral clusters
 319   TH2F       * fhMassDispEtaNLocMax2Ebin[4] ;           //! Mass of 2 cells local maxima, vs M02, 4 E bins, neutral clusters
 320   TH2F       * fhMassDispEtaNLocMaxNEbin[4] ;           //! Mass of >2 cells local maxima, vs M02, 4 E bins, neutral clusters
 321
 322   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 particle types
 323   TH2F       * fhMassDispPhiNLocMax2[7][2]  ;           //! Mass of 2 cells local maxima, vs M02, for E > 8 GeV,  1-6 for different MC particle types
 324   TH2F       * fhMassDispPhiNLocMaxN[7][2]  ;           //! Mass of >2 cells local maxima, vs M02, for E > 8 GeV, 1-6 for different MC particle types
 325
 326   TH2F       * fhMassDispPhiNLocMax1Ebin[4] ;           //! Mass of 2 highest energy cells when 1 local max, vs M02, 4 E bins, neutral clusters
 327   TH2F       * fhMassDispPhiNLocMax2Ebin[4] ;           //! Mass of 2 cells local maxima, vs M02, 4 E bins, neutral clusters
 328   TH2F       * fhMassDispPhiNLocMaxNEbin[4] ;           //! Mass of >2 cells local maxima, vs M02, 4 E bins, neutral clusters
 329
 330   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 particle types
 331   TH2F       * fhMassDispAsyNLocMax2[7][2]  ;           //! Mass of 2 cells local maxima, vs M02, for E > 8 GeV,  1-6 for different MC particle types
 332   TH2F       * fhMassDispAsyNLocMaxN[7][2]  ;           //! Mass of >2 cells local maxima, vs M02, for E > 8 GeV, 1-6 for different MC particle types
 333
 334   TH2F       * fhMassDispAsyNLocMax1Ebin[4] ;           //! Mass of 2 highest energy cells when 1 local max, vs M02, 4 E bins, neutral clusters
 335   TH2F       * fhMassDispAsyNLocMax2Ebin[4] ;           //! Mass of 2 cells local maxima, vs M02, 4 E bins, neutral clusters
 336   TH2F       * fhMassDispAsyNLocMaxNEbin[4] ;           //! Mass of >2 cells local maxima, vs M02, 4 E bins, neutral clusters
 337
 338   TH2F       * fhNLocMax      [7][2] ;                  //! Number of maxima in cluster vs E, 1-6 for different MC particle types
 339   TH2F       * fhNLocMaxM02Cut[7][2] ;                  //! Number of maxima in cluster vs E, 1-6 for different MC particle types, after SS cut
 340   TH2F       * fhNLocMaxIdPi0 [7][2] ;                  //! Number of maxima in cluster vs E, 1-6 for different MC particle types, after pi0 selection
 341
 342   TH2F       * fhSplitClusterENLocMax[7][2] ;           //! Number of maxima in cluster vs E of splitted clusters, 1-6 for different MC particle types
 343   TH2F       * fhSplitClusterEPi0NLocMax[7][2] ;        //! Number of maxima in cluster vs E of splitted clusters when cluster id as pi0, 1-6 for different MC particle types
 344
 345   TH2F       * fhLM1NLocMax      [7][2] ;               //! Split cluster 1 E distribution vs Number of maxima in cluster vs E, 1-6 for different MC particle types
 346   TH2F       * fhLM1NLocMaxM02Cut[7][2] ;               //! Split cluster 1 E distribution vs Number of maxima in cluster vs E, 1-6 for different MC particle types, after SS cut
 347   TH2F       * fhLM1NLocMaxIdPi0 [7][2] ;               //! Split cluster 1 E distribution vs Number of maxima in cluster vs E, 1-6 for different MC particle types, pi0 selection
 348
 349   TH2F       * fhLM2NLocMax      [7][2] ;               //! Split cluster 2 E distribution vs Number of maxima in cluster vs E, 1-6 for different MC particle types
 350   TH2F       * fhLM2NLocMaxM02Cut[7][2] ;               //! Split cluster 2 E distribution vs Number of maxima in cluster vs E, 1-6 for different MC particle types, after SS cut
 351   TH2F       * fhLM2NLocMaxIdPi0 [7][2] ;               //! Split cluster 2 E distribution vs Number of maxima in cluster vs E, 1-6 for different MC particle types, pi0 selection
 352
 353   TH2F       * fhM02NLocMax1  [7][2] ;                  //! M02 vs E for N max in cluster = 1, 1-6 for different MC particle types
 354   TH2F       * fhM02NLocMax2  [7][2] ;                  //! M02 vs E for N max in cluster = 2, 1-6 for different MC particle types
 355   TH2F       * fhM02NLocMaxN  [7][2] ;                  //! M02 vs E for N max in cluster > 2, 1-6 for different MC particle types
 356
 357   TH2F       * fhMCAsymM02NLocMax1MCPi0Ebin[4] ;        //! M02 vs decay asymmetry for N max in cluster = 1, for 4 energy bins
 358   TH2F       * fhMCAsymM02NLocMax2MCPi0Ebin[4] ;        //! M02 vs decay asymmetry for N max in cluster = 2, for 4 energy bins
 359   TH2F       * fhMCAsymM02NLocMaxNMCPi0Ebin[4] ;        //! M02 vs decay asymmetry for N max in cluster > 2, for 4 energy bins
 360
 361   TH2F       * fhMCGenFracNLocMax1[7][2] ;              //! E generated particle / E reconstructed vs E reconstructed for N max in cluster = 1, 1-6 for different MC particle types
 362   TH2F       * fhMCGenFracNLocMax2[7][2] ;              //! E generated particle / E reconstructed vs E reconstructed for N max in cluster = 2, 1-6 for different MC particle types
 363   TH2F       * fhMCGenFracNLocMaxN[7][2] ;              //! E generated particle / E reconstructed vs E reconstructed for N max in cluster > 2, 1-6 for different MC particle types
 364
 365   TH2F       * fhMCGenFracNLocMax1NoOverlap[7][2] ;     //! E generated particle / E reconstructed vs E reconstructed for N max in cluster = 1, 1-6 for different MC particle types, no overlap found
 366   TH2F       * fhMCGenFracNLocMax2NoOverlap[7][2] ;     //! E generated particle / E reconstructed vs E reconstructed for N max in cluster = 2, 1-6 for different MC particle types, no overlap found
 367   TH2F       * fhMCGenFracNLocMaxNNoOverlap[7][2] ;     //! E generated particle / E reconstructed vs E reconstructed for N max in cluster > 2, 1-6 for different MC particle types, no overlap found
 368
 369   TH2F       * fhMCGenFracAfterCutsNLocMax1MCPi0 ;      //! E generated particle / E reconstructed vs E reconstructed for N max in cluster = 1, MCPi0 after M02 and asymmetry cut
 370   TH2F       * fhMCGenFracAfterCutsNLocMax2MCPi0 ;      //! E generated particle / E reconstructed vs E reconstructed for N max in cluster = 2, MCPi0, after M02 and asymmetry cut
 371   TH2F       * fhMCGenFracAfterCutsNLocMaxNMCPi0 ;      //! E generated particle / E reconstructed vs E reconstructed for N max in cluster > 2, MCPi0, after M02 and asymmetry cut
 372
 373   TH2F       * fhMCGenSplitEFracNLocMax1[7][2] ;        //! E generated particle / E1+E2 reconstructed vs E reconstructed for N max in cluster = 1, 1-6 for different MC particle types
 374   TH2F       * fhMCGenSplitEFracNLocMax2[7][2] ;        //! E generated particle / E1+E2 reconstructed vs E reconstructed for N max in cluster = 2, 1-6 for different MC particle types
 375   TH2F       * fhMCGenSplitEFracNLocMaxN[7][2] ;        //! E generated particle / E1+E2 reconstructed vs E reconstructed for N max in cluster > 2, 1-6 for different MC particle types
 376
 377   TH2F       * fhMCGenSplitEFracNLocMax1NoOverlap[7][2];//! E generated particle / E1+E2 reconstructed vs E reconstructed for N max in cluster = 1, 1-6 for different MC particle types, no overlap
 378   TH2F       * fhMCGenSplitEFracNLocMax2NoOverlap[7][2];//! E generated particle / E1+E2 reconstructed vs E reconstructed for N max in cluster = 2, 1-6 for different MC particle types, no overlap
 379   TH2F       * fhMCGenSplitEFracNLocMaxNNoOverlap[7][2];//! E generated particle / E1+E2 reconstructed vs E reconstructed for N max in cluster > 2, 1-6 for different MC particle types, no overlap
 380
 381   TH2F       * fhMCGenSplitEFracAfterCutsNLocMax1MCPi0; //! E generated particle / E1+E2 reconstructed vs E reconstructed for N max in cluster = 1, 1-6 for different MC particle types
 382   TH2F       * fhMCGenSplitEFracAfterCutsNLocMax2MCPi0; //! E generated particle / E1+E2 reconstructed vs E reconstructed for N max in cluster = 2, 1-6 for different MC particle types
 383   TH2F       * fhMCGenSplitEFracAfterCutsNLocMaxNMCPi0; //! E generated particle / E1+E2 reconstructed vs E reconstructed for N max in cluster > 2, 1-6 for different MC particle types
 384
 385   TH2F       * fhMCGenEFracvsSplitEFracNLocMax1[7][2] ; //! E generated particle / E reconstructed vs E1+E2 reconstructed / E reconstructed for N max in cluster = 1, MC pi0
 386   TH2F       * fhMCGenEFracvsSplitEFracNLocMax2[7][2] ; //! E generated particle / E reconstructed vs E1+E2 reconstructed / E reconstructed for N max in cluster = 2, MC pi0
 387   TH2F       * fhMCGenEFracvsSplitEFracNLocMaxN[7][2] ; //! E generated particle / E reconstructed vs E1+E2 reconstructed / E reconstructed for N max in cluster > 2, MC pi0
 388
 389   TH2F       * fhMCGenEvsSplitENLocMax1[7][2] ;         //! E generated particle vs E1+E2 for N max in cluster = 1, 1-6 for different MC particle types
 390   TH2F       * fhMCGenEvsSplitENLocMax2[7][2] ;         //! E generated particle vs E1+E2 for N max in cluster = 2, 1-6 for different MC particle types
 391   TH2F       * fhMCGenEvsSplitENLocMaxN[7][2] ;         //! E generated particle vs E1+E2 for N max in cluster > 2, 1-6 for different MC particle types
 392
 393   TH2F       * fhMCGenFracNLocMaxEbin[7][4] ;           //! NLM vs E generated particle / E reconstructed vs E reconstructed 1-6 for different MC particle types, not matched to track
 394   TH2F       * fhMCGenFracNLocMaxEbinMatched[7][4] ;    //! NLM vs E generated particle / E reconstructed vs E reconstructed 1-6 for different MC particle types, matched to track
 395
 396   TH2F       * fhM02MCGenFracNLocMax1Ebin[7][4] ;       //! M02 vs E generated particle / E reconstructed vs E reconstructed for N max in cluster = 1, 1-6 for different MC particle types, not track matched
 397   TH2F       * fhM02MCGenFracNLocMax2Ebin[7][4] ;       //! M02 vs E generated particle / E reconstructed vs E reconstructed for N max in cluster = 2, 1-6 for different MC particle types, not track matched
 398   TH2F       * fhM02MCGenFracNLocMaxNEbin[7][4] ;       //! M02 vs E generated particle / E reconstructed vs E reconstructed for N max in cluster > 2, 1-6 for different MC particle types, not track matched
 399
 400   TH2F       * fhMassMCGenFracNLocMax1Ebin[7][4] ;      //! Mass vs E generated particle / E reconstructed vs E reconstructed for N max in cluster = 1, 1-6 for different MC particle types, not track matched
 401   TH2F       * fhMassMCGenFracNLocMax2Ebin[7][4] ;      //! Mass vs E generated particle / E reconstructed vs E reconstructed for N max in cluster = 2, 1-6 for different MC particle types, not track matched
 402   TH2F       * fhMassMCGenFracNLocMaxNEbin[7][4] ;      //! Mass vs E generated particle / E reconstructed vs E reconstructed for N max in cluster > 2, 1-6 for different MC particle types, not track matched
 403
 404   TH2F       * fhNCellNLocMax1[7][2] ;                  //! n cells in cluster vs E for N max in cluster = 1, 1-6 for different MC particle types
 405   TH2F       * fhNCellNLocMax2[7][2] ;                  //! n cells in cluster vs E for N max in cluster = 2, 1-6 for different MC particle types
 406   TH2F       * fhNCellNLocMaxN[7][2] ;                  //! n cells in cluster vs E for N max in cluster > 2, 1-6 for different MC particle types
 407
 408   TH2F       * fhNCellMassEHighNLocMax1MCPi0 ;          //! n cells in cluster vs mass for high energy clusters,  for N max in cluster = 1, for MC pi0
 409   TH2F       * fhNCellM02EHighNLocMax1MCPi0  ;          //! n cells in cluster vs m02  for high energy clusters,  for N max in cluster = 1, for MC pi0
 410   TH2F       * fhNCellMassELowNLocMax1MCPi0  ;          //! n cells in cluster vs mass for low  energy clusters,  for N max in cluster = 1, for MC pi0
 411   TH2F       * fhNCellM02ELowNLocMax1MCPi0   ;          //! n cells in cluster vs m02  for low  energy clusters,  for N max in cluster = 1, for MC pi0
 412
 413   TH2F       * fhNCellMassEHighNLocMax2MCPi0 ;          //! n cells in cluster vs mass for high energy clusters,  for N max in cluster = 2, for MC pi0
 414   TH2F       * fhNCellM02EHighNLocMax2MCPi0  ;          //! n cells in cluster vs m02  for high energy clusters,  for N max in cluster = 2, for MC pi0
 415   TH2F       * fhNCellMassELowNLocMax2MCPi0  ;          //! n cells in cluster vs mass for low  energy clusters,  for N max in cluster = 2, for MC pi0
 416   TH2F       * fhNCellM02ELowNLocMax2MCPi0   ;          //! n cells in cluster vs m02  for low  energy clusters,  for N max in cluster = 2, for MC pi0
 417
 418   TH2F       * fhNCellMassEHighNLocMaxNMCPi0 ;          //! n cells in cluster vs mass for high energy clusters,  for N max in cluster > 2, for MC pi0
 419   TH2F       * fhNCellM02EHighNLocMaxNMCPi0  ;          //! n cells in cluster vs m02  for high energy clusters,  for N max in cluster > 2, for MC pi0
 420   TH2F       * fhNCellMassELowNLocMaxNMCPi0  ;          //! n cells in cluster vs mass for low  energy clusters,  for N max in cluster > 2, for MC pi0
 421   TH2F       * fhNCellM02ELowNLocMaxNMCPi0   ;          //! n cells in cluster vs m02  for low  energy clusters,  for N max in cluster > 2, for MC pi0
 422
 423   TH2F       * fhM02Pi0NLocMax1[7][2] ;                 //! M02 for Mass around pi0, N Local Maxima = 1
 424   TH2F       * fhM02EtaNLocMax1[7][2] ;                 //! M02 for Mass around eta, N Local Maxima = 1
 425   TH2F       * fhM02ConNLocMax1[7][2] ;                 //! M02 for Mass around close to 0, N Local Maxima = 1
 426
 427   TH2F       * fhM02Pi0NLocMax2[7][2] ;                 //! M02 for Mass around pi0, N Local Maxima = 2
 428   TH2F       * fhM02EtaNLocMax2[7][2] ;                 //! M02 for Mass around eta, N Local Maxima = 2
 429   TH2F       * fhM02ConNLocMax2[7][2] ;                 //! M02 for Mass around close to 0, N Local Maxima = 2
 430
 431   TH2F       * fhM02Pi0NLocMaxN[7][2] ;                 //! M02 for Mass around pi0, N Local Maxima > 2
 432   TH2F       * fhM02EtaNLocMaxN[7][2] ;                 //! M02 for Mass around eta, N Local Maxima > 2
 433   TH2F       * fhM02ConNLocMaxN[7][2] ;                 //! M02 for Mass around close to 0, N Local Maxima > 2
 434
 435   TH2F       * fhMassPi0NLocMax1[7][2] ;                //! Mass for selected pi0, N Local Maxima = 1
 436   TH2F       * fhMassEtaNLocMax1[7][2] ;                //! Mass for selected around eta, N Local Maxima = 1
 437   TH2F       * fhMassConNLocMax1[7][2] ;                //! Mass for selected around close to 0, N Local Maxima = 1
 438
 439   TH2F       * fhMassPi0NLocMax2[7][2] ;                //! Mass for selected around pi0, N Local Maxima = 2
 440   TH2F       * fhMassEtaNLocMax2[7][2] ;                //! Mass for selected around eta, N Local Maxima = 2
 441   TH2F       * fhMassConNLocMax2[7][2] ;                //! Mass for selected around close to 0, N Local Maxima = 2
 442
 443   TH2F       * fhMassPi0NLocMaxN[7][2] ;                //! Mass for selected around pi0, N Local Maxima > 2
 444   TH2F       * fhMassEtaNLocMaxN[7][2] ;                //! Mass for selected around eta, N Local Maxima > 2
 445   TH2F       * fhMassConNLocMaxN[7][2] ;                //! Mass for selected around close to 0, N Local Maxima > 2
 446
 447   TH2F       * fhNCellPi0NLocMax1[7][2] ;               //! n cells for selected around pi0, N Local Maxima = 1
 448   TH2F       * fhNCellEtaNLocMax1[7][2] ;               //! n cells for selected around eta, N Local Maxima = 1
 449   TH2F       * fhNCellPi0NLocMax2[7][2] ;               //! n cells for selected around pi0, N Local Maxima = 2
 450   TH2F       * fhNCellEtaNLocMax2[7][2] ;               //! n cells for selected around eta, N Local Maxima = 2
 451   TH2F       * fhNCellPi0NLocMaxN[7][2] ;               //! n cells for selected around pi0, N Local Maxima > 2
 452   TH2F       * fhNCellEtaNLocMaxN[7][2] ;               //! n cells for selected around eta, N Local Maxima > 2
 453
 454   TH2F       * fhMassAfterCutsNLocMax1[7][2] ;          //! Mass after M02, asymmetry cuts for MC part, N Local Maxima = 1
 455   TH2F       * fhMassAfterCutsNLocMax2[7][2] ;          //! Mass after M02, asymmetry cuts for MC part, N Local Maxima = 2
 456   TH2F       * fhMassAfterCutsNLocMaxN[7][2] ;          //! Mass after M02, asymmetry cuts for MC part, N Local Maxima > 2
 457
 458   TH2F       * fhMassSplitEAfterCutsNLocMax1[7][2]  ;   //! Split Inv Mass vs E1+E2, NLM=1, after M02, asymmetry cuts, different MC particle types, track matching  on/off
 459   TH2F       * fhMassSplitEAfterCutsNLocMax2[7][2]  ;   //! Split Inv Mass vs E1+E2, NLM=2, after M02, asymmetry cuts, different MC particle types, track matching  on/off
 460   TH2F       * fhMassSplitEAfterCutsNLocMaxN[7][2]  ;   //! Split Inv Mass vs E1+E2, NLM>2, after M02, asymmetry cuts, different MC particle types, track matching  on/off
 461
 462   TH2F       * fhMassSplitEPi0NLocMax1[7][2]  ;         //! Split Inv Mass vs E1+E2, NLM=1, after pi0 selection, different MC particle types, track matching  on/off
 463   TH2F       * fhMassSplitEPi0NLocMax2[7][2]  ;         //! Split Inv Mass vs E1+E2, NLM=2, after pi0 selection, different MC particle types, track matching  on/off
 464   TH2F       * fhMassSplitEPi0NLocMaxN[7][2]  ;         //! Split Inv Mass vs E1+E2, NLM>2, after pi0 selection, different MC particle types, track matching  on/off
 465
 466   TH2F       * fhAsyPi0NLocMax1[7][2] ;                 //! Asy for Mass around pi0, N Local Maxima = 1
 467   TH2F       * fhAsyEtaNLocMax1[7][2] ;                 //! Asy for Mass around eta, N Local Maxima = 1
 468   TH2F       * fhAsyConNLocMax1[7][2] ;                 //! Asy for Mass around close to 0, N Local Maxima = 1
 469
 470   TH2F       * fhAsyPi0NLocMax2[7][2] ;                 //! Asy for Mass around pi0, N Local Maxima = 2
 471   TH2F       * fhAsyEtaNLocMax2[7][2] ;                 //! Asy for Mass around eta, N Local Maxima = 2
 472   TH2F       * fhAsyConNLocMax2[7][2] ;                 //! Asy for Mass around close to 0, N Local Maxima = 2
 473
 474   TH2F       * fhAsyPi0NLocMaxN[7][2] ;                 //! Asy for Mass around pi0, N Local Maxima > 2
 475   TH2F       * fhAsyEtaNLocMaxN[7][2] ;                 //! Asy for Mass around eta, N Local Maxima > 2
 476   TH2F       * fhAsyConNLocMaxN[7][2] ;                 //! Asy for Mass around close to 0, N Local Maxima > 2
 477
 478   TH2F       * fhSplitEFractionNLocMax1[7][2] ;         //! sum of splitted cluster energy / cluster energy for N Local Maxima = 1
 479   TH2F       * fhSplitEFractionNLocMax2[7][2] ;         //! sum of splitted cluster energy / cluster energy for N Local Maxima = 2
 480   TH2F       * fhSplitEFractionNLocMaxN[7][2] ;         //! sum of splitted cluster energy / cluster energy for N Local Maxima > 2
 481
 482   TH2F       * fhSplitEFractionAfterCutsNLocMax1[7][2] ; //! sum of splitted cluster energy / cluster energy for N Local Maxima = 1, after M02 and asymmetry cut
 483   TH2F       * fhSplitEFractionAfterCutsNLocMax2[7][2] ; //! sum of splitted cluster energy / cluster energy for N Local Maxima = 2, after M02 and asymmetry cut
 484   TH2F       * fhSplitEFractionAfterCutsNLocMaxN[7][2] ; //! sum of splitted cluster energy / cluster energy for N Local Maxima > 2, after M02 and asymmetry cut
 485
 486   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
 487   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
 488   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
 489
 490   TH2F       * fhAnglePairNLocMax1[7][2] ;              //! pair opening angle vs E
 491   TH2F       * fhAnglePairNLocMax2[7][2] ;              //! pair opening angle vs E
 492   TH2F       * fhAnglePairNLocMaxN[7][2] ;              //! pair opening angle vs E
 493
 494   TH2F       * fhAnglePairAfterCutsNLocMax1[7][2] ;     //! pair opening angle vs E
 495   TH2F       * fhAnglePairAfterCutsNLocMax2[7][2] ;     //! pair opening angle vs E
 496   TH2F       * fhAnglePairAfterCutsNLocMaxN[7][2] ;     //! pair opening angle vs E
 497
 498   TH2F       * fhAnglePairPi0NLocMax1[7][2] ;           //! pair opening angle vs E
 499   TH2F       * fhAnglePairPi0NLocMax2[7][2] ;           //! pair opening angle vs E
 500   TH2F       * fhAnglePairPi0NLocMaxN[7][2] ;           //! pair opening angle vs E
 501
 502   TH2F       * fhAnglePairMassNLocMax1[7][2] ;          //! pair opening angle vs Mass for E > 7 GeV
 503   TH2F       * fhAnglePairMassNLocMax2[7][2] ;          //! pair opening angle vs Mass for E > 7 GeV
 504   TH2F       * fhAnglePairMassNLocMaxN[7][2] ;          //! pair opening angle vs Mass for E > 7 GeV
 505
 506   TH2F       * fhAnglePairM02NLocMax1[7][2] ;           //! pair opening angle vs M02 for E > 7 GeV
 507   TH2F       * fhAnglePairM02NLocMax2[7][2] ;           //! pair opening angle vs M02 for E > 7 GeV
 508   TH2F       * fhAnglePairM02NLocMaxN[7][2] ;           //! pair opening angle vs M02 for E > 7 GeV
 509
 510   TH2F       * fhAnglePairPrimPi0RecoNLocMax1;          //! pair opening angle pi0 generated/reconstructed vs E
 511   TH2F       * fhAnglePairPrimPi0RecoNLocMax2;          //! pair opening angle pi0 generated/reconstructed vs E
 512   TH2F       * fhAnglePairPrimPi0RecoNLocMaxN;          //! pair opening angle pi0 generated/reconstructed vs E
 513
 514   TH2F       * fhAnglePairPrimPi0vsRecoNLocMax1;        //! pair opening angle pi0 generated vs reconstructed
 515   TH2F       * fhAnglePairPrimPi0vsRecoNLocMax2;        //! pair opening angle pi0 generated vs reconstructed
 516   TH2F       * fhAnglePairPrimPi0vsRecoNLocMaxN;        //! pair opening angle pi0 generated vs reconstructed
 517
 518   TH2F       * fhAnglePairOverM02NLocMax1[7][2];        //! pair opening angle / m02 vs E, NLM=1
 519   TH2F       * fhAnglePairOverM02NLocMax2[7][2];        //! pair opening angle / m02 vs E, NLM=2
 520   TH2F       * fhAnglePairOverM02NLocMaxN[7][2];        //! pair opening angle / m02 vs E, NLM=N
 521
 522   TH2F       * fhAnglePairOverM02NLocMax1Overlap0[7][2];//! pair opening angle / m02 vs E, NLM=1
 523   TH2F       * fhAnglePairOverM02NLocMax2Overlap0[7][2];//! pair opening angle / m02 vs E, NLM=2
 524   TH2F       * fhAnglePairOverM02NLocMaxNOverlap0[7][2];//! pair opening angle / m02 vs E, NLM=N
 525
 526   TH2F       * fhAnglePairPrimPi0OverM02NLocMax1;       //! pair opening angle / m02 vs E, NLM=1, prim pi0
 527   TH2F       * fhAnglePairPrimPi0OverM02NLocMax2;       //! pair opening angle / m02 vs E, NLM=2, prim pi0
 528   TH2F       * fhAnglePairPrimPi0OverM02NLocMaxN;       //! pair opening angle / m02 vs E, NLM=N, prim pi0
 529
 530   TH2F       * fhArmNLocMax1[7][4]  ;                   //! Armenteros of 2 highest energy cells when 1 local max vs E, 1-6 for different MC particle types
 531   TH2F       * fhArmNLocMax2[7][4]  ;                   //! Armenteros of 2 cells local maxima vs E,  1-6 for different MC particle types
 532   TH2F       * fhArmNLocMaxN[7][4]  ;                   //! Armenteros of >2 cells local maxima vs E, 1-6 for different MC particle types
 533
 534   TH2F       * fhArmAfterCutsNLocMax1[7][4] ;           //! Armenteros after M02, asymmetry cuts for pi0, N Local Maxima = 1
 535   TH2F       * fhArmAfterCutsNLocMax2[7][4] ;           //! Armenteros after M02, asymmetry cuts for pi0, N Local Maxima = 2
 536   TH2F       * fhArmAfterCutsNLocMaxN[7][4] ;           //! Armenteros after M02, asymmetry cuts for pi0, N Local Maxima > 2
 537
 538   TH2F       * fhArmPi0NLocMax1[7][4] ;                 //! Armenteros for selected pi0, N Local Maxima = 1
 539   TH2F       * fhArmPi0NLocMax2[7][4] ;                 //! Armenteros for selected pi0, N Local Maxima = 2
 540   TH2F       * fhArmPi0NLocMaxN[7][4] ;                 //! Armenteros for selected pi0, N Local Maxima > 2
 541
 542   TH2F       * fhCosThStarNLocMax1[7][2] ;              //! cos(theta^star) vs E, NLM=1
 543   TH2F       * fhCosThStarNLocMax2[7][2] ;              //! cos(theta^star) vs E, NLM=2
 544   TH2F       * fhCosThStarNLocMaxN[7][2] ;              //! cos(theta^star) vs E, NLM>2
 545
 546   TH2F       * fhCosThStarAfterCutsNLocMax1[7][2] ;     //! cos(theta^star) vs E, after M02, asymmetry cuts, NLM=1
 547   TH2F       * fhCosThStarAfterCutsNLocMax2[7][2] ;     //! cos(theta^star) vs E, after M02, asymmetry cuts, NLM=2
 548   TH2F       * fhCosThStarAfterCutsNLocMaxN[7][2] ;     //! cos(theta^star) vs E, after M02, asymmetry cuts, NLM>2
 549
 550   TH2F       * fhCosThStarPi0NLocMax1[7][2] ;           //! cos(theta^star) vs E, after M02, asymmetry and pi0 mass cuts, NLM=1
 551   TH2F       * fhCosThStarPi0NLocMax2[7][2] ;           //! cos(theta^star) vs E, after M02, asymmetry and pi0 mass cuts, NLM=2
 552   TH2F       * fhCosThStarPi0NLocMaxN[7][2] ;           //! cos(theta^star) vs E, after M02, asymmetry and pi0 mass cuts, NLM>2
 553
 554   TH2F       * fhTrackMatchedDEtaNLocMax1[7] ;          //! Eta distance between track and cluster vs cluster E, 1 local maximum
 555   TH2F       * fhTrackMatchedDPhiNLocMax1[7] ;          //! Phi distance between track and cluster vs cluster E, 1 local maximum
 556   TH2F       * fhTrackMatchedDEtaNLocMax2[7] ;          //! Eta distance between track and cluster vs cluster E, 2 local maximum
 557   TH2F       * fhTrackMatchedDPhiNLocMax2[7] ;          //! Phi distance between track and cluster vs cluster E, 2 local maximum
 558   TH2F       * fhTrackMatchedDEtaNLocMaxN[7] ;          //! Eta distance between track and cluster vs cluster E, more than 2 local maximum
 559   TH2F       * fhTrackMatchedDPhiNLocMaxN[7] ;          //! Phi distance between track and cluster vs cluster E, more than 2 local maximum
 560
 561   TH2F       * fhTrackMatchedDEtaNLocMax1Pos[7] ;       //! Eta distance between track and cluster vs cluster E, 1 local maximum
 562   TH2F       * fhTrackMatchedDPhiNLocMax1Pos[7] ;       //! Phi distance between track and cluster vs cluster E, 1 local maximum
 563   TH2F       * fhTrackMatchedDEtaNLocMax2Pos[7] ;       //! Eta distance between track and cluster vs cluster E, 2 local maximum
 564   TH2F       * fhTrackMatchedDPhiNLocMax2Pos[7] ;       //! Phi distance between track and cluster vs cluster E, 2 local maximum
 565   TH2F       * fhTrackMatchedDEtaNLocMaxNPos[7] ;       //! Eta distance between track and cluster vs cluster E, more than 2 local maximum
 566   TH2F       * fhTrackMatchedDPhiNLocMaxNPos[7] ;       //! Phi distance between track and cluster vs cluster E, more than 2 local maximum
 567
 568   TH2F       * fhTrackMatchedDEtaNLocMax1Neg[7] ;       //! Eta distance between track and cluster vs cluster E, 1 local maximum
 569   TH2F       * fhTrackMatchedDPhiNLocMax1Neg[7] ;       //! Phi distance between track and cluster vs cluster E, 1 local maximum
 570   TH2F       * fhTrackMatchedDEtaNLocMax2Neg[7] ;       //! Eta distance between track and cluster vs cluster E, 2 local maximum
 571   TH2F       * fhTrackMatchedDPhiNLocMax2Neg[7] ;       //! Phi distance between track and cluster vs cluster E, 2 local maximum
 572   TH2F       * fhTrackMatchedDEtaNLocMaxNNeg[7] ;       //! Eta distance between track and cluster vs cluster E, more than 2 local maximum
 573   TH2F       * fhTrackMatchedDPhiNLocMaxNNeg[7] ;       //! Phi distance between track and cluster vs cluster E, more than 2 local maximum
 574
 575   TH2F       * fhCentralityPi0NLocMax1 ;                //! Centrality for selected pi0, N Local Maxima = 1
 576   TH2F       * fhCentralityEtaNLocMax1 ;                //! Centrality for selected eta, N Local Maxima = 1
 577   TH2F       * fhCentralityPi0NLocMax2 ;                //! Centrality for selected pi0, N Local Maxima = 2
 578   TH2F       * fhCentralityEtaNLocMax2 ;                //! Centrality for selected eta, N Local Maxima = 2
 579   TH2F       * fhCentralityPi0NLocMaxN ;                //! Centrality for selected pi0, N Local Maxima > 2
 580   TH2F       * fhCentralityEtaNLocMaxN ;                //! Centrality for selected eta, N Local Maxima > 2
 581
 582   TH2F       * fhEventPlanePi0NLocMax1 ;                //! Event plane for selected pi0, N Local Maxima = 1
 583   TH2F       * fhEventPlaneEtaNLocMax1 ;                //! Event plane for selected eta, N Local Maxima = 1
 584   TH2F       * fhEventPlanePi0NLocMax2 ;                //! Event plane for selected pi0, N Local Maxima = 2
 585   TH2F       * fhEventPlaneEtaNLocMax2 ;                //! Event plane for selected eta, N Local Maxima = 2
 586   TH2F       * fhEventPlanePi0NLocMaxN ;                //! Event plane for selected pi0, N Local Maxima > 2
 587   TH2F       * fhEventPlaneEtaNLocMaxN ;                //! Event plane for selected eta, N Local Maxima > 2
 588
 589   TH2F       * fhClusterEtaPhiNLocMax1 ;                //! Eta vs Phi of clusters with N Local Maxima = 1, E > 8 GeV
 590   TH2F       * fhClusterEtaPhiNLocMax2 ;                //! Eta vs Phi of clusters with N Local Maxima = 2, E > 8 GeV
 591   TH2F       * fhClusterEtaPhiNLocMaxN ;                //! Eta vs Phi of clusters with N Local Maxima > 2, E > 8 GeV
 592   TH2F       * fhPi0EtaPhiNLocMax1 ;                    //! Eta vs Phi of pi0's with N Local Maxima = 1, E > 8 GeV
 593   TH2F       * fhPi0EtaPhiNLocMax2 ;                    //! Eta vs Phi of pi0's with N Local Maxima = 2, E > 8 GeV
 594   TH2F       * fhPi0EtaPhiNLocMaxN ;                    //! Eta vs Phi of pi0's with N Local Maxima > N, E > 8 GeV
 595   TH2F       * fhEtaEtaPhiNLocMax1 ;                    //! Eta vs Phi of eta's with N Local Maxima = 1, E > 8 GeV
 596   TH2F       * fhEtaEtaPhiNLocMax2 ;                    //! Eta vs Phi of eta's with N Local Maxima = 2, E > 8 GeV
 597   TH2F       * fhEtaEtaPhiNLocMaxN ;                    //! Eta vs Phi of eta's with N Local Maxima > N, E > 8 GeV
 598
 599   TH2F       * fhPi0CellE[3] ;                          //! pi0's energy vs cluster cell energy with NLM = 1, = 2, > 2
 600   TH2F       * fhPi0CellEFrac[3] ;                      //! pi0's energy vs cluster cell energy fraction with NLM = 1, = 2, > 2
 601   TH2F       * fhPi0CellLogEFrac[3] ;                   //! pi0's energy vs cluster log cell energy fraction with NLM = 1, = 2, > 2
 602   TH2F       * fhPi0CellEMaxEMax2Frac   [3];            //! pi0's energy vs fraction of 2 main maxima energy with NLM = 1, = 2, > 2
 603   TH2F       * fhPi0CellEMaxClusterFrac [3];            //! pi0's energy vs energy fraction of main   LM and cluster energy with NLM = 1, = 2, > 2
 604   TH2F       * fhPi0CellEMax2ClusterFrac[3];            //! pi0's energy vs energy fraction of second LM and cluster energy with NLM = 1, = 2, > 2
 605   TH2F       * fhPi0CellEMaxFrac [3];                   //! pi0's energy vs energy fraction of main LM and cluster cell energy with NLM = 1, = 2, > 2
 606   TH2F       * fhPi0CellEMax2Frac [3];                  //! pi0's energy vs energy fraction of second LM and cluster cell energy with NLM = 1, = 2, > 2
 607
 608   TH2F       * fhM02WeightPi0[3][20] ;                  //! M02 for selected pi0 with different weight, with NLM = 1, = 2, > 2
 609   TH2F       * fhM02ECellCutPi0[3][20] ;                //! M02 for selected pi0 with different cut on cell energy, with NLM = 1, = 2, > 2
 610
 611   TH2F       * fhPi0EPairDiffTimeNLM1;                  //! E vs Pair of clusters time difference vs E, for selected pi0, NLM=1
 612   TH2F       * fhPi0EPairDiffTimeNLM2;                  //! E vs Pair of clusters time difference vs E, for selected pi0, NLM=2
 613   TH2F       * fhPi0EPairDiffTimeNLMN;                  //! E vs Pair of clusters time difference vs E, for selected pi0, NLM>2
 614   TH2F       * fhEtaEPairDiffTimeNLM1;                  //! E vs Pair of clusters time difference vs E, for selected eta, NLM=1
 615   TH2F       * fhEtaEPairDiffTimeNLM2;                  //! E vs Pair of clusters time difference vs E, for selected eta, NLM=2
 616   TH2F       * fhEtaEPairDiffTimeNLMN;                  //! E vs Pair of clusters time difference vs E, for selected eta, NLM>2
 617
 618   TH2F       * fhMCEM02Overlap0[3][7];                  //! E vs M02 for different MC origin, no other MC particles contributes, neutral cluster
 619   TH2F       * fhMCEM02Overlap1[3][7];                  //! E vs M02 for different MC origin, 1  other MC particles contributes, neutral cluster
 620   TH2F       * fhMCEM02OverlapN[3][7];                  //! E vs M02 for different MC origin, N  other MC particles contributes, neutral cluster
 621   TH2F       * fhMCEM02Overlap0Match[3][7];             //! E vs M02 for different MC origin, no other MC particles contributes, charged cluster
 622   TH2F       * fhMCEM02Overlap1Match[3][7];             //! E vs M02 for different MC origin, 1  other MC particles contributes, charged cluster
 623   TH2F       * fhMCEM02OverlapNMatch[3][7];             //! E vs M02 for different MC origin, N  other MC particles contributes, charged cluster
 624
 625   TH2F       * fhMCEMassOverlap0[3][7];                 //! E vs Mass for different MC origin, no other MC particles contributes, neutral cluster
 626   TH2F       * fhMCEMassOverlap1[3][7];                 //! E vs Mass for different MC origin, 1  other MC particles contributes, neutral cluster
 627   TH2F       * fhMCEMassOverlapN[3][7];                 //! E vs Mass for different MC origin, N  other MC particles contributes, neutral cluster
 628   TH2F       * fhMCEMassOverlap0Match[3][7];            //! E vs Mass for different MC origin, no other MC particles contributes, charged cluster
 629   TH2F       * fhMCEMassOverlap1Match[3][7];            //! E vs Mass for different MC origin, 1  other MC particles contributes, charged cluster
 630   TH2F       * fhMCEMassOverlapNMatch[3][7];            //! E vs Mass for different MC origin, N  other MC particles contributes, charged cluster
 631
 632   TH2F       * fhMCESplitEFracOverlap0[3][7];           //! E vs sum of splitted cluster energy / cluster energy for different MC origin, no other MC particles contributes, neutral cluster
 633   TH2F       * fhMCESplitEFracOverlap1[3][7];           //! E vs sum of splitted cluster energy / cluster energy for different MC origin, 1  other MC particles contributes, neutral cluster
 634   TH2F       * fhMCESplitEFracOverlapN[3][7];           //! E vs sum of splitted cluster energy / cluster energy for different MC origin, N  other MC particles contributes, neutral cluster
 635   TH2F       * fhMCESplitEFracOverlap0Match[3][7];      //! E vs sum of splitted cluster energy / cluster energy for different MC origin, no other MC particles contributes, charged cluster
 636   TH2F       * fhMCESplitEFracOverlap1Match[3][7];      //! E vs sum of splitted cluster energy / cluster energy for different MC origin, 1  other MC particles contributes, charged cluster
 637   TH2F       * fhMCESplitEFracOverlapNMatch[3][7];      //! E vs sum of splitted cluster energy / cluster energy for different MC origin, N  other MC particles contributes, charged cluster
 638
 639   TH2F       * fhMCEAsymOverlap0[3][7];                 //! E vs sum of splitted cluster energy asymmetry for different MC origin, no other MC particles contributes, neutral cluster
 640   TH2F       * fhMCEAsymOverlap1[3][7];                 //! E vs sum of splitted cluster energy asymmetry for different MC origin, 1  other MC particles contributes, neutral cluster
 641   TH2F       * fhMCEAsymOverlapN[3][7];                 //! E vs sum of splitted cluster energy asymmetry for different MC origin, N  other MC particles contributes, neutral cluster
 642   TH2F       * fhMCEAsymOverlap0Match[3][7];            //! E vs sum of splitted cluster energy asymmetry for different MC origin, no other MC particles contributes, charged cluster
 643   TH2F       * fhMCEAsymOverlap1Match[3][7];            //! E vs sum of splitted cluster energy asymmetry for different MC origin, 1  other MC particles contributes, charged cluster
 644   TH2F       * fhMCEAsymOverlapNMatch[3][7];            //! E vs sum of splitted cluster energy asymmetry for different MC origin, N  other MC particles contributes, charged cluster
 645
 646   TH2F       * fhMCENCellOverlap0[3][7];                //! E vs sum of splitted cluster energy asymmetry for different MC origin, no other MC particles contributes, neutral cluster
 647   TH2F       * fhMCENCellOverlap1[3][7];                //! E vs sum of splitted cluster energy asymmetry for different MC origin, 1  other MC particles contributes, neutral cluster
 648   TH2F       * fhMCENCellOverlapN[3][7];                //! E vs sum of splitted cluster energy asymmetry for different MC origin, N  other MC particles contributes, neutral cluster
 649   TH2F       * fhMCENCellOverlap0Match[3][7];           //! E vs sum of splitted cluster energy asymmetry for different MC origin, no other MC particles contributes, charged cluster
 650   TH2F       * fhMCENCellOverlap1Match[3][7];           //! E vs sum of splitted cluster energy asymmetry for different MC origin, 1  other MC particles contributes, charged cluster
 651   TH2F       * fhMCENCellOverlapNMatch[3][7];           //! E vs sum of splitted cluster energy asymmetry for different MC origin, N  other MC particles contributes, charged cluster
 652
 653   TH2F       * fhMCEEpriOverlap0[3][7];                 //! E reco vs primary for different MC origin, no other MC particles contributes, neutral cluster
 654   TH2F       * fhMCEEpriOverlap1[3][7];                 //! E reco vs primary for different MC origin, 1  other MC particles contributes, neutral cluster
 655   TH2F       * fhMCEEpriOverlapN[3][7];                 //! E reco vs primary for different MC origin, N  other MC particles contributes, neutral cluster
 656   TH2F       * fhMCEEpriOverlap0Match[3][7];            //! E reco vs primary for different MC origin, no other MC particles contributes, charged cluster
 657   TH2F       * fhMCEEpriOverlap1Match[3][7];            //! E reco vs primary for different MC origin, 1  other MC particles contributes, charged cluster
 658   TH2F       * fhMCEEpriOverlapNMatch[3][7];            //! E reco vs primary for different MC origin, N  other MC particles contributes, charged cluster
 659
 660   TH2F       * fhMCEEpriOverlap0IdPi0[3][7];            //! E reco vs primary for different MC origin, no other MC particles contributes, neutral cluster, neutral clusters id as pi0
 661   TH2F       * fhMCEEpriOverlap1IdPi0[3][7];            //! E reco vs primary for different MC origin, 1  other MC particles contributes, neutral cluster, neutral clusters id as pi0
 662   TH2F       * fhMCEEpriOverlapNIdPi0[3][7];            //! E reco vs primary for different MC origin, 1  other MC particles contributes, neutral cluster, neutral clusters is as pi0
 663
 664   TH2F       * fhMCPi0MassM02Overlap0[3][4];            //! MC Pi0 M02 vs Mass for different MC origin, no other MC particles contributes, neutral cluster, 4 E bins
 665   TH2F       * fhMCPi0MassM02Overlap1[3][4];            //! MC Pi0 M02 vs Mass for different MC origin, 1  other MC particles contributes, neutral cluster, 4 E bins
 666   TH2F       * fhMCPi0MassM02OverlapN[3][4];            //! MC Pi0 M02 vs Mass for different MC origin, N  other MC particles contributes, neutral cluster, 4 E bins
 667   TH2F       * fhMCPi0MassM02Overlap0Match[3][4];       //! MC Pi0 M02 vs Mass for different MC origin, no other MC particles contributes, charged cluster, 4 E bins
 668   TH2F       * fhMCPi0MassM02Overlap1Match[3][4];       //! MC Pi0 M02 vs Mass for different MC origin, 1  other MC particles contributes, charged cluster, 4 E bins
 669   TH2F       * fhMCPi0MassM02OverlapNMatch[3][4];       //! MC Pi0 M02 vs Mass for different MC origin, N  other MC particles contributes, charged cluster, 4 E bins
 670
 671   TH2F       * fhMCENOverlaps[3][7];                    //! E vs number of Overlaps in MC, neutral cluster
 672   TH2F       * fhMCENOverlapsMatch[3][7];               //! E vs number of Overlaps in MC, charged cluster
 673
 674   TH2F       * fhMCPi0HighNLMPair;                      //! E vs NLM when cluster originated in pi0 merging and highest energy local maxima correspond to 2 photons
 675   TH2F       * fhMCPi0LowNLMPair;                       //! E vs NLM when cluster originated in pi0 merging and a pair of local maxima except highest energy correspond to 2 photons
 676   TH2F       * fhMCPi0AnyNLMPair;                       //! E vs NLM when cluster originated in pi0 merging and a both highest energy pairs and other pairs correspond to 2 photons
 677   TH2F       * fhMCPi0NoneNLMPair;                      //! E vs NLM when cluster originated in pi0 merging and a both no NLM corresponds to the photons
 678   // No match between highest energy local maxima and highest energy MC particle
 679   TH2F       * fhMCPi0HighNLMPairNoMCMatch;             //! E vs NLM when cluster originated in pi0 merging and highest energy local maxima correspond to 2 photons
 680   TH2F       * fhMCPi0LowNLMPairNoMCMatch;              //! E vs NLM when cluster originated in pi0 merging and a pair of local maxima except highest energy correspond to 2 photons
 681   TH2F       * fhMCPi0AnyNLMPairNoMCMatch;              //! E vs NLM when cluster originated in pi0 merging and a both highest energy pairs and other pairs correspond to 2 photons
 682   TH2F       * fhMCPi0NoneNLMPairNoMCMatch;             //! E vs NLM when cluster originated in pi0 merging and a both no NLM corresponds to the photons
 683
 684   TH2F       * fhMCPi0HighNLMPairOverlap;              //! E vs NLM when cluster originated in pi0 merging and highest energy local maxima correspond to 2 photons, overlap
 685   TH2F       * fhMCPi0LowNLMPairOverlap;               //! E vs NLM when cluster originated in pi0 merging and a pair of local maxima except highest energy correspond to 2 photons, overlap
 686   TH2F       * fhMCPi0AnyNLMPairOverlap;               //! E vs NLM when cluster originated in pi0 merging and a both highest energy pairs and other pairs correspond to 2 photons, overlap
 687   TH2F       * fhMCPi0NoneNLMPairOverlap;              //! E vs NLM when cluster originated in pi0 merging and a both no NLM corresponds to the photons, overlap
 688   // No match between highest energy local maxima and highest energy MC particle
 689   TH2F       * fhMCPi0HighNLMPairNoMCMatchOverlap;     //! E vs NLM when cluster originated in pi0 merging and highest energy local maxima correspond to 2 photons, overlap
 690   TH2F       * fhMCPi0LowNLMPairNoMCMatchOverlap;      //! E vs NLM when cluster originated in pi0 merging and a pair of local maxima except highest energy correspond to 2 photons, overlap
 691   TH2F       * fhMCPi0AnyNLMPairNoMCMatchOverlap;      //! E vs NLM when cluster originated in pi0 merging and a both highest energy pairs and other pairs correspond to 2 photons, overlap
 692   TH2F       * fhMCPi0NoneNLMPairNoMCMatchOverlap;     //! E vs NLM when cluster originated in pi0 merging and a both no NLM corresponds to the photons, overlap
 693
 694   TH2F       * fhMCPi0DecayPhotonHitHighLM;             //! E vs NLM when cluster originated in pi0 merging and MC photon decay hit the cell local maxima
 695   TH2F       * fhMCPi0DecayPhotonAdjHighLM;             //! E vs NLM when cluster originated in pi0 merging and MC photon decay hit the adjacent cell local maxima
 696   TH2F       * fhMCPi0DecayPhotonHitOtherLM;            //! E vs NLM when cluster originated in pi0 merging and MC photon decay hit the cell local maximas, not high
 697   TH2F       * fhMCPi0DecayPhotonAdjOtherLM;            //! E vs NLM when cluster originated in pi0 merging and MC photon decay do not hit the adjacent cell local maximas, not high
 698   TH2F       * fhMCPi0DecayPhotonAdjacent;              //! E vs NLM when cluster originated in pi0 merging and MC photon decay hit adjacen cells, not 2 LM
 699   TH2F       * fhMCPi0DecayPhotonHitNoLM;               //! E vs NLM when cluster originated in pi0 merging and MC photon decay do not hit the cell local maximas
 700
 701   TH2F       * fhMCPi0DecayPhotonHitHighLMOverlap;      //! E vs NLM when cluster originated in pi0 merging and MC photon decay hit the cell local maxima, overlap
 702   TH2F       * fhMCPi0DecayPhotonAdjHighLMOverlap;      //! E vs NLM when cluster originated in pi0 merging and MC photon decay hit the adjacent cell local maxima, overlap
 703   TH2F       * fhMCPi0DecayPhotonHitOtherLMOverlap;     //! E vs NLM when cluster originated in pi0 merging and MC photon decay hit the cell local maximas, not high, overlap
 704   TH2F       * fhMCPi0DecayPhotonAdjOtherLMOverlap;     //! E vs NLM when cluster originated in pi0 merging and MC photon decay do not hit the adjacent cell local maximas, not high, overlap
 705   TH2F       * fhMCPi0DecayPhotonAdjacentOverlap;       //! E vs NLM when cluster originated in pi0 merging and MC photon decay hit adjacen cells, not 2 LM, overlap
 706   TH2F       * fhMCPi0DecayPhotonHitNoLMOverlap;        //! E vs NLM when cluster originated in pi0 merging and MC photon decay do not hit the cell local maximas, overlap
 707
 708   TH2F       * fhMCPi0DecayPhotonHitHighLMDiffELM1[3];             //! E vs Ephoton-Esplit cluster when cluster originated in pi0 merging and MC photon decay hit the cell local maxima
 709   TH2F       * fhMCPi0DecayPhotonAdjHighLMDiffELM1[3];             //! E vs Ephoton-Esplit cluster when cluster originated in pi0 merging and MC photon decay hit the adjacent cell local maxima
 710   TH2F       * fhMCPi0DecayPhotonHitOtherLMDiffELM1[3];            //! E vs Ephoton-Esplit when cluster originated in pi0 merging and MC photon decay hit the cell local maximas, not high
 711   TH2F       * fhMCPi0DecayPhotonAdjOtherLMDiffELM1[3];            //! E vs Ephoton-Esplit when cluster originated in pi0 merging and MC photon decay do not hit the adjacent cell local maximas, not high
 712
 713   TH2F       * fhMCPi0DecayPhotonHitHighLMOverlapDiffELM1[3];      //! E vs Ephoton-Esplit cluster when cluster originated in pi0 merging and MC photon decay hit the cell local maxima
 714   TH2F       * fhMCPi0DecayPhotonAdjHighLMOverlapDiffELM1[3];      //! E vs Ephoton-Esplit when cluster originated in pi0 merging and MC photon decay hit the adjacent cell local maxima, overlap
 715   TH2F       * fhMCPi0DecayPhotonHitOtherLMOverlapDiffELM1[3];     //! E vs Ephoton-Esplit when cluster originated in pi0 merging and MC photon decay hit the cell local maximas, not high, overlap
 716   TH2F       * fhMCPi0DecayPhotonAdjOtherLMOverlapDiffELM1[3];     //! E vs Ephoton-Esplit when cluster originated in pi0 merging and MC photon decay do not hit the adjacent cell local maximas, not high, overlap
 717
 718   TH2F       * fhMCPi0DecayPhotonHitHighLMDiffELM2[3];             //! E vs Ephoton-Esplit when cluster originated in pi0 merging and MC photon decay hit the cell local maxima
 719   TH2F       * fhMCPi0DecayPhotonAdjHighLMDiffELM2[3];             //! E vs Ephoton-Esplit when cluster originated in pi0 merging and MC photon decay hit the adjacent cell local maxima
 720   TH2F       * fhMCPi0DecayPhotonHitOtherLMDiffELM2[3];            //! E vs Ephoton-Esplit when cluster originated in pi0 merging and MC photon decay hit the cell local maximas, not high
 721   TH2F       * fhMCPi0DecayPhotonAdjOtherLMDiffELM2[3];            //! E vs Ephoton-Esplit when cluster originated in pi0 merging and MC photon decay do not hit the adjacent cell local maximas, not high
 722
 723   TH2F       * fhMCPi0DecayPhotonHitHighLMOverlapDiffELM2[3];      //! E vs Ephoton-Esplit cluster when cluster originated in pi0 merging and MC photon decay hit the cell local maxima
 724   TH2F       * fhMCPi0DecayPhotonAdjHighLMOverlapDiffELM2[3];      //! E vs Ephoton-Esplit when cluster originated in pi0 merging and MC photon decay hit the adjacent cell local maxima, overlap
 725   TH2F       * fhMCPi0DecayPhotonHitOtherLMOverlapDiffELM2[3];     //! E vs Ephoton-Esplit when cluster originated in pi0 merging and MC photon decay hit the cell local maximas, not high, overlap
 726   TH2F       * fhMCPi0DecayPhotonAdjOtherLMOverlapDiffELM2[3];     //! E vs Ephoton-Esplit when cluster originated in pi0 merging and MC photon decay do not hit the adjacent cell local maximas, not high, overlap
 727
 728
 729   TH2F       * fhMCPi0DecayPhotonHitHighLMDiffELM1vsELM1[3];             //! E vs Ephoton-Esplit cluster when cluster originated in pi0 merging and MC photon decay hit the cell local maxima
 730   TH2F       * fhMCPi0DecayPhotonAdjHighLMDiffELM1vsELM1[3];             //! E vs Ephoton-Esplit cluster when cluster originated in pi0 merging and MC photon decay hit the adjacent cell local maxima
 731   TH2F       * fhMCPi0DecayPhotonHitOtherLMDiffELM1vsELM1[3];            //! E vs Ephoton-Esplit when cluster originated in pi0 merging and MC photon decay hit the cell local maximas, not high
 732   TH2F       * fhMCPi0DecayPhotonAdjOtherLMDiffELM1vsELM1[3];            //! E vs Ephoton-Esplit when cluster originated in pi0 merging and MC photon decay do not hit the adjacent cell local maximas, not high
 733
 734   TH2F       * fhMCPi0DecayPhotonHitHighLMOverlapDiffELM1vsELM1[3];      //! E vs Ephoton-Esplit cluster when cluster originated in pi0 merging and MC photon decay hit the cell local maxima
 735   TH2F       * fhMCPi0DecayPhotonAdjHighLMOverlapDiffELM1vsELM1[3];      //! E vs Ephoton-Esplit when cluster originated in pi0 merging and MC photon decay hit the adjacent cell local maxima, overlap
 736   TH2F       * fhMCPi0DecayPhotonHitOtherLMOverlapDiffELM1vsELM1[3];     //! E vs Ephoton-Esplit when cluster originated in pi0 merging and MC photon decay hit the cell local maximas, not high, overlap
 737   TH2F       * fhMCPi0DecayPhotonAdjOtherLMOverlapDiffELM1vsELM1[3];     //! E vs Ephoton-Esplit when cluster originated in pi0 merging and MC photon decay do not hit the adjacent cell local maximas, not high, overlap
 738
 739   TH2F       * fhMCPi0DecayPhotonHitHighLMDiffELM2vsELM2[3];             //! E vs Ephoton-Esplit when cluster originated in pi0 merging and MC photon decay hit the cell local maxima
 740   TH2F       * fhMCPi0DecayPhotonAdjHighLMDiffELM2vsELM2[3];             //! E vs Ephoton-Esplit when cluster originated in pi0 merging and MC photon decay hit the adjacent cell local maxima
 741   TH2F       * fhMCPi0DecayPhotonHitOtherLMDiffELM2vsELM2[3];            //! E vs Ephoton-Esplit when cluster originated in pi0 merging and MC photon decay hit the cell local maximas, not high
 742   TH2F       * fhMCPi0DecayPhotonAdjOtherLMDiffELM2vsELM2[3];            //! E vs Ephoton-Esplit when cluster originated in pi0 merging and MC photon decay do not hit the adjacent cell local maximas, not high
 743
 744   TH2F       * fhMCPi0DecayPhotonHitHighLMOverlapDiffELM2vsELM2[3];      //! E vs Ephoton-Esplit cluster when cluster originated in pi0 merging and MC photon decay hit the cell local maxima
 745   TH2F       * fhMCPi0DecayPhotonAdjHighLMOverlapDiffELM2vsELM2[3];      //! E vs Ephoton-Esplit when cluster originated in pi0 merging and MC photon decay hit the adjacent cell local maxima, overlap
 746   TH2F       * fhMCPi0DecayPhotonHitOtherLMOverlapDiffELM2vsELM2[3];     //! E vs Ephoton-Esplit when cluster originated in pi0 merging and MC photon decay hit the cell local maximas, not high, overlap
 747   TH2F       * fhMCPi0DecayPhotonAdjOtherLMOverlapDiffELM2vsELM2[3];     //! E vs Ephoton-Esplit when cluster originated in pi0 merging and MC photon decay do not hit the adjacent cell local maximas, not high, overlap
 748
 749
 750   TH2F       * fhMCPi0DecayPhotonHitHighLMMass[3];                 //! E vs Mass when cluster originated in pi0 merging and MC photon decay hit the cell local maxima
 751   TH2F       * fhMCPi0DecayPhotonAdjHighLMMass[3];                 //! E vs Mass when cluster originated in pi0 merging and MC photon decay hit the adjacent cell local maxima
 752   TH2F       * fhMCPi0DecayPhotonHitOtherLMMass[3];                //! E vs Mass when cluster originated in pi0 merging and MC photon decay hit the cell local maximas, not high
 753   TH2F       * fhMCPi0DecayPhotonAdjOtherLMMass[3];                //! E vs Mass when cluster originated in pi0 merging and MC photon decay do not hit the adjacent cell local maximas, not high
 754   TH2F       * fhMCPi0DecayPhotonAdjacentMass[3];                  //! E vs Mass when cluster originated in pi0 merging and MC photon decay hit adjacen cells, not 2 LM
 755   TH2F       * fhMCPi0DecayPhotonHitNoLMMass[3];                   //! E vs Mass when cluster originated in pi0 merging and MC photon decay do not hit the cell local maximas
 756
 757   TH2F       * fhMCPi0DecayPhotonHitHighLMOverlapMass[3];          //! E vs Mass when cluster originated in pi0 merging and MC photon decay hit the cell local maxima, overlap
 758   TH2F       * fhMCPi0DecayPhotonAdjHighLMOverlapMass[3];          //! E vs Mass when cluster originated in pi0 merging and MC photon decay hit the adjacent cell local maxima, overlap
 759   TH2F       * fhMCPi0DecayPhotonHitOtherLMOverlapMass[3];         //! E vs Mass when cluster originated in pi0 merging and MC photon decay hit the cell local maximas, not high, overlap
 760   TH2F       * fhMCPi0DecayPhotonAdjOtherLMOverlapMass[3];         //! E vs Mass when cluster originated in pi0 merging and MC photon decay do not hit the adjacent cell local maximas, not high, overlap
 761   TH2F       * fhMCPi0DecayPhotonAdjacentOverlapMass[3];           //! E vs Mass when cluster originated in pi0 merging and MC photon decay hit adjacen cells, not 2 LM, overlap
 762   TH2F       * fhMCPi0DecayPhotonHitNoLMOverlapMass[3];            //! E vs Mass when cluster originated in pi0 merging and MC photon decay do not hit the cell local maximas, overlap
 763
 764
 765   TH2F       * fhMCEOverlapType;                        //! what particles overlap with pi0, neutral clusters
 766   TH2F       * fhMCEOverlapTypeMatch;                   //! what particles overlap with pi0, charged clusters
 767
 768   TH2F       * fhMassBadDistClose[3];                   //! split mass of clusters with second LM close to bad channel
 769   TH2F       * fhM02BadDistClose[3];                    //! m02 of clusters with second LM close to bad channel
 770   TH2F       * fhMassOnBorder[3];                       //! split mass of clusters with second LM on EMCAL border
 771   TH2F       * fhM02OnBorder[3];                        //! m02 of clusters with second LM close to EMCAL border
 772
 773
 774   TH2F       * fhNLocMaxDiffCut   [5][5]   [2] ;        //! Number of maxima for different values of min Loc Max value and min difference between cells, matched/unmatched with tracks
 775   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 between cells, matched/unmatched with tracks
 776   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 between cells, matched/unmatched with tracks
 777
 778   TH2F       * fhNLocMaxDiffCutPi0   [5][5]   [2] ;     //! Number of maxima for different values of min Loc Max value and min difference between cells, matched/unmatched with tracks, cluster selected as pi0
 779   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 between cells, matched/unmatched with tracks, cluster selected as pi0
 780   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 between cells, matched/unmatched with tracks
 781
 782   AliAnaInsideClusterInvariantMass(              const AliAnaInsideClusterInvariantMass & split) ; // cpy ctor
 783   AliAnaInsideClusterInvariantMass & operator = (const AliAnaInsideClusterInvariantMass & split) ; // cpy assignment
 784
 785   ClassDef(AliAnaInsideClusterInvariantMass,30)
 786
 787 } ;
 788
 789 #endif //ALIANAINSIDECLUSTERINVARIANTMASS_H
 790
 791
 792