1 |
+ |
|
2 |
|
//-------------------------------------------------------------------------------------------------- |
3 |
|
// $Id $ |
4 |
|
// |
13 |
|
#include "MitAna/DataTree/interface/PFJetFwd.h" |
14 |
|
#include "MitAna/DataTree/interface/VertexFwd.h" |
15 |
|
#include "MitAna/DataTree/interface/TrackFwd.h" |
16 |
+ |
#include "MitAna/DataTree/interface/MuonCol.h" |
17 |
+ |
#include "MitAna/DataTree/interface/PFTauCol.h" |
18 |
+ |
#include "MitAna/DataTree/interface/ElectronCol.h" |
19 |
+ |
#include "MitAna/DataTree/interface/PhotonCol.h" |
20 |
|
#include "MitAna/DataTree/interface/Met.h" |
21 |
|
#include "MitAna/DataTree/interface/PFMet.h" |
22 |
|
#include "MitAna/DataTree/interface/PFMetCol.h" |
27 |
|
#include "MitCommon/MathTools/interface/MathUtils.h" |
28 |
|
|
29 |
|
#include "MitPhysics/Utils/interface/RecoilTools.h" |
30 |
+ |
#include "MitPhysics/Utils/interface/MetLeptonTools.h" |
31 |
|
#include "CondFormats/EgammaObjects/interface/GBRForest.h" |
32 |
|
|
33 |
|
class TRandom3; |
57 |
|
TString iPhiWeights ="$CMSSW_BASE/src/MitPhysics/data/gbrmetphi.root", |
58 |
|
TString iCovU1Weights ="$CMSSW_BASE/src/MitPhysics/data/gbrcovu1_52.root", |
59 |
|
TString iCovU2Weights ="$CMSSW_BASE/src/MitPhysics/data/gbrcovu2_52.root", |
60 |
< |
MVAMet::MVAType iType=kBaseline); |
60 |
> |
JetIDMVA::MVAType iType=JetIDMVA::kBaseline); |
61 |
|
|
62 |
|
Bool_t IsInitialized() const { return fIsInitialized; } |
63 |
+ |
Float_t* getVals(); |
64 |
|
Double_t evaluatePhi(); |
65 |
|
Double_t evaluateU1(); |
66 |
|
Double_t evaluateCovU1(); |
91 |
|
Float_t iNAllJet, |
92 |
|
Float_t iNPV ); |
93 |
|
|
94 |
< |
Met GetMet( Bool_t iPhi,Float_t iPtVis,Float_t iPhiVis,Float_t iSumEtVis, |
94 |
> |
Met GetMet( Bool_t iPhi, |
95 |
> |
Float_t iPtVis,Float_t iPhiVis,Float_t iSumEtVis, |
96 |
> |
Float_t iPtQ ,Float_t iPhiQ ,Float_t iSumEtQ, |
97 |
|
const PFMet *iMet , |
98 |
|
const PFCandidateCol *iCands, |
99 |
|
const Vertex *iVertex,const VertexCol *iVertices, |
103 |
|
int iNPV, |
104 |
|
Bool_t printDebug=false); |
105 |
|
|
106 |
< |
Met GetMet( Bool_t iPhi,Float_t iPtVis,Float_t iPhiVis,Float_t iSumEtVis, |
106 |
> |
Met GetMet( Bool_t iPhi, |
107 |
> |
Float_t iPtVis,Float_t iPhiVis,Float_t iSumEtVis, |
108 |
> |
Float_t iPtQ ,Float_t iPhiQ ,Float_t iSumEtQ, |
109 |
|
const PFMet *iMet , |
110 |
|
const PFCandidateCol *iCands, |
111 |
|
const Vertex *iVertex,const VertexCol *iVertices,Double_t iRho, |
114 |
|
Bool_t printDebug=false); |
115 |
|
|
116 |
|
Met GetMet( Bool_t iPhi, |
117 |
< |
Float_t iPt1,Float_t iPhi1,Float_t iEta1, |
118 |
< |
Float_t iPt2,Float_t iPhi2,Float_t iEta2, |
117 |
> |
Float_t iPt1,Float_t iPhi1,Float_t iEta1,Float_t iChargedFrac1, |
118 |
> |
Float_t iPt2,Float_t iPhi2,Float_t iEta2,Float_t iChargedFrac2, |
119 |
|
const PFMet *iMet , |
120 |
|
const PFCandidateCol *iCands, |
121 |
|
const Vertex *iVertex,const VertexCol *iVertices, |
126 |
|
Bool_t printDebug=false); |
127 |
|
|
128 |
|
Met GetMet( Bool_t iPhi, |
129 |
< |
Float_t iPt1,Float_t iPhi1,Float_t iEta1, |
130 |
< |
Float_t iPt2,Float_t iPhi2,Float_t iEta2, |
129 |
> |
Float_t iPt1,Float_t iPhi1,Float_t iEta1,Float_t iChargedFrac1, |
130 |
> |
Float_t iPt2,Float_t iPhi2,Float_t iEta2,Float_t iChargedFrac2, |
131 |
|
const PFMet *iMet , |
132 |
|
const PFCandidateCol *iCands, |
133 |
|
const Vertex *iVertex,const VertexCol *iVertices,Double_t iRho, |
135 |
|
int iNPV, |
136 |
|
Bool_t printDebug=false); |
137 |
|
|
138 |
+ |
Met GetMet( const MuonCol *iMuons,const ElectronCol *iElectrons,const PFTauCol *iTaus, |
139 |
+ |
const PFCandidateCol *iCands,const PFJetCol *iJets,const Vertex *iPV,const VertexCol *iVertices,const PFMetCol *iPFMet, |
140 |
+ |
FactorizedJetCorrector *iJetCorrector,const PileupEnergyDensityCol* iPUEnergyDensity); |
141 |
+ |
|
142 |
+ |
Met GetMet(const PhotonCol *iPhotons, |
143 |
+ |
const PFCandidateCol *iCands,const PFJetCol *iJets,const Vertex *iPV,const VertexCol *iVertices,const PFMetCol *iPFMet, |
144 |
+ |
FactorizedJetCorrector *iJetCorrector,const PileupEnergyDensityCol* iPUEnergyDensity); |
145 |
+ |
|
146 |
|
TMatrixD* GetMetCovariance() { return fCov; } |
147 |
|
Float_t GetSignificance () { return fSignificance;} |
148 |
|
Float_t GetUncertainty () { return fUncertainty;} |
149 |
|
RecoilTools *fRecoilTools; |
150 |
+ |
|
151 |
+ |
Int_t fNPhotons; |
152 |
|
|
153 |
|
protected: |
154 |
|
TString fPhiMethodName; |
156 |
|
TString fCovU1MethodName; |
157 |
|
TString fCovU2MethodName; |
158 |
|
Bool_t fIsInitialized; |
159 |
< |
MVAType fType; |
159 |
> |
JetIDMVA::MVAType fType; |
160 |
|
Bool_t f42; |
161 |
|
|
162 |
|
Float_t fSumEt ; |
199 |
|
Float_t fUncertainty; |
200 |
|
TMatrixD *fCov; |
201 |
|
|
202 |
+ |
Int_t fNMuons; |
203 |
+ |
Int_t fNElectrons; |
204 |
+ |
Int_t fNTaus; |
205 |
+ |
|
206 |
+ |
MetLeptonTools *fMetLeptonTools; |
207 |
|
//TMVA::Reader* fPhiReader; |
208 |
|
//TMVA::Reader* fU1Reader; |
209 |
|
ClassDef(MVAMet,0) |