[Rivet-svn] r2924 - trunk/src/Analyses

[blackhole at projects.hepforge.org]

Author: hoeth
Date: Thu Feb 10 14:18:49 2011
New Revision: 2924

Log:
Fix pT cuts for dphi plots at 900 GeV in ATLAS_2010_S8894728 track
based UE analysis, since they differ from the 7 TeV cuts.

Modified:
   trunk/src/Analyses/ATLAS_2010_S8894728.cc

Modified: trunk/src/Analyses/ATLAS_2010_S8894728.cc
==============================================================================
--- trunk/src/Analyses/ATLAS_2010_S8894728.cc	Thu Feb 10 14:11:31 2011	(r2923)
+++ trunk/src/Analyses/ATLAS_2010_S8894728.cc	Thu Feb 10 14:18:49 2011	(r2924)
@@ -217,21 +217,29 @@
       // The values tabulated in the note are for an (undefined) signed Delta(phi) rather than
       // |Delta(phi)| and so differ by a factor of 2: we have to actually norm for angular range = 2pi
       const size_t nbins = binEdges(13,1,1).size() - 1;
+      std::vector<double> ptcut;
+      if (fuzzyEquals(sqrtS(), 900*GeV)) {
+        ptcut += 1.0; ptcut += 1.5; ptcut += 2.0; ptcut += 2.5;
+      }
+      else if (fuzzyEquals(sqrtS(), 7*TeV)) {
+        ptcut += 1.0; ptcut += 2.0; ptcut += 3.0; ptcut += 5.0;
+      }
+      assert(ptcut.size == 4);
       for (size_t i = 0; i < nbins; ++i) {
         // First Nch
         const double binmean_num = hist_num_dphi_500.binMean(i);
         const double binvalue_num = hist_num_dphi_500.binHeight(i)/hist_num_dphi_500.axis().binWidth(i)/10.0;
-        if (pTlead/GeV >= 1.0) _hist_N_vs_dPhi_1_500->fill(binmean_num, binvalue_num, weight);
-        if (pTlead/GeV >= 2.0) _hist_N_vs_dPhi_2_500->fill(binmean_num, binvalue_num, weight);
-        if (pTlead/GeV >= 3.0) _hist_N_vs_dPhi_3_500->fill(binmean_num, binvalue_num, weight);
-        if (pTlead/GeV >= 5.0) _hist_N_vs_dPhi_5_500->fill(binmean_num, binvalue_num, weight);
+        if (pTlead/GeV >= ptcut[0]) _hist_N_vs_dPhi_1_500->fill(binmean_num, binvalue_num, weight);
+        if (pTlead/GeV >= ptcut[1]) _hist_N_vs_dPhi_2_500->fill(binmean_num, binvalue_num, weight);
+        if (pTlead/GeV >= ptcut[2]) _hist_N_vs_dPhi_3_500->fill(binmean_num, binvalue_num, weight);
+        if (pTlead/GeV >= ptcut[3]) _hist_N_vs_dPhi_5_500->fill(binmean_num, binvalue_num, weight);
         // Then pT
         const double binmean_pt = hist_pt_dphi_500.binMean(i);
         const double binvalue_pt = hist_pt_dphi_500.binHeight(i)/hist_pt_dphi_500.axis().binWidth(i)/10.0;
-        if (pTlead/GeV >= 1.0) _hist_pT_vs_dPhi_1_500->fill(binmean_pt, binvalue_pt, weight);
-        if (pTlead/GeV >= 2.0) _hist_pT_vs_dPhi_2_500->fill(binmean_pt, binvalue_pt, weight);
-        if (pTlead/GeV >= 3.0) _hist_pT_vs_dPhi_3_500->fill(binmean_pt, binvalue_pt, weight);
-        if (pTlead/GeV >= 5.0) _hist_pT_vs_dPhi_5_500->fill(binmean_pt, binvalue_pt, weight);
+        if (pTlead/GeV >= ptcut[0]) _hist_pT_vs_dPhi_1_500->fill(binmean_pt, binvalue_pt, weight);
+        if (pTlead/GeV >= ptcut[1]) _hist_pT_vs_dPhi_2_500->fill(binmean_pt, binvalue_pt, weight);
+        if (pTlead/GeV >= ptcut[2]) _hist_pT_vs_dPhi_3_500->fill(binmean_pt, binvalue_pt, weight);
+        if (pTlead/GeV >= ptcut[3]) _hist_pT_vs_dPhi_5_500->fill(binmean_pt, binvalue_pt, weight);
       }