fHistNtrVsNchMCPhysicalPrimary(0),
fHistNtrCorrVsNchMCPhysicalPrimary(0),
fHistGenPrimaryParticlesInelGt0(0),
+ fHistNchMCVsNchMCPrimaryVsNchMCPhysicalPrimary(0),
fHistNtrCorrEvSel(0),
fHistNtrCorrEvWithCand(0),
fHistNtrCorrEvWithD(0),
fHistNtrVsNchMCPhysicalPrimary(0),
fHistNtrCorrVsNchMCPhysicalPrimary(0),
fHistGenPrimaryParticlesInelGt0(0),
+ fHistNchMCVsNchMCPrimaryVsNchMCPhysicalPrimary(0),
fHistNtrCorrEvSel(0),
fHistNtrCorrEvWithCand(0),
fHistNtrCorrEvWithD(0),
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();
fOutput->Add(fHistNtrVsNchMCPhysicalPrimary);
fOutput->Add(fHistNtrCorrVsNchMCPhysicalPrimary);
fOutput->Add(fHistGenPrimaryParticlesInelGt0);
+ fOutput->Add(fHistNchMCVsNchMCPrimaryVsNchMCPhysicalPrimary);
fHistNEvents = new TH1F("fHistNEvents", "number of events ",11,-0.5,10.5);
fHistNtrVsNchMCPhysicalPrimary->Fill(nChargedMCPhysicalPrimary,countTreta1);
fHistNtrCorrVsNchMCPhysicalPrimary->Fill(nChargedMCPhysicalPrimary,countTreta1corr);
+ fHistNchMCVsNchMCPrimaryVsNchMCPhysicalPrimary->Fill(nChargedMC,nChargedMCPrimary,nChargedMCPhysicalPrimary);
}
Int_t nCand = arrayCand->GetEntriesFast();
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
Int_t fPdgMeson; // pdg code of analyzed meson
- ClassDef(AliAnalysisTaskSEDvsMultiplicity,3); // D vs. mult task
+ ClassDef(AliAnalysisTaskSEDvsMultiplicity,4); // D vs. mult task
};
#endif