Add T0 and SPD pile up and SPD background flag to counters (Cynthia)

diff --git a/PWGPP/MUON/lite/AliAnalysisTaskMuonQA.cxx b/PWGPP/MUON/lite/AliAnalysisTaskMuonQA.cxx
index 6ecf76b4a06bcf652235b21bd7caa0de4dd06f87..676f5d38db1b5976d295d90d377b4ec040b6e71e 100644 (file)
--- a/PWGPP/MUON/lite/AliAnalysisTaskMuonQA.cxx
+++ b/PWGPP/MUON/lite/AliAnalysisTaskMuonQA.cxx
@@ -43,6 +43,8 @@
 #include "AliAnalysisManager.h"
 #include "AliAnalysisTaskMuonQA.h"
 #include "AliCounterCollection.h"
+#include "AliPhysicsSelection.h"
+#include "AliTriggerAnalysis.h"
 
 // PWG includes
 #include "AliAnalysisMuonUtility.h"
@@ -361,7 +363,11 @@ void AliAnalysisTaskMuonQA::UserCreateOutputObjects()
   fTrackCounters->AddRubric("pt", "low/high/any");
   fTrackCounters->AddRubric("acc", "in/out");
   fTrackCounters->AddRubric("tagTrack", "beamGas/good");
+  fTrackCounters->AddRubric("t0pileup", "yes/no");
+  fTrackCounters->AddRubric("bgID", "yes/no");
+  fTrackCounters->AddRubric("spdpileup", "yes/no");
   fTrackCounters->Init();
+
   
   // initialize event counters
   fEventCounters = new AliCounterCollection("eventCounters");
@@ -371,6 +377,9 @@ void AliAnalysisTaskMuonQA::UserCreateOutputObjects()
   fEventCounters->AddRubric("selected", "yes/no");
   fEventCounters->AddRubric("triggerRO", "good/bad");
   fEventCounters->AddRubric("v0mult", "low/int/high/any");
+  fEventCounters->AddRubric("t0pileup", "yes/no");
+  fEventCounters->AddRubric("spdpileup", "yes/no");
+  fEventCounters->AddRubric("bgID", "yes/no");
   fEventCounters->Init();
   
   // Post data at least once per task to ensure data synchronisation (required for merging)
@@ -384,19 +393,32 @@ void AliAnalysisTaskMuonQA::UserCreateOutputObjects()
 void AliAnalysisTaskMuonQA::UserExec(Option_t *)
 {
   /// Called for each event
-  
+
   AliESDEvent* fESD = dynamic_cast<AliESDEvent*>(InputEvent());
   if (!fESD) {
     Printf("ERROR: fESD not available");
     return;
   }
-  
+
+  //Flag for T0Pileup, SPDPileup and bgID (SPD cluster vs tracket)
+  Bool_t t0PileUp = kFALSE, spdPileUp = kFALSE, bgID = kFALSE;
+  spdPileUp = fESD->IsPileupFromSPDInMultBins();
+  AliPhysicsSelection *physicsSelection = (AliPhysicsSelection*)((AliInputEventHandler*)(AliAnalysisManager::GetAnalysisManager()->GetInputEventHandler()))->GetEventSelection();
+  AliTriggerAnalysis * triggerAnalysis = (AliTriggerAnalysis*)physicsSelection->GetTriggerAnalysis();
+  if ( physicsSelection && triggerAnalysis ) {
+    t0PileUp = triggerAnalysis->EvaluateTrigger(fESD, (AliTriggerAnalysis::Trigger) (AliTriggerAnalysis::kOfflineFlag | AliTriggerAnalysis::kT0Pileup));
+    bgID = triggerAnalysis->EvaluateTrigger(fESD,  (AliTriggerAnalysis::Trigger) (AliTriggerAnalysis::kSPDClsVsTrkBG | AliTriggerAnalysis::kOfflineFlag)); 
+  }
+
   UInt_t geomAccMask = ( AliMuonTrackCuts::kMuEta | AliMuonTrackCuts::kMuThetaAbs );
   
   // check physics selection
   UInt_t triggerWord = (fInputHandler) ? fInputHandler->IsEventSelected() : 0;
   Bool_t isPhysicsSelected = (triggerWord != 0);
   TString selected = isPhysicsSelected ? "selected:yes" : "selected:no";
+  selected += t0PileUp ? "/t0pileup:yes" : "/t0pileup:no";
+  selected += bgID ? "/bgID:yes" : "/bgID:no";
+  selected += spdPileUp ? "/spdpileup:yes" : "/spdpileup:no";
 
   // fill muon trigger cases
   for ( Int_t idx=0; idx<fMuonTrigIndex.GetSize(); idx++ ) {

diff --git a/PWGPP/MUON/lite/PlotMuonQA.C b/PWGPP/MUON/lite/PlotMuonQA.C
index 56cdb9b9f9a9e8e2a0fa8d38f537d9bb86589c63..f1bce36a418d43eccd0828df22b326d5236ae709 100644 (file)
--- a/PWGPP/MUON/lite/PlotMuonQA.C
+++ b/PWGPP/MUON/lite/PlotMuonQA.C
@@ -238,6 +238,8 @@ void PlotMuonQA(const char* baseDir, const char* runList = 0x0, const char * tri
 
   TH1* hBNoPS[centBinMax][kNMaxTriggers];
   TH1* hBWithPS[centBinMax][kNMaxTriggers];
+  TH1* hBWithPST0Flag[centBinMax][kNMaxTriggers]; 
+  TH1* hBWithPST0SPDFlag[centBinMax][kNMaxTriggers]; 
   TH1* hB[centBinMax][kNMaxTriggers]={};
   TH1* hTriggerB[centBinMax][kNMaxTriggers], *hTrackerB[centBinMax][kNMaxTriggers], *hMatchedB[centBinMax][kNMaxTriggers], *hAllTracksB[centBinMax][kNMaxTriggers], *hMatchedLowPtB[centBinMax][kNMaxTriggers], *hMatchedHighPtB[centBinMax][kNMaxTriggers];
   TH1* hMatchedLowPtBNoPS[centBinMax][kNMaxTriggers], *hMatchedHighPtBNoPS[centBinMax][kNMaxTriggers];
@@ -276,6 +278,16 @@ void PlotMuonQA(const char* baseDir, const char* runList = 0x0, const char * tri
       histoName = histoNameBase;
       histoName += "BWithPS";
       hBWithPS[centBin][i] = (TH1*) ProcessHisto(eventCounters, "run", selection, histoName, "", "Trigger content w/ Phys. Sel.", colorInd->At(i));
+    // Histo trigger with Phys. Sel. and T0 pile up not flagged
+      selection = selectionCent; selection += Form("trigger:%s/%s/selected:yes/t0pileup:no", triggerName.Data(), selectRuns.Data());
+      histoName = histoNameBase;
+      histoName += "BWithPST0Flag";
+      hBWithPST0Flag[centBin][i] = (TH1*) ProcessHisto(eventCounters, "run", selection, histoName, "", "Trigger content w/ Phys. Sel. and no pile up from T0 flag", colorInd->At(i));
+    // Histo trigger with Phys. Sel. and T0 + SPD pile up not flagged
+      selection = selectionCent; selection += Form("trigger:%s/%s/selected:yes/t0pileup:no/spdpileup:no", triggerName.Data(), selectRuns.Data());
+      histoName = histoNameBase;
+      histoName += "BWithPST0SPDFlag";
+      hBWithPST0SPDFlag[centBin][i] = (TH1*) ProcessHisto(eventCounters, "run", selection, histoName, "", "Trigger content w/ Phys. Sel. and no pile up from T0 and SPD flag", colorInd->At(i));
       // Histo trigger : Phys. Sel.  is selected or not depending on the macro arguments
       selection = selectionCent; selection += Form("trigger:%s/%s/%s", triggerName.Data(), selectRuns.Data(), select.Data());
       histoName = histoNameBase;
@@ -360,6 +372,8 @@ void PlotMuonQA(const char* baseDir, const char* runList = 0x0, const char * tri
        
   Int_t NumOfBNoPS[centBinMax][kNMaxTriggers];
   Int_t NumOfBWithPS[centBinMax][kNMaxTriggers];
+  Int_t NumOfBWithPST0Flag[centBinMax][kNMaxTriggers];
+  Int_t NumOfBWithPST0SPDFlag[centBinMax][kNMaxTriggers];
   Int_t NumOfACNoPS[centBinMax][kNMaxTriggers];
   Int_t NumOfENoPS[centBinMax][kNMaxTriggers];
   Int_t NumOfACWithPS[centBinMax][kNMaxTriggers];
@@ -369,6 +383,8 @@ void PlotMuonQA(const char* baseDir, const char* runList = 0x0, const char * tri
     for(Int_t i = 0; i < triggersB->GetEntriesFast(); i++){
       NumOfBNoPS[centBin][i] = hBNoPS[centBin][i]->Integral();
       NumOfBWithPS[centBin][i] = hBWithPS[centBin][i]->Integral();
+      NumOfBWithPST0Flag[centBin][i] = hBWithPST0Flag[centBin][i]->Integral();
+      NumOfBWithPST0SPDFlag[centBin][i] = hBWithPST0Flag[centBin][i]->Integral();
       NumOfACNoPS[centBin][i] = hACNoPS[centBin][i]->Integral();
       NumOfENoPS[centBin][i] = hENoPS[centBin][i]->Integral();
       NumOfACWithPS[centBin][i] = hACWithPS[centBin][i]->Integral();
@@ -438,6 +454,33 @@ void PlotMuonQA(const char* baseDir, const char* runList = 0x0, const char * tri
   cPhysSelCut->Write();
   rootFileOut->cd();
        
+  cout<<"//==================================================================================="<<endl;
+  cout<<"// new canvas with effect from T0 pile-up flag for each trigger and background trigger "<<endl;
+  cout<<"//==================================================================================="<<endl;
+       
+  TCanvas *cPhysSelCutT0Flag = 0x0;
+  cPhysSelCutT0Flag = ProcessCanvasPhysSelCut(triggersB, triggersAC, triggersE, hBWithPS[centBin], hACNoPS[centBin], hENoPS[centBin], hBWithPST0Flag[centBin], "PhysSelCutOnCollTriggerT0Flag", colorInd,isHeavyIon);
+  cPhysSelCutT0Flag->Draw();
+  cPhysSelCutT0Flag->Print(OutFileNamePDF.Data());
+  dirGlobal->cd();
+  cPhysSelCutT0Flag->Write();
+  rootFileOut->cd();
+       
+
+ cout<<"//==================================================================================="<<endl;
+  cout<<"// new canvas with effect from T0 + SPD pile-up flag for each trigger and background trigger "<<endl;
+  cout<<"//==================================================================================="<<endl;
+       
+  TCanvas *cPhysSelCutT0SPDFlag = 0x0;
+  cPhysSelCutT0SPDFlag = ProcessCanvasPhysSelCut(triggersB, triggersAC, triggersE, hBWithPS[centBin], hACNoPS[centBin], hENoPS[centBin], hBWithPST0SPDFlag[centBin], "PhysSelCutOnCollTriggerT0Flag", colorInd,isHeavyIon);
+  cPhysSelCutT0SPDFlag->Draw();
+  cPhysSelCutT0SPDFlag->Print(OutFileNamePDF.Data());
+  dirGlobal->cd();
+  cPhysSelCutT0SPDFlag->Write();
+  rootFileOut->cd();
+
+
+
   cout<<"//==================================================================================="<<endl;
   cout<<"// new canvas with effect from physics selection for each trigger and centrality bin (only in PbPb) "<<endl;
   cout<<"//==================================================================================="<<endl;
@@ -707,9 +750,9 @@ void PlotMuonQA(const char* baseDir, const char* runList = 0x0, const char * tri
     } 
   }
   
-  
-  //return;
-        
+  //temporary
+  //    return;
+
   
   //--------------------------------------------- //
   //        monitor quantities run per run        //
@@ -767,6 +810,7 @@ void PlotMuonQA(const char* baseDir, const char* runList = 0x0, const char * tri
     isAlienFile = 1;
     alienBaseDir.ReplaceAll("alien://","");
   }
+
   cout<<"============================================================"<<endl;
   cout<< "Monitoring quantities run per run: loop over "<<runs->GetEntriesFast()<<" runs."<<endl;
   cout<<"  directory = "<<alienBaseDir.Data()<<endl;