add centrality dependent histograms when checking the origin of the particle

diff --git a/PWGGA/CaloTrackCorrelations/AliAnaInsideClusterInvariantMass.cxx b/PWGGA/CaloTrackCorrelations/AliAnaInsideClusterInvariantMass.cxx
index 14302f127a67e1763c89642647bd58b95c2afeb8..d144082ac30d8d67260c05501bc12671b6122c98 100755 (executable)
--- a/PWGGA/CaloTrackCorrelations/AliAnaInsideClusterInvariantMass.cxx
+++ b/PWGGA/CaloTrackCorrelations/AliAnaInsideClusterInvariantMass.cxx
@@ -68,9 +68,7 @@ AliAnaInsideClusterInvariantMass::AliAnaInsideClusterInvariantMass() :
   fhMCGenSplitEFracAfterCutsNLocMax1MCPi0(0),
   fhMCGenSplitEFracAfterCutsNLocMax2MCPi0(0),
   fhMCGenSplitEFracAfterCutsNLocMaxNMCPi0(0),
-  fhCentralityPi0NLocMax1(0), fhCentralityEtaNLocMax1(0),
-  fhCentralityPi0NLocMax2(0), fhCentralityEtaNLocMax2(0),
-  fhCentralityPi0NLocMaxN(0), fhCentralityEtaNLocMaxN(0),
+
   fhEventPlanePi0NLocMax1(0), fhEventPlaneEtaNLocMax1(0),
   fhEventPlanePi0NLocMax2(0), fhEventPlaneEtaNLocMax2(0),
   fhEventPlanePi0NLocMaxN(0), fhEventPlaneEtaNLocMaxN(0)
@@ -168,6 +166,15 @@ AliAnaInsideClusterInvariantMass::AliAnaInsideClusterInvariantMass() :
       fhSplitEFractionAfterCutsNLocMax1[i][j] = 0 ;
       fhSplitEFractionAfterCutsNLocMax2[i][j] = 0 ;
       fhSplitEFractionAfterCutsNLocMaxN[i][j] = 0 ;
+      
+      fhCentralityPi0NLocMax1[i][j] = 0 ;
+      fhCentralityEtaNLocMax1[i][j] = 0 ;
+
+      fhCentralityPi0NLocMax2[i][j] = 0 ;
+      fhCentralityEtaNLocMax2[i][j] = 0 ;
+
+      fhCentralityPi0NLocMaxN[i][j] = 0 ;
+      fhCentralityEtaNLocMaxN[i][j] = 0 ;      
     }
     
     for(Int_t jj = 0; jj < 4; jj++)
@@ -304,84 +311,6 @@ TList * AliAnaInsideClusterInvariantMass::GetCreateOutputObjects()
   
   TString sMatched[] = {"","Matched"};
   
-  
-  // E vs centrality
-  
-  fhCentralityPi0NLocMax1  = new TH2F("hCentralityPi0NLocMax1","E vs Centrality, selected pi0 cluster with NLM=1",
-                                nptbins,ptmin,ptmax,100,0,100);
-  fhCentralityPi0NLocMax1->SetYTitle("Centrality");
-  fhCentralityPi0NLocMax1->SetXTitle("E (GeV)");
-  outputContainer->Add(fhCentralityPi0NLocMax1) ;
-  
-  fhCentralityPi0NLocMax2  = new TH2F("hCentralityPi0NLocMax2","E vs Centrality, selected pi0 cluster with NLM=2",
-                                nptbins,ptmin,ptmax,100,0,100);
-  fhCentralityPi0NLocMax2->SetYTitle("Centrality");
-  fhCentralityPi0NLocMax2->SetXTitle("E (GeV)");
-  outputContainer->Add(fhCentralityPi0NLocMax2) ;
-  
-  fhCentralityPi0NLocMaxN  = new TH2F("hCentralityPi0NLocMaxN","E vs Centrality, selected pi0 cluster with NLM>1",
-                                nptbins,ptmin,ptmax,100,0,100);
-  fhCentralityPi0NLocMaxN->SetYTitle("Centrality");
-  fhCentralityPi0NLocMaxN->SetXTitle("E (GeV)");
-  outputContainer->Add(fhCentralityPi0NLocMaxN) ;
-  
-  fhCentralityEtaNLocMax1  = new TH2F("hCentralityEtaNLocMax1","E vs Centrality, selected pi0 cluster with NLM=1",
-                                nptbins,ptmin,ptmax,100,0,100);
-  fhCentralityEtaNLocMax1->SetYTitle("Centrality");
-  fhCentralityEtaNLocMax1->SetXTitle("E (GeV)");
-  outputContainer->Add(fhCentralityEtaNLocMax1) ;
-  
-  fhCentralityEtaNLocMax2  = new TH2F("hCentralityEtaNLocMax2","E vs Centrality, selected pi0 cluster with NLM=2",
-                                nptbins,ptmin,ptmax,100,0,100);
-  fhCentralityEtaNLocMax2->SetYTitle("Centrality");
-  fhCentralityEtaNLocMax2->SetXTitle("E (GeV)");
-  outputContainer->Add(fhCentralityEtaNLocMax2) ;
-  
-  fhCentralityEtaNLocMaxN  = new TH2F("hCentralityEtaNLocMaxN","E vs Centrality, selected pi0 cluster with NLM>1",
-                                nptbins,ptmin,ptmax,100,0,100);
-  fhCentralityEtaNLocMaxN->SetYTitle("Centrality");
-  fhCentralityEtaNLocMaxN->SetXTitle("E (GeV)");
-  outputContainer->Add(fhCentralityEtaNLocMaxN) ;
-  
-  // E vs Event plane angle
-  
-  fhEventPlanePi0NLocMax1  = new TH2F("hEventPlanePi0NLocMax1","E vs Event Plane Angle, selected pi0 cluster with NLM=1",
-                                      nptbins,ptmin,ptmax,100,0,TMath::Pi());
-  fhEventPlanePi0NLocMax1->SetYTitle("Event Plane Angle (rad)");
-  fhEventPlanePi0NLocMax1->SetXTitle("E (GeV)");
-  outputContainer->Add(fhEventPlanePi0NLocMax1) ;
-  
-  fhEventPlanePi0NLocMax2  = new TH2F("hEventPlanePi0NLocMax2","E vs Event Plane Angle, selected pi0 cluster with NLM=2",
-                                      nptbins,ptmin,ptmax,100,0,TMath::Pi());
-  fhEventPlanePi0NLocMax2->SetYTitle("Event Plane Angle (rad)");
-  fhEventPlanePi0NLocMax2->SetXTitle("E (GeV)");
-  outputContainer->Add(fhEventPlanePi0NLocMax2) ;
-  
-  fhEventPlanePi0NLocMaxN  = new TH2F("hEventPlanePi0NLocMaxN","E vs Event Plane Angle, selected pi0 cluster with NLM>1",
-                                      nptbins,ptmin,ptmax,100,0,TMath::Pi());
-  fhEventPlanePi0NLocMaxN->SetYTitle("Event Plane Angle (rad)");
-  fhEventPlanePi0NLocMaxN->SetXTitle("E (GeV)");
-  outputContainer->Add(fhEventPlanePi0NLocMaxN) ;
-  
-  fhEventPlaneEtaNLocMax1  = new TH2F("hEventPlaneEtaNLocMax1","E vs Event Plane Angle, selected pi0 cluster with NLM=1",
-                                      nptbins,ptmin,ptmax,100,0,TMath::Pi());
-  fhEventPlaneEtaNLocMax1->SetYTitle("Event Plane Angle (rad)");
-  fhEventPlaneEtaNLocMax1->SetXTitle("E (GeV)");
-  outputContainer->Add(fhEventPlaneEtaNLocMax1) ;
-  
-  fhEventPlaneEtaNLocMax2  = new TH2F("hEventPlaneEtaNLocMax2","E vs Event Plane Angle, selected pi0 cluster with NLM=2",
-                                      nptbins,ptmin,ptmax,100,0,TMath::Pi());
-  fhEventPlaneEtaNLocMax2->SetYTitle("Event Plane Angle (rad)");
-  fhEventPlaneEtaNLocMax2->SetXTitle("E (GeV)");
-  outputContainer->Add(fhEventPlaneEtaNLocMax2) ;
-  
-  fhEventPlaneEtaNLocMaxN  = new TH2F("hEventPlaneEtaNLocMaxN","E vs Event Plane Angle, selected pi0 cluster with NLM>1",
-                                      nptbins,ptmin,ptmax,100,0,TMath::Pi());
-  fhEventPlaneEtaNLocMaxN->SetYTitle("Event Plane Angle (rad)");
-  fhEventPlaneEtaNLocMaxN->SetXTitle("E (GeV)");
-  outputContainer->Add(fhEventPlaneEtaNLocMaxN) ;
-
-  
   for(Int_t i = 0; i < n; i++)
   {
     for(Int_t j = 0; j < 2; j++)
@@ -739,21 +668,67 @@ TList * AliAnaInsideClusterInvariantMass::GetCreateOutputObjects()
         outputContainer->Add(fhNCellNLocMaxN[i][j]) ;
       }
       
+      
+      // E vs centrality
+      
+      fhCentralityPi0NLocMax1[i][j]  = new TH2F(Form("hCentralityPi0NLocMax1%s%s",pname[i].Data(),sMatched[j].Data()),
+                                                Form("E vs Centrality, selected pi0 cluster with NLM=1, %s",ptype[i].Data()),
+                                                nptbins,ptmin,ptmax,100,0,100);
+      fhCentralityPi0NLocMax1[i][j]->SetYTitle("Centrality");
+      fhCentralityPi0NLocMax1[i][j]->SetXTitle("E (GeV)");
+      outputContainer->Add(fhCentralityPi0NLocMax1[i][j]) ;
+      
+      fhCentralityPi0NLocMax2[i][j]  = new TH2F(Form("hCentralityPi0NLocMax2%s%s",pname[i].Data(),sMatched[j].Data()),
+                                                Form("E vs Centrality, selected pi0 cluster with NLM=2, %s",ptype[i].Data()),
+                                                nptbins,ptmin,ptmax,100,0,100);
+      fhCentralityPi0NLocMax2[i][j]->SetYTitle("Centrality");
+      fhCentralityPi0NLocMax2[i][j]->SetXTitle("E (GeV)");
+      outputContainer->Add(fhCentralityPi0NLocMax2[i][j]) ;
+      
+      fhCentralityPi0NLocMaxN[i][j]  = new TH2F(Form("hCentralityPi0NLocMaxN%s%s",pname[i].Data(),sMatched[j].Data()),
+                                                Form("E vs Centrality, selected pi0 cluster with NLM>1, %s",ptype[i].Data()),
+                                                nptbins,ptmin,ptmax,100,0,100);
+      fhCentralityPi0NLocMaxN[i][j]->SetYTitle("Centrality");
+      fhCentralityPi0NLocMaxN[i][j]->SetXTitle("E (GeV)");
+      outputContainer->Add(fhCentralityPi0NLocMaxN[i][j]) ;
+      
+      fhCentralityEtaNLocMax1[i][j]  = new TH2F(Form("hCentralityEtaNLocMax1%s%s",pname[i].Data(),sMatched[j].Data()),
+                                                Form("E vs Centrality, selected pi0 cluster with NLM=1, %s",ptype[i].Data()),
+                                                nptbins,ptmin,ptmax,100,0,100);
+      fhCentralityEtaNLocMax1[i][j]->SetYTitle("Centrality");
+      fhCentralityEtaNLocMax1[i][j]->SetXTitle("E (GeV)");
+      outputContainer->Add(fhCentralityEtaNLocMax1[i][j]) ;
+      
+      fhCentralityEtaNLocMax2[i][j]  = new TH2F(Form("hCentralityEtaNLocMax2%s%s",pname[i].Data(),sMatched[j].Data()),
+                                                Form("E vs Centrality, selected pi0 cluster with NLM=2, %s",ptype[i].Data()),
+                                                nptbins,ptmin,ptmax,100,0,100);
+      fhCentralityEtaNLocMax2[i][j]->SetYTitle("Centrality");
+      fhCentralityEtaNLocMax2[i][j]->SetXTitle("E (GeV)");
+      outputContainer->Add(fhCentralityEtaNLocMax2[i][j]) ;
+      
+      fhCentralityEtaNLocMaxN[i][j]  = new TH2F(Form("hCentralityEtaNLocMaxN%s%s",pname[i].Data(),sMatched[j].Data()),
+                                                Form("E vs Centrality, selected pi0 cluster with NLM>1, %s",ptype[i].Data()),
+                                                nptbins,ptmin,ptmax,100,0,100);
+      fhCentralityEtaNLocMaxN[i][j]->SetYTitle("Centrality");
+      fhCentralityEtaNLocMaxN[i][j]->SetXTitle("E (GeV)");
+      outputContainer->Add(fhCentralityEtaNLocMaxN[i][j]) ;
+      
+      
       fhM02Pi0NLocMax1[i][j]     = new TH2F(Form("hM02Pi0NLocMax1%s%s",pname[i].Data(),sMatched[j].Data()),
                                            Form("#lambda_{0}^{2} vs E for mass range [%2.2f-%2.2f] GeV/c^{2} %s, for N Local max = 1",
                                                 GetCaloPID()->GetPi0MinMass(),GetCaloPID()->GetPi0MaxMass(),ptype[i].Data()),
-                                           nptbins,ptmin,ptmax,ssbins,ssmin,ssmax); 
+                                           nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
       fhM02Pi0NLocMax1[i][j]   ->SetYTitle("#lambda_{0}^{2}");
       fhM02Pi0NLocMax1[i][j]   ->SetXTitle("E (GeV)");
-      outputContainer->Add(fhM02Pi0NLocMax1[i][j]) ; 
+      outputContainer->Add(fhM02Pi0NLocMax1[i][j]) ;
       
       fhM02EtaNLocMax1[i][j]     = new TH2F(Form("hM02EtaNLocMax1%s%s",pname[i].Data(),sMatched[j].Data()),
                                            Form("#lambda_{0}^{2} vs E for mass range [%2.2f-%2.2f] GeV/c^{2}, %s, for N Local max = 1",
                                                 GetCaloPID()->GetEtaMinMass(),GetCaloPID()->GetEtaMaxMass(),ptype[i].Data()),
-                                           nptbins,ptmin,ptmax,ssbins,ssmin,ssmax); 
+                                           nptbins,ptmin,ptmax,ssbins,ssmin,ssmax);
       fhM02EtaNLocMax1[i][j]   ->SetYTitle("#lambda_{0}^{2}");
       fhM02EtaNLocMax1[i][j]   ->SetXTitle("E (GeV)");
-      outputContainer->Add(fhM02EtaNLocMax1[i][j]) ; 
+      outputContainer->Add(fhM02EtaNLocMax1[i][j]) ;
       
       fhM02ConNLocMax1[i][j]    = new TH2F(Form("hM02ConNLocMax1%s%s",pname[i].Data(),sMatched[j].Data()),
                                           Form("#lambda_{0}^{2} vs E for mass range [%2.2f-%2.2f] GeV/c^{2}, %s, for N Local max = 1",
@@ -1046,6 +1021,45 @@ TList * AliAnaInsideClusterInvariantMass::GetCreateOutputObjects()
       }
   } // MC particle list
  
+  // E vs Event plane angle
+  
+  fhEventPlanePi0NLocMax1  = new TH2F("hEventPlanePi0NLocMax1","E vs Event Plane Angle, selected pi0 cluster with NLM=1",
+                                      nptbins,ptmin,ptmax,100,0,TMath::Pi());
+  fhEventPlanePi0NLocMax1->SetYTitle("Event Plane Angle (rad)");
+  fhEventPlanePi0NLocMax1->SetXTitle("E (GeV)");
+  outputContainer->Add(fhEventPlanePi0NLocMax1) ;
+  
+  fhEventPlanePi0NLocMax2  = new TH2F("hEventPlanePi0NLocMax2","E vs Event Plane Angle, selected pi0 cluster with NLM=2",
+                                      nptbins,ptmin,ptmax,100,0,TMath::Pi());
+  fhEventPlanePi0NLocMax2->SetYTitle("Event Plane Angle (rad)");
+  fhEventPlanePi0NLocMax2->SetXTitle("E (GeV)");
+  outputContainer->Add(fhEventPlanePi0NLocMax2) ;
+  
+  fhEventPlanePi0NLocMaxN  = new TH2F("hEventPlanePi0NLocMaxN","E vs Event Plane Angle, selected pi0 cluster with NLM>1",
+                                      nptbins,ptmin,ptmax,100,0,TMath::Pi());
+  fhEventPlanePi0NLocMaxN->SetYTitle("Event Plane Angle (rad)");
+  fhEventPlanePi0NLocMaxN->SetXTitle("E (GeV)");
+  outputContainer->Add(fhEventPlanePi0NLocMaxN) ;
+  
+  fhEventPlaneEtaNLocMax1  = new TH2F("hEventPlaneEtaNLocMax1","E vs Event Plane Angle, selected pi0 cluster with NLM=1",
+                                      nptbins,ptmin,ptmax,100,0,TMath::Pi());
+  fhEventPlaneEtaNLocMax1->SetYTitle("Event Plane Angle (rad)");
+  fhEventPlaneEtaNLocMax1->SetXTitle("E (GeV)");
+  outputContainer->Add(fhEventPlaneEtaNLocMax1) ;
+  
+  fhEventPlaneEtaNLocMax2  = new TH2F("hEventPlaneEtaNLocMax2","E vs Event Plane Angle, selected pi0 cluster with NLM=2",
+                                      nptbins,ptmin,ptmax,100,0,TMath::Pi());
+  fhEventPlaneEtaNLocMax2->SetYTitle("Event Plane Angle (rad)");
+  fhEventPlaneEtaNLocMax2->SetXTitle("E (GeV)");
+  outputContainer->Add(fhEventPlaneEtaNLocMax2) ;
+  
+  fhEventPlaneEtaNLocMaxN  = new TH2F("hEventPlaneEtaNLocMaxN","E vs Event Plane Angle, selected pi0 cluster with NLM>1",
+                                      nptbins,ptmin,ptmax,100,0,TMath::Pi());
+  fhEventPlaneEtaNLocMaxN->SetYTitle("Event Plane Angle (rad)");
+  fhEventPlaneEtaNLocMaxN->SetXTitle("E (GeV)");
+  outputContainer->Add(fhEventPlaneEtaNLocMaxN) ;
+  
+  
   for(Int_t i = 0; i < 4; i++)
   {  
     fhMassM02NLocMax1Ebin[i]  = new TH2F(Form("hMassM02NLocMax1Ebin%d",i),
@@ -1413,27 +1427,27 @@ TList * AliAnaInsideClusterInvariantMass::GetCreateOutputObjects()
   }
   
   for(Int_t j = 0; j < 2; j++)
-  {  
-  fhSplitEFractionvsAsyNLocMax1[j]     = new TH2F(Form("hSplitEFractionvsAsyNLocMax1%s",sMatched[j].Data()),
-                                                Form("(E1+E2)/E_{cluster} vs (E_{split1}-E_{split2})/(E_{split1}+E_{split2}) for N max  = 1, E>12, %s",sMatched[j].Data()),
-                                                100,-1,1,120,0,1.2); 
-  fhSplitEFractionvsAsyNLocMax1[j]   ->SetXTitle("(E_{split1}-E_{split2})/(E_{split1}+E_{split2})");
-  fhSplitEFractionvsAsyNLocMax1[j]   ->SetYTitle("(E_{split1}+E_{split2})/E_{cluster}");
-  outputContainer->Add(fhSplitEFractionvsAsyNLocMax1[j]) ; 
-  
-  fhSplitEFractionvsAsyNLocMax2[j]     = new TH2F(Form("hSplitEFractionvsAsyNLocMax2%s",sMatched[j].Data()),
-                                                Form("(E1+E2)/E_{cluster} vs (E_{split1}-E_{split2})/(E_{split1}+E_{split2}) for N max  = 2,E>12, %s",sMatched[j].Data()),
-                                                100,-1,1,120,0,1.2); 
-  fhSplitEFractionvsAsyNLocMax2[j]   ->SetXTitle("(E_{split1}-E_{split2})/(E_{split1}+E_{split2})");
-  fhSplitEFractionvsAsyNLocMax2[j]   ->SetYTitle("(E_{split1}+E_{split2})/E_{cluster}");
-  outputContainer->Add(fhSplitEFractionvsAsyNLocMax2[j]) ; 
-  
-  fhSplitEFractionvsAsyNLocMaxN[j]    = new TH2F(Form("hSplitEFractionvsAsyNLocMaxN%s",sMatched[j].Data()),
-                                               Form("(E1+E2)/E_{cluster} vs (E_{split1}-E_{split2})/(E_{split1}+E_{split2}) for N max  > 2, E>12, %s",sMatched[j].Data()),
-                                               100,-1,1,120,0,1.2); 
-  fhSplitEFractionvsAsyNLocMaxN[j]   ->SetXTitle("(E_{split1}-E_{split2})/(E_{split1}+E_{split2})");
-  fhSplitEFractionvsAsyNLocMaxN[j]   ->SetYTitle("(E_{split1}+E_{split2})/E_{cluster}");
-  outputContainer->Add(fhSplitEFractionvsAsyNLocMaxN[j]) ; 
+  {
+    fhSplitEFractionvsAsyNLocMax1[j]     = new TH2F(Form("hSplitEFractionvsAsyNLocMax1%s",sMatched[j].Data()),
+                                                    Form("(E1+E2)/E_{cluster} vs (E_{split1}-E_{split2})/(E_{split1}+E_{split2}) for N max  = 1, E>12, %s",sMatched[j].Data()),
+                                                    100,-1,1,120,0,1.2); 
+    fhSplitEFractionvsAsyNLocMax1[j]   ->SetXTitle("(E_{split1}-E_{split2})/(E_{split1}+E_{split2})");
+    fhSplitEFractionvsAsyNLocMax1[j]   ->SetYTitle("(E_{split1}+E_{split2})/E_{cluster}");
+    outputContainer->Add(fhSplitEFractionvsAsyNLocMax1[j]) ; 
+    
+    fhSplitEFractionvsAsyNLocMax2[j]     = new TH2F(Form("hSplitEFractionvsAsyNLocMax2%s",sMatched[j].Data()),
+                                                    Form("(E1+E2)/E_{cluster} vs (E_{split1}-E_{split2})/(E_{split1}+E_{split2}) for N max  = 2,E>12, %s",sMatched[j].Data()),
+                                                    100,-1,1,120,0,1.2); 
+    fhSplitEFractionvsAsyNLocMax2[j]   ->SetXTitle("(E_{split1}-E_{split2})/(E_{split1}+E_{split2})");
+    fhSplitEFractionvsAsyNLocMax2[j]   ->SetYTitle("(E_{split1}+E_{split2})/E_{cluster}");
+    outputContainer->Add(fhSplitEFractionvsAsyNLocMax2[j]) ; 
+    
+    fhSplitEFractionvsAsyNLocMaxN[j]    = new TH2F(Form("hSplitEFractionvsAsyNLocMaxN%s",sMatched[j].Data()),
+                                                   Form("(E1+E2)/E_{cluster} vs (E_{split1}-E_{split2})/(E_{split1}+E_{split2}) for N max  > 2, E>12, %s",sMatched[j].Data()),
+                                                   100,-1,1,120,0,1.2); 
+    fhSplitEFractionvsAsyNLocMaxN[j]   ->SetXTitle("(E_{split1}-E_{split2})/(E_{split1}+E_{split2})");
+    fhSplitEFractionvsAsyNLocMaxN[j]   ->SetYTitle("(E_{split1}+E_{split2})/E_{cluster}");
+    outputContainer->Add(fhSplitEFractionvsAsyNLocMaxN[j]) ; 
   }
    
   
@@ -1503,7 +1517,7 @@ void  AliAnaInsideClusterInvariantMass::MakeAnalysisFillHistograms()
     cells = GetEMCALCells();
   }
   
-  const Float_t ecut = 12.; // Fixed cut for some histograms
+  const Float_t ecut = 8.; // Fixed cut for some histograms
   
   if(!pl || !cells) 
   {
@@ -1908,11 +1922,9 @@ void  AliAnaInsideClusterInvariantMass::MakeAnalysisFillHistograms()
       else if(pidTag==AliCaloPID::kPi0   ) { fhM02Pi0NLocMax1[0][matched]->Fill(en,l0); fhMassPi0NLocMax1[0][matched]->Fill(en,mass);  fhAsyPi0NLocMax1[0][matched]->Fill(en,asym); }
       else if(pidTag==AliCaloPID::kEta)    { fhM02EtaNLocMax1[0][matched]->Fill(en,l0); fhMassEtaNLocMax1[0][matched]->Fill(en,mass);  fhAsyEtaNLocMax1[0][matched]->Fill(en,asym); }
       
-      if(!matched)
-      {
-        if     (pidTag==AliCaloPID::kPi0) { fhCentralityPi0NLocMax1->Fill(en,cent) ; fhEventPlanePi0NLocMax1->Fill(en,evp) ; }
-        else if(pidTag==AliCaloPID::kEta) { fhCentralityEtaNLocMax1->Fill(en,cent) ; fhEventPlaneEtaNLocMax1->Fill(en,evp) ; }
-      }
+      if     (pidTag==AliCaloPID::kPi0) { fhCentralityPi0NLocMax1[0][matched]->Fill(en,cent) ; if(!matched)fhEventPlanePi0NLocMax1->Fill(en,evp) ; }
+      else if(pidTag==AliCaloPID::kEta) { fhCentralityEtaNLocMax1[0][matched]->Fill(en,cent) ; if(!matched)fhEventPlaneEtaNLocMax1->Fill(en,evp) ; }
+      
     }
     else if(nMax==2) 
     {
@@ -1954,11 +1966,8 @@ void  AliAnaInsideClusterInvariantMass::MakeAnalysisFillHistograms()
       else if(pidTag==AliCaloPID::kPi0   ) { fhM02Pi0NLocMax2[0][matched]->Fill(en,l0); fhMassPi0NLocMax2[0][matched]->Fill(en,mass);  fhAsyPi0NLocMax2[0][matched]->Fill(en,asym); }        
       else if(pidTag==AliCaloPID::kEta)    { fhM02EtaNLocMax2[0][matched]->Fill(en,l0); fhMassEtaNLocMax2[0][matched]->Fill(en,mass);  fhAsyEtaNLocMax2[0][matched]->Fill(en,asym); }
       
-      if(!matched)
-      {
-        if     (pidTag==AliCaloPID::kPi0) { fhCentralityPi0NLocMax2->Fill(en,cent) ; fhEventPlanePi0NLocMax2->Fill(en,evp) ; }
-        else if(pidTag==AliCaloPID::kEta) { fhCentralityEtaNLocMax2->Fill(en,cent) ; fhEventPlaneEtaNLocMax2->Fill(en,evp) ; }
-      }      
+      if     (pidTag==AliCaloPID::kPi0) { fhCentralityPi0NLocMax2[0][matched]->Fill(en,cent) ; if(!matched)fhEventPlanePi0NLocMax2->Fill(en,evp) ; }
+      else if(pidTag==AliCaloPID::kEta) { fhCentralityEtaNLocMax2[0][matched]->Fill(en,cent) ; if(!matched)fhEventPlaneEtaNLocMax2->Fill(en,evp) ; }
     }
     else if(nMax >2) 
     {
@@ -2000,11 +2009,8 @@ void  AliAnaInsideClusterInvariantMass::MakeAnalysisFillHistograms()
       else if(pidTag==AliCaloPID::kPi0   ) { fhM02Pi0NLocMaxN[0][matched]->Fill(en,l0); fhMassPi0NLocMaxN[0][matched]->Fill(en,mass);  fhAsyPi0NLocMaxN[0][matched]->Fill(en,asym); }
       else if(pidTag==AliCaloPID::kEta)    { fhM02EtaNLocMaxN[0][matched]->Fill(en,l0); fhMassEtaNLocMaxN[0][matched]->Fill(en,mass);  fhAsyEtaNLocMaxN[0][matched]->Fill(en,asym); }
       
-      if(!matched)
-      {
-        if     (pidTag==AliCaloPID::kPi0) { fhCentralityPi0NLocMaxN->Fill(en,cent) ; fhEventPlanePi0NLocMaxN->Fill(en,evp) ; }
-        else if(pidTag==AliCaloPID::kEta) { fhCentralityEtaNLocMaxN->Fill(en,cent) ; fhEventPlaneEtaNLocMaxN->Fill(en,evp) ; }
-      }
+      if     (pidTag==AliCaloPID::kPi0) { fhCentralityPi0NLocMaxN[0][matched]->Fill(en,cent) ; if(!matched)fhEventPlanePi0NLocMaxN->Fill(en,evp) ; }
+      else if(pidTag==AliCaloPID::kEta) { fhCentralityEtaNLocMaxN[0][matched]->Fill(en,cent) ; if(!matched)fhEventPlaneEtaNLocMaxN->Fill(en,evp) ; }
       
     }
     
@@ -2025,8 +2031,11 @@ void  AliAnaInsideClusterInvariantMass::MakeAnalysisFillHistograms()
 
         if     (pidTag==AliCaloPID::kPhoton) { fhM02ConNLocMax1[mcindex][matched]->Fill(en,l0); fhMassConNLocMax1[mcindex][matched]->Fill(en,mass); fhAsyConNLocMax1[mcindex][matched]->Fill(en,asym); }
         else if(pidTag==AliCaloPID::kPi0   ) { fhM02Pi0NLocMax1[mcindex][matched]->Fill(en,l0); fhMassPi0NLocMax1[mcindex][matched]->Fill(en,mass); fhAsyPi0NLocMax1[mcindex][matched]->Fill(en,asym); }
-        else if(pidTag==AliCaloPID::kEta   ) { fhM02EtaNLocMax1[mcindex][matched]->Fill(en,l0); fhMassEtaNLocMax1[mcindex][matched]->Fill(en,mass); fhAsyEtaNLocMax1[mcindex][matched]->Fill(en,asym); } 
-      }  
+        else if(pidTag==AliCaloPID::kEta   ) { fhM02EtaNLocMax1[mcindex][matched]->Fill(en,l0); fhMassEtaNLocMax1[mcindex][matched]->Fill(en,mass); fhAsyEtaNLocMax1[mcindex][matched]->Fill(en,asym); }
+        
+        if     (pidTag==AliCaloPID::kPi0) fhCentralityPi0NLocMax1[mcindex][matched]->Fill(en,cent) ;
+        else if(pidTag==AliCaloPID::kEta) fhCentralityEtaNLocMax1[mcindex][matched]->Fill(en,cent) ;
+      }
       else if(nMax==2) 
       {
         fhMassNLocMax2[mcindex][matched]->Fill(en,mass);
@@ -2042,6 +2051,9 @@ void  AliAnaInsideClusterInvariantMass::MakeAnalysisFillHistograms()
         if     (pidTag==AliCaloPID::kPhoton) { fhM02ConNLocMax2[mcindex][matched]->Fill(en,l0); fhMassConNLocMax2[mcindex][matched]->Fill(en,mass); fhAsyConNLocMax2[mcindex][matched]->Fill(en,asym); }
         else if(pidTag==AliCaloPID::kPi0   ) { fhM02Pi0NLocMax2[mcindex][matched]->Fill(en,l0); fhMassPi0NLocMax2[mcindex][matched]->Fill(en,mass); fhAsyPi0NLocMax2[mcindex][matched]->Fill(en,asym); }
         else if(pidTag==AliCaloPID::kEta   ) { fhM02EtaNLocMax2[mcindex][matched]->Fill(en,l0); fhMassEtaNLocMax2[mcindex][matched]->Fill(en,mass); fhAsyEtaNLocMax2[mcindex][matched]->Fill(en,asym); } 
+       
+        if     (pidTag==AliCaloPID::kPi0) fhCentralityPi0NLocMax2[mcindex][matched]->Fill(en,cent) ;
+        else if(pidTag==AliCaloPID::kEta) fhCentralityEtaNLocMax2[mcindex][matched]->Fill(en,cent) ;
         
       }
       else if(nMax >2) 
@@ -2058,7 +2070,10 @@ void  AliAnaInsideClusterInvariantMass::MakeAnalysisFillHistograms()
         
         if     (pidTag==AliCaloPID::kPhoton) { fhM02ConNLocMaxN[mcindex][matched]->Fill(en,l0); fhMassConNLocMaxN[mcindex][matched]->Fill(en,mass); fhAsyConNLocMaxN[mcindex][matched]->Fill(en,asym); }
         else if(pidTag==AliCaloPID::kPi0   ) { fhM02Pi0NLocMaxN[mcindex][matched]->Fill(en,l0); fhMassPi0NLocMaxN[mcindex][matched]->Fill(en,mass); fhAsyPi0NLocMaxN[mcindex][matched]->Fill(en,asym); }
-        else if(pidTag==AliCaloPID::kEta   ) { fhM02EtaNLocMaxN[mcindex][matched]->Fill(en,l0); fhMassEtaNLocMaxN[mcindex][matched]->Fill(en,mass); fhAsyEtaNLocMaxN[mcindex][matched]->Fill(en,asym); } 
+        else if(pidTag==AliCaloPID::kEta   ) { fhM02EtaNLocMaxN[mcindex][matched]->Fill(en,l0); fhMassEtaNLocMaxN[mcindex][matched]->Fill(en,mass); fhAsyEtaNLocMaxN[mcindex][matched]->Fill(en,asym); }
+        
+        if     (pidTag==AliCaloPID::kPi0) fhCentralityPi0NLocMaxN[mcindex][matched]->Fill(en,cent) ;
+        else if(pidTag==AliCaloPID::kEta) fhCentralityEtaNLocMaxN[mcindex][matched]->Fill(en,cent) ;
       }
       
     }//Work with MC truth first

diff --git a/PWGGA/CaloTrackCorrelations/AliAnaInsideClusterInvariantMass.h b/PWGGA/CaloTrackCorrelations/AliAnaInsideClusterInvariantMass.h
index 53fee334738d04a88319e06c7f5691e78af068b7..a660507d9a35d8dd71437e36232cbc36db7512ab 100755 (executable)
--- a/PWGGA/CaloTrackCorrelations/AliAnaInsideClusterInvariantMass.h
+++ b/PWGGA/CaloTrackCorrelations/AliAnaInsideClusterInvariantMass.h
@@ -263,12 +263,12 @@ class AliAnaInsideClusterInvariantMass : public AliAnaCaloTrackCorrBaseClass {
   TH2F       * fhTrackMatchedDEtaNLocMaxN[8] ;          //! Eta distance between track and cluster vs cluster E, more than 2 local maximum
   TH2F       * fhTrackMatchedDPhiNLocMaxN[8] ;          //! Phi distance between track and cluster vs cluster E, more than 2 local maximum
   
-  TH2F       * fhCentralityPi0NLocMax1 ;                //! Centrality for selected pi0, N Local Maxima = 1
-  TH2F       * fhCentralityEtaNLocMax1 ;                //! Centrality for selected eta, N Local Maxima = 1
-  TH2F       * fhCentralityPi0NLocMax2 ;                //! Centrality for selected pi0, N Local Maxima = 2
-  TH2F       * fhCentralityEtaNLocMax2 ;                //! Centrality for selected eta, N Local Maxima = 2
-  TH2F       * fhCentralityPi0NLocMaxN ;                //! Centrality for selected pi0, N Local Maxima > 2
-  TH2F       * fhCentralityEtaNLocMaxN ;                //! Centrality for selected eta, N Local Maxima > 2
+  TH2F       * fhCentralityPi0NLocMax1[8][2] ;          //! Centrality for selected pi0, N Local Maxima = 1
+  TH2F       * fhCentralityEtaNLocMax1[8][2] ;          //! Centrality for selected eta, N Local Maxima = 1
+  TH2F       * fhCentralityPi0NLocMax2[8][2] ;          //! Centrality for selected pi0, N Local Maxima = 2
+  TH2F       * fhCentralityEtaNLocMax2[8][2] ;          //! Centrality for selected eta, N Local Maxima = 2
+  TH2F       * fhCentralityPi0NLocMaxN[8][2] ;          //! Centrality for selected pi0, N Local Maxima > 2
+  TH2F       * fhCentralityEtaNLocMaxN[8][2] ;          //! Centrality for selected eta, N Local Maxima > 2
 
   TH2F       * fhEventPlanePi0NLocMax1 ;                //! Event plane for selected pi0, N Local Maxima = 1
   TH2F       * fhEventPlaneEtaNLocMax1 ;                //! Event plane for selected eta, N Local Maxima = 1

diff --git a/PWGGA/CaloTrackCorrelations/AliAnaPi0EbE.cxx b/PWGGA/CaloTrackCorrelations/AliAnaPi0EbE.cxx
index f4e8e04d027d3e96241d6285fbe6e7b5e88e3944..1828d95d3ed72318652f9347ec4262a7cc2de8b3 100755 (executable)
--- a/PWGGA/CaloTrackCorrelations/AliAnaPi0EbE.cxx
+++ b/PWGGA/CaloTrackCorrelations/AliAnaPi0EbE.cxx
@@ -60,7 +60,7 @@ AliAnaPi0EbE::AliAnaPi0EbE() :
     fhPt(0),                       fhE(0),                    
     fhEEta(0),                     fhEPhi(0),                    fhEtaPhi(0),
     fhPtCentrality(),              fhPtEventPlane(0),
-    fhPtReject(0),                 fhEReject(0),                    
+    fhPtReject(0),                 fhEReject(0),
     fhEEtaReject(0),               fhEPhiReject(0),              fhEtaPhiReject(0),
     fhMass(0),                     fhAsymmetry(0), 
     fhSelectedMass(0),             fhSelectedAsymmetry(0),
@@ -77,7 +77,8 @@ AliAnaPi0EbE::AliAnaPi0EbE() :
     // MC histos
     fhMCE(),                       fhMCPt(),        
     fhMCPhi(),                     fhMCEta(),
-    fhMCEReject(),                 fhMCPtReject(),       
+    fhMCEReject(),                 fhMCPtReject(),
+    fhMCPtCentrality(),            
     fhMCPi0PtGenRecoFraction(0),   fhMCEtaPtGenRecoFraction(0),
     fhMCPi0DecayPt(0),             fhMCPi0DecayPtFraction(0),      
     fhMCEtaDecayPt(0),             fhMCEtaDecayPtFraction(0),
@@ -106,6 +107,8 @@ AliAnaPi0EbE::AliAnaPi0EbE() :
     fhMCPt             [i] = 0;
     fhMCPhi            [i] = 0;                   
     fhMCEta            [i] = 0;
+    fhMCPtCentrality   [i] = 0;
+    
     fhEMCLambda0       [i] = 0;
     fhEMCLambda0NoTRD  [i] = 0;
     fhEMCLambda0FracMaxCellCut[i]= 0;
@@ -1158,6 +1161,15 @@ TList *  AliAnaPi0EbE::GetCreateOutputObjects()
         fhMCPt[i]->SetXTitle("p_{T} (GeV/c)");
         outputContainer->Add(fhMCPt[i]) ; 
         
+        fhMCPtCentrality[i]  = new TH2F
+        (Form("hPtCentrality_MC%s",pname[i].Data()),
+         Form("Identified as #pi^{0} (#eta), cluster from %s",
+              ptype[i].Data()),
+         nptbins,ptmin,ptmax, 100,0,100);
+        fhMCPtCentrality[i]->SetYTitle("centrality");
+        fhMCPtCentrality[i]->SetXTitle("p_{T} (GeV/c)");
+        outputContainer->Add(fhMCPtCentrality[i]) ;
+        
         if(fAnaType == kSSCalo)
         {
           fhMCEReject[i]  = new TH1F
@@ -2236,6 +2248,8 @@ void  AliAnaPi0EbE::MakeAnalysisFillHistograms()
       fhMCPhi[mcIndex] ->Fill(pt,phi);
       fhMCEta[mcIndex] ->Fill(pt,eta);
       
+      fhMCPtCentrality[mcIndex]->Fill(pt,cen);
+
       if((mcIndex==kmcPhoton || mcIndex==kmcPi0 || mcIndex==kmcEta) && fAnaType==kSSCalo)
       {
         Float_t efracMC = 0;

diff --git a/PWGGA/CaloTrackCorrelations/AliAnaPi0EbE.h b/PWGGA/CaloTrackCorrelations/AliAnaPi0EbE.h
index b28210c341cae0f398b7bb9f5b1a57e286d17c47..134b18ad1086a56fe5817fe21551c6d002ec26aa 100755 (executable)
--- a/PWGGA/CaloTrackCorrelations/AliAnaPi0EbE.h
+++ b/PWGGA/CaloTrackCorrelations/AliAnaPi0EbE.h
@@ -215,6 +215,9 @@ class AliAnaPi0EbE : public AliAnaCaloTrackCorrBaseClass {
   TH1F         * fhMCEReject[6];              //! Number of rejected as pi0 vs E coming from X
   TH1F         * fhMCPtReject[6];             //! Number of rejected as pi0 vs Pt coming from X
 
+  TH2F         * fhMCPtCentrality[6] ;        //! centrality  vs pi0/eta pT  coming from X
+
+  
   TH2F         * fhMCPi0PtGenRecoFraction;    //! SS id, clusters id as pi0 (eta), coming from 2 photon, pi0 primary, pt vs E prim pi0 / E reco
   TH2F         * fhMCEtaPtGenRecoFraction;    //! SS id, clusters id as pi0 (eta), coming from 2 photon, eta primary, pt vs E prim eta / E reco  
   TH1F         * fhMCPi0DecayPt;              //! SS id, clusters id as pi0 (eta), coming from 1 photon, pi0 decay primary, pt