[Rivet-svn] rivet: add 7TeV CMS photon routine from DG

[Rivet Mercurial]

details:   https://rivet.hepforge.org/hg/rivet/rev/f1b8ed04d8d4
branches:  
changeset: 6094:f1b8ed04d8d4
user:      = Jon Butterworth <j.butterworth at cern.ch>
date:      Thu Oct 26 10:08:36 2017 +0200
description:
add 7TeV CMS photon routine from DG

diffs (truncated from 387 to 50 lines):

--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/analyses/pluginCMS/CMS_2014_I1266056.cc	Thu Oct 26 10:08:36 2017 +0200
@@ -0,0 +1,142 @@
+// -*- C++ -*-
+#include "Rivet/Analysis.hh"
+#include "Rivet/Projections/FinalState.hh"
+#include "Rivet/Projections/LeadingParticlesFinalState.hh"
+#include "Rivet/Projections/VetoedFinalState.hh"
+#include "Rivet/Projections/FastJets.hh"
+
+namespace Rivet {
+
+
+  /// @brief Measurement of gamma + jets + X triple differential cross-sections
+  ///
+  /// @author David Grellscheid
+  class CMS_2014_I1266056 : public Analysis {
+  public:
+
+    // Constructor
+    CMS_2014_I1266056()
+      : Analysis("CMS_2014_I1266056")
+    {    }
+
+
+    // Book histograms and initialise projections before the run
+    void init() {
+      // Final state
+      FinalState fs(-3, 3);
+      addProjection(fs, "FS");
+
+      // Leading photon
+      LeadingParticlesFinalState photonfs(FinalState(-2.5, 2.5, 40.0*GeV));
+      photonfs.addParticleId(PID::PHOTON);
+      addProjection(photonfs, "LeadingPhoton");
+
+      // FS excluding the leading photon
+      VetoedFinalState vfs(fs);
+      vfs.addVetoOnThisFinalState(photonfs);
+      addProjection(vfs, "JetFS");
+
+      // Jets
+      FastJets jetpro(vfs, FastJets::ANTIKT, 0.5);
+      //jetpro.useInvisibles();
+      addProjection(jetpro, "Jets");
+
+      _h_phverycentral_jetcentral = bookHisto1D(1, 1, 1);
+      _h_phcentral_jetcentral     = bookHisto1D(2, 1, 1);
+      _h_phforward_jetcentral     = bookHisto1D(3, 1, 1);
+      _h_phveryforward_jetcentral = bookHisto1D(4, 1, 1);