[Rivet-svn] rivet: 4mom shouldnt be pointers cos they get reordered behind t...

[Rivet Mercurial]

details:   https://rivet.hepforge.org/hg/rivet/rev/ce40a6dcf214
branches:  release-2-6-x
changeset: 6246:ce40a6dcf214
user:      Jon Butterworth <j.butterworth at cern.ch>
date:      Wed Feb 28 23:10:34 2018 +0100
description:
4mom shouldnt be pointers cos they get reordered behind the scenes and mT is then always zero.

diffs (20 lines):

--- a/analyses/pluginCMS/CMS_2016_I1491953.cc	Tue Feb 27 22:34:33 2018 +0000
+++ b/analyses/pluginCMS/CMS_2016_I1491953.cc	Wed Feb 28 23:10:34 2018 +0100
@@ -107,9 +107,14 @@
       const WFinder& wfinder_mu = applyProjection<WFinder>(event, "WFinder_mu");
       if (wfinder_mu.bosons().size() != 1) vetoEvent;
 
-      const FourMomentum& lepton0 = wfinder_mu.constituentLeptons()[0].momentum();
-      const FourMomentum& neutrino = wfinder_mu.constituentNeutrinos()[0].momentum();
-      double WmT = sqrt( 2 * lepton0.pT() * neutrino.pT() * (1 - cos(deltaPhi(lepton0, neutrino))) );
+      //const FourMomentum& lepton0 = wfinder_mu.constituentLeptons()[0].momentum();
+      //const FourMomentum& neutrino = wfinder_mu.constituentNeutrinos()[0].momentum();
+      //double WmT = sqrt( 2 * lepton0.pT() * neutrino.pT() * (1 - cos(deltaPhi(lepton0, neutrino))) );
+
+      const FourMomentum lepton0 = wfinder_mu.constituentLepton().momentum();
+      const FourMomentum neutrino = wfinder_mu.constituentNeutrino().momentum();
+      double WmT = wfinder_mu.mT();
+
       if (WmT < 50.0*GeV) vetoEvent;
       if (lepton0.abseta() > 2.1 || lepton0.pT() < 25.0*GeV) vetoEvent;