[Rivet-svn] rivet: 4 new changesets

[Rivet Mercurial]

details:   https://rivet.hepforge.org/hg/rivet/rev/e032e92e6149
branches:  
changeset: 6095:e032e92e6149
user:      = Jon Butterworth <j.butterworth at cern.ch>
date:      Thu Oct 26 10:45:53 2017 +0200
description:
restore analysis description - but in the right place

details:   https://rivet.hepforge.org/hg/rivet/rev/e03f0e8fcce6
branches:  
changeset: 6096:e03f0e8fcce6
user:      = Jon Butterworth <j.butterworth at cern.ch>
date:      Thu Oct 26 11:29:18 2017 +0200
description:
remove specific requirements on neutrino flavour and 4-mom

details:   https://rivet.hepforge.org/hg/rivet/rev/ade365f19a31
branches:  
changeset: 6097:ade365f19a31
user:      = Jon Butterworth <j.butterworth at cern.ch>
date:      Thu Oct 26 11:33:14 2017 +0200
description:
use prompt leptons only

details:   https://rivet.hepforge.org/hg/rivet/rev/72b5b960cb6e
branches:  
changeset: 6098:72b5b960cb6e
user:      = Jon Butterworth <j.butterworth at cern.ch>
date:      Thu Oct 26 11:54:30 2017 +0200
description:
use femtobarns

diffs (truncated from 210 to 50 lines):

--- a/analyses/pluginATLAS/ATLAS_2012_I1203852.cc	Thu Oct 26 10:08:36 2017 +0200
+++ b/analyses/pluginATLAS/ATLAS_2012_I1203852.cc	Thu Oct 26 11:54:30 2017 +0200
@@ -12,6 +12,8 @@
 #include "Rivet/Projections/MissingMomentum.hh"
 #include "Rivet/Projections/InvMassFinalState.hh"
 
+#define ZMASS 91.1876 // GeV
+
 namespace Rivet {
 
 
@@ -26,7 +28,7 @@
 
 
   /// 4l to ZZ assignment -- algorithm
-  void identifyZstates(Zstate& Z1, Zstate& Z2, const Particles& leptons_sel4l, const double ZMASS) {
+  void identifyZstates(Zstate& Z1, Zstate& Z2, const Particles& leptons_sel4l) {
 
     /////////////////////////////////////////////////////////////////////////////
     /// ZZ->4l pairing
@@ -94,16 +96,14 @@
     {    }
 
     void init() {
+
+      // NB Missing ET is not required to be neutrinos
       FinalState fs(-5.0, 5.0, 0.0*GeV);
 
       // Final states to form Z bosons
       vids.push_back(make_pair(PID::ELECTRON, PID::POSITRON));
       vids.push_back(make_pair(PID::MUON, PID::ANTIMUON));
 
-      vnuids.push_back(make_pair(PID::NU_E, PID::NU_EBAR) );
-      vnuids.push_back(make_pair(PID::NU_MU, PID::NU_MUBAR) );
-      vnuids.push_back(make_pair(PID::NU_TAU, PID::NU_TAUBAR) );
-
       IdentifiedFinalState Photon(fs);
       Photon.acceptIdPair(PID::PHOTON);
 
@@ -133,11 +133,13 @@
 
 
       /// Get all neutrinos. These will not be used to form jets.
-      /// We'll use the highest 2 pT neutrinos to calculate the MET
       IdentifiedFinalState neutrino_fs(Cuts::abseta < 4.5);
       neutrino_fs.acceptNeutrinos();
       declare(neutrino_fs, "NEUTRINO_FS");
 
+      // Calculate missing ET from the visible final state, not by requiring neutrinos
+      addProjection(MissingMomentum(Cuts::abseta < 4.5), "MISSING");