allow parameterization of pile-up rejection via spd method

diff --git a/PWG/CaloTrackCorrBase/AliCaloTrackReader.cxx b/PWG/CaloTrackCorrBase/AliCaloTrackReader.cxx
index 54dea3abed372eb1e7fc6cac4a9b2413001f87a4..f8befc36293db89e2a2786ef550db6aa85e09737 100755 (executable)
--- a/PWG/CaloTrackCorrBase/AliCaloTrackReader.cxx
+++ b/PWG/CaloTrackCorrBase/AliCaloTrackReader.cxx
@@ -433,6 +433,9 @@ void AliCaloTrackReader::InitParameters()
 
   fImportGeometryFromFile = kFALSE;
   
+  fPileUpParam[0] = 3   ; fPileUpParam[1] = 0.8 ;
+  fPileUpParam[2] = 3.0 ; fPileUpParam[3] = 2.0 ; fPileUpParam[4] = 5.0;
+  
 }
 
 //________________________________________________________
@@ -601,9 +604,12 @@ Bool_t AliCaloTrackReader::FillInputEvent(const Int_t iEntry,
   {
     if(!fCaloFilterPatch)
     {
-      //Do not analyze events with pileup
-      Bool_t bPileup = fInputEvent->IsPileupFromSPD(3, 0.8, 3., 2., 5.); //Default values, if not it does not compile
-      //Bool_t bPileup = event->IsPileupFromSPD(); 
+      // Do not analyze events with pileup
+      // Default values: (3, 0.8, 3., 2., 5.)
+      Bool_t bPileup = fInputEvent->IsPileupFromSPD((Int_t) fPileUpParam[0] , fPileUpParam[1] , 
+                                                    fPileUpParam[2] ,fPileUpParam[3] , fPileUpParam[4] ); 
+      //IsPileupFromSPDInMultBins() // method to try
+      //printf("pile-up %d, %d, %2.2f, %2.2f, %2.2f, %2.2f\n",bPileup, (Int_t) fPileUpParam[0], fPileUpParam[1], fPileUpParam[2], fPileUpParam[3], fPileUpParam[4]);
       if(bPileup) return kFALSE;
       
       if(fDoV0ANDEventSelection)

diff --git a/PWG/CaloTrackCorrBase/AliCaloTrackReader.h b/PWG/CaloTrackCorrBase/AliCaloTrackReader.h
index 8b02694ceb28ba32fcfa63cf29a0a97e3127bd86..802de1befb63213c2346115c4a12b9c9b03c2e0a 100755 (executable)
--- a/PWG/CaloTrackCorrBase/AliCaloTrackReader.h
+++ b/PWG/CaloTrackCorrBase/AliCaloTrackReader.h
@@ -233,6 +233,8 @@ public:
   void             SwitchOffPrimaryVertexSelection()       { fUseEventsWithPrimaryVertex = kFALSE ; }
   Bool_t           IsPrimaryVertexSelectionDone()    const { return fUseEventsWithPrimaryVertex   ; } 
   
+  void             SetPileUpParam(Int_t i, Double_t param) { fPileUpParam[i] = param   ; }
+  
   // Track selection
   ULong_t          GetTrackStatus()                  const { return fTrackStatus       ; }
   void             SetTrackStatus(ULong_t bit)             { fTrackStatus = bit        ; }             
@@ -487,7 +489,11 @@ public:
   Bool_t           fDoV0ANDEventSelection;       // Select events depending on V0, fDoEventSelection should be on
   Bool_t           fUseEventsWithPrimaryVertex ; // Select events with primary vertex
   AliTriggerAnalysis* fTriggerAnalysis;          // Access to trigger selection algorithm for V0AND calculation
-  
+  Double_t         fPileUpParam[5];              // Parameters to pass to method IsPileupFromSPD: Int_t minContributors, 
+                                                 // Double_t minZdist, 
+                                                 // Double_t nSigmaZdist, 
+                                                 // Double_t nSigmaDiamXY, 
+                                                 // Double_t nSigmaDiamZ)
   //Centrality/Event plane
   TString          fCentralityClass;        // Name of selected centrality class     
   Int_t            fCentralityOpt;          // Option for the returned value of the centrality, possible options 5, 10, 100
@@ -500,7 +506,7 @@ public:
   AliCaloTrackReader(              const AliCaloTrackReader & r) ; // cpy ctor
   AliCaloTrackReader & operator = (const AliCaloTrackReader & r) ; // cpy assignment
   
-  ClassDef(AliCaloTrackReader,43)
+  ClassDef(AliCaloTrackReader,44)
   
 } ;