|
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
|
// -*- C++ -*-
// Package: METReco
// Class: MET
//
/**\class MET
Description: This is the fundemental class for missing transverse
momentum (a.k.a MET). The class inherits from RecoCandidate, and so
is stored in the event as a candidate. The actual MET information is
contained in the RecoCandidate LorentzVector while supplimentary
information is stored as varibles in the MET class itself (such as
SumET). As there may have been more than one correction applied to
the missing transverse momentum, a vector of corrections to the
missing px and the missing py is maintained so that one can always
recover the uncorrected MET by applying the negative of each
correction.
*/
// Original authors: Michael Schmitt, Richard Cavanaugh The University of Florida
// changes by: Freya Blekman, Cornell University
//
//____________________________________________________________________________||
#ifndef METRECO_MET_H
#define METRECO_MET_H
//____________________________________________________________________________||
#include "DataFormats/RecoCandidate/interface/RecoCandidate.h"
#include "DataFormats/METReco/interface/CorrMETData.h"
#include <cmath>
#include <vector>
#include <cstring>
#include <Math/SMatrix.h>
#include <Math/SVector.h>
//____________________________________________________________________________||
namespace reco {
typedef ROOT::Math::SMatrix<double, 2> METCovMatrix;
class MET : public RecoCandidate {
public:
MET();
MET(const LorentzVector& p4_, const Point& vtx_, bool isWeighted = false);
MET(double sumet_, const LorentzVector& p4_, const Point& vtx_, bool isWeighted = false);
MET(double sumet_,
const std::vector<CorrMETData>& corr_,
const LorentzVector& p4_,
const Point& vtx_,
bool isWeighted = false);
MET* clone() const override;
//________________________________________________________________________||
//scalar sum of transverse energy over all objects
double sumEt() const { return sumet; }
//MET Significance = MET / std::sqrt(SumET)
double mEtSig() const { return (sumet ? (this->et() / std::sqrt(sumet)) : (0.0)); }
//real MET significance
double significance() const;
//longitudinal component of the vector sum of energy over all object
//(useful for data quality monitoring)
double e_longitudinal() const { return elongit; }
//________________________________________________________________________||
//Define different methods for the corrections to individual MET elements
std::vector<double> dmEx() const;
std::vector<double> dmEy() const;
std::vector<double> dsumEt() const;
std::vector<CorrMETData> mEtCorr() const { return corr; }
//________________________________________________________________________||
void setSignificanceMatrix(const reco::METCovMatrix& matrix);
reco::METCovMatrix getSignificanceMatrix(void) const;
/// Set boolean if weights were applied by algorithm (e.g. PUPPI weights)
void setIsWeighted(bool isWeighted) { mIsWeighted = isWeighted; }
/// boolean if weights were applied by algorithm (e.g. PUPPI weights)
int isWeighted() const { return mIsWeighted; }
private:
bool overlap(const Candidate&) const override;
double sumet;
double elongit;
// bookkeeping for the significance
double signif_dxx;
double signif_dyy;
double signif_dyx;
double signif_dxy;
std::vector<CorrMETData> corr;
bool mIsWeighted;
};
} // namespace reco
//____________________________________________________________________________||
#endif // METRECO_MET_H
|