[Rivet-svn] rivet: Cleanup, fix and mark D0_1996_S3324664 jet decorrelation ...

[Rivet Mercurial]

details:   https://rivet.hepforge.org/hg/rivet/rev/79dcaf6395c9
branches:  release-2-5-x
changeset: 5390:79dcaf6395c9
user:      Holger Schulz <holger.schulz at durham.ac.uk>
date:      Fri Aug 05 16:59:01 2016 +0100
description:
Cleanup, fix and mark D0_1996_S3324664 jet decorrelation analysis validated

diffs (truncated from 65 to 50 lines):

--- a/ChangeLog	Fri Aug 05 14:03:13 2016 +0100
+++ b/ChangeLog	Fri Aug 05 16:59:01 2016 +0100
@@ -10,6 +10,10 @@
 	algorithm, I tried with the implementation of fastjet
 	fastjet/D0RunIConePlugin.hh but that really does not help.
 
+	* D0_1996_S3324664 --- fix normalisations, sorting jets properly now,
+	cleanup and mark validated.
+
+
 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_S3324664.info	Fri Aug 05 14:03:13 2016 +0100
+++ b/data/anainfo/D0_1996_S3324664.info	Fri Aug 05 16:59:01 2016 +0100
@@ -4,7 +4,7 @@
 Experiment: D0
 Collider: Tevatron Run 1
 SpiresID: 3324664
-Status: UNVALIDATED
+Status: VALIDATED
 Authors:
  - Frank Siegert <frank.siegert at cern.ch>
 References:
@@ -21,6 +21,8 @@
   First measurement of the azimuthal decorrelation between jets with
   pseudorapidity separation up to five units. The data were accumulated using
   the D0 detector during Tevatron Run 1 at $\sqrt{s} = 1.8\;\TeV$.
+  Requires di-jet events with at least one jet above 50 GeV and two jets
+  above 20 GeV.
 BibKey: Abachi:1996et
 BibTeX: '@Article{Abachi:1996et,
      author    = "Abachi, S. and others",
--- a/src/Analyses/D0_1996_S3324664.cc	Fri Aug 05 14:03:13 2016 +0100
+++ b/src/Analyses/D0_1996_S3324664.cc	Fri Aug 05 16:59:01 2016 +0100
@@ -38,11 +38,7 @@
 
     void analyze(const Event& event) {
       const double weight = event.weight();
-
-      Jets jets;
-      foreach (const Jet& jet, apply<FastJets>(event, "ConeJets").jets(Cuts::pT > 20*GeV)) {
-        if (jet.abseta() < 3.0) jets.push_back(jet);
-      }
+      Jets jets = apply<FastJets>(event, "ConeJets").jets(Cuts::Et > 20*GeV && Cuts::abseta<3, cmpMomByEt);
 
       if (jets.size() < 2) vetoEvent;
 
@@ -75,13 +71,11 @@