|
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
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
|
// This is the header file "FactorizedJetCorrectorCalculator.h". This is the interface for the
// class FactorizedJetCorrectorCalculator.
// Author: Konstantinos Kousouris, Philipp Schieferdecker
// Email: kkousour@fnal.gov, philipp.schieferdecker@cern.ch
#ifndef CondFormats_JetMETObjects_FactorizedJetCorrectorCalculator_h
#define CondFormats_JetMETObjects_FactorizedJetCorrectorCalculator_h
#include <vector>
#include <string>
#include "TLorentzVector.h"
class SimpleJetCorrector;
class JetCorrectorParameters;
class FactorizedJetCorrectorCalculator {
public:
class VariableValues {
public:
friend class FactorizedJetCorrector;
friend class FactorizedJetCorrectorCalculator;
VariableValues();
void setNPV(int fNPV);
void setJetEta(float fEta);
void setJetPt(float fPt);
void setJetE(float fE);
void setJetPhi(float fE);
void setJetEMF(float fEMF);
void setJetA(float fA);
void setRho(float fRho);
void setJPTrawP4(const TLorentzVector& fJPTrawP4);
void setJPTrawOff(float fJPTrawOff);
void setLepPx(float fLepPx);
void setLepPy(float fLepPy);
void setLepPz(float fLepPz);
void setAddLepToJet(bool fAddLepToJet);
void reset();
private:
//---- Member Data ---------
int mNPV;
float mJetE;
float mJetEta;
float mJetPt;
float mJetPhi;
float mJetEMF;
float mJetA;
float mRho;
float mJPTrawE;
float mJPTrawEt;
float mJPTrawPt;
float mJPTrawEta;
float mJPTrawOff;
float mLepPx;
float mLepPy;
float mLepPz;
bool mAddLepToJet;
bool mIsNPVset;
bool mIsJetEset;
bool mIsJetPtset;
bool mIsJetPhiset;
bool mIsJetEtaset;
bool mIsJetEMFset;
bool mIsJetAset;
bool mIsRhoset;
bool mIsJPTrawP4set;
bool mIsJPTrawOFFset;
bool mIsLepPxset;
bool mIsLepPyset;
bool mIsLepPzset;
bool mIsAddLepToJetset;
};
enum VarTypes {
kJetPt,
kJetEta,
kJetPhi,
kJetE,
kJetEMF,
kRelLepPt,
kPtRel,
kNPV,
kJetA,
kRho,
kJPTrawE,
kJPTrawEt,
kJPTrawPt,
kJPTrawEta,
kJPTrawOff
};
enum LevelTypes { kL1, kL2, kL3, kL4, kL5, kL6, kL7, kL1fj, kL1JPT, kL2L3Res };
FactorizedJetCorrectorCalculator();
FactorizedJetCorrectorCalculator(const std::string& fLevels,
const std::string& fTags,
const std::string& fOptions = "");
FactorizedJetCorrectorCalculator(const std::vector<JetCorrectorParameters>& fParameters);
FactorizedJetCorrectorCalculator(const FactorizedJetCorrectorCalculator&) = delete;
FactorizedJetCorrectorCalculator& operator=(const FactorizedJetCorrectorCalculator&) = delete;
~FactorizedJetCorrectorCalculator();
float getCorrection(VariableValues&) const;
std::vector<float> getSubCorrections(VariableValues&) const;
private:
//---- Member Functions ----
float getLepPt(const VariableValues&) const;
float getRelLepPt(const VariableValues&) const;
float getPtRel(const VariableValues&) const;
std::string parseOption(const std::string& ss, const std::string& type) const;
std::string removeSpaces(const std::string& ss) const;
std::vector<std::string> parseLevels(const std::string& ss) const;
void initCorrectors(const std::string& fLevels, const std::string& fFiles, const std::string& fOptions);
void checkConsistency(const std::vector<std::string>& fLevels, const std::vector<std::string>& fTags);
std::vector<float> fillVector(const std::vector<VarTypes>& fVarTypes, const VariableValues&) const;
std::vector<VarTypes> mapping(const std::vector<std::string>& fNames) const;
//---- Member Data ---------
std::vector<LevelTypes> mLevels;
std::vector<std::vector<VarTypes> > mParTypes, mBinTypes;
std::vector<SimpleJetCorrector const*> mCorrectors;
};
#endif
|