Page MenuHomeHEPForge
Files F7879810

No OneTemporary
Actions

View Options

diff --git a/src/Core/Jet.cc b/src/Core/Jet.cc
--- a/src/Core/Jet.cc
+++ b/src/Core/Jet.cc
@@ -1,154 +1,125 @@
#include "Rivet/Jet.hh"
#include "Rivet/Tools/Logging.hh"
#include "Rivet/Tools/ParticleIdUtils.hh"
#include "Rivet/Tools/RivetBoost.hh"
#include "Rivet/ParticleName.hh"
namespace Rivet {
Jet& Jet::setState(const vector<Particle>& particles, const FourMomentum& pjet) {
setParticles(particles);
setMomentum(pjet);
return *this;
}
- // Jet& Jet::setState(const vector<FourMomentum>& momenta, const FourMomentum& pjet) {
- // setParticles(momenta);
- // setMomentum(pjet);
- // return *this;
- // }
-
-
Jet& Jet::setMomentum(const FourMomentum& momentum) {
_momentum = momentum;
return *this;
}
- Jet& Jet::setParticles(const vector<Particle>& particles) {
- _particles = particles;
- // foreach (const Particle& p, particles) {
- // _momenta.push_back(p.momentum());
- // }
- return *this;
- }
-
-
- // Jet& Jet::setParticles(const vector<FourMomentum>& momenta) {
- // _momenta = momenta;
- // return *this;
- // }
-
-
- // Jet& Jet::addParticle(const FourMomentum& particle) {
- // _momenta.push_back(particle);
- // return *this;
- // }
-
-
- // Jet& Jet::addParticle(const Particle& particle) {
- // _particles.push_back(particle);
- // _momenta.push_back(particle.momentum());
- // return *this;
- // }
+ Jet& Jet::setParticles(const vector<Particle>& particles) {
+ _particles = particles;
+ return *this;
+ }
bool Jet::containsParticle(const Particle& particle) const {
const int barcode = particle.genParticle()->barcode();
foreach (const Particle& p, particles()) {
if (p.genParticle()->barcode() == barcode) return true;
}
return false;
}
bool Jet::containsParticleId(PdgId pid) const {
foreach (const Particle& p, particles()) {
if (p.pdgId() == pid) return true;
}
return false;
}
bool Jet::containsParticleId(const vector<PdgId>& pids) const {
foreach (const Particle& p, particles()) {
foreach (PdgId pid, pids) {
if (p.pdgId() == pid) return true;
}
}
return false;
}
/// @todo Jet::containsMatch(Matcher m) { ... if m(pid) return true; ... }
double Jet::neutralEnergy() const {
double e_neutral = 0.0;
foreach (const Particle& p, particles()) {
const PdgId pid = p.pdgId();
if (PID::threeCharge(pid) == 0) {
- e_neutral += p.momentum().E();
+ e_neutral += p.E();
}
}
return e_neutral;
}
double Jet::hadronicEnergy() const {
double e_hadr = 0.0;
foreach (const Particle& p, particles()) {
const PdgId pid = p.pdgId();
if (PID::isHadron(pid)) {
- e_hadr += p.momentum().E();
+ e_hadr += p.E();
}
}
return e_hadr;
}
bool Jet::containsCharm() const {
foreach (const Particle& p, particles()) {
const PdgId pid = p.pdgId();
if (abs(pid) == PID::CQUARK) return true;
if (PID::isHadron(pid) && PID::hasCharm(pid)) return true;
HepMC::GenVertex* gv = p.genParticle()->production_vertex();
if (gv) {
foreach (const GenParticle* pi, Rivet::particles(gv, HepMC::ancestors)) {
const PdgId pid2 = pi->pdg_id();
if (PID::isHadron(pid2) && PID::hasCharm(pid2)) return true;
}
}
}
return false;
}
bool Jet::containsBottom() const {
foreach (const Particle& p, particles()) {
const PdgId pid = p.pdgId();
if (abs(pid) == PID::BQUARK) return true;
if (PID::isHadron(pid) && PID::hasBottom(pid)) return true;
HepMC::GenVertex* gv = p.genParticle()->production_vertex();
if (gv) {
foreach (const GenParticle* pi, Rivet::particles(gv, HepMC::ancestors)) {
const PdgId pid2 = pi->pdg_id();
if (PID::isHadron(pid2) && PID::hasBottom(pid2)) return true;
}
}
}
return false;
}
Jet& Jet::clear() {
//_momenta.clear();
_particles.clear();
_momentum = FourMomentum();
return *this;
}
}

File Metadata

Mime Type
text/x-diff
Expires
Tue, Nov 19, 9:07 PM (1 d, 45 m)
Storage Engine
blob
Storage Format
Raw Data
Storage Handle
3806202
Default Alt Text
(4 KB)

Event Timeline