[Rivet-svn] rivet: Enable Nsubjettiness build

Rivet Mercurial rivet at projects.hepforge.org
Thu May 24 11:45:02 BST 2018


details:   https://rivet.hepforge.org/hg/rivet/rev/b7bc0ddb9a6f
branches:  release-2-6-x
changeset: 6318:b7bc0ddb9a6f
user:      David Grellscheid <david.grellscheid at durham.ac.uk>
date:      Thu May 24 11:37:03 2018 +0100
description:
Enable Nsubjettiness build

diffs (truncated from 8141 to 50 lines):

--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/analyses/pluginATLAS/ATLAS_2017_I1614149.cc	Thu May 24 11:37:03 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 "Rivet/Tools/Nsubjettiness/Njettiness.hh"
+#include "Rivet/Tools/Nsubjettiness/Nsubjettiness.hh"
+#include "Rivet/Tools/Nsubjettiness/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");
+


More information about the Rivet-svn mailing list