[Rivet-svn] r3400 - trunk/src/Analyses

blackhole at projects.hepforge.org blackhole at projects.hepforge.org
Thu Sep 29 13:28:19 BST 2011


Author: hoeth
Date: Thu Sep 29 13:28:18 2011
New Revision: 3400

Log:
use units

Modified:
   trunk/src/Analyses/MC_DIJET.cc

Modified: trunk/src/Analyses/MC_DIJET.cc
==============================================================================
--- trunk/src/Analyses/MC_DIJET.cc	Thu Sep 29 08:12:39 2011	(r3399)
+++ trunk/src/Analyses/MC_DIJET.cc	Thu Sep 29 13:28:18 2011	(r3400)
@@ -47,7 +47,7 @@
 
     void analyze(const Event& event) {
       const FastJets& fastjets = applyProjection<FastJets>(event, "Jets");
-      const Jets jets = fastjets.jetsByPt(20.);
+      const Jets jets = fastjets.jetsByPt(20.*GeV);
       const double weight = event.weight();
 
       if (jets.size() < 2 || jets.size() >= 3) vetoEvent;
@@ -59,7 +59,7 @@
       _hist_jetdphi->fill(angle , weight);
       _hist_jetdeta->fill(prapidity, weight);
 
-      foreach(Jet j, fastjets.jetsByPt(20*GeV)) {
+      foreach(Jet j, fastjets.jetsByPt(20.*GeV)) {
         _hist_jetpt->fill(j.momentum().pT(), weight);
         _hist_jetptlog->fill(log(j.momentum().pT()), weight);
         _hist_jetphi->fill(j.momentum().azimuthalAngle(), weight);
@@ -70,7 +70,7 @@
       // const FastJets& cfastjets = applyProjection<FastJets>(event, "ChargedJets");
       double meanpt(0), rmspt(0);
       /// @todo Add jets
-      // foreach(Jet cj, cfastjets.jetsByPt(20.)){
+      // foreach(Jet cj, cfastjets.jetsByPt(20.*GeV)){
       _hist_chargemultiplicity->fill(cfs.particles().size(), weight);
       foreach(Particle cp, cfs.particles()) {
         meanpt= meanpt + cp.momentum().pT();


More information about the Rivet-svn mailing list