Fix in histogram binning

[u/mrichter/AliRoot.git] / PWGHF / vertexingHF / AliAnalysisTaskSEHFQA.cxx

diff --git a/PWGHF/vertexingHF/AliAnalysisTaskSEHFQA.cxx b/PWGHF/vertexingHF/AliAnalysisTaskSEHFQA.cxx
index ce6d322c08454be4390ed5f21dc53549281c7c1a..5924880de3999d88d6340223d35e89a14d9b61bf 100644 (file)
--- a/PWGHF/vertexingHF/AliAnalysisTaskSEHFQA.cxx
+++ b/PWGHF/vertexingHF/AliAnalysisTaskSEHFQA.cxx
@@ -805,7 +805,7 @@ void AliAnalysisTaskSEHFQA::UserCreateOutputObjects()
     hWhichVertSelEv->GetXaxis()->SetBinLabel(3,"SPD-3D");
     hWhichVertSelEv->GetXaxis()->SetBinLabel(4,"SPD-z");
 
-    TH2F* hTrigCent=new TH2F("hTrigCent","Centrality vs. Trigger types",12,-1.5,10.5,12,-10,110);
+    TH2F* hTrigCent=new TH2F("hTrigCent","Centrality vs. Trigger types",14,-1.5,12.5,12,-10,110);
     hTrigCent->GetXaxis()->SetBinLabel(1,"All");
     hTrigCent->GetXaxis()->SetBinLabel(2,"kAny");
     hTrigCent->GetXaxis()->SetBinLabel(3,"kMB");
@@ -817,9 +817,11 @@ void AliAnalysisTaskSEHFQA::UserCreateOutputObjects()
     hTrigCent->GetXaxis()->SetBinLabel(9,"kEMCJET+GAMMA");
     hTrigCent->GetXaxis()->SetBinLabel(10,"Muons");
     hTrigCent->GetXaxis()->SetBinLabel(11,"PHOS");
-    hTrigCent->GetXaxis()->SetBinLabel(12,"Others");
+    hTrigCent->GetXaxis()->SetBinLabel(12,"TRD");
+    hTrigCent->GetXaxis()->SetBinLabel(13,"HighMult");
+    hTrigCent->GetXaxis()->SetBinLabel(14,"Others");
 
-    TH2F* hTrigMul=new TH2F("hTrigMul","Multiplicity vs. Trigger types",12,-1.5,10.5,100,0.,10000.);
+    TH2F* hTrigMul=new TH2F("hTrigMul","Multiplicity vs. Trigger types",14,-1.5,12.5,100,0.,10000.);
     hTrigMul->GetXaxis()->SetBinLabel(1,"All");
     hTrigMul->GetXaxis()->SetBinLabel(2,"kAny");
     hTrigMul->GetXaxis()->SetBinLabel(3,"kMB");
@@ -831,9 +833,11 @@ void AliAnalysisTaskSEHFQA::UserCreateOutputObjects()
     hTrigMul->GetXaxis()->SetBinLabel(9,"kEMCJET+GAMMA");
     hTrigMul->GetXaxis()->SetBinLabel(10,"Muons");
     hTrigMul->GetXaxis()->SetBinLabel(11,"PHOS");
-    hTrigMul->GetXaxis()->SetBinLabel(12,"Others");
+    hTrigMul->GetXaxis()->SetBinLabel(12,"TRD");
+    hTrigMul->GetXaxis()->SetBinLabel(13,"HighMult");
+    hTrigMul->GetXaxis()->SetBinLabel(14,"Others");
 
-    TH2F* hTrigCentSel=new TH2F("hTrigCentSel","Trigger types",12,-1.5,10.5,12,-10,110);
+    TH2F* hTrigCentSel=new TH2F("hTrigCentSel","Trigger types",14,-1.5,12.5,12,-10,110);
     hTrigCentSel->GetXaxis()->SetBinLabel(1,"All");
     hTrigCentSel->GetXaxis()->SetBinLabel(2,"kAny");
     hTrigCentSel->GetXaxis()->SetBinLabel(3,"kMB");
@@ -845,21 +849,24 @@ void AliAnalysisTaskSEHFQA::UserCreateOutputObjects()
     hTrigCentSel->GetXaxis()->SetBinLabel(9,"kEMCJET+GAMMA");
     hTrigCentSel->GetXaxis()->SetBinLabel(10,"Muons");
     hTrigCentSel->GetXaxis()->SetBinLabel(11,"PHOS");
-    hTrigCentSel->GetXaxis()->SetBinLabel(12,"Others");
+    hTrigCentSel->GetXaxis()->SetBinLabel(12,"TRD");
+    hTrigCentSel->GetXaxis()->SetBinLabel(13,"HighMult");
+    hTrigCentSel->GetXaxis()->SetBinLabel(14,"Others");
 
     AliCounterCollection *trigCounter=new AliCounterCollection("trigCounter");
     trigCounter->AddRubric("run",500000);
     trigCounter->AddRubric("triggerType","All/Any/MB/Cent/SemiCent/EMCAL/MUON/NoPhysSelMUON/NoPhysSelEvNot7/NoPhysSelCMUP1/NoPhysSelMB/NoPhysSelCent/NoPhysSelSemiCent/CINT7");
     trigCounter->Init();
 
-    TH1F* hWhyEvRejected=new TH1F("hWhyEvRejected", "Why Event rejected",6,-0.5,5.5);
+    TH1F* hWhyEvRejected=new TH1F("hWhyEvRejected", "Why Event rejected",7,-1.5,5.5);
 
-    hWhyEvRejected->GetXaxis()->SetBinLabel(1,"pileup");
-    hWhyEvRejected->GetXaxis()->SetBinLabel(2,"centrality");
-    hWhyEvRejected->GetXaxis()->SetBinLabel(3,"Vertex (more reasons)");
-    hWhyEvRejected->GetXaxis()->SetBinLabel(4,"trigger");
-    hWhyEvRejected->GetXaxis()->SetBinLabel(5,"z vertex out");
-    hWhyEvRejected->GetXaxis()->SetBinLabel(6,"physics sel");
+    hWhyEvRejected->GetXaxis()->SetBinLabel(1,"N events");
+    hWhyEvRejected->GetXaxis()->SetBinLabel(2,"pileup");
+    hWhyEvRejected->GetXaxis()->SetBinLabel(3,"centrality");
+    hWhyEvRejected->GetXaxis()->SetBinLabel(4,"Vertex not found");
+    hWhyEvRejected->GetXaxis()->SetBinLabel(5,"trigger");
+    hWhyEvRejected->GetXaxis()->SetBinLabel(6,"z vertex out of 10 cm");
+    hWhyEvRejected->GetXaxis()->SetBinLabel(7,"physics sel");
 
 
     fOutputEvSelection->Add(evselection);
@@ -1292,11 +1299,18 @@ void AliAnalysisTaskSEHFQA::UserExec(Option_t */*option*/)
       hTrigC->Fill(9.,centrality);
       hTrigM->Fill(9.,multiplicity);
     }
-    if(evSelMask & (AliVEvent::kDG5 | AliVEvent::kZED)){
+    if(evSelMask & (AliVEvent::kTRD)){
       hTrigC->Fill(10.,centrality);
       hTrigM->Fill(10.,multiplicity);
     }
-
+    if(evSelMask & (AliVEvent::kHighMult)){
+      hTrigC->Fill(11.,centrality);
+      hTrigM->Fill(11.,multiplicity);
+    }
+    if(evSelMask & (AliVEvent::kDG5 | AliVEvent::kZED)){
+      hTrigC->Fill(12.,centrality);
+      hTrigM->Fill(12.,multiplicity);
+    }
   }
   
 
@@ -1334,35 +1348,36 @@ void AliAnalysisTaskSEHFQA::UserExec(Option_t */*option*/)
   TH1F* hWhyEvRejected=0x0;
   if(fOnOff[3]){
     hWhyEvRejected=(TH1F*)fOutputEvSelection->FindObject("hWhyEvRejected");
+     if(hWhyEvRejected) hWhyEvRejected->Fill(-1); 
   }
 
   //select event
   if(!fCuts->IsEventSelected(aod)) {
     evSelected=kFALSE;
     if(fCuts->IsEventRejectedDueToPileupSPD()) {
-      if(hWhyEvRejected) hWhyEvRejected->Fill(1); 
+      if(hWhyEvRejected) hWhyEvRejected->Fill(0); 
       evselByPileup=kFALSE;
     }// rejected for pileup
     if(fCuts->IsEventRejectedDueToCentrality()) {
-      if(hWhyEvRejected) hWhyEvRejected->Fill(2); 
+      if(hWhyEvRejected) hWhyEvRejected->Fill(1); 
       evSelbyCentrality=kFALSE; //rejected by centrality
     }
     if(fCuts->IsEventRejectedDueToNotRecoVertex() ||
-       fCuts->IsEventRejectedDueToVertexContributors() ||
-       fCuts->IsEventRejectedDueToZVertexOutsideFiducialRegion()){ 
+       fCuts->IsEventRejectedDueToVertexContributors()){ 
       evSelByVertex=kFALSE; 
-      if(hWhyEvRejected) hWhyEvRejected->Fill(3);
+      if(hWhyEvRejected) hWhyEvRejected->Fill(2);
     }
     if(fCuts->IsEventRejectedDueToTrigger()){
-      if(hWhyEvRejected) hWhyEvRejected->Fill(4);
+      if(hWhyEvRejected) hWhyEvRejected->Fill(3);
     }
     if(fCuts->IsEventRejectedDueToZVertexOutsideFiducialRegion()) {
+      evSelByVertex=kFALSE; 
       if(fOnOff[3]) ((AliCounterCollection*)fOutputEvSelection->FindObject("evselection"))->Count(Form("evnonsel:zvtx/Run:%d",runNumber)); 
-      if(hWhyEvRejected) hWhyEvRejected->Fill(5);
+      if(hWhyEvRejected) hWhyEvRejected->Fill(4);
     }
     if(fCuts->IsEventRejectedDuePhysicsSelection()) { 
       evSelByPS=kFALSE;
-      if(hWhyEvRejected) hWhyEvRejected->Fill(6); 
+      if(hWhyEvRejected) hWhyEvRejected->Fill(5); 
     }
   }
   if(evSelected && fOnOff[3]){
@@ -1379,7 +1394,9 @@ void AliAnalysisTaskSEHFQA::UserExec(Option_t */*option*/)
     if(evSelMask & (AliVEvent::kEMCEJE | AliVEvent::kEMCEGA)) hTrigS->Fill(7.,centrality);
     if(evSelMask & (((AliVEvent::kCMUS5 | AliVEvent::kMUSH7) | (AliVEvent::kMUL7 | AliVEvent::kMUU7)) |  (AliVEvent::kMUS7 | AliVEvent::kMUON))) hTrigS->Fill(8.,centrality);
     if(evSelMask & (AliVEvent::kPHI1 | AliVEvent::kPHI7)) hTrigS->Fill(9.,centrality);
-    if(evSelMask & (AliVEvent::kDG5 | AliVEvent::kZED)) hTrigS->Fill(10.,centrality);
+    if(evSelMask & (AliVEvent::kTRD)) hTrigS->Fill(10.,centrality);
+    if(evSelMask & (AliVEvent::kHighMult)) hTrigS->Fill(11.,centrality);
+    if(evSelMask & (AliVEvent::kDG5 | AliVEvent::kZED)) hTrigS->Fill(12.,centrality);
   }
   
   if(evSelected || (!evSelbyCentrality && evSelByVertex && evselByPileup && evSelByPS)){ //events selected or not selected because of centrality