[Rivet-svn] rivet: More cuts in an attempt to validate but still stuff is mi...

[Rivet Mercurial]

details:   https://rivet.hepforge.org/hg/rivet/rev/0def9462a8ae
branches:  release-2-5-x
changeset: 5403:0def9462a8ae
user:      Holger Schulz <holger.schulz at durham.ac.uk>
date:      Fri Aug 12 11:15:05 2016 +0100
description:
More cuts in an attempt to validate but still stuff is missing --- needs proper c.m. cuts maybe

diffs (truncated from 77 to 50 lines):

--- a/ChangeLog	Thu Aug 11 15:52:30 2016 +0100
+++ b/ChangeLog	Fri Aug 12 11:15:05 2016 +0100
@@ -1,3 +1,12 @@
+2016-08-12  Holger Schulz  <holger.schulz at cern.ch>
+
+        * Implemented a few more cuts in prompt photon analysis
+	(CDF_1993_S2742446) but to no avail, the rise of the data towards
+	larger costheta values cannot be reproduced --- maybe this is a
+	candidate for more scrutiny and using the boosting machinery such that
+	the c.m. cuts can be done in a non-approximate way
+
+        
 2016-08-11  Holger Schulz  <holger.schulz at cern.ch>
         
 	* Rename CDF_2009_S8383952 to CDF_2009_I856131 due to invalid Spires
--- a/data/anainfo/CDF_1993_S2742446.info	Thu Aug 11 15:52:30 2016 +0100
+++ b/data/anainfo/CDF_1993_S2742446.info	Fri Aug 12 11:15:05 2016 +0100
@@ -23,7 +23,10 @@
   Data taken with the Collider Detector at Fermilab (CDF) during the 1988-1989
   run of the Tevatron are used to measure the distribution of the center-of-mass
   (rest frame of the initial state partons) angle between isolated prompt
-  photons and the beam direction.
+  photons and the beam direction. 
+  WARNING --- Implemented is the simplified c.m. definition
+  given in the paper. The rise of the data and the MC predictions presented in the 
+  latter could not be reproduced.
 BibKey: Abe:1993cv
 BibTeX: '@Article{Abe:1993cv,
      author    = "Abe, F. and others",
--- a/data/plotinfo/CDF_1993_S2742446.plot	Thu Aug 11 15:52:30 2016 +0100
+++ b/data/plotinfo/CDF_1993_S2742446.plot	Fri Aug 12 11:15:05 2016 +0100
@@ -2,5 +2,6 @@
 Title=Angular distribution of prompt photons
 XLabel=$|\cos\theta^*|=|\tanh((\eta_j-\eta_\gamma)/2.0)|$
 YLabel=$\mathrm{d}N/\mathrm{d}\cos\theta^*$
+YMax=7
 LogY=0
 # END PLOT
--- a/src/Analyses/CDF_1993_S2742446.cc	Thu Aug 11 15:52:30 2016 +0100
+++ b/src/Analyses/CDF_1993_S2742446.cc	Fri Aug 12 11:15:05 2016 +0100
@@ -58,13 +58,34 @@
       }
       if (Etsum > 2*GeV) vetoEvent;
 
+      FourMomentum jetsum;
+      Jets jets = apply<FastJets>(event, "Jets").jets(Cuts::pT > 10*GeV, cmpMomByPt);
+  
+      // Require at least one jet with pT>10 GeV
+      if (jets.size()==0) vetoEvent;
+