For multiplicity plot
authorrbailhac <rbailhac@f7af4fe6-9843-0410-8265-dc069ae4e863>
Wed, 18 Jul 2012 16:47:16 +0000 (16:47 +0000)
committerrbailhac <rbailhac@f7af4fe6-9843-0410-8265-dc069ae4e863>
Wed, 18 Jul 2012 16:47:16 +0000 (16:47 +0000)
PWGHF/hfe/AliAnalysisTaskHFEFlow.cxx
PWGHF/hfe/AliAnalysisTaskHFEFlow.h

index 675acad..baccaab 100644 (file)
@@ -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
index 0aaa6fc..6f86ac9 100644 (file)
@@ -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