trigger differentiated QA histos
authorkleinb <kleinb@f7af4fe6-9843-0410-8265-dc069ae4e863>
Mon, 30 Apr 2012 12:13:22 +0000 (12:13 +0000)
committerkleinb <kleinb@f7af4fe6-9843-0410-8265-dc069ae4e863>
Mon, 30 Apr 2012 12:13:22 +0000 (12:13 +0000)
PWGJE/AliAnalysisTaskJetServices.cxx
PWGJE/AliAnalysisTaskJetServices.h
PWGJE/AliAnalysisTaskJetSpectrum2.cxx

index fe7e3f8..ad765ac 100644 (file)
@@ -114,8 +114,8 @@ AliAnalysisTaskJetServices::AliAnalysisTaskJetServices():
   fh1EventCutInfoESD(0),
   fh1CentralityESD(0),
   fh1Centrality(0),
-  fh1ReducedTrigger(0),
   fh1RP(0),
+  fh2ReducedTrigger(0),
   fh2CentralityTriggerESD(0),
   fh2CentralityTrigger(0),
   fh2TriggerCount(0x0),
@@ -187,8 +187,8 @@ AliAnalysisTaskJetServices::AliAnalysisTaskJetServices(const char* name):
   fh1EventCutInfoESD(0),
   fh1CentralityESD(0),
   fh1Centrality(0),
-  fh1ReducedTrigger(0),
   fh1RP(0),
+  fh2ReducedTrigger(0),
   fh2CentralityTriggerESD(0),
   fh2CentralityTrigger(0),
   fh2TriggerCount(0x0),
@@ -300,8 +300,8 @@ void AliAnalysisTaskJetServices::UserCreateOutputObjects()
   fh1RP = new TH1F("fh1RP","RP;#Psi",440, -1.*TMath::Pi(), 2.*TMath::Pi());
   fHistList->Add(fh1RP);
 
-  fh1ReducedTrigger = new TH1F("fh1ReducedTrigger","red trigger;red trigger",1<<fNTrigger,-0.5,(1<<fNTrigger)-0.5);
-  fHistList->Add(fh1ReducedTrigger);
+  fh2ReducedTrigger = new TH2F("fh2ReducedTrigger","red trigger;red trigger;centrality",1<<fNTrigger,-0.5,(1<<fNTrigger)-0.5,100,-0.5,99.5);
+  fHistList->Add(fh2ReducedTrigger);
 
   fh2CentralityTriggerESD = new TH2F("fh2CentralityTriggerESD",";cent;trigger no",103,-1,102,fNTrigger,-0.5,fNTrigger-0.5);
   fHistList->Add(fh2CentralityTriggerESD);
@@ -599,7 +599,7 @@ void AliAnalysisTaskJetServices::UserExec(Option_t */*option*/)
            ir |= (1<<it);
          }
        }
-       fh1ReducedTrigger->Fill(ir);
+       fh2ReducedTrigger->Fill(ir,cent);
       }
     }
   }
index 376fde8..4ffc606 100644 (file)
@@ -157,8 +157,8 @@ class AliAnalysisTaskJetServices : public AliAnalysisTaskSE
     TH1F*         fh1EventCutInfoESD;  //! Masks that satisfy fSelectionInfo
     TH1F*         fh1CentralityESD;    //! centrality 
     TH1F*         fh1Centrality;       //! centrality 
-    TH1F*         fh1ReducedTrigger;   //! reduced trigger count 
     TH1F*         fh1RP;               //! RP distribution
+    TH2F*         fh2ReducedTrigger;   //! reduced trigger count 
     TH2F*         fh2CentralityTriggerESD;    //! centrality 
     TH2F*         fh2CentralityTrigger;       //! centrality 
     TH2F*         fh2TriggerCount;     //! number of fire triggers in each case
index c36e1ac..eecd865 100644 (file)
@@ -504,21 +504,21 @@ void AliAnalysisTaskJetSpectrum2::UserCreateOutputObjects()
     fHistList->Add(fhnJetPt[ij]);
     
     // Bins:  Jet number: pTJet, cent, eta, phi, Area.   total bins = 9.72 M
-    const Int_t nBinsSparse2 = 6;
-    Int_t nBins2[nBinsSparse2] = {     kMaxJets+1, 25,   5,  18,             360, 10};
+    const Int_t nBinsSparse2 = 7;
+    Int_t nBins2[nBinsSparse2] = {     kMaxJets+1, 25,   5,  18,             180, 10,fNTrigger};
     if(cJetBranch.Contains("RandomCone")){
       nBins2[5] = 1;
     }
-    const Double_t xmin2[nBinsSparse2]  = {        -0.5,  0,   0,-0.9,              0,  0.};
-    const Double_t xmax2[nBinsSparse2]  = {kMaxJets+0.5,250, 100, 0.9, 2.*TMath::Pi(),1.0};
+    const Double_t xmin2[nBinsSparse2]  = {        -0.5,  0,   0,-0.9,              0,  0.,-0.5};
+    const Double_t xmax2[nBinsSparse2]  = {kMaxJets+0.5,250, 100, 0.9, 2.*TMath::Pi(),1.0,fNTrigger-0.5};
     fhnJetPtQA[ij] = new THnSparseF(Form("fhnJetPtQA%s",cAdd.Data()),";jet number;p_{T,jet};cent;#eta;#phi;area",nBinsSparse2,nBins2,xmin2,xmax2);
     fHistList->Add(fhnJetPtQA[ij]);
     
     // Bins:track number  pTtrack, cent, mult, RP.   total bins = 224 k
-    const Int_t nBinsSparse3 = 5;
-    const Int_t nBins3[nBinsSparse3] = {       2,    100,     10,  20,    fNRPBins};
-    const Double_t xmin3[nBinsSparse3]  = { -0.5,     0,   0,      0,        -0.5};
-    const Double_t xmax3[nBinsSparse3]  = { 1.5,    200, 100,   4000,fNRPBins-0.5};  
+    const Int_t nBinsSparse3 = 6;
+    const Int_t nBins3[nBinsSparse3] = {       2,    100,     10,  20,    fNRPBins,fNTrigger};
+    const Double_t xmin3[nBinsSparse3]  = { -0.5,      0,   0,      0,        -0.5,-0.5};
+    const Double_t xmax3[nBinsSparse3]  = { 1.5,     200, 100,   4000,fNRPBins-0.5,fNTrigger-0.5};  
     
       // change the binning ot the pT axis:
     Double_t *xPt3 = new Double_t[nBins3[1]+1];