set minimum pT for dicing (M. Verweij)
authorkleinb <kleinb@f7af4fe6-9843-0410-8265-dc069ae4e863>
Thu, 21 Jun 2012 08:29:23 +0000 (08:29 +0000)
committerkleinb <kleinb@f7af4fe6-9843-0410-8265-dc069ae4e863>
Thu, 21 Jun 2012 08:29:23 +0000 (08:29 +0000)
JETAN/AliAnalysisTaskJetCluster.cxx
JETAN/AliAnalysisTaskJetCluster.h
JETAN/DEV/AliAnalysisTaskJetCluster.cxx
JETAN/DEV/AliAnalysisTaskJetCluster.h

index 3e1bf13..26f7c57 100644 (file)
@@ -138,6 +138,7 @@ AliAnalysisTaskJetCluster::AliAnalysisTaskJetCluster():
   fhEffH3(0x0),
   fUseTrPtResolutionSmearing(kFALSE),
   fUseDiceEfficiency(kFALSE),
+  fDiceEfficiencyMinPt(0.),
   fUseTrPtResolutionFromOADB(kFALSE),
   fUseTrEfficiencyFromOADB(kFALSE),
   fPathTrPtResolution(""),
@@ -275,6 +276,7 @@ AliAnalysisTaskJetCluster::AliAnalysisTaskJetCluster(const char* name):
   fhEffH3(0x0),
   fUseTrPtResolutionSmearing(kFALSE),
   fUseDiceEfficiency(kFALSE),
+  fDiceEfficiencyMinPt(0.),
   fUseTrPtResolutionFromOADB(kFALSE),
   fUseTrEfficiencyFromOADB(kFALSE),
   fPathTrPtResolution(""),
@@ -926,7 +928,7 @@ void AliAnalysisTaskJetCluster::UserExec(Option_t */*option*/)
       
       Double_t sumEff = eff[0]+eff[1]+eff[2];
       fp1Efficiency->Fill(vp->Pt(),sumEff);
-      if(rnd>sumEff) continue;
+      if(rnd>sumEff && pT > fDiceEfficiencyMinPt) continue;
 
       if(fUseTrPtResolutionSmearing) {
        //Smear momentum of generated particle
index f055d3c..14b6c55 100644 (file)
@@ -50,7 +50,6 @@ class AliAnalysisTaskJetCluster : public AliAnalysisTaskSE
     virtual Bool_t Notify();
 
     
-
     virtual void SetAODTrackInput(Bool_t b){fUseAODTrackInput = b;}
     virtual void SetAODMCInput(Bool_t b){fUseAODMCInput = b;}
     virtual void SetEventSelection(Bool_t b){fEventSelection = b;}
@@ -89,6 +88,7 @@ class AliAnalysisTaskJetCluster : public AliAnalysisTaskSE
     virtual void SetChangeEfficiencyFraction(Double_t p) {fChangeEfficiencyFraction = p;}
     virtual void SetSmearResolution(Bool_t b){fUseTrPtResolutionSmearing = b;} 
     virtual void SetDiceEfficiency(Bool_t b){fUseDiceEfficiency = b;} 
+    virtual void SetDiceEfficiencyMinPt(Double_t pt) {fDiceEfficiencyMinPt = pt;}
     virtual void SetMomentumResolutionHybrid(TProfile *p1, TProfile *p2, TProfile *p3);
     virtual void SetEfficiencyHybrid(TH1 *h1, TH1 *h2, TH1 *h3);
 
@@ -186,6 +186,7 @@ class AliAnalysisTaskJetCluster : public AliAnalysisTaskSE
     TH1      *fhEffH3;        // Efficiency for Spectra Hybrid Category 3
     Bool_t    fUseTrPtResolutionSmearing;  // Apply momentum smearing on track level
     Bool_t    fUseDiceEfficiency;          // Apply efficiency on track level by dicing
+    Double_t  fDiceEfficiencyMinPt;        // Only do efficiency dicing for tracks above this pt
     Bool_t fUseTrPtResolutionFromOADB;     // Load track pt resolution root file from OADB path
     Bool_t fUseTrEfficiencyFromOADB;       // Load tracking efficiency root file from OADB path
     TString fPathTrPtResolution;           // OADB path to root file
index 3e1bf13..26f7c57 100644 (file)
@@ -138,6 +138,7 @@ AliAnalysisTaskJetCluster::AliAnalysisTaskJetCluster():
   fhEffH3(0x0),
   fUseTrPtResolutionSmearing(kFALSE),
   fUseDiceEfficiency(kFALSE),
+  fDiceEfficiencyMinPt(0.),
   fUseTrPtResolutionFromOADB(kFALSE),
   fUseTrEfficiencyFromOADB(kFALSE),
   fPathTrPtResolution(""),
@@ -275,6 +276,7 @@ AliAnalysisTaskJetCluster::AliAnalysisTaskJetCluster(const char* name):
   fhEffH3(0x0),
   fUseTrPtResolutionSmearing(kFALSE),
   fUseDiceEfficiency(kFALSE),
+  fDiceEfficiencyMinPt(0.),
   fUseTrPtResolutionFromOADB(kFALSE),
   fUseTrEfficiencyFromOADB(kFALSE),
   fPathTrPtResolution(""),
@@ -926,7 +928,7 @@ void AliAnalysisTaskJetCluster::UserExec(Option_t */*option*/)
       
       Double_t sumEff = eff[0]+eff[1]+eff[2];
       fp1Efficiency->Fill(vp->Pt(),sumEff);
-      if(rnd>sumEff) continue;
+      if(rnd>sumEff && pT > fDiceEfficiencyMinPt) continue;
 
       if(fUseTrPtResolutionSmearing) {
        //Smear momentum of generated particle
index f055d3c..14b6c55 100644 (file)
@@ -50,7 +50,6 @@ class AliAnalysisTaskJetCluster : public AliAnalysisTaskSE
     virtual Bool_t Notify();
 
     
-
     virtual void SetAODTrackInput(Bool_t b){fUseAODTrackInput = b;}
     virtual void SetAODMCInput(Bool_t b){fUseAODMCInput = b;}
     virtual void SetEventSelection(Bool_t b){fEventSelection = b;}
@@ -89,6 +88,7 @@ class AliAnalysisTaskJetCluster : public AliAnalysisTaskSE
     virtual void SetChangeEfficiencyFraction(Double_t p) {fChangeEfficiencyFraction = p;}
     virtual void SetSmearResolution(Bool_t b){fUseTrPtResolutionSmearing = b;} 
     virtual void SetDiceEfficiency(Bool_t b){fUseDiceEfficiency = b;} 
+    virtual void SetDiceEfficiencyMinPt(Double_t pt) {fDiceEfficiencyMinPt = pt;}
     virtual void SetMomentumResolutionHybrid(TProfile *p1, TProfile *p2, TProfile *p3);
     virtual void SetEfficiencyHybrid(TH1 *h1, TH1 *h2, TH1 *h3);
 
@@ -186,6 +186,7 @@ class AliAnalysisTaskJetCluster : public AliAnalysisTaskSE
     TH1      *fhEffH3;        // Efficiency for Spectra Hybrid Category 3
     Bool_t    fUseTrPtResolutionSmearing;  // Apply momentum smearing on track level
     Bool_t    fUseDiceEfficiency;          // Apply efficiency on track level by dicing
+    Double_t  fDiceEfficiencyMinPt;        // Only do efficiency dicing for tracks above this pt
     Bool_t fUseTrPtResolutionFromOADB;     // Load track pt resolution root file from OADB path
     Bool_t fUseTrEfficiencyFromOADB;       // Load tracking efficiency root file from OADB path
     TString fPathTrPtResolution;           // OADB path to root file