For multiplicity plot

diff --git a/PWGHF/hfe/AliAnalysisTaskHFEFlow.cxx b/PWGHF/hfe/AliAnalysisTaskHFEFlow.cxx
index 675acad963644344975f8adb70afda7b3cf150b0..baccaabb245b56ae58bb6d4169ed364db34ab740 100644 (file)
--- a/PWGHF/hfe/AliAnalysisTaskHFEFlow.cxx
+++ b/PWGHF/hfe/AliAnalysisTaskHFEFlow.cxx
@@ -133,7 +133,7 @@ AliAnalysisTaskHFEFlow::AliAnalysisTaskHFEFlow() :
   fArraytrack(NULL),\r
   fCounterPoolBackground(0),\r
   fHFEVZEROEventPlane(0x0),\r
-  fHistEV(0),\r
+  fMultEv(0),\r
   fEventPlane(0x0),\r
   fEventPlaneaftersubtraction(0x0),\r
   fFractionContamination(0x0),\r
@@ -239,7 +239,7 @@ AliAnalysisTaskHFEFlow:: AliAnalysisTaskHFEFlow(const char *name) :
   fArraytrack(NULL),\r
   fCounterPoolBackground(0),\r
   fHFEVZEROEventPlane(0x0),\r
-  fHistEV(0),\r
+  fMultEv(0),\r
   fEventPlane(0x0),\r
   fEventPlaneaftersubtraction(0x0),\r
   fFractionContamination(0x0),\r
@@ -367,7 +367,7 @@ AliAnalysisTaskHFEFlow::AliAnalysisTaskHFEFlow(const AliAnalysisTaskHFEFlow &ref
   fArraytrack(NULL),\r
   fCounterPoolBackground(ref.fCounterPoolBackground),\r
   fHFEVZEROEventPlane(NULL),\r
-  fHistEV(NULL),\r
+  fMultEv(NULL),\r
   fEventPlane(NULL),\r
   fEventPlaneaftersubtraction(NULL),\r
   fFractionContamination(NULL),\r
@@ -772,7 +772,7 @@ void AliAnalysisTaskHFEFlow::UserCreateOutputObjects()
   // Minimum histos\r
 \r
   // Histos\r
-  fHistEV = new TH2D("fHistEV", "events", 3, 0, 3, 3, 0,3);\r
+  fMultEv = new TH1D("fMultEv", "", 5000, 0, 10000.);\r
   \r
   AliDebug(2,"AliAnalysisTaskHFEFlow: histev");\r
 \r
@@ -1125,7 +1125,7 @@ void AliAnalysisTaskHFEFlow::UserCreateOutputObjects()
   // Add to the list\r
   //******************************\r
 \r
-  fListHist->Add(fHistEV);\r
+  fListHist->Add(fMultEv);\r
   fListHist->Add(fEventPlane);\r
   fListHist->Add(fFractionContamination);\r
   fListHist->Add(fCosRes);\r
@@ -1344,7 +1344,7 @@ void AliAnalysisTaskHFEFlow::UserExec(Option_t */*option*/)
   fPIDTOFOnly->SetPIDResponse(pidResponse);\r
   fPIDBackground->SetPIDResponse(pidResponse);\r
 \r
-  fHistEV->Fill(binctt,0.0);\r
+  //fMultEv->Fill(binctt,0.0);\r
  \r
 \r
   //////////////////\r
@@ -1356,7 +1356,7 @@ void AliAnalysisTaskHFEFlow::UserExec(Option_t */*option*/)
     return;\r
   }\r
 \r
-  fHistEV->Fill(binctt,1.0);\r
+  //fMultEv->Fill(binctt,1.0);\r
 \r
   ////////////////////////////////////  \r
   // First method event plane\r
@@ -1537,7 +1537,7 @@ void AliAnalysisTaskHFEFlow::UserExec(Option_t */*option*/)
 \r
   if(eventplanedefined) {\r
     \r
-    fHistEV->Fill(binctt,2.0);\r
+    //fMultEv->Fill(binctt,2.0);\r
     \r
     // Fill\r
     valuensparsea[0] = eventPlaneV0A;\r
@@ -1684,6 +1684,7 @@ void AliAnalysisTaskHFEFlow::UserExec(Option_t */*option*/)
   //////////////////////////\r
   // Loop over track\r
   //////////////////////////\r
+  Double_t nboftrackTPCcut = 0.0;\r
   \r
   for(Int_t k = 0; k < nbtracks; k++){\r
       \r
@@ -1710,6 +1711,7 @@ void AliAnalysisTaskHFEFlow::UserExec(Option_t */*option*/)
       // RecKine: ITSTPC cuts  \r
       if(!fHFECuts->CheckParticleCuts(AliHFEcuts::kStepRecKineITSTPC + AliHFEcuts::kNcutStepsMCTrack, (TObject *)track)) continue;\r
       if(fMonitorTrackCuts) fTrackingCuts->Fill(&valuetrackingcuts[0]);\r
+      nboftrackTPCcut++;\r
       \r
       // Reject kink mother\r
       if(fAODAnalysis) {\r
@@ -2037,6 +2039,8 @@ void AliAnalysisTaskHFEFlow::UserExec(Option_t */*option*/)
       fArraytrack = NULL;\r
     }\r
   }\r
+\r
+  if((cntr>20) && (cntr<40)) fMultEv->Fill(nboftrackTPCcut);\r
   \r
   PostData(1, fListHist);\r
 \r

diff --git a/PWGHF/hfe/AliAnalysisTaskHFEFlow.h b/PWGHF/hfe/AliAnalysisTaskHFEFlow.h
index 0aaa6fc3aba3182d5c6729a6b92223ebcc79efea..6f86ac910c9e5bdfffa0370d469856f19296abff 100644 (file)
--- a/PWGHF/hfe/AliAnalysisTaskHFEFlow.h
+++ b/PWGHF/hfe/AliAnalysisTaskHFEFlow.h
@@ -221,7 +221,7 @@ private:
   AliHFEVZEROEventPlane *fHFEVZEROEventPlane; // VZERO event plane calibrated\r
   \r
   // Histos\r
-  TH2D *fHistEV;               //! Number of events\r
+  TH1D *fMultEv;               //! Multiplicity in the event\r
   \r
   // A Event plane as function of phiepa, phiepb, phiepc, phiepd centrality \r
   // a V0A, b V0C, c TPC,\r