22 |
|
#include "MitCommon/MathTools/interface/MathUtils.h" |
23 |
|
#include "MitPhysics/Utils/interface/MuonTools.h" |
24 |
|
|
25 |
+ |
// for Rho definitons |
26 |
+ |
#include "MitPhysics/Utils/interface/RhoUtilities.h" |
27 |
+ |
|
28 |
|
class TRandom3; |
29 |
|
namespace TMVA { |
30 |
|
class Reader; |
43 |
|
kV8, |
44 |
|
kIDIsoCombinedDetIso, |
45 |
|
kIsoRingsV0, |
46 |
< |
kIDIsoCombinedIsoRingsV0 |
46 |
> |
kIDIsoCombinedIsoRingsV0, |
47 |
> |
kIDV0, |
48 |
> |
kIsoDeltaR |
49 |
|
}; |
50 |
|
|
51 |
|
|
52 |
|
void Initialize( std::string methodName, |
53 |
|
std::string weightsfile, |
54 |
< |
MuonIDMVA::MVAType type); |
54 |
> |
MuonIDMVA::MVAType type, |
55 |
> |
RhoUtilities::RhoType theRhoType = RhoUtilities::DEFAULT); |
56 |
|
void Initialize( std::string methodName, |
57 |
|
MuonIDMVA::MVAType type, |
58 |
|
Bool_t useBinnedVersion, |
59 |
< |
std::vector<std::string> weightsfiles ); |
59 |
> |
std::vector<std::string> weightsfiles, |
60 |
> |
RhoUtilities::RhoType theRhoType = RhoUtilities::DEFAULT); |
61 |
|
void Initialize(TString methodName, |
62 |
|
TString Subdet0Pt10To14p5Weights , |
63 |
|
TString Subdet1Pt10To14p5Weights , |
65 |
|
TString Subdet1Pt14p5To20Weights, |
66 |
|
TString Subdet0Pt20ToInfWeights, |
67 |
|
TString Subdet1Pt20ToInfWeights, |
68 |
< |
MuonIDMVA::MVAType type); |
68 |
> |
MuonIDMVA::MVAType type, |
69 |
> |
RhoUtilities::RhoType theRhoType = RhoUtilities::DEFAULT); |
70 |
|
|
71 |
|
Bool_t IsInitialized() const { return fIsInitialized; } |
72 |
|
UInt_t GetMVABin(double eta,double pt, |
73 |
< |
Bool_t isGlobal, Bool_t isTrackerMuon ) const; |
73 |
> |
Bool_t isGlobal = kTRUE, Bool_t isTrackerMuon = kTRUE ) const; |
74 |
|
Double_t MVAValue(const Muon *mu, const Vertex *vertex, MuonTools *fMuonTools, |
75 |
|
const PFCandidateCol *PFCands, |
76 |
|
const PileupEnergyDensityCol *PileupEnergyDensity, |
134 |
|
Double_t MuHadIso05OverPt, |
135 |
|
Bool_t printDebug = kFALSE |
136 |
|
); |
137 |
< |
|
137 |
> |
Double_t MVAValue_IsoRings( Double_t MuPt, |
138 |
> |
Double_t MuEta, |
139 |
> |
Double_t ChargedIso_DR0p0To0p1, |
140 |
> |
Double_t ChargedIso_DR0p1To0p2, |
141 |
> |
Double_t ChargedIso_DR0p2To0p3, |
142 |
> |
Double_t ChargedIso_DR0p3To0p4, |
143 |
> |
Double_t ChargedIso_DR0p4To0p5, |
144 |
> |
Double_t GammaIso_DR0p0To0p1, |
145 |
> |
Double_t GammaIso_DR0p1To0p2, |
146 |
> |
Double_t GammaIso_DR0p2To0p3, |
147 |
> |
Double_t GammaIso_DR0p3To0p4, |
148 |
> |
Double_t GammaIso_DR0p4To0p5, |
149 |
> |
Double_t NeutralHadronIso_DR0p0To0p1, |
150 |
> |
Double_t NeutralHadronIso_DR0p1To0p2, |
151 |
> |
Double_t NeutralHadronIso_DR0p2To0p3, |
152 |
> |
Double_t NeutralHadronIso_DR0p3To0p4, |
153 |
> |
Double_t NeutralHadronIso_DR0p4To0p5, |
154 |
> |
Bool_t printDebug = kFALSE); |
155 |
> |
Double_t MVAValue_ID( Double_t MuPt, |
156 |
> |
Double_t MuEta, |
157 |
> |
Bool_t MuIsGlobal, |
158 |
> |
Bool_t MuIsTracker, |
159 |
> |
Double_t MuTkNchi2, |
160 |
> |
Double_t MuGlobalNchi2, |
161 |
> |
Double_t MuNValidHits, |
162 |
> |
Double_t MuNTrackerHits, |
163 |
> |
Double_t MuNPixelHits, |
164 |
> |
Double_t MuNMatches, |
165 |
> |
Double_t MuTrkKink, |
166 |
> |
Double_t MuSegmentCompatibility, |
167 |
> |
Double_t MuCaloCompatibility, |
168 |
> |
Double_t MuHadEnergy, |
169 |
> |
Double_t MuEmEnergy, |
170 |
> |
Double_t MuHadS9Energy, |
171 |
> |
Double_t MuEmS9Energy, |
172 |
> |
Bool_t printDebug = kFALSE); |
173 |
|
|
174 |
|
protected: |
175 |
|
std::vector<TMVA::Reader*> fTMVAReader; |
178 |
|
MVAType fMVAType; |
179 |
|
Bool_t fUseBinnedVersion; |
180 |
|
UInt_t fNMVABins; |
181 |
< |
|
181 |
> |
RhoUtilities::RhoType fTheRhoType; |
182 |
> |
|
183 |
|
Float_t fMVAVar_MuPt; |
184 |
|
Float_t fMVAVar_MuEta; |
185 |
|
Float_t fMVAVar_MuTkNchi2; |
234 |
|
Float_t fMVAVar_NeutralHadronIso_DR0p4To0p5; |
235 |
|
Float_t fMVAVar_NeutralHadronIso_DR0p5To0p7; |
236 |
|
|
237 |
+ |
// isolation variables II |
238 |
+ |
Float_t fMVAVar_MuRelIsoPFCharged; |
239 |
+ |
Float_t fMVAVar_MuRelIsoPFNeutral; |
240 |
+ |
Float_t fMVAVar_MuRelIsoPFPhotons; |
241 |
+ |
Float_t fMVAVar_MuDeltaRMean; |
242 |
+ |
Float_t fMVAVar_MuDeltaRSum; |
243 |
+ |
Float_t fMVAVar_MuDensity; |
244 |
+ |
|
245 |
+ |
|
246 |
|
ClassDef(MuonIDMVA, 0) // Muon MVA |
247 |
|
}; |
248 |
|
} |