Add histograms to check effect of M02 cut on split cluster asymmetry

diff --git a/PWGGA/CaloTrackCorrelations/AliAnaInsideClusterInvariantMass.cxx b/PWGGA/CaloTrackCorrelations/AliAnaInsideClusterInvariantMass.cxx
index b812e5dd108f42af64f5b8af141099a8632b420e..da78c2dd82fec54c3aa13fac905558598581acc6 100755 (executable)
--- a/PWGGA/CaloTrackCorrelations/AliAnaInsideClusterInvariantMass.cxx
+++ b/PWGGA/CaloTrackCorrelations/AliAnaInsideClusterInvariantMass.cxx
@@ -60,6 +60,7 @@ AliAnaInsideClusterInvariantMass::AliAnaInsideClusterInvariantMass() :
   fFillSSExtraHisto(kFALSE),
   fFillMCFractionHisto(kFALSE),
   fhMassM02CutNLocMax1(0),    fhMassM02CutNLocMax2(0),    fhMassM02CutNLocMaxN(0),
+  fhAsymM02CutNLocMax1(0),    fhAsymM02CutNLocMax2(0),    fhAsymM02CutNLocMaxN(0),
   fhMassSplitECutNLocMax1(0), fhMassSplitECutNLocMax2(0), fhMassSplitECutNLocMaxN(0),
   fhMassAsyCutNLocMax1(0),    fhMassAsyCutNLocMax2(0),    fhMassAsyCutNLocMaxN(0)
 {
@@ -312,6 +313,21 @@ TList * AliAnaInsideClusterInvariantMass::GetCreateOutputObjects()
   fhMassM02CutNLocMaxN->SetXTitle("E (GeV)");
   outputContainer->Add(fhMassM02CutNLocMaxN) ;   
   
+  fhAsymM02CutNLocMax1  = new TH2F("hAsymM02CutNLocMax1","Asymmetry of NLM=1  vs cluster Energy, M02Cut", nptbins,ptmin,ptmax,200,-1,1);
+  fhAsymM02CutNLocMax1->SetYTitle("(E_{1}-E_{2})/(E_{1}+E_{2})");
+  fhAsymM02CutNLocMax1->SetXTitle("E (GeV)");
+  outputContainer->Add(fhAsymM02CutNLocMax1) ;
+  
+  fhAsymM02CutNLocMax2  = new TH2F("hAsymM02CutNLocMax2","Asymmetry of NLM=2  vs cluster Energy, M02Cut", nptbins,ptmin,ptmax,200,-1,1);
+  fhAsymM02CutNLocMax2->SetYTitle("(E_{1}-E_{2})/(E_{1}+E_{2})");
+  fhAsymM02CutNLocMax2->SetXTitle("E (GeV)");
+  outputContainer->Add(fhAsymM02CutNLocMax2) ;
+
+  fhAsymM02CutNLocMaxN  = new TH2F("hAsymM02CutNLocMaxN","Asymmetry of NLM>2  vs cluster Energy, M02Cut", nptbins,ptmin,ptmax,200,-1,1);
+  fhAsymM02CutNLocMaxN->SetYTitle("(E_{1}-E_{2})/(E_{1}+E_{2})");
+  fhAsymM02CutNLocMaxN->SetXTitle("E (GeV)");
+  outputContainer->Add(fhAsymM02CutNLocMaxN) ;
+  
   fhMassAsyCutNLocMax1  = new TH2F("hMassAsyCutNLocMax1","Invariant mass of splitted cluster with NLM=1 vs E, with |A|>0.8",
                                    nptbins,ptmin,ptmax,mbins,mmin,mmax); 
   fhMassAsyCutNLocMax1->SetYTitle("M (GeV/c^{2})");
@@ -1621,6 +1637,7 @@ void  AliAnaInsideClusterInvariantMass::MakeAnalysisFillHistograms()
         if(GetCaloPID()->IsInMergedM02Range(en,l0,nMax))
         {
           fhMassM02CutNLocMax1->Fill(en,mass);
+          fhAsymM02CutNLocMax1->Fill(en,asym );
           if(GetCaloPID()->IsInPi0SplitAsymmetryRange(en,asym,nMax)) fhMassAsyCutNLocMax1->Fill(en,mass);
         }
       }
@@ -1647,6 +1664,7 @@ void  AliAnaInsideClusterInvariantMass::MakeAnalysisFillHistograms()
         if(GetCaloPID()->IsInMergedM02Range(en,l0,nMax))
         {
           fhMassM02CutNLocMax2->Fill(en,mass);
+          fhAsymM02CutNLocMax2->Fill(en,asym );
           if(GetCaloPID()->IsInPi0SplitAsymmetryRange(en,asym,nMax)) fhMassAsyCutNLocMax2->Fill(en,mass);
         }
       }
@@ -1673,6 +1691,7 @@ void  AliAnaInsideClusterInvariantMass::MakeAnalysisFillHistograms()
         if(GetCaloPID()->IsInMergedM02Range(en,l0,nMax))
         {
           fhMassM02CutNLocMaxN->Fill(en,mass);
+          fhAsymM02CutNLocMaxN->Fill(en,asym );
           if(GetCaloPID()->IsInPi0SplitAsymmetryRange(en,asym,nMax)) fhMassAsyCutNLocMaxN->Fill(en,mass);
         }
       }

diff --git a/PWGGA/CaloTrackCorrelations/AliAnaInsideClusterInvariantMass.h b/PWGGA/CaloTrackCorrelations/AliAnaInsideClusterInvariantMass.h
index d3524833eb76c4d54b024586605de8b41cbf3fb1..7caccc8d21f02904d9256161b998f4274ff52498 100755 (executable)
--- a/PWGGA/CaloTrackCorrelations/AliAnaInsideClusterInvariantMass.h
+++ b/PWGGA/CaloTrackCorrelations/AliAnaInsideClusterInvariantMass.h
@@ -94,17 +94,21 @@ class AliAnaInsideClusterInvariantMass : public AliAnaCaloTrackCorrBaseClass {
   TH2F       * fhSplitEFractionvsAsyNLocMax2[2] ;       //! sum of splitted cluster energy / cluster energy for N Local Maxima = 2 vs |A|
   TH2F       * fhSplitEFractionvsAsyNLocMaxN[2] ;       //! sum of splitted cluster energy / cluster energy for N Local Maxima > 2 vs |A|  
   
-  TH2F       * fhMassM02CutNLocMax1  ;                  //! M02(E) selection, not matched, Mass of 2 highest energy cells when 1 local max vs E, 1-6 for different MC particle types 
-  TH2F       * fhMassM02CutNLocMax2  ;                  //! M02(E) selection, not matched, Mass of 2 cells local maxima vs E,  1-6 for different MC particle types
-  TH2F       * fhMassM02CutNLocMaxN  ;                  //! M02(E) selection, not matched, Mass of >2 cells local maxima vs E, 1-6 for different MC particle types  
+  TH2F       * fhMassM02CutNLocMax1  ;                  //! M02(E) selection, not matched, Mass of split clusters, NLM = 1
+  TH2F       * fhMassM02CutNLocMax2  ;                  //! M02(E) selection, not matched, Mass of split clusters, NLM = 1
+  TH2F       * fhMassM02CutNLocMaxN  ;                  //! M02(E) selection, not matched, Mass of split clusters, NLM > 2
 
-  TH2F       * fhMassSplitECutNLocMax1 ;                //! 85% of split energy, not matched, Mass of 2 highest energy cells when 1 local max vs E, 1-6 for different MC particle types 
-  TH2F       * fhMassSplitECutNLocMax2 ;                //! 85% of split energy, not matched, Mass of 2 cells local maxima vs E,  1-6 for different MC particle types
-  TH2F       * fhMassSplitECutNLocMaxN ;                //! 85% of split energy, not matched, Mass of >2 cells local maxima vs E, 1-6 for different MC particle types    
+  TH2F       * fhAsymM02CutNLocMax1  ;                  //! M02(E) selection, not matched, energy asymmetry of split clusters, NLM = 1
+  TH2F       * fhAsymM02CutNLocMax2  ;                  //! M02(E) selection, not matched, energy asymmetry of split clusters, NLM = 2
+  TH2F       * fhAsymM02CutNLocMaxN  ;                  //! M02(E) selection, not matched, energy asymmetry of split clusters, NLM > 2
   
-  TH2F       * fhMassAsyCutNLocMax1  ;                  //! |A|>0.8 selection, not matched, Mass of 2 highest energy cells when 1 local max vs E, 1-6 for different MC particle types 
-  TH2F       * fhMassAsyCutNLocMax2  ;                  //! |A|>0.8 selection, not matched, Mass of 2 cells local maxima vs E,  1-6 for different MC particle types
-  TH2F       * fhMassAsyCutNLocMaxN  ;                  //! |A|>0.8 selection, not matched, Mass of >2 cells local maxima vs E, 1-6 for different MC particle types  
+  TH2F       * fhMassSplitECutNLocMax1 ;                //! 85% of split energy, not matched, Mass of split clusters, NLM = 1
+  TH2F       * fhMassSplitECutNLocMax2 ;                //! 85% of split energy, not matched, Mass of split clusters, NLM = 1
+  TH2F       * fhMassSplitECutNLocMaxN ;                //! 85% of split energy, not matched, Mass of split clusters, NLM > 2    
+  
+  TH2F       * fhMassAsyCutNLocMax1  ;                  //! |A|>0.8 selection, Mass of split clusters, NLM = 1
+  TH2F       * fhMassAsyCutNLocMax2  ;                  //! |A|>0.8 selection, Mass of split clusters, NLM = 1
+  TH2F       * fhMassAsyCutNLocMaxN  ;                  //! |A|>0.8 selection, Mass of split clusters, NLM > 2
   
   TH2F       * fhMassM02NLocMax1[7][2]  ;               //! Mass of 2 highest energy cells when 1 local max, vs M02, for E > 7 GeV, 1-6 for different MC particle types 
   TH2F       * fhMassM02NLocMax2[7][2]  ;               //! Mass of 2 cells local maxima, vs M02, for E > 7 GeV,  1-6 for different MC particle types
@@ -242,7 +246,7 @@ class AliAnaInsideClusterInvariantMass : public AliAnaCaloTrackCorrBaseClass {
   AliAnaInsideClusterInvariantMass(              const AliAnaInsideClusterInvariantMass & split) ; // cpy ctor
   AliAnaInsideClusterInvariantMass & operator = (const AliAnaInsideClusterInvariantMass & split) ; // cpy assignment
   
-  ClassDef(AliAnaInsideClusterInvariantMass,17)
+  ClassDef(AliAnaInsideClusterInvariantMass,18)
   
 } ;