Adding histos to look at MB trigger efficiency

diff --git a/PWGLF/totEt/AliAnalysisHadEtMonteCarlo.cxx b/PWGLF/totEt/AliAnalysisHadEtMonteCarlo.cxx
index a1981287120e523d471b9e7a34d2033ff4321d1e..7ec0bb7910d345c02f9d56f5d7be939218aaffc5 100644 (file)
--- a/PWGLF/totEt/AliAnalysisHadEtMonteCarlo.cxx
+++ b/PWGLF/totEt/AliAnalysisHadEtMonteCarlo.cxx
@@ -58,6 +58,7 @@ AliAnalysisHadEtMonteCarlo::AliAnalysisHadEtMonteCarlo():AliAnalysisHadEt()
                                                        ,fBaryonEnhancement(0)
                                                        ,fUseRecoPt(0)
                                                        ,kIsOfflineV0AND(0)
+                                                       ,kIsOfflineMB(0)
                                                        ,kDoTriggerChecks(0)
                                                        ,kDoTriggerChecksOnly(0)
                                                        ,fPtSmearer(0)
@@ -99,7 +100,7 @@ Int_t AliAnalysisHadEtMonteCarlo::AnalyseEvent(AliVEvent* ev,AliVEvent* ev2)
   AnalyseEvent(ev);
   if(kDoTriggerChecksOnly) return 1;//If we are only doing trigger checks, don't bother with all of the reconstructed stuff
   //for PID
-  if(kDoTriggerChecks && !kIsOfflineV0AND){return 1;}//In this case we are just after trigger efficiencies and don't care about the ET reconstructed.
+  if(kDoTriggerChecks && (!kIsOfflineV0AND ||!kIsOfflineMB ) ){return 1;}//In this case we are just after trigger efficiencies and don't care about the ET reconstructed.
   AliESDpid *pID = new AliESDpid();//This is identified as a memory leak in valgrind but I delete this object so I think it may be a problem with AliESDpid.
 
   //=============================================
@@ -1671,6 +1672,12 @@ Int_t AliAnalysisHadEtMonteCarlo::AnalyseEvent(AliVEvent* ev)
        FillHisto1D("NEventsNDV0AND",0.5,1);
        FillHisto1D("SimPiKPEtNDV0AND",fSimPiKPEt,1.0);
       }
+      if(kIsOfflineMB){
+       FillHisto1D("SimHadEtNDMB",fSimHadEt,1.0);
+       FillHisto1D("SimTotEtNDMB",fSimHadEt,1.0);
+       FillHisto1D("NEventsNDMB",0.5,1);
+       FillHisto1D("SimPiKPEtNDMB",fSimPiKPEt,1.0);
+      }
     }
     if(fDataSet!=20100 && AliPWG0Helper::GetEventProcessType(mcEvent->Header()) == AliPWG0Helper::kSD){
       FillHisto1D("SimHadEtSD",fSimHadEt,1.0);
@@ -1683,6 +1690,12 @@ Int_t AliAnalysisHadEtMonteCarlo::AnalyseEvent(AliVEvent* ev)
        FillHisto1D("NEventsSDV0AND",0.5,1);
        FillHisto1D("SimPiKPEtSDV0AND",fSimPiKPEt,1.0);
       }
+      if(kIsOfflineMB){
+       FillHisto1D("SimHadEtSDMB",fSimHadEt,1.0);
+       FillHisto1D("SimTotEtSDMB",fSimHadEt,1.0);
+       FillHisto1D("NEventsSDMB",0.5,1);
+       FillHisto1D("SimPiKPEtSDMB",fSimPiKPEt,1.0);
+      }
     }
     if(fDataSet!=20100 && AliPWG0Helper::GetEventProcessType(mcEvent->Header()) == AliPWG0Helper::kDD){
       FillHisto1D("SimHadEtDD",fSimHadEt,1.0);
@@ -1695,6 +1708,12 @@ Int_t AliAnalysisHadEtMonteCarlo::AnalyseEvent(AliVEvent* ev)
        FillHisto1D("NEventsDDV0AND",0.5,1);
        FillHisto1D("SimPiKPEtDDV0AND",fSimPiKPEt,1.0);
       }
+      if(kIsOfflineMB){
+       FillHisto1D("SimHadEtDDMB",fSimHadEt,1.0);
+       FillHisto1D("SimTotEtDDMB",fSimHadEt,1.0);
+       FillHisto1D("NEventsDDMB",0.5,1);
+       FillHisto1D("SimPiKPEtDDMB",fSimPiKPEt,1.0);
+      }
     }
     if(fCentBin != -1){//if we have Pb+Pb and a centrality bin was found
       if(fSimTotEt>0.0) FillHisto1D(Form("SimTotEtCB%i",fCentBin),fSimTotEt,1.0);
@@ -2104,20 +2123,29 @@ void AliAnalysisHadEtMonteCarlo::CreateHistograms(){
   CreateHisto1D("SimPiKPEtNDV0AND","Simulated #pi,K,p E_{T}","Simulated #pi,K,p E_{T}","Number of events",nbinsEt*4,minEtPiKP,maxEtPiKP);
   CreateHisto1D("SimPiKPEtDDV0AND","Simulated #pi,K,p E_{T}","Simulated #pi,K,p E_{T}","Number of events",nbinsEt*4,minEtPiKP,maxEtPiKP);
   CreateHisto1D("SimPiKPEtSDV0AND","Simulated #pi,K,p E_{T}","Simulated #pi,K,p E_{T}","Number of events",nbinsEt*4,minEtPiKP,maxEtPiKP);
+  CreateHisto1D("SimPiKPEtNDMB","Simulated #pi,K,p E_{T}","Simulated #pi,K,p E_{T}","Number of events",nbinsEt*4,minEtPiKP,maxEtPiKP);
+  CreateHisto1D("SimPiKPEtDDMB","Simulated #pi,K,p E_{T}","Simulated #pi,K,p E_{T}","Number of events",nbinsEt*4,minEtPiKP,maxEtPiKP);
+  CreateHisto1D("SimPiKPEtSDMB","Simulated #pi,K,p E_{T}","Simulated #pi,K,p E_{T}","Number of events",nbinsEt*4,minEtPiKP,maxEtPiKP);
   CreateHisto1D("SimTotEt","Simulated Total E_{T}","Simulated Total E_{T}","Number of events",nbinsEt*4,minEt,maxEt);
   CreateHisto1D("SimHadEt","Simulated Hadronic E_{T}","Simulated Hadronic E_{T}","Number of events",nbinsEt*4,minEt,maxEt);
   CreateHisto1D("SimTotEtND","Simulated Total E_{T}","Simulated Total E_{T} for non-diffractive events","Number of events",nbinsEt*4,minEt,maxEt);
   CreateHisto1D("SimHadEtND","Simulated Hadronic E_{T}","Simulated Hadronic E_{T} for non-diffractive events","Number of events",nbinsEt*4,minEt,maxEt);
   CreateHisto1D("SimTotEtNDV0AND","Simulated Total E_{T}","Simulated Total E_{T} for non-diffractive events","Number of events",nbinsEt*4,minEt,maxEt);
   CreateHisto1D("SimHadEtNDV0AND","Simulated Hadronic E_{T}","Simulated Hadronic E_{T} for non-diffractive events","Number of events",nbinsEt*4,minEt,maxEt);
+  CreateHisto1D("SimTotEtNDMB","Simulated Total E_{T}","Simulated Total E_{T} for non-diffractive events","Number of events",nbinsEt*4,minEt,maxEt);
+  CreateHisto1D("SimHadEtNDMB","Simulated Hadronic E_{T}","Simulated Hadronic E_{T} for non-diffractive events","Number of events",nbinsEt*4,minEt,maxEt);
   CreateHisto1D("SimTotEtSD","Simulated Total E_{T}","Simulated Total E_{T} for singly diffractive events","Number of events",nbinsEt*4,minEt,maxEt);
   CreateHisto1D("SimHadEtSD","Simulated Hadronic E_{T}","Simulated Hadronic E_{T} for singly diffractive events","Number of events",nbinsEt*4,minEt,maxEt);
   CreateHisto1D("SimTotEtSDV0AND","Simulated Total E_{T}","Simulated Total E_{T} for singly diffractive events","Number of events",nbinsEt*4,minEt,maxEt);
   CreateHisto1D("SimHadEtSDV0AND","Simulated Hadronic E_{T}","Simulated Hadronic E_{T} for singly diffractive events","Number of events",nbinsEt*4,minEt,maxEt);
+  CreateHisto1D("SimTotEtSDMB","Simulated Total E_{T}","Simulated Total E_{T} for singly diffractive events","Number of events",nbinsEt*4,minEt,maxEt);
+  CreateHisto1D("SimHadEtSDMB","Simulated Hadronic E_{T}","Simulated Hadronic E_{T} for singly diffractive events","Number of events",nbinsEt*4,minEt,maxEt);
   CreateHisto1D("SimTotEtDD","Simulated Total E_{T}","Simulated Total E_{T} for doubly diffractive events","Number of events",nbinsEt*4,minEt,maxEt);
   CreateHisto1D("SimHadEtDD","Simulated Hadronic E_{T}","Simulated Hadronic E_{T} for doubly diffractive events","Number of events",nbinsEt*4,minEt,maxEt);
   CreateHisto1D("SimTotEtDDV0AND","Simulated Total E_{T}","Simulated Total E_{T} for doubly diffractive events","Number of events",nbinsEt*4,minEt,maxEt);
   CreateHisto1D("SimHadEtDDV0AND","Simulated Hadronic E_{T}","Simulated Hadronic E_{T} for doubly diffractive events","Number of events",nbinsEt*4,minEt,maxEt);
+  CreateHisto1D("SimTotEtDDMB","Simulated Total E_{T}","Simulated Total E_{T} for doubly diffractive events","Number of events",nbinsEt*4,minEt,maxEt);
+  CreateHisto1D("SimHadEtDDMB","Simulated Hadronic E_{T}","Simulated Hadronic E_{T} for doubly diffractive events","Number of events",nbinsEt*4,minEt,maxEt);
   if(fDataSet==20100){
     Int_t width = 5;
     if(fNCentBins<21) width = 10;
@@ -2313,6 +2341,9 @@ void AliAnalysisHadEtMonteCarlo::CreateHistograms(){
   CreateIntHisto1D("NEventsSDV0AND","Number of events","number of singly diffractive events","Number of events",1,0,1);
   CreateIntHisto1D("NEventsDDV0AND","Number of events","number of doubly diffractive events","Number of events",1,0,1);
   CreateIntHisto1D("NEventsNDV0AND","Number of events","number of non-diffractive events","Number of events",1,0,1);
+  CreateIntHisto1D("NEventsSDMB","Number of events","number of singly diffractive events","Number of events",1,0,1);
+  CreateIntHisto1D("NEventsDDMB","Number of events","number of doubly diffractive events","Number of events",1,0,1);
+  CreateIntHisto1D("NEventsNDMB","Number of events","number of non-diffractive events","Number of events",1,0,1);
   if( !fRunLightweight){
     CreateResolutionPtHisto2D("presolutionTPC","p resolution","p^{rec}","(p^{sim}-p^{rec})/p^{rec}");
     CreateResolutionPtHisto2D("pTresolutionTPC","p_{T} resolution","p_{T}^{rec}","(p_{T}^{sim}-p_{T}^{rec})/p_{T}^{rec}");

diff --git a/PWGLF/totEt/AliAnalysisHadEtMonteCarlo.h b/PWGLF/totEt/AliAnalysisHadEtMonteCarlo.h
index aca44f0789b651202a90c478629b27fcc948492e..9328feb0d505807c4c0acb1ece103ef578906d3f 100644 (file)
--- a/PWGLF/totEt/AliAnalysisHadEtMonteCarlo.h
+++ b/PWGLF/totEt/AliAnalysisHadEtMonteCarlo.h
@@ -125,6 +125,7 @@ public:
     void DoTriggerChecks(){kDoTriggerChecks=kTRUE;}
     void DoTriggerChecksOnly(){kDoTriggerChecksOnly=kTRUE;}
     void SetIsOfflineV0AND(Bool_t val){kIsOfflineV0AND = val;}
+    void SetIsOfflineMB(Bool_t val){kIsOfflineMB = val;}
  protected:
 
  private:
@@ -150,6 +151,7 @@ public:
     Bool_t fBaryonEnhancement;//Turns on and off baryon enhancement
     Bool_t fUseRecoPt;//Switch to decide whether to use reconstructed or simulated pT for reconstructed tracks
     Bool_t kIsOfflineV0AND;//Boolean to keep track of whether or not this matches the offline trigger
+    Bool_t kIsOfflineMB;//Boolean to keep track of whether or not this matches the offline trigger
     Bool_t kDoTriggerChecks;//Boolean to keep track of whether or not I want to let the physics selection do the work
     Bool_t kDoTriggerChecksOnly;//Boolean to keep track of whether or not I want to let the physics selection do the work
 

diff --git a/PWGLF/totEt/AliAnalysisTaskHadEt.cxx b/PWGLF/totEt/AliAnalysisTaskHadEt.cxx
index 15a1c7d36d25810f0c3e131b4b1afd42fe4314e8..e3fc58f9917513c9cd990cff0e712d41477366e4 100644 (file)
--- a/PWGLF/totEt/AliAnalysisTaskHadEt.cxx
+++ b/PWGLF/totEt/AliAnalysisTaskHadEt.cxx
@@ -42,6 +42,7 @@ ClassImp(AliAnalysisTaskHadEt)
        ,fMCAnalysis(0)
        ,fIsSim(isMc)
        ,kIsOfflineV0AND(0)
+       ,kIsOfflineMB(0)
 {
     // Constructor
   fMCConfigFile = mcConfigFile;
@@ -179,8 +180,10 @@ if (!fESDEvent) {
  AliTriggerAnalysis *fTriggerAnalysis = new AliTriggerAnalysis();
 
   kIsOfflineV0AND = fTriggerAnalysis->IsOfflineTriggerFired(fESDEvent, AliTriggerAnalysis::kV0AND);  
+  kIsOfflineMB = fTriggerAnalysis->IsOfflineTriggerFired(fESDEvent, AliTriggerAnalysis::kMB1);  
   fRecAnalysis->SetIsOfflineV0AND(kIsOfflineV0AND);
   fMCAnalysis->SetIsOfflineV0AND(kIsOfflineV0AND);
+  fMCAnalysis->SetIsOfflineMB(kIsOfflineMB);
 
   Int_t eventtype =    AliPWG0Helper::kInvalidProcess;
   if(fIsSim && fRecAnalysis->DataSet()!=20100) eventtype = (Int_t) AliPWG0Helper::GetEventProcessType(MCEvent()->Header());

diff --git a/PWGLF/totEt/AliAnalysisTaskHadEt.h b/PWGLF/totEt/AliAnalysisTaskHadEt.h
index 41f5b9c4ee595185fd8cf619d6f95d426d66a36b..7ec47d2a2aa620088c3bbfd1b10301cf21c95419 100644 (file)
--- a/PWGLF/totEt/AliAnalysisTaskHadEt.h
+++ b/PWGLF/totEt/AliAnalysisTaskHadEt.h
@@ -39,6 +39,7 @@ private:
     AliAnalysisHadEtMonteCarlo *fMCAnalysis; // MC\r
     Bool_t fIsSim;//Boolean to keep track of whether or not this is running on simulations\r
     Bool_t kIsOfflineV0AND;//Boolean to keep track of whether or not this matches the offline trigger\r
+    Bool_t kIsOfflineMB;//Boolean to keep track of whether or not this matches the offline trigger\r
     ClassDef(AliAnalysisTaskHadEt, 2); // example of analysis\r
 };\r
 \r