|
[Rivet-svn] r3400 - trunk/src/Analysesblackhole at projects.hepforge.org blackhole at projects.hepforge.orgThu 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 |