rivet is hosted by Hepforge, IPPP Durham
Rivet 4.0.0
SmearedMET.hh
1// -*- C++ -*-
2#ifndef RIVET_SmearedMET_HH
3#define RIVET_SmearedMET_HH
4
5#include "Rivet/Projection.hh"
6#include "Rivet/Projections/METFinder.hh"
7#include "Rivet/Projections/MissingMomentum.hh"
8#include "Rivet/Tools/SmearingFunctions.hh"
9#include <functional>
10
11namespace Rivet {
12
13
15 class SmearedMET : public METFinder {
16 public:
17
20
24 template <typename V2VFN>
25 SmearedMET(const MissingMomentum& mm, const V2VFN& metSmearFn)
26 : _metSmearFn(metSmearFn)
27 {
28 setName("SmearedMET");
29 declare(mm, "TruthMET");
30 }
31
33 template <typename V2VFN>
34 SmearedMET(const V2VFN& metSmearFn, const Cut& cut)
35 : _metSmearFn(metSmearFn)
36 {
37 setName("SmearedMET");
38 declare(MissingMomentum(cut), "TruthMET");
39 }
40
41
44
46
48 using Projection::operator =;
49
50
52 CmpState compare(const Projection& p) const {
53 const SmearedMET& other = dynamic_cast<const SmearedMET&>(p);
54 if (get_address(_metSmearFn) == 0) return cmp((size_t)this, (size_t)&p);
55 MSG_TRACE("Smear hashes = " << get_address(_metSmearFn) << "," << get_address(other._metSmearFn));
56 return mkPCmp(other, "TruthMET") || cmp(get_address(_metSmearFn), get_address(other._metSmearFn));
57 }
58
59
61 void project(const Event& e) {
62 const auto& mm = apply<MissingMomentum>(e, "TruthMET");
63 _vet = mm.vectorEt();
64 if (_metSmearFn) _vet = _metSmearFn(_vet, mm.scalarEt()); //< smearing
65 }
66
67
72
77 const Vector3& vectorPt() const { return vectorEt(); }
78
80
81
86
90 const Vector3& vectorEt() const { return _vet; }
91
93
94
96 void reset() { }
97
98
99 protected:
100
101 Vector3 _vet;
102
104 std::function<Vector3(const Vector3&, double)> _metSmearFn;
105
106 };
107
108
109}
110
111#endif
Representation of a HepMC event, and enabler of Projection caching.
Definition Event.hh:22
Interface for projections that find missing transverse energy/momentum.
Definition METFinder.hh:11
Calculate missing , etc. as complements to the total visible momentum.
Definition MissingMomentum.hh:22
const PROJ & declare(const PROJ &proj, const std::string &name) const
Register a contained projection (user-facing version)
Definition ProjectionApplier.hh:175
Base class for all Rivet projections.
Definition Projection.hh:29
Cmp< Projection > mkPCmp(const Projection &otherparent, const std::string &pname) const
void setName(const std::string &name)
Used by derived classes to set their name.
Definition Projection.hh:148
Wrapper projection for smearing missing (transverse) energy/momentum with detector resolutions.
Definition SmearedMET.hh:15
SmearedMET(const MissingMomentum &mm, const V2VFN &metSmearFn)
Constructor from a MissingMomentum projection and a smearing function.
Definition SmearedMET.hh:25
void reset()
Reset the projection. Smearing functions will be unchanged.
Definition SmearedMET.hh:96
const Vector3 & vectorPt() const
Definition SmearedMET.hh:77
const Vector3 & vectorEt() const
Definition SmearedMET.hh:90
RIVET_DEFAULT_PROJ_CLONE(SmearedMET)
Clone on the heap.
void project(const Event &e)
Perform the MET finding & smearing calculation.
Definition SmearedMET.hh:61
SmearedMET(const V2VFN &metSmearFn, const Cut &cut)
Constructor from a Cut (on the particles used to determine missing momentum) and a smearing function.
Definition SmearedMET.hh:34
CmpState compare(const Projection &p) const
Compare to another SmearedMET.
Definition SmearedMET.hh:52
Three-dimensional specialisation of Vector.
Definition Vector3.hh:40
#define MSG_TRACE(x)
Lowest-level, most verbose messaging, using MSG_LVL.
Definition Logging.hh:180
Definition MC_CENT_PPB_Projections.hh:10
Cmp< T > cmp(const T &t1, const T &t2)
Global helper function for easy creation of Cmp objects.
Definition Cmp.hh:255
uintptr_t get_address(std::function< T(U...)> f)
Get a function pointer / hash integer from an std::function.
Definition RivetSTL.hh:216