Adding a track candle for p-A correlations and comparison with CMS.

diff --git a/ANALYSIS/AliCentralitySelectionTask.cxx b/ANALYSIS/AliCentralitySelectionTask.cxx
index 93153c0bb4a4efda2351b37fcbb37565e4e0f03a..6b5f13b4a5fd3eabf20f8526d59f86067d34a0e2 100644 (file)
--- a/ANALYSIS/AliCentralitySelectionTask.cxx
+++ b/ANALYSIS/AliCentralitySelectionTask.cxx
@@ -106,6 +106,9 @@ AliAnalysisTaskSE(),
   fV0MZDCEcalOutlierPar0(0),   
   fV0MZDCEcalOutlierPar1(0),   
   fTrackCuts(0),
+  fEsdTrackCuts(0),
+  fEsdTrackCutsExtra1(0),
+  fEsdTrackCutsExtra2(0),
   fZVCut(10),
   fOutliersCut(5),
   fQuality(999),
@@ -295,6 +298,9 @@ AliCentralitySelectionTask::AliCentralitySelectionTask(const char *name):
   fV0MZDCEcalOutlierPar0(0),   
   fV0MZDCEcalOutlierPar1(0),   
   fTrackCuts(0),
+  fEsdTrackCuts(0),
+  fEsdTrackCutsExtra1(0),
+  fEsdTrackCutsExtra2(0),
   fZVCut(10),
   fOutliersCut(5),
   fQuality(999),
@@ -494,6 +500,9 @@ AliCentralitySelectionTask::AliCentralitySelectionTask(const AliCentralitySelect
   fV0MZDCEcalOutlierPar0(ana.fV0MZDCEcalOutlierPar0),   
   fV0MZDCEcalOutlierPar1(ana.fV0MZDCEcalOutlierPar1),   
   fTrackCuts(ana.fTrackCuts),
+  fEsdTrackCuts(ana.fEsdTrackCuts),
+  fEsdTrackCutsExtra1(ana.fEsdTrackCutsExtra1),
+  fEsdTrackCutsExtra2(ana.fEsdTrackCutsExtra2),
   fZVCut(ana.fZVCut),
   fOutliersCut(ana.fOutliersCut),
   fQuality(ana.fQuality),
@@ -655,6 +664,9 @@ AliCentralitySelectionTask::~AliCentralitySelectionTask()
   // Destructor  
   if (fOutputList && !AliAnalysisManager::GetAnalysisManager()->IsProofMode()) delete fOutputList;
   if (fTrackCuts) delete fTrackCuts;
+  if (fEsdTrackCuts) delete fEsdTrackCuts;
+  if (fEsdTrackCutsExtra1) delete fEsdTrackCutsExtra1;
+  if (fEsdTrackCutsExtra2) delete fEsdTrackCutsExtra2;
 }  
 
 //________________________________________________________________________
@@ -884,6 +896,14 @@ void AliCentralitySelectionTask::UserCreateOutputObjects()
   }
   
   fTrackCuts = AliESDtrackCuts::GetStandardTPCOnlyTrackCuts();
+  fEsdTrackCuts = AliESDtrackCuts::GetStandardITSTPCTrackCuts2011();
+  fEsdTrackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD,AliESDtrackCuts::kOff);
+  // Add SPD requirement
+  fEsdTrackCutsExtra1 = new AliESDtrackCuts("SPD", "Require 1 cluster in SPD");
+  fEsdTrackCutsExtra1->SetClusterRequirementITS(AliESDtrackCuts::kSPD,AliESDtrackCuts::kAny);
+  // Add SDD requirement
+  fEsdTrackCutsExtra2 = new AliESDtrackCuts("SDD", "Require 1 cluster in first layer SDD");
+  fEsdTrackCutsExtra2->SetClusterRequirementITS(AliESDtrackCuts::kSDD,AliESDtrackCuts::kFirst);  
 }
 
 //________________________________________________________________________
@@ -1041,11 +1061,20 @@ void AliCentralitySelectionTask::UserExec(Option_t */*option*/)
     // ***** CB info (tracklets, clusters, chips)
     //nTracks    = event->GetNumberOfTracks();     
     nTracks    = fTrackCuts ? (Short_t)fTrackCuts->GetReferenceMultiplicity(esd,kTRUE):-1;
+
     Short_t nTrTPCcandle = 0;
     for (Int_t iTracks = 0; iTracks < esd->GetNumberOfTracks(); iTracks++) {
+
       AliESDtrack* track = esd->GetTrack(iTracks);
       if (!track) continue;
-      if (track->Pt() > 0.5 && TMath::Abs(track->Eta()) < 0.8) nTrTPCcandle++;
+      
+      if (! fEsdTrackCuts->IsSelected(track) )continue;
+      
+      if (fEsdTrackCutsExtra1 && fEsdTrackCutsExtra2 &&
+                 !fEsdTrackCutsExtra1->IsSelected(track) &&
+                 !fEsdTrackCutsExtra2->IsSelected(track)) continue;
+      
+      if (track->Pt() > 0.4 && TMath::Abs(track->Eta()) < 0.9) nTrTPCcandle++;
     } 
     multCND = nTrTPCcandle;
 

diff --git a/ANALYSIS/AliCentralitySelectionTask.h b/ANALYSIS/AliCentralitySelectionTask.h
index 445fc2c234ef7545a42fe7b73be4aaba83c01082..80d3cc2898b20176d6736f89a02aab281a78a9f5 100644 (file)
--- a/ANALYSIS/AliCentralitySelectionTask.h
+++ b/ANALYSIS/AliCentralitySelectionTask.h
@@ -76,6 +76,9 @@ class AliCentralitySelectionTask : public AliAnalysisTaskSE {
   Float_t  fV0MZDCEcalOutlierPar1;   // outliers parameter
 
   AliESDtrackCuts* fTrackCuts;  //! optional track cuts
+  AliESDtrackCuts* fEsdTrackCuts;  //! optional track cuts
+  AliESDtrackCuts* fEsdTrackCutsExtra1;  //! optional track cuts
+  AliESDtrackCuts* fEsdTrackCutsExtra2;  //! optional track cuts
 
   Float_t  fZVCut;              //! z-vertex cut (in cm)
   Float_t  fOutliersCut;        //! outliers cut (in n-sigma)
@@ -241,7 +244,7 @@ class AliCentralitySelectionTask : public AliAnalysisTaskSE {
   TH1F *fHOutVertex ;           //control histogram for vertex SPD
   TH1F *fHOutVertexT0 ;         //control histogram for vertex T0
 
-  ClassDef(AliCentralitySelectionTask, 23); 
+  ClassDef(AliCentralitySelectionTask, 24); 
 };
 
 #endif