[Rivet-svn] r1717 - in trunk: data/anainfo src/Analyses

[blackhole at projects.hepforge.org]

Author: fsiegert
Date: Thu Jul 23 23:14:51 2009
New Revision: 1717

Log:
D0_2009_S8320160:
 - add cut on y_boost
 - use FS, not CFS for jets
 - go back to normalize() for the histograms, because that's what the
paper actually does.
 - set status to UNVALIDATED, because I have to find out, whether the
cut on y_boost is done after selecting the two hardest jets (i.e. veto
the event) or before (i.e. selecting two different hardest jets). The
paper seems to suggest the former, but I can't imagine that's what an
experiment can do.

Modified:
   trunk/data/anainfo/D0_2009_S8320160.info
   trunk/src/Analyses/D0_2009_S8320160.cc

Modified: trunk/data/anainfo/D0_2009_S8320160.info
==============================================================================
--- trunk/data/anainfo/D0_2009_S8320160.info	Thu Jul 23 21:55:35 2009	(r1716)
+++ trunk/data/anainfo/D0_2009_S8320160.info	Thu Jul 23 23:14:51 2009	(r1717)
@@ -4,7 +4,7 @@
 Experiment: D0
 Collider: Tevatron Run 2
 SpiresID: 8320160
-Status: VALIDATED
+Status: UNVALIDATED
 Authors:
  - Frank Siegert <frank.siegert at durham.ac.uk>
 References:

Modified: trunk/src/Analyses/D0_2009_S8320160.cc
==============================================================================
--- trunk/src/Analyses/D0_2009_S8320160.cc	Thu Jul 23 21:55:35 2009	(r1716)
+++ trunk/src/Analyses/D0_2009_S8320160.cc	Thu Jul 23 23:14:51 2009	(r1717)
@@ -14,8 +14,8 @@
   {
     setBeams(PROTON, ANTIPROTON);
     
-    ChargedFinalState cfs;
-    FastJets conefinder(cfs, FastJets::D0ILCONE, 0.7);
+    FinalState fs;
+    FastJets conefinder(fs, FastJets::D0ILCONE, 0.7);
     addProjection(conefinder, "ConeFinder");
   } 
 
@@ -40,17 +40,19 @@
   void D0_2009_S8320160::analyze(const Event & e) {
     double weight = e.weight();
 
-    const JetAlg& jetpro = applyProjection<JetAlg>(e, "ConeFinder");
-    const Jets& jets = jetpro.jetsByPt();
+    const Jets& jets = applyProjection<JetAlg>(e, "ConeFinder").jetsByPt();
     
-    if(jets.size()<2) {
-      vetoEvent;
-    }
+    if(jets.size()<2) vetoEvent;
     
     FourMomentum j0(jets[0].momentum());
     FourMomentum j1(jets[1].momentum());
+    double y0 = j0.rapidity();
+    double y1 = j1.rapidity();
+    
+    if (fabs(y0+y1)>2) vetoEvent;
+    
     double mjj = FourMomentum(j0+j1).mass();
-    double chi = exp(fabs(j0.rapidity()-j1.rapidity()));
+    double chi = exp(fabs(y0-y1));
     _h_chi_dijet.fill(mjj, chi, weight);
   }
 
@@ -59,7 +61,7 @@
   // Finalize
   void D0_2009_S8320160::finalize() {
     foreach (AIDA::IHistogram1D* hist, _h_chi_dijet.getHistograms()) {
-      scale(hist, 1.0/sumOfWeights());
+      normalize(hist);
     }
   }