add pile-up event counter based on SPD

diff --git a/PWG/CaloTrackCorrBase/AliAnaCaloTrackCorrMaker.cxx b/PWG/CaloTrackCorrBase/AliAnaCaloTrackCorrMaker.cxx
index c2a6217a5f9169c9091a5a880ba1cedc679a1965..1d61a8db8ee1394f4295322ff078e51a7b669dda 100755 (executable)
--- a/PWG/CaloTrackCorrBase/AliAnaCaloTrackCorrMaker.cxx
+++ b/PWG/CaloTrackCorrBase/AliAnaCaloTrackCorrMaker.cxx
@@ -47,7 +47,8 @@ fOutputContainer(new TList ), fAnalysisContainer(new TList ),
 fMakeHisto(kFALSE),           fMakeAOD(kFALSE), 
 fAnaDebug(0),                 fCuts(new TList), 
 fScaleFactor(-1),
-fhNEvents(0),                 fhZVertex(0),                 
+fhNEvents(0),                 fhNPileUpEvents(0),
+fhZVertex(0),                 
 fhTrackMult(0),               fhCentrality(0),
 fhEventPlaneAngle(0),
 fhNMergedFiles(0),            fhScaleFactor(0)
@@ -68,8 +69,11 @@ fOutputContainer(new TList()), fAnalysisContainer(new TList()),
 fMakeHisto(maker.fMakeHisto),  fMakeAOD(maker.fMakeAOD),
 fAnaDebug(maker.fAnaDebug),    fCuts(new TList()),
 fScaleFactor(maker.fScaleFactor),
-fhNEvents(maker.fhNEvents),    fhZVertex(maker.fhZVertex),    
-fhTrackMult(maker.fhTrackMult),fhCentrality(maker.fhCentrality),
+fhNEvents(maker.fhNEvents), 
+fhNPileUpEvents(maker.fhNPileUpEvents),
+fhZVertex(maker.fhZVertex),    
+fhTrackMult(maker.fhTrackMult),
+fhCentrality(maker.fhCentrality),
 fhEventPlaneAngle(maker.fhEventPlaneAngle),
 fhNMergedFiles(maker.fhNMergedFiles),          
 fhScaleFactor(maker.fhScaleFactor)
@@ -175,6 +179,12 @@ TList *AliAnaCaloTrackCorrMaker::GetOutputContainer()
   fhNEvents->SetYTitle("# events");
   fOutputContainer->Add(fhNEvents);
   
+  fhNPileUpEvents      = new TH1I("hNPileUpEvents",   "Number of events considered as pile-up", 2 , 0 , 2  ) ;
+  fhNPileUpEvents->SetYTitle("# events");
+  fhNPileUpEvents->GetXaxis()->SetBinLabel(1 ,"SPD");
+  fhNPileUpEvents->GetXaxis()->SetBinLabel(2 ,"Multi SPD");
+  fOutputContainer->Add(fhNPileUpEvents);
+  
   fhZVertex      = new TH1F("hZVertex", " Z vertex distribution"   , 200 , -50 , 50  ) ;
   fhZVertex->SetXTitle("v_{z} (cm)");
   fOutputContainer->Add(fhZVertex);
@@ -414,16 +424,18 @@ void AliAnaCaloTrackCorrMaker::ProcessEvent(const Int_t iEntry,
   // Event control histograms
  
   fhNEvents        ->Fill(0); // Number of events analyzed
+  if(fReader->IsPileUpFromSPD()) 
+    fhNPileUpEvents->Fill(0.5);
+  if(fReader->GetInputEvent()->IsPileupFromSPDInMultBins()) 
+    fhNPileUpEvents->Fill(1.5);
   fhTrackMult      ->Fill(fReader->GetTrackMultiplicity()); 
   fhCentrality     ->Fill(fReader->GetEventCentrality  ());
   fhEventPlaneAngle->Fill(fReader->GetEventPlaneAngle  ());
-  
-  
+
   Double_t v[3];
   fReader->GetInputEvent()->GetPrimaryVertex()->GetXYZ(v) ;
   fhZVertex->Fill(v[2]);
   
-  
   //printf(">>>>>>>>>> AFTER >>>>>>>>>>>\n");
   //gObjectTable->Print();
        

diff --git a/PWG/CaloTrackCorrBase/AliAnaCaloTrackCorrMaker.h b/PWG/CaloTrackCorrBase/AliAnaCaloTrackCorrMaker.h
index 8670dd5c2d48c7e0cccd831d98b04977d41b8a74..e5746fc9128541c4f13da985ec20c30f260d1150 100755 (executable)
--- a/PWG/CaloTrackCorrBase/AliAnaCaloTrackCorrMaker.h
+++ b/PWG/CaloTrackCorrBase/AliAnaCaloTrackCorrMaker.h
@@ -92,6 +92,7 @@ class AliAnaCaloTrackCorrMaker : public TObject {
 
   // Control histograms
   TH1I *   fhNEvents;           //! Number of events counter histogram
+  TH1I *   fhNPileUpEvents;     //! N events pasing pile up cut
   TH1F *   fhZVertex;           //! Vertex of accepted event
   TH1I *   fhTrackMult;         //! Number of tracks per event histogram
   TH1F *   fhCentrality;        //! Histogram with centrality bins

diff --git a/PWG/CaloTrackCorrBase/AliCaloTrackReader.cxx b/PWG/CaloTrackCorrBase/AliCaloTrackReader.cxx
index 02b0b9257b1b9032485f7dc5e58b8c9fe41dd482..f2207ce2d96e2e21ff971bc4686310701e1f724f 100755 (executable)
--- a/PWG/CaloTrackCorrBase/AliCaloTrackReader.cxx
+++ b/PWG/CaloTrackCorrBase/AliCaloTrackReader.cxx
@@ -446,6 +446,8 @@ Bool_t AliCaloTrackReader::IsPileUpFromSPD() const
   // Default values: (3, 0.8, 3., 2., 5.)
   return fInputEvent->IsPileupFromSPD((Int_t) fPileUpParam[0] , fPileUpParam[1] , 
                                               fPileUpParam[2] , fPileUpParam[3] , fPileUpParam[4] ); 
+  //printf("Param : %d, %2.2f, %2.2f, %2.2f, %2.2f\n",(Int_t) fPileUpParam[0], fPileUpParam[1], fPileUpParam[2], fPileUpParam[3], fPileUpParam[4]);
+
 }
 
 //________________________________________________________
@@ -610,6 +612,8 @@ Bool_t AliCaloTrackReader::FillInputEvent(const Int_t iEntry,
         TMath::Abs(fVertex[0][2] ) < 1.e-6    ) return kFALSE;
   }
   
+  //printf("Reader : IsPileUp %d, Multi %d\n",IsPileUpFromSPD(),fInputEvent->IsPileupFromSPDInMultBins());
+  
   if(fDoEventSelection)
   {
     if(!fCaloFilterPatch)