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

[blackhole at projects.hepforge.org]

Author: buckley
Date: Tue Sep 29 17:00:10 2009
New Revision: 1840

Log:
Fixed dphi distributions: only copy temp histo enties to long-term profiles if jet1 pT condition is met (otherwise lots of zero-fill temp profiles drag the average down) and get rid of the bin width factor (as the temp histos are also profiles).

Modified:
   trunk/src/Analyses/CDF_2001_S4751469.cc

Modified: trunk/src/Analyses/CDF_2001_S4751469.cc
==============================================================================
--- trunk/src/Analyses/CDF_2001_S4751469.cc	Tue Sep 29 16:39:47 2009	(r1839)
+++ trunk/src/Analyses/CDF_2001_S4751469.cc	Tue Sep 29 17:00:10 2009	(r1840)
@@ -181,19 +181,20 @@
 
       // Update the "proper" dphi profile histograms
       for (int i= 0; i < 50; i++) {
-        // pT > 2 GeV
-        _numvsDeltaPhi2->fill(hist_num_dphi_2.binMean(i), hist_num_dphi_2.binHeight(i)*hist_num_dphi_2.axis().binWidth(i));
-        _pTvsDeltaPhi2->fill(hist_pt_dphi_2.binMean(i), hist_pt_dphi_2.binHeight(i)*hist_pt_dphi_2.axis().binWidth(i)); 
-        // pT > 5 GeV
-        /// @todo Normalisation?
-        _numvsDeltaPhi5->fill(hist_num_dphi_5.binMean(i), hist_num_dphi_5.binHeight(i)*hist_num_dphi_5.axis().binWidth(i));
-        _pTvsDeltaPhi5->fill(hist_pt_dphi_5.binMean(i), hist_pt_dphi_5.binHeight(i)*hist_pt_dphi_5.axis().binWidth(i));
-        // pT > 30 GeV
-        /// @todo Normalisation?
-        _numvsDeltaPhi30->fill(hist_num_dphi_30.binMean(i), hist_num_dphi_30.binHeight(i)*hist_num_dphi_30.axis().binWidth(i));
-        _pTvsDeltaPhi30->fill(hist_pt_dphi_30.binMean(i), hist_pt_dphi_30.binHeight(i)*hist_pt_dphi_30.axis().binWidth(i));
+        if (ptLead/GeV > 2.0) {
+          _numvsDeltaPhi2->fill(hist_num_dphi_2.binMean(i), hist_num_dphi_2.binHeight(i));
+          _pTvsDeltaPhi2->fill(hist_pt_dphi_2.binMean(i), hist_pt_dphi_2.binHeight(i)); 
+        }
+        if (ptLead/GeV > 5.0) {
+          _numvsDeltaPhi5->fill(hist_num_dphi_5.binMean(i), hist_num_dphi_5.binHeight(i));
+          _pTvsDeltaPhi5->fill(hist_pt_dphi_5.binMean(i), hist_pt_dphi_5.binHeight(i));
+        }
+        if (ptLead/GeV > 30.0) {
+          _numvsDeltaPhi30->fill(hist_num_dphi_30.binMean(i), hist_num_dphi_30.binHeight(i));
+          _pTvsDeltaPhi30->fill(hist_pt_dphi_30.binMean(i), hist_pt_dphi_30.binHeight(i));
+        }
       }
-
+      
       // Log some event details about pT
       getLog() << Log::DEBUG 
                << "pT [lead; twd, away, trans] = ["