Page MenuHomeHEPForge
Files F10275334

No OneTemporary
Actions

View Options

diff --git a/include/Rivet/Projections/ZFinder.hh b/include/Rivet/Projections/ZFinder.hh
--- a/include/Rivet/Projections/ZFinder.hh
+++ b/include/Rivet/Projections/ZFinder.hh
@@ -1,138 +1,124 @@
// -*- C++ -*-
#ifndef RIVET_ZFinder_HH
#define RIVET_ZFinder_HH
#include "Rivet/Projections/FinalState.hh"
#include "Rivet/Projections/DressedLeptons.hh"
#include "Rivet/Projections/VetoedFinalState.hh"
namespace Rivet {
/// @brief Convenience finder of leptonically decaying Zs
///
/// Chain together different projections as convenience for finding Z's
/// from two leptons in the final state, including photon clustering.
///
/// @todo Alias then rename as Dileptons
class ZFinder : public ParticleFinder {
public:
enum ChargedLeptons { PROMPTCHLEPTONS=0, ALLCHLEPTONS };
enum ClusterPhotons { NOCLUSTER=0, CLUSTERNODECAY=1, CLUSTERALL };
enum PhotonTracking { NOTRACK=0, TRACK=1 };
/// @name Constructors
//@{
/// Constructor taking cuts object
/// @param inputfs Input final state
/// @param cuts lepton cuts
/// @param pid type of the leptons
/// @param minmass,maxmass mass window
/// @param dRmax maximum dR of photons around leptons to take into account
/// for Z reconstruction (only relevant if one of the following are true)
/// @param clusterPhotons whether such photons are supposed to be
/// clustered to the lepton objects and thus Z mom
/// @param trackPhotons whether such photons should be added to _theParticles
/// (cf. _trackPhotons)
ZFinder(const FinalState& inputfs,
const Cut& cuts,
PdgId pid,
double minmass, double maxmass,
double dRmax=0.1,
ChargedLeptons chLeptons=PROMPTCHLEPTONS,
ClusterPhotons clusterPhotons=CLUSTERNODECAY,
PhotonTracking trackPhotons=NOTRACK,
double masstarget=91.2*GeV);
/// Backward-compatible constructor with implicit chLeptons mode = PROMPTCHLEPTONS
/// @deprecated Remove this and always use the constructor with chLeptons argument.
ZFinder(const FinalState& inputfs,
const Cut& cuts,
PdgId pid,
double minmass, double maxmass,
double dRmax,
ClusterPhotons clusterPhotons,
PhotonTracking trackPhotons=NOTRACK,
double masstarget=91.2*GeV)
: ZFinder(inputfs, cuts, pid, minmass, maxmass,
dRmax, PROMPTCHLEPTONS, clusterPhotons, trackPhotons, masstarget)
{ }
/// Clone on the heap.
DEFAULT_RIVET_PROJ_CLONE(ZFinder);
//@}
/// Access to the found bosons
///
/// @note Currently either 0 or 1 boson can be found.
const Particles& bosons() const { return particles(); }
/// Access to the found boson (assuming it exists).
const Particle& boson() const { return bosons().front(); }
/// Access to the Z constituent clustered leptons
///
/// For example, to make more fine-grained cuts on the clustered leptons.
/// The positive charge constituent is first in the list (if not empty), and
/// the negative one second.
- Particles constituentLeptons() const;
- Particles constituents() const { return constituentLeptons(); }
-
-
- /// Access to the Z constituent clustered leptons, sorted by a comparison functor
- ///
- /// Unlike the no-arg version, this returns by value (i.e. is less efficient)
- Particles constituentLeptons(const ParticleSorter& cmp) const;
- Particles constituents(const ParticleSorter& cmp) const { return constituentLeptons(cmp); }
-
-
- // /// Access to all DressedLeptons in the fiducial region.
- // ///
- // /// This includes those DressedLeptons that could not
- // /// be paired up with any other DressedLepton to form a Z candidate.
- // const vector<DressedLepton>& allLeptons() const { return _allLeptons; }
+ const Particles & constituentLeptons() const;
+ const Particles & constituents() const { return constituentLeptons(); }
/// Access to the particles other than the Z leptons and clustered photons
///
/// Useful for e.g. input to a jet finder
const VetoedFinalState& remainingFinalState() const;
protected:
/// Apply the projection on the supplied event.
void project(const Event& e);
/// Compare projections.
int compare(const Projection& p) const;
public:
/// Clear the projection
void clear() { _theParticles.clear(); }
private:
/// Mass cuts to apply to clustered leptons (cf. InvMassFinalState)
double _minmass, _maxmass, _masstarget;
/// Switch for tracking of photons (whether to include them in the Z particle)
/// This is relevant when the clustered photons need to be excluded from e.g. a jet finder
PhotonTracking _trackPhotons;
/// Lepton flavour
PdgId _pid;
};
}
#endif
diff --git a/src/Projections/ZFinder.cc b/src/Projections/ZFinder.cc
--- a/src/Projections/ZFinder.cc
+++ b/src/Projections/ZFinder.cc
@@ -1,127 +1,118 @@
// -*- C++ -*-
#include "Rivet/Projections/ZFinder.hh"
#include "Rivet/Projections/PromptFinalState.hh"
#include "Rivet/Projections/InvMassFinalState.hh"
#include "Rivet/Projections/VetoedFinalState.hh"
namespace Rivet {
ZFinder::ZFinder(const FinalState& inputfs,
const Cut & fsCut,
PdgId pid,
double minmass, double maxmass,
double dRmax,
ChargedLeptons chLeptons,
ClusterPhotons clusterPhotons,
PhotonTracking trackPhotons,
double masstarget)
{
setName("ZFinder");
_minmass = minmass;
_maxmass = maxmass;
_masstarget = masstarget;
_pid = abs(pid);
_trackPhotons = trackPhotons;
// Identify bare leptons for dressing
// Bit of a code nightmare -- FS projection copy constructors don't work?
/// @todo Fix FS copy constructors!!
if (chLeptons == PROMPTCHLEPTONS) {
PromptFinalState inputfs_prompt(inputfs);
IdentifiedFinalState bareleptons = IdentifiedFinalState(inputfs_prompt);
bareleptons.acceptIdPair(_pid);
declare(bareleptons, "BareLeptons");
} else {
IdentifiedFinalState bareleptons = IdentifiedFinalState(inputfs);
bareleptons.acceptIdPair(_pid);
declare(bareleptons, "BareLeptons");
}
// Dress the bare leptons
const bool doClustering = (clusterPhotons != NOCLUSTER);
const bool useDecayPhotons = (clusterPhotons == CLUSTERALL);
DressedLeptons leptons(inputfs, get<FinalState>("BareLeptons"), (doClustering ? dRmax : -1.0), fsCut, useDecayPhotons);
addProjection(leptons, "DressedLeptons");
// Identify the non-Z part of the event
VetoedFinalState remainingFS;
remainingFS.addVetoOnThisFinalState(*this);
addProjection(remainingFS, "RFS");
}
/////////////////////////////////////////////////////
- Particles ZFinder::constituentLeptons() const {
- /// @note Ugly -- if there's no boson, the user should already have
- /// stopped. Being forced to return by value because of an edge case which
- /// is anyway wrong = yuck
- if (empty()) return Particles();
+ const Particles & ZFinder::constituentLeptons() const {
+ static const Particles none;
+ if (empty()) return none;
return boson().constituents();
- // return boson().constituents(isChargedLepton);
}
-
- Particles ZFinder::constituentLeptons(const ParticleSorter& cmp) const {
- return sortBy(constituentLeptons(), cmp);
- }
-
-
const VetoedFinalState& ZFinder::remainingFinalState() const {
return getProjection<VetoedFinalState>("RFS");
}
int ZFinder::compare(const Projection& p) const {
PCmp LCcmp = mkNamedPCmp(p, "DressedLeptons");
if (LCcmp != EQUIVALENT) return LCcmp;
const ZFinder& other = dynamic_cast<const ZFinder&>(p);
return (cmp(_minmass, other._minmass) ||
cmp(_maxmass, other._maxmass) ||
cmp(_pid, other._pid) ||
cmp(_trackPhotons, other._trackPhotons));
}
void ZFinder::project(const Event& e) {
clear();
// Get leptons and find an acceptable invariant mass OSSF pair
const DressedLeptons& leptons = applyProjection<DressedLeptons>(e, "DressedLeptons");
InvMassFinalState imfs({_pid, -_pid}, _minmass, _maxmass, _masstarget);
imfs.calc(leptons.particles());
if (imfs.particlePairs().empty()) {
MSG_TRACE("No acceptable inv-mass lepton/antilepton pairs found");
return;
}
// Assemble a pseudo-Z particle
const ParticlePair& Zconstituents = imfs.particlePairs().front();
const Particle& p1(Zconstituents.first), p2(Zconstituents.second);
const FourMomentum pZ = p1.momentum() + p2.momentum();
assert(p1.charge3() + p2.charge3() == 0);
Particle z(PID::Z0BOSON, pZ);
MSG_DEBUG(z << " reconstructed from: " << p1 << " + " << p2);
// Add (dressed) lepton constituents to the Z (skipping photons if requested)
// Keep the DressedLeptons found by the ZFinder
const Particle& l1 = p1.charge() > 0 ? p1 : p2;
const Particle& l2 = p2.charge() < 0 ? p2 : p1;
MSG_TRACE("l1 = " << l1.constituents());
MSG_TRACE("l2 = " << l2.constituents());
z.addConstituent(_trackPhotons == TRACK ? l1 : l1.constituents().front());
z.addConstituent(_trackPhotons == TRACK ? l2 : l2.constituents().front());
MSG_DEBUG("Number of stored raw Z constituents = " << z.rawConstituents().size() << " " << z.rawConstituents());
// Register the completed Z
_theParticles.push_back(z);
}
}

File Metadata

Mime Type
text/x-diff
Expires
Fri, Apr 4, 9:26 PM (1 m, 50 s)
Storage Engine
blob
Storage Format
Raw Data
Storage Handle
4737372
Default Alt Text
(9 KB)

Event Timeline