Added centrality histo for accepted events

diff --git a/PWGLF/RESONANCES/AliRsnMiniAnalysisTask.cxx b/PWGLF/RESONANCES/AliRsnMiniAnalysisTask.cxx
index be893b6313220d59d7ff22bce13acede2b22f1d3..e9bad3b0596e0ed74fdd2a72149102213462c383 100644 (file)
--- a/PWGLF/RESONANCES/AliRsnMiniAnalysisTask.cxx
+++ b/PWGLF/RESONANCES/AliRsnMiniAnalysisTask.cxx
@@ -57,6 +57,7 @@ AliRsnMiniAnalysisTask::AliRsnMiniAnalysisTask() :
    fHistograms("AliRsnMiniOutput", 0),
    fValues("AliRsnMiniValue", 0),
    fHEventStat(0x0),
+   fHAEventsVsMulti(0x0),
    fEventCuts(0x0),
    fTrackCuts(0),
    fRsnEvent(),
@@ -89,6 +90,7 @@ AliRsnMiniAnalysisTask::AliRsnMiniAnalysisTask(const char *name, Bool_t useMC) :
    fHistograms("AliRsnMiniOutput", 0),
    fValues("AliRsnMiniValue", 0),
    fHEventStat(0x0),
+   fHAEventsVsMulti(0x0),
    fEventCuts(0x0),
    fTrackCuts(0),
    fRsnEvent(),
@@ -126,6 +128,7 @@ AliRsnMiniAnalysisTask::AliRsnMiniAnalysisTask(const AliRsnMiniAnalysisTask &cop
    fHistograms(copy.fHistograms),
    fValues(copy.fValues),
    fHEventStat(0x0),
+   fHAEventsVsMulti(0x0),
    fEventCuts(copy.fEventCuts),
    fTrackCuts(copy.fTrackCuts),
    fRsnEvent(),
@@ -166,6 +169,8 @@ AliRsnMiniAnalysisTask &AliRsnMiniAnalysisTask::operator=(const AliRsnMiniAnalys
    fMaxDiffAngle = copy.fMaxDiffAngle;
    fHistograms = copy.fHistograms;
    fValues = copy.fValues;
+   fHEventStat = copy.fHEventStat;
+   fHAEventsVsMulti = copy.fHAEventsVsMulti;
    fEventCuts = copy.fEventCuts;
    fTrackCuts = copy.fTrackCuts;
    fTriggerAna = copy.fTriggerAna;
@@ -249,6 +254,13 @@ void AliRsnMiniAnalysisTask::UserCreateOutputObjects()
    fHEventStat->GetXaxis()->SetBinLabel(4, "Accepted");
    fOutput->Add(fHEventStat);
 
+   if (fUseCentrality)
+      fHAEventsVsMulti = new TH1F("hAEventsVsMulti", "Accepted events vs Centrality", 100, 0, 100.0);
+   else
+      fHAEventsVsMulti = new TH1F("hAEventsVsMulti", "Accepted events vs Multiplicity",1000, 0, 1000.0);
+   
+   fOutput->Add(fHAEventsVsMulti);
+   
    TIter next(&fTrackCuts);
    AliRsnCutSet *cs;
    while ((cs = (AliRsnCutSet *) next())) {
@@ -670,6 +682,8 @@ Char_t AliRsnMiniAnalysisTask::CheckCurrentEvent()
    AliDebugClass(2, Form("Stats: %s", msg.Data()));
    if (isSelected) {
       fHEventStat->Fill(3.1);
+      Double_t multi = ComputeCentrality((output == 'E'));
+      fHAEventsVsMulti->Fill(multi);
       return output;
    } else {
       return 0;

diff --git a/PWGLF/RESONANCES/AliRsnMiniAnalysisTask.h b/PWGLF/RESONANCES/AliRsnMiniAnalysisTask.h
index 0f9be88d5d334925189ba26937a09cf63d7fb048..46b9033b3d94cfb11544d8ec59090c70704dbf91 100644 (file)
--- a/PWGLF/RESONANCES/AliRsnMiniAnalysisTask.h
+++ b/PWGLF/RESONANCES/AliRsnMiniAnalysisTask.h
@@ -89,6 +89,7 @@ private:
    TClonesArray         fHistograms;      //  list of histogram definitions
    TClonesArray         fValues;          //  list of values to be computed
    TH1F                *fHEventStat;      //  histogram of event statistics
+   TH1F                *fHAEventsVsMulti; //  histogram of event statistics
 
    AliRsnCutSet        *fEventCuts;       //  cuts on events
    TObjArray            fTrackCuts;       //  list of single track cuts
@@ -100,7 +101,7 @@ private:
    Bool_t               fBigOutput;       // flag if open file for output list
    Int_t                fMixPrintRefresh; // how often info in mixing part is printed
 
-   ClassDef(AliRsnMiniAnalysisTask, 4);   // AliRsnMiniAnalysisTask
+   ClassDef(AliRsnMiniAnalysisTask, 5);   // AliRsnMiniAnalysisTask
 };
 
 inline Int_t AliRsnMiniAnalysisTask::CreateValue(AliRsnMiniValue::EType type, Bool_t useMC)