New histogram for D vs. multiplicity studies (Aamer)

diff --git a/PWGHF/vertexingHF/AliAnalysisTaskSEDvsMultiplicity.cxx b/PWGHF/vertexingHF/AliAnalysisTaskSEDvsMultiplicity.cxx
index 869f660ba3bfed399e31a05ffe07339e5ef24ac6..5b73c6ebbf8e348a6625892589fe472440ec5296 100644 (file)
--- a/PWGHF/vertexingHF/AliAnalysisTaskSEDvsMultiplicity.cxx
+++ b/PWGHF/vertexingHF/AliAnalysisTaskSEDvsMultiplicity.cxx
@@ -69,6 +69,7 @@ AliAnalysisTaskSE(),
   fHistNtrVsNchMCPhysicalPrimary(0),
   fHistNtrCorrVsNchMCPhysicalPrimary(0),
   fHistGenPrimaryParticlesInelGt0(0),
+  fHistNchMCVsNchMCPrimaryVsNchMCPhysicalPrimary(0),
   fHistNtrCorrEvSel(0),
   fHistNtrCorrEvWithCand(0),
   fHistNtrCorrEvWithD(0),
@@ -117,6 +118,7 @@ AliAnalysisTaskSEDvsMultiplicity::AliAnalysisTaskSEDvsMultiplicity(const char *n
   fHistNtrVsNchMCPhysicalPrimary(0),
   fHistNtrCorrVsNchMCPhysicalPrimary(0),
   fHistGenPrimaryParticlesInelGt0(0),
+  fHistNchMCVsNchMCPrimaryVsNchMCPhysicalPrimary(0),
   fHistNtrCorrEvSel(0),
   fHistNtrCorrEvWithCand(0),
   fHistNtrCorrEvWithD(0),
@@ -268,6 +270,8 @@ void AliAnalysisTaskSEDvsMultiplicity::UserCreateOutputObjects()
   
   fHistGenPrimaryParticlesInelGt0 = new TH1F("hGenPrimaryParticlesInelGt0","Multiplcity of generated charged particles ; Nparticles ; Entries",200,-0.5,199.5);
 
+  fHistNchMCVsNchMCPrimaryVsNchMCPhysicalPrimary = new TH3F("fHistNchMCVsNchMCPrimaryVsNchMCPhysicalPrimary", "MC: Nch (Physical Primary) vs Nch (Primary) vs Nch (Generated); Nch (Generated); Nch (Primary); Nch (Physical Primary)",200,0.,200.,200,0.,200.,200,0.,200.);
+
   fHistNtrCorrEvSel->Sumw2();
   fHistNtrCorrEvWithCand->Sumw2();
   fHistNtrCorrEvWithD->Sumw2();
@@ -287,6 +291,7 @@ void AliAnalysisTaskSEDvsMultiplicity::UserCreateOutputObjects()
   fOutput->Add(fHistNtrVsNchMCPhysicalPrimary);
   fOutput->Add(fHistNtrCorrVsNchMCPhysicalPrimary);
   fOutput->Add(fHistGenPrimaryParticlesInelGt0);
+  fOutput->Add(fHistNchMCVsNchMCPrimaryVsNchMCPhysicalPrimary);
 
   
   fHistNEvents = new TH1F("fHistNEvents", "number of events ",11,-0.5,10.5);
@@ -482,6 +487,7 @@ void AliAnalysisTaskSEDvsMultiplicity::UserExec(Option_t */*option*/)
     fHistNtrVsNchMCPhysicalPrimary->Fill(nChargedMCPhysicalPrimary,countTreta1);
     fHistNtrCorrVsNchMCPhysicalPrimary->Fill(nChargedMCPhysicalPrimary,countTreta1corr);
 
+    fHistNchMCVsNchMCPrimaryVsNchMCPhysicalPrimary->Fill(nChargedMC,nChargedMCPrimary,nChargedMCPhysicalPrimary);
   }
   
   Int_t nCand = arrayCand->GetEntriesFast(); 

diff --git a/PWGHF/vertexingHF/AliAnalysisTaskSEDvsMultiplicity.h b/PWGHF/vertexingHF/AliAnalysisTaskSEDvsMultiplicity.h
index 59bef6dca1cb6c88efa914c8e76340724d578269..be601b34a6088efda59a99fbf07f4bf2969d40e7 100644 (file)
--- a/PWGHF/vertexingHF/AliAnalysisTaskSEDvsMultiplicity.h
+++ b/PWGHF/vertexingHF/AliAnalysisTaskSEDvsMultiplicity.h
@@ -114,6 +114,7 @@ class AliAnalysisTaskSEDvsMultiplicity : public AliAnalysisTaskSE
   TH2F* fHistNtrVsNchMCPhysicalPrimary; //!  hist of ntracklets vs Nch (Physical Primary)
   TH2F* fHistNtrCorrVsNchMCPhysicalPrimary; //!  hist of ntracklets vs Nch (Physical Primary)
   TH1F* fHistGenPrimaryParticlesInelGt0; //!hist. of geenrated multiplcity
+  TH3F* fHistNchMCVsNchMCPrimaryVsNchMCPhysicalPrimary; //! hist of Nch (generated) vs Nch (Primary) vs Nch (Physical Primary) 
   
   TH1F* fHistNtrCorrEvSel; //! hist. of ntracklets for selected events
   TH1F* fHistNtrCorrEvWithCand; //! hist. of ntracklets for evnts with a candidate
@@ -150,7 +151,7 @@ class AliAnalysisTaskSEDvsMultiplicity : public AliAnalysisTaskSE
   Int_t fPdgMeson;   // pdg code of analyzed meson
 
   
-  ClassDef(AliAnalysisTaskSEDvsMultiplicity,3); // D vs. mult task
+  ClassDef(AliAnalysisTaskSEDvsMultiplicity,4); // D vs. mult task
 };
 
 #endif