[Rivet-svn] rivet: Fix a few things and mark D0_1996_S3214044 jet analysis v...

[Rivet Mercurial]

details:   https://rivet.hepforge.org/hg/rivet/rev/ed82e65e6f77
branches:  release-2-5-x
changeset: 5389:ed82e65e6f77
user:      Holger Schulz <holger.schulz at durham.ac.uk>
date:      Fri Aug 05 14:03:13 2016 +0100
description:
Fix a few things and mark D0_1996_S3214044 jet analysis validated

diffs (truncated from 57 to 50 lines):

--- a/ChangeLog	Fri Aug 05 11:09:38 2016 +0100
+++ b/ChangeLog	Fri Aug 05 14:03:13 2016 +0100
@@ -5,6 +5,11 @@
 
 	* Mark CDF_2001_S4563131 validated
 
+	* D0_1996_S3214044 --- cut on jet Et rather than pT, fix filling of
+	costheta and theta plots, mark validated. Concerning the jet
+	algorithm, I tried with the implementation of fastjet
+	fastjet/D0RunIConePlugin.hh but that really does not help.
+
 2016-08-04  Holger Schulz  <holger.schulz at cern.ch>
 
 	* Use Jet mass and energy smearing in CDF_1996_S310 ... jet properties
--- a/data/anainfo/D0_1996_S3214044.info	Fri Aug 05 11:09:38 2016 +0100
+++ b/data/anainfo/D0_1996_S3214044.info	Fri Aug 05 14:03:13 2016 +0100
@@ -4,7 +4,7 @@
 Experiment: D0
 Collider: Tevatron Run 1
 SpiresID: 3214044
-Status: UNVALIDATED
+Status: VALIDATED
 Authors:
  - Frank Siegert <frank.siegert at cern.ch>
 References:
--- a/src/Analyses/D0_1996_S3214044.cc	Fri Aug 05 11:09:38 2016 +0100
+++ b/src/Analyses/D0_1996_S3214044.cc	Fri Aug 05 14:03:13 2016 +0100
@@ -26,7 +26,8 @@
     void init() {
       const FinalState fs;
       declare(fs, "FS");
-      /// @todo Use correct jet algorithm
+      /// @todo Use correct jet algorithm --- tried FJ3 D0RunICone but does
+      // not look as good as the Run2 cone alg used here
       declare(FastJets(fs, FastJets::D0ILCONE, 0.7), "ConeJets");
 
       _h_3j_x3 = bookHisto1D(1, 1, 1);
@@ -67,7 +68,7 @@
       const double weight = event.weight();
 
       Jets jets_in = apply<FastJets>(event, "ConeJets")
-        .jets(cmpMomByEt, Cuts::pT > 20*GeV && Cuts::abseta < 3);
+        .jets(Cuts::Et > 20*GeV && Cuts::abseta < 3, cmpMomByEt);
 
       Jets jets_isolated;
       for (size_t i = 0; i < jets_in.size(); ++i) {
@@ -204,8 +205,8 @@
       _h_4j_mu45->fill(_safeMass(FourMomentum(p4+p5))/sqrts, weight);
       _h_4j_mu46->fill(_safeMass(FourMomentum(p4+p6))/sqrts, weight);
       _h_4j_mu56->fill(_safeMass(FourMomentum(p5+p6))/sqrts, weight);