New histogram for D vs. multiplicity studies (Aamer)
authorprino <prino@f7af4fe6-9843-0410-8265-dc069ae4e863>
Mon, 11 Jun 2012 16:44:36 +0000 (16:44 +0000)
committerprino <prino@f7af4fe6-9843-0410-8265-dc069ae4e863>
Mon, 11 Jun 2012 16:44:36 +0000 (16:44 +0000)
PWGHF/vertexingHF/AliAnalysisTaskSEDvsMultiplicity.cxx
PWGHF/vertexingHF/AliAnalysisTaskSEDvsMultiplicity.h

index 869f660..5b73c6e 100644 (file)
@@ -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(); 
index 59bef6d..be601b3 100644 (file)
@@ -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