add MC histogram for decay pi0 with lost companion

diff --git a/PWGGA/CaloTrackCorrelations/AliAnaParticleIsolation.cxx b/PWGGA/CaloTrackCorrelations/AliAnaParticleIsolation.cxx
index b1748a95c5d9abbf6097ae4681f14011c74e3517..d11da923a06ead54668ad52892fca8f015532589 100755 (executable)
--- a/PWGGA/CaloTrackCorrelations/AliAnaParticleIsolation.cxx
+++ b/PWGGA/CaloTrackCorrelations/AliAnaParticleIsolation.cxx
@@ -1245,6 +1245,9 @@ void AliAnaParticleIsolation::FillTrackMatchingShowerShapeControlHistograms(AliA
       if(GetMCAnalysisUtils()->CheckTagBit(mcTag,AliMCAnalysisUtils::kMCPhoton))
         fhPtLambda0MC[kmcPhoton][isolated]->Fill(pt,m02);
       
+      if(GetMCAnalysisUtils()->CheckTagBit(mcTag,AliMCAnalysisUtils::kMCDecayPairInCalo) && mcIndex==kmcPi0Decay )
+        fhPtLambda0MC[kmcPi0DecayLostPair][isolated]->Fill(pt,m02);
+      
       fhPtLambda0MC[mcIndex][isolated]->Fill(pt,m02);
     }
     
@@ -1475,12 +1478,12 @@ TList *  AliAnaParticleIsolation::GetCreateOutputObjects()
   
   // MC histograms title and name
   TString mcPartType[] = { "#gamma", "#gamma_{prompt}", "#gamma_{fragmentation}",
-    "#pi^{0} (merged #gamma)","#gamma_{#pi decay}",
+    "#pi^{0} (merged #gamma)","#gamma_{#pi decay}","#gamma_{#pi decay} lost companion",
     "#gamma_{#eta decay}","#gamma_{other decay}",
     "e^{#pm}","hadrons?"} ;
   
   TString mcPartName[] = { "Photon","PhotonPrompt","PhotonFrag",
-    "Pi0","Pi0Decay","EtaDecay","OtherDecay",
+    "Pi0","Pi0Decay","Pi0DecayLostPair","EtaDecay","OtherDecay",
     "Electron","Hadron"} ;
   
   // Primary MC histograms title and name
@@ -1519,7 +1522,6 @@ TList *  AliAnaParticleIsolation::GetCreateOutputObjects()
     
     for(Int_t imc = 0; imc < fgkNmcTypes; imc++)
     {
-      
       fhPtNoIsoMC[imc]  = new TH1F(Form("hPtNoIsoMC%s",mcPartName[imc].Data()),
                                    Form("#it{p}_{T} of NOT isolated %s, %s",mcPartType[imc].Data(),parTitle.Data()),
                                    nptbins,ptmin,ptmax);
@@ -3743,6 +3745,13 @@ void  AliAnaParticleIsolation::MakeAnalysisFillHistograms()
           fhEtaIsoMC[kmcPhoton]->Fill(pt,eta);
         }
         
+        if(GetMCAnalysisUtils()->CheckTagBit(mcTag,AliMCAnalysisUtils::kMCDecayPairInCalo) && mcIndex==kmcPi0Decay )
+        {
+          fhPtIsoMC [kmcPi0DecayLostPair]->Fill(pt);
+          fhPhiIsoMC[kmcPi0DecayLostPair]->Fill(pt,phi);
+          fhEtaIsoMC[kmcPi0DecayLostPair]->Fill(pt,eta);
+        }
+        
         fhPtIsoMC [mcIndex]->Fill(pt);
         fhPhiIsoMC[mcIndex]->Fill(pt,phi);
         fhEtaIsoMC[mcIndex]->Fill(pt,eta);
@@ -3763,6 +3772,9 @@ void  AliAnaParticleIsolation::MakeAnalysisFillHistograms()
               if(GetMCAnalysisUtils()->CheckTagBit(mcTag,AliMCAnalysisUtils::kMCPhoton))
                 fhPtDecayIsoMC[ibit][kmcPhoton]->Fill(pt);
               
+              if(GetMCAnalysisUtils()->CheckTagBit(mcTag,AliMCAnalysisUtils::kMCDecayPairInCalo) && mcIndex==kmcPi0Decay )
+                fhPtDecayIsoMC[ibit][kmcPi0DecayLostPair]->Fill(pt);
+              
               fhPtDecayIsoMC[ibit][mcIndex]->Fill(pt);
             }
           } // bit ok
@@ -3808,6 +3820,9 @@ void  AliAnaParticleIsolation::MakeAnalysisFillHistograms()
         if( GetMCAnalysisUtils()->CheckTagBit(mcTag,AliMCAnalysisUtils::kMCPhoton) )
           fhPtNoIsoMC[kmcPhoton]->Fill(pt);
         
+        if(GetMCAnalysisUtils()->CheckTagBit(mcTag,AliMCAnalysisUtils::kMCDecayPairInCalo) && mcIndex==kmcPi0Decay )
+          fhPtNoIsoMC[kmcPi0DecayLostPair]->Fill(pt);
+
         fhPtNoIsoMC[mcIndex]->Fill(pt);
       }
       
@@ -3826,6 +3841,9 @@ void  AliAnaParticleIsolation::MakeAnalysisFillHistograms()
               if(GetMCAnalysisUtils()->CheckTagBit(mcTag,AliMCAnalysisUtils::kMCPhoton))
                 fhPtDecayNoIsoMC[ibit][kmcPhoton]->Fill(pt);
               
+              if(GetMCAnalysisUtils()->CheckTagBit(mcTag,AliMCAnalysisUtils::kMCDecayPairInCalo) && mcIndex==kmcPi0Decay )
+                fhPtDecayNoIsoMC[ibit][kmcPi0DecayLostPair]->Fill(pt);
+              
               fhPtDecayNoIsoMC[ibit][mcIndex]->Fill(pt);
             }
           } // bit ok
@@ -4354,6 +4372,9 @@ void  AliAnaParticleIsolation::MakeSeveralICAnalysis(AliAODPWG4ParticleCorrelati
     if(GetMCAnalysisUtils()->CheckTagBit(tag,AliMCAnalysisUtils::kMCPhoton))
       fhPtNoIsoMC[kmcPhoton]->Fill(ptC);
     
+    if(GetMCAnalysisUtils()->CheckTagBit(tag,AliMCAnalysisUtils::kMCDecayPairInCalo) && mcIndex==kmcPi0Decay )
+      fhPtNoIsoMC[kmcPi0DecayLostPair]->Fill(ptC);
+    
     fhPtNoIsoMC[mcIndex]->Fill(ptC);
   }
   
@@ -4371,6 +4392,9 @@ void  AliAnaParticleIsolation::MakeSeveralICAnalysis(AliAODPWG4ParticleCorrelati
         {
           if(GetMCAnalysisUtils()->CheckTagBit(tag,AliMCAnalysisUtils::kMCPhoton))
             fhPtDecayNoIsoMC[ibit][kmcPhoton]->Fill(ptC);
+
+          if(GetMCAnalysisUtils()->CheckTagBit(tag,AliMCAnalysisUtils::kMCDecayPairInCalo) && mcIndex==kmcPi0Decay )
+            fhPtDecayNoIsoMC[ibit][kmcPi0DecayLostPair]->Fill(ptC);
           
           fhPtDecayNoIsoMC[ibit][mcIndex]->Fill(ptC);
         }
@@ -4481,6 +4505,9 @@ void  AliAnaParticleIsolation::MakeSeveralICAnalysis(AliAODPWG4ParticleCorrelati
       if( GetMCAnalysisUtils()->CheckTagBit(tag,AliMCAnalysisUtils::kMCPhoton))
         fhSumPtLeadingPtMC[kmcPhoton][icone]->Fill(ptC,coneptsum) ;
       
+      if(GetMCAnalysisUtils()->CheckTagBit(tag,AliMCAnalysisUtils::kMCDecayPairInCalo) && mcIndex==kmcPi0Decay )
+        fhSumPtLeadingPtMC[kmcPi0DecayLostPair][icone]->Fill(ptC,coneptsum) ;
+      
       fhSumPtLeadingPtMC[mcIndex][icone]->Fill(ptC,coneptsum) ;
     }
     
@@ -4535,6 +4562,9 @@ void  AliAnaParticleIsolation::MakeSeveralICAnalysis(AliAODPWG4ParticleCorrelati
           if( GetMCAnalysisUtils()->CheckTagBit(tag,AliMCAnalysisUtils::kMCPhoton) )
             fhPtThresIsolatedMC[kmcPhoton][icone][ipt]->Fill(ptC) ;
           
+          if(GetMCAnalysisUtils()->CheckTagBit(tag,AliMCAnalysisUtils::kMCDecayPairInCalo) && mcIndex==kmcPi0Decay )
+            fhPtThresIsolatedMC[kmcPi0DecayLostPair][icone][ipt]->Fill(ptC) ;
+
           fhPtThresIsolatedMC[mcIndex][icone][ipt]->Fill(ptC) ;
           
         }
@@ -4563,6 +4593,9 @@ void  AliAnaParticleIsolation::MakeSeveralICAnalysis(AliAODPWG4ParticleCorrelati
           if( GetMCAnalysisUtils()->CheckTagBit(tag,AliMCAnalysisUtils::kMCPhoton))
             fhPtFracIsolatedMC[kmcPhoton][icone][ipt]->Fill(ptC) ;
           
+          if(GetMCAnalysisUtils()->CheckTagBit(tag,AliMCAnalysisUtils::kMCDecayPairInCalo) && mcIndex==kmcPi0Decay )
+            fhPtFracIsolatedMC[kmcPi0DecayLostPair][icone][ipt]->Fill(ptC) ;
+          
           fhPtFracIsolatedMC[mcIndex][icone][ipt]->Fill(ptC) ;
         }
       }
@@ -4593,6 +4626,9 @@ void  AliAnaParticleIsolation::MakeSeveralICAnalysis(AliAODPWG4ParticleCorrelati
           if( GetMCAnalysisUtils()->CheckTagBit(tag,AliMCAnalysisUtils::kMCPhoton))
             fhSumPtIsolatedMC[kmcPhoton][icone][ipt]->Fill(ptC) ;
           
+          if(GetMCAnalysisUtils()->CheckTagBit(tag,AliMCAnalysisUtils::kMCDecayPairInCalo) && mcIndex==kmcPi0Decay )
+            fhSumPtIsolatedMC[kmcPi0DecayLostPair][icone][ipt]->Fill(ptC) ;
+          
           fhSumPtIsolatedMC[mcIndex][icone][ipt]->Fill(ptC) ;
         }
       }

diff --git a/PWGGA/CaloTrackCorrelations/AliAnaParticleIsolation.h b/PWGGA/CaloTrackCorrelations/AliAnaParticleIsolation.h
index 761fee203cefc57b87c731f39cb90d5e3acf8fbd..2995ac81e015ed4b76450cce394e10f44e8a3403 100755 (executable)
--- a/PWGGA/CaloTrackCorrelations/AliAnaParticleIsolation.h
+++ b/PWGGA/CaloTrackCorrelations/AliAnaParticleIsolation.h
@@ -156,10 +156,11 @@ class AliAnaParticleIsolation : public AliAnaCaloTrackCorrBaseClass {
   static const Int_t fgkNmcPrimTypes = 7;
   
   // For histograms in arrays, index in the array, corresponding to any particle origin
-  enum mcTypes     { kmcPhoton   = 0, kmcPrompt   = 1, kmcFragment = 2,
-                     kmcPi0      = 3, kmcPi0Decay = 4, kmcEtaDecay = 5, kmcOtherDecay = 6,
-                     kmcElectron = 7, kmcHadron   = 8                                     } ;
-  static const Int_t fgkNmcTypes = 9;
+  enum mcTypes     { kmcPhoton   = 0, kmcPrompt     = 1, kmcFragment         = 2,
+                     kmcPi0      = 3, kmcPi0Decay   = 4, kmcPi0DecayLostPair = 5,
+                     kmcEtaDecay = 6, kmcOtherDecay = 7,
+                     kmcElectron = 8, kmcHadron     = 9                          } ;
+  static const Int_t fgkNmcTypes = 10;
 
  private: