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

[blackhole at projects.hepforge.org]

Author: fsiegert
Date: Thu Feb 10 08:29:21 2011
New Revision: 2920

Log:
Bugfixes in ATLAS inclusive+dijet XS analysis. Now the MC isn't off by a factor of 0.5 (inclusive jet pT) or two orders of magnitude anymore ;-)

Modified:
   trunk/src/Analyses/ATLAS_2010_S8817804.cc

Modified: trunk/src/Analyses/ATLAS_2010_S8817804.cc
==============================================================================
--- trunk/src/Analyses/ATLAS_2010_S8817804.cc	Wed Feb  9 15:47:24 2011	(r2919)
+++ trunk/src/Analyses/ATLAS_2010_S8817804.cc	Thu Feb 10 08:29:21 2011	(r2920)
@@ -17,7 +17,6 @@
     {
       setBeams(PROTON, PROTON);
       setNeedsCrossSection(true);
-      _sumW = 0;
     }
 
 
@@ -82,7 +81,6 @@
           MSG_DEBUG("Could not find two suitable leading jets");
           continue;
         }
-        _sumW += evt.weight();
 
         const double rap1 = leadjets[0].rapidity();
         const double rap2 = leadjets[1].rapidity();
@@ -100,18 +98,16 @@
 
     void finalize() {
       for (size_t alg = 0; alg < 2; ++alg) {
-        _pThistos[alg].scale(crossSectionPerEvent()/picobarn, this);
-        _massVsY[alg].scale(crossSection()/_sumW/picobarn, this);
-        _chiVsMass[alg].scale(crossSection()/_sumW/picobarn, this);
+        // factor 0.5 needed because it is differential in dy and not d|y|
+        _pThistos[alg].scale(0.5*crossSectionPerEvent()/picobarn, this);
+        _massVsY[alg].scale(crossSectionPerEvent()/picobarn, this);
+        _chiVsMass[alg].scale(crossSectionPerEvent()/picobarn, this);
       }
     }
 
 
   private:
 
-    /// Counter for weights passing the dijet requirement cut
-    double _sumW;
-
     /// The inclusive pT spectrum for akt6 and akt4 jets (array index is jet type from enum above)
     BinnedHistogram<double> _pThistos[2];