[Rivet-svn] rivet: atlas top routine - requires nsubjettiness

[Rivet Mercurial]

details:   https://rivet.hepforge.org/hg/rivet/rev/16851ce753b0
branches:  release-2-6-x
changeset: 6313:16851ce753b0
user:      Jon Butterworth <j.butterworth at cern.ch>
date:      Thu May 24 10:14:54 2018 +0100
description:
atlas top routine - requires nsubjettiness

diffs (truncated from 805 to 50 lines):

--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/analyses/pluginATLAS/ATLAS_2017_I1614149.cc.waiting-for-nsubjettiness	Thu May 24 10:14:54 2018 +0100
@@ -0,0 +1,353 @@
+// -*- C++ -*
+#include "Rivet/Analysis.hh"
+#include "Rivet/Projections/FinalState.hh"
+#include "Rivet/Projections/VetoedFinalState.hh"
+#include "Rivet/Projections/IdentifiedFinalState.hh"
+#include "Rivet/Projections/PromptFinalState.hh"
+#include "Rivet/Projections/DressedLeptons.hh"
+#include "Rivet/Projections/FastJets.hh"
+#include "Rivet/Projections/MissingMomentum.hh"
+#include "fastjet/tools/Filter.hh"
+
+// substructure includes included in fjcontrib-1.021 (http://fastjet.hepforge.org/contrib/)
+#include "fastjet/contrib/Njettiness.hh"
+#include "fastjet/contrib/Nsubjettiness.hh"
+#include "fastjet/contrib/NjettinessPlugin.hh"
+
+namespace Rivet {
+  class ATLAS_2017_I1614149 : public Analysis {
+    public:
+
+      /// Constructor
+      ///@brief: Resolved and boosted ttbar l+jets cross sections at 13 TeV
+      DEFAULT_RIVET_ANALYSIS_CTOR(ATLAS_2017_I1614149);
+
+      void init() {
+        // Eta ranges
+	      Cut eta_full = (Cuts::abseta < 5.0);
+        Cut lep_cuts = (Cuts::abseta < 2.5) && (Cuts::pT > 25*GeV);
+
+        // All final state particles
+        FinalState fs(eta_full);
+
+	      IdentifiedFinalState all_photons(fs);
+        all_photons.acceptIdPair(PID::PHOTON);
+
+	      // Get photons to dress leptons
+        IdentifiedFinalState ph_id(fs);
+        ph_id.acceptIdPair(PID::PHOTON);
+
+        // Projection to find the electrons
+        IdentifiedFinalState el_id(fs);
+        el_id.acceptIdPair(PID::ELECTRON);
+
+	      PromptFinalState photons(ph_id);
+        photons.acceptTauDecays(true);
+	      addProjection(photons, "photons");
+