Added histograms for decay particles (Martin)

diff --git a/PWG4/GammaConv/AliAnalysisTaskGammaConversion.cxx b/PWG4/GammaConv/AliAnalysisTaskGammaConversion.cxx
index 5ec78d5ef98de0a9a79ad8a46ead069942ab880d..c0730a1bd4e8d9a702641ec035581e2061733016 100644 (file)
--- a/PWG4/GammaConv/AliAnalysisTaskGammaConversion.cxx
+++ b/PWG4/GammaConv/AliAnalysisTaskGammaConversion.cxx
@@ -777,7 +777,31 @@ void AliAnalysisTaskGammaConversion::ProcessMCData(){
       if(particle->GetMother(0) >= fStack->GetNprimary()){
        continue; // the gamma has a mother, and it is not a primary particle
       }
-                       
+               
+      if(particle->GetMother(0) >-1){
+        fHistograms->FillHistogram("MC_DecayAllGamma_Pt", particle->Pt()); // All
+        switch(fStack->Particle(particle->GetMother(0))->GetPdgCode()){
+        case 111: // Pi0
+           fHistograms->FillHistogram("MC_DecayPi0Gamma_Pt", particle->Pt());
+           break;
+        case 113: // Rho0
+           fHistograms->FillHistogram("MC_DecayRho0Gamma_Pt", particle->Pt());
+           break;
+        case 221: // Eta
+           fHistograms->FillHistogram("MC_DecayEtaGamma_Pt", particle->Pt());
+           break;
+        case 223: // Omega
+           fHistograms->FillHistogram("MC_DecayOmegaGamma_Pt", particle->Pt());
+           break;
+        case 310: // K_s0
+           fHistograms->FillHistogram("MC_DecayK0sGamma_Pt", particle->Pt());
+           break;
+        case 331: // Eta'
+           fHistograms->FillHistogram("MC_DecayEtapGamma_Pt", particle->Pt());
+           break;
+        }
+      }
+       
       fHistograms->FillHistogram("MC_allGamma_Energy", particle->Energy());
       fHistograms->FillHistogram("MC_allGamma_Pt", particle->Pt());
       fHistograms->FillHistogram("MC_allGamma_Eta", particle->Eta());

diff --git a/PWG4/macros/ConfigGammaConversion.C b/PWG4/macros/ConfigGammaConversion.C
index f188714656047b83da70380eec636c59ca8e8c48..2deb431ba20e2426ff0dea76cdb523a9cced256c 100644 (file)
--- a/PWG4/macros/ConfigGammaConversion.C
+++ b/PWG4/macros/ConfigGammaConversion.C
@@ -150,6 +150,9 @@ Bool_t kGCForceAOD = kFALSE;  // Call AliAnalysisManager::SetFillAOD(kTRUE) ever
 
 Bool_t kGCplotEventQuality              = kTRUE;
 
+
+// MC GAMMA DECAY PLOTS
+Bool_t kGCplotMCGammaFromDecay          = kTRUE;
 // NEUTRAL MESON PLOTS
 Bool_t kGCplotMCConversionR             = kTRUE;
 Bool_t kGCplotMCConversionZR            = kTRUE;
@@ -1776,6 +1779,18 @@ void AddHistograms(AliGammaConversionHistograms *histograms){
 
   }// end kGCrunChic
        
+
+  //---------------------------------------------- Gamma from Decay ----------------------------------------------------
+  if(kGCplotMCGammaFromDecay == kTRUE && kGCdoMCTruth == kTRUE){
+     histograms->AddHistogram("MC_DecayPi0Gamma_Pt" ,"" , kGCnXBinsPt, kGCfirstXBinPt, kGClastXBinPt, "", "");  
+     histograms->AddHistogram("MC_DecayRho0Gamma_Pt" ,"" , kGCnXBinsPt, kGCfirstXBinPt, kGClastXBinPt, "", "");
+     histograms->AddHistogram("MC_DecayEtaGamma_Pt" ,"" , kGCnXBinsPt, kGCfirstXBinPt, kGClastXBinPt, "", "");
+     histograms->AddHistogram("MC_DecayOmegaGamma_Pt" ,"" , kGCnXBinsPt, kGCfirstXBinPt, kGClastXBinPt, "", "");
+     histograms->AddHistogram("MC_DecayK0sGamma_Pt" ,"" , kGCnXBinsPt, kGCfirstXBinPt, kGClastXBinPt, "", "");
+     histograms->AddHistogram("MC_DecayEtapGamma_Pt" ,"" , kGCnXBinsPt, kGCfirstXBinPt, kGClastXBinPt, "", "");
+     histograms->AddHistogram("MC_DecayAllGamma_Pt" ,"" , kGCnXBinsPt, kGCfirstXBinPt, kGClastXBinPt, "", "");
+  }
+
   //---------------------------------------------- Neutral Meson ---------------------------------------------------------
   if(kGCrunNeutralMeson){