CMS 3D CMS Logo
EGRegressionModifierV1.cc
Go to the documentation of this file.
1 #include "CondFormats/DataRecord/interface/GBRWrapperRcd.h"
2 #include "CondFormats/GBRForest/interface/GBRForest.h"
3 #include "DataFormats/EcalDetId/interface/EBDetId.h"
4 #include "DataFormats/EcalDetId/interface/EEDetId.h"
5 #include "RecoEgamma/EgammaTools/plugins/EGRegressionModifierHelpers.h"
6 #include "CommonTools/CandAlgos/interface/ModifyObjectValueBase.h"
7 #include "DataFormats/Common/interface/ValueMap.h"
8 #include "DataFormats/EgammaCandidates/interface/GsfElectron.h"
9 #include "DataFormats/EgammaCandidates/interface/Photon.h"
10 #include "FWCore/Utilities/interface/EDGetToken.h"
11 #include "FWCore/Utilities/interface/InputTag.h"
12 #include "RecoEgamma/EgammaTools/interface/EcalClusterLocal.h"
13 #include "Geometry/Records/interface/CaloGeometryRecord.h"
14 #include "Geometry/CaloGeometry/interface/CaloGeometry.h"
15 
16 #include <vdt/vdtMath.h>
17 
18 class EGRegressionModifierV1 : public ModifyObjectValueBase {
19 public:
20  EGRegressionModifierV1(const edm::ParameterSet& conf, edm::ConsumesCollector& cc);
21 
22  void setEvent(const edm::Event&) final;
23  void setEventContent(const edm::EventSetup&) final;
24 
25  void modifyObject(reco::GsfElectron&) const final;
26  void modifyObject(reco::Photon&) const final;
27 
28  // just calls reco versions
29  void modifyObject(pat::Electron& ele) const final { modifyObject(static_cast<reco::GsfElectron&>(ele)); }
30  void modifyObject(pat::Photon& pho) const final { modifyObject(static_cast<reco::Photon&>(pho)); }
31 
32 private:
33  EGRegressionModifierCondTokens eleCond50nsTokens_;
34  EGRegressionModifierCondTokens phoCond50nsTokens_;
35  EGRegressionModifierCondTokens eleCond25nsTokens_;
36  EGRegressionModifierCondTokens phoCond25nsTokens_;
37 
38  edm::ESGetToken<GBRForest, GBRWrapperRcd> condNamesWeight50nsToken_;
39  edm::ESGetToken<GBRForest, GBRWrapperRcd> condNamesWeight25nsToken_;
40 
41  const bool autoDetectBunchSpacing_;
42  int bunchspacing_;
43  edm::EDGetTokenT<unsigned int> bunchSpacingToken_;
44  float rhoValue_;
45  edm::EDGetTokenT<double> rhoToken_;
46  int nVtx_;
47  edm::EDGetTokenT<reco::VertexCollection> vtxToken_;
48  edm::ESGetToken<CaloGeometry, CaloGeometryRecord> caloGeomToken_;
49  CaloGeometry const* caloGeom_;
50  const bool applyExtraHighEnergyProtection_;
51 
52  std::vector<const GBRForestD*> phoForestsMean_;
53  std::vector<const GBRForestD*> phoForestsSigma_;
54  std::vector<const GBRForestD*> eleForestsMean_;
55  std::vector<const GBRForestD*> eleForestsSigma_;
56  const GBRForest* epForest_;
57 };
58 
59 DEFINE_EDM_PLUGIN(ModifyObjectValueFactory, EGRegressionModifierV1, "EGRegressionModifierV1");
60 
61 EGRegressionModifierV1::EGRegressionModifierV1(const edm::ParameterSet& conf, edm::ConsumesCollector& cc)
62  : ModifyObjectValueBase(conf),
63  eleCond50nsTokens_{conf.getParameterSet("electron_config"), "regressionKey_50ns", "uncertaintyKey_50ns", cc},
64  phoCond50nsTokens_{conf.getParameterSet("photon_config"), "regressionKey_50ns", "uncertaintyKey_50ns", cc},
65  eleCond25nsTokens_{conf.getParameterSet("electron_config"), "regressionKey_25ns", "uncertaintyKey_25ns", cc},
66  phoCond25nsTokens_{conf.getParameterSet("photon_config"), "regressionKey_25ns", "uncertaintyKey_25ns", cc},
67  autoDetectBunchSpacing_(conf.getParameter<bool>("autoDetectBunchSpacing")),
68  bunchspacing_(autoDetectBunchSpacing_ ? 450 : conf.getParameter<int>("manualBunchSpacing")),
69  rhoToken_(cc.consumes(conf.getParameter<edm::InputTag>("rhoCollection"))),
70  vtxToken_(cc.consumes(conf.getParameter<edm::InputTag>("vertexCollection"))),
71  caloGeomToken_{cc.esConsumes()},
72  applyExtraHighEnergyProtection_(conf.getParameter<bool>("applyExtraHighEnergyProtection")) {
73  if (autoDetectBunchSpacing_)
74  bunchSpacingToken_ = cc.consumes(conf.getParameter<edm::InputTag>("bunchSpacingTag"));
75 
76  auto const& electrons = conf.getParameterSet("electron_config");
77  condNamesWeight50nsToken_ = cc.esConsumes(electrons.getParameter<edm::ESInputTag>("combinationKey_50ns"));
78  condNamesWeight25nsToken_ = cc.esConsumes(electrons.getParameter<edm::ESInputTag>("combinationKey_25ns"));
79 }
80 
81 void EGRegressionModifierV1::setEvent(const edm::Event& evt) {
82  if (autoDetectBunchSpacing_) {
83  bunchspacing_ = evt.get(bunchSpacingToken_);
84  }
85  rhoValue_ = evt.get(rhoToken_);
86  nVtx_ = evt.get(vtxToken_).size();
87 }
88 
89 void EGRegressionModifierV1::setEventContent(const edm::EventSetup& evs) {
90  caloGeom_ = &evs.getData(caloGeomToken_);
91 
92  phoForestsMean_ = retrieveGBRForests(evs, (bunchspacing_ == 25) ? phoCond25nsTokens_.mean : phoCond50nsTokens_.mean);
93  phoForestsSigma_ =
94  retrieveGBRForests(evs, (bunchspacing_ == 25) ? phoCond25nsTokens_.sigma : phoCond50nsTokens_.sigma);
95 
96  eleForestsMean_ = retrieveGBRForests(evs, (bunchspacing_ == 25) ? eleCond25nsTokens_.mean : eleCond50nsTokens_.mean);
97  eleForestsSigma_ =
98  retrieveGBRForests(evs, (bunchspacing_ == 25) ? eleCond25nsTokens_.sigma : eleCond50nsTokens_.sigma);
99 
100  epForest_ = &evs.getData((bunchspacing_ == 25) ? condNamesWeight25nsToken_ : condNamesWeight50nsToken_);
101 }
102 
103 void EGRegressionModifierV1::modifyObject(reco::GsfElectron& ele) const {
104  // regression calculation needs no additional valuemaps
105 
106  const reco::SuperClusterRef& superClus = ele.superCluster();
107  const edm::Ptr<reco::CaloCluster>& theseed = superClus->seed();
108  const int numberOfClusters = superClus->clusters().size();
109  const bool missing_clusters = !superClus->clusters()[numberOfClusters - 1].isAvailable();
110 
111  if (missing_clusters)
112  return; // do not apply corrections in case of missing info (slimmed MiniAOD electrons)
113 
114  std::array<float, 33> eval;
115  const double rawEnergy = superClus->rawEnergy();
116  const auto& ess = ele.showerShape();
117 
118  // SET INPUTS
119  eval[0] = nVtx_;
120  eval[1] = rawEnergy;
121  eval[2] = superClus->eta();
122  eval[3] = superClus->phi();
123  eval[4] = superClus->etaWidth();
124  eval[5] = superClus->phiWidth();
125  eval[6] = ess.r9;
126  eval[7] = theseed->energy() / rawEnergy;
127  eval[8] = ess.eMax / rawEnergy;
128  eval[9] = ess.e2nd / rawEnergy;
129  eval[10] = (ess.eLeft + ess.eRight != 0.f ? (ess.eLeft - ess.eRight) / (ess.eLeft + ess.eRight) : 0.f);
130  eval[11] = (ess.eTop + ess.eBottom != 0.f ? (ess.eTop - ess.eBottom) / (ess.eTop + ess.eBottom) : 0.f);
131  eval[12] = ess.sigmaIetaIeta;
132  eval[13] = ess.sigmaIetaIphi;
133  eval[14] = ess.sigmaIphiIphi;
134  eval[15] = std::max(0, numberOfClusters - 1);
135 
136  // calculate sub-cluster variables
137  std::vector<float> clusterRawEnergy;
138  clusterRawEnergy.resize(std::max(3, numberOfClusters), 0);
139  std::vector<float> clusterDEtaToSeed;
140  clusterDEtaToSeed.resize(std::max(3, numberOfClusters), 0);
141  std::vector<float> clusterDPhiToSeed;
142  clusterDPhiToSeed.resize(std::max(3, numberOfClusters), 0);
143  float clusterMaxDR = 999.;
144  float clusterMaxDRDPhi = 999.;
145  float clusterMaxDRDEta = 999.;
146  float clusterMaxDRRawEnergy = 0.;
147 
148  size_t iclus = 0;
149  float maxDR = 0;
150  // loop over all clusters that aren't the seed
151  for (auto const& pclus : superClus->clusters()) {
152  if (theseed == pclus)
153  continue;
154  clusterRawEnergy[iclus] = pclus->energy();
155  clusterDPhiToSeed[iclus] = reco::deltaPhi(pclus->phi(), theseed->phi());
156  clusterDEtaToSeed[iclus] = pclus->eta() - theseed->eta();
157 
158  // find cluster with max dR
159  const auto the_dr = reco::deltaR(*pclus, *theseed);
160  if (the_dr > maxDR) {
161  maxDR = the_dr;
162  clusterMaxDR = maxDR;
163  clusterMaxDRDPhi = clusterDPhiToSeed[iclus];
164  clusterMaxDRDEta = clusterDEtaToSeed[iclus];
165  clusterMaxDRRawEnergy = clusterRawEnergy[iclus];
166  }
167  ++iclus;
168  }
169 
170  eval[16] = clusterMaxDR;
171  eval[17] = clusterMaxDRDPhi;
172  eval[18] = clusterMaxDRDEta;
173  eval[19] = clusterMaxDRRawEnergy / rawEnergy;
174  eval[20] = clusterRawEnergy[0] / rawEnergy;
175  eval[21] = clusterRawEnergy[1] / rawEnergy;
176  eval[22] = clusterRawEnergy[2] / rawEnergy;
177  eval[23] = clusterDPhiToSeed[0];
178  eval[24] = clusterDPhiToSeed[1];
179  eval[25] = clusterDPhiToSeed[2];
180  eval[26] = clusterDEtaToSeed[0];
181  eval[27] = clusterDEtaToSeed[1];
182  eval[28] = clusterDEtaToSeed[2];
183 
184  // calculate coordinate variables
185  const bool isEB = ele.isEB();
186  float dummy;
187  int iPhi;
188  int iEta;
189  float cryPhi;
190  float cryEta;
191  if (ele.isEB())
192  egammaTools::localEcalClusterCoordsEB(*theseed, *caloGeom_, cryEta, cryPhi, iEta, iPhi, dummy, dummy);
193  else
194  egammaTools::localEcalClusterCoordsEE(*theseed, *caloGeom_, cryEta, cryPhi, iEta, iPhi, dummy, dummy);
195 
196  if (isEB) {
197  eval[29] = cryEta;
198  eval[30] = cryPhi;
199  eval[31] = iEta;
200  eval[32] = iPhi;
201  } else {
202  eval[29] = superClus->preshowerEnergy() / superClus->rawEnergy();
203  }
204 
205  //magic numbers for MINUIT-like transformation of BDT output onto limited range
206  //(These should be stored inside the conditions object in the future as well)
207  constexpr double meanlimlow = 0.2;
208  constexpr double meanlimhigh = 2.0;
209  constexpr double meanoffset = meanlimlow + 0.5 * (meanlimhigh - meanlimlow);
210  constexpr double meanscale = 0.5 * (meanlimhigh - meanlimlow);
211 
212  constexpr double sigmalimlow = 0.0002;
213  constexpr double sigmalimhigh = 0.5;
214  constexpr double sigmaoffset = sigmalimlow + 0.5 * (sigmalimhigh - sigmalimlow);
215  constexpr double sigmascale = 0.5 * (sigmalimhigh - sigmalimlow);
216 
217  const int coridx = isEB ? 0 : 1;
218 
219  //these are the actual BDT responses
220  double rawmean = eleForestsMean_[coridx]->GetResponse(eval.data());
221  double rawsigma = eleForestsSigma_[coridx]->GetResponse(eval.data());
222 
223  //apply transformation to limited output range (matching the training)
224  double mean = meanoffset + meanscale * vdt::fast_sin(rawmean);
225  double sigma = sigmaoffset + sigmascale * vdt::fast_sin(rawsigma);
226 
227  //regression target is ln(Etrue/Eraw)
228  //so corrected energy is ecor=exp(mean)*e, uncertainty is exp(mean)*eraw*sigma=ecor*sigma
229  double ecor = mean * (eval[1]);
230  if (!isEB)
231  ecor = mean * (eval[1] + superClus->preshowerEnergy());
232  const double sigmacor = sigma * ecor;
233 
234  ele.setCorrectedEcalEnergy(ecor);
235  ele.setCorrectedEcalEnergyError(sigmacor);
236 
237  // E-p combination
238  std::array<float, 11> eval_ep;
239 
240  const float ep = ele.trackMomentumAtVtx().R();
241  const float tot_energy = superClus->rawEnergy() + superClus->preshowerEnergy();
242  const float momentumError = ele.trackMomentumError();
243  const float trkMomentumRelError = ele.trackMomentumError() / ep;
244  const float eOverP = tot_energy * mean / ep;
245  eval_ep[0] = tot_energy * mean;
246  eval_ep[1] = sigma / mean;
247  eval_ep[2] = ep;
248  eval_ep[3] = trkMomentumRelError;
249  eval_ep[4] = sigma / mean / trkMomentumRelError;
250  eval_ep[5] = tot_energy * mean / ep;
251  eval_ep[6] = tot_energy * mean / ep * sqrt(sigma / mean * sigma / mean + trkMomentumRelError * trkMomentumRelError);
252  eval_ep[7] = ele.ecalDriven();
253  eval_ep[8] = ele.trackerDrivenSeed();
254  eval_ep[9] = int(ele.classification()); //eleClass;
255  eval_ep[10] = isEB;
256 
257  // CODE FOR FUTURE SEMI_PARAMETRIC
258  //double rawweight = ep_forestsMean_[coridx]->GetResponse(eval_ep.data());
260  //double weight = meanoffset + meanscale*vdt::fast_sin(rawweight);
262 
263  // CODE FOR STANDARD BDT
264  double weight = 0.;
265  if (eOverP > 0.025 && std::abs(ep - ecor) < 15. * std::sqrt(momentumError * momentumError + sigmacor * sigmacor) &&
266  (!applyExtraHighEnergyProtection_ || ((momentumError < 10. * ep) || (ecor < 200.)))) {
267  // protection against crazy track measurement
268  weight = std::clamp(epForest_->GetResponse(eval_ep.data()), 0., 1.);
269  }
270 
271  double combinedMomentum = weight * ele.trackMomentumAtVtx().R() + (1. - weight) * ecor;
272  double combinedMomentumError = sqrt(weight * weight * ele.trackMomentumError() * ele.trackMomentumError() +
273  (1. - weight) * (1. - weight) * sigmacor * sigmacor);
274 
275  math::XYZTLorentzVector oldMomentum = ele.p4();
276  math::XYZTLorentzVector newMomentum(oldMomentum.x() * combinedMomentum / oldMomentum.t(),
277  oldMomentum.y() * combinedMomentum / oldMomentum.t(),
278  oldMomentum.z() * combinedMomentum / oldMomentum.t(),
279  combinedMomentum);
280 
281  ele.correctMomentum(newMomentum, ele.trackMomentumError(), combinedMomentumError);
282 }
283 
284 void EGRegressionModifierV1::modifyObject(reco::Photon& pho) const {
285  // regression calculation needs no additional valuemaps
286 
287  std::array<float, 35> eval;
288  const reco::SuperClusterRef& superClus = pho.superCluster();
289  const edm::Ptr<reco::CaloCluster>& theseed = superClus->seed();
290 
291  const int numberOfClusters = superClus->clusters().size();
292  const bool missing_clusters = !superClus->clusters()[numberOfClusters - 1].isAvailable();
293 
294  if (missing_clusters)
295  return; // do not apply corrections in case of missing info (slimmed MiniAOD electrons)
296 
297  const double rawEnergy = superClus->rawEnergy();
298  const auto& ess = pho.showerShapeVariables();
299 
300  // SET INPUTS
301  eval[0] = rawEnergy;
302  eval[1] = pho.r9();
303  eval[2] = superClus->etaWidth();
304  eval[3] = superClus->phiWidth();
305  eval[4] = std::max(0, numberOfClusters - 1);
306  eval[5] = pho.hadronicOverEm();
307  eval[6] = rhoValue_;
308  eval[7] = nVtx_;
309  eval[8] = theseed->eta() - superClus->position().Eta();
310  eval[9] = reco::deltaPhi(theseed->phi(), superClus->position().Phi());
311  eval[10] = theseed->energy() / rawEnergy;
312  eval[11] = ess.e3x3 / ess.e5x5;
313  eval[12] = ess.sigmaIetaIeta;
314  eval[13] = ess.sigmaIphiIphi;
315  eval[14] = ess.sigmaIetaIphi / (ess.sigmaIphiIphi * ess.sigmaIetaIeta);
316  eval[15] = ess.maxEnergyXtal / ess.e5x5;
317  eval[16] = ess.e2nd / ess.e5x5;
318  eval[17] = ess.eTop / ess.e5x5;
319  eval[18] = ess.eBottom / ess.e5x5;
320  eval[19] = ess.eLeft / ess.e5x5;
321  eval[20] = ess.eRight / ess.e5x5;
322  eval[21] = ess.e2x5Max / ess.e5x5;
323  eval[22] = ess.e2x5Left / ess.e5x5;
324  eval[23] = ess.e2x5Right / ess.e5x5;
325  eval[24] = ess.e2x5Top / ess.e5x5;
326  eval[25] = ess.e2x5Bottom / ess.e5x5;
327 
328  const bool isEB = pho.isEB();
329  if (isEB) {
330  EBDetId ebseedid(theseed->seed());
331  eval[26] = pho.e5x5() / theseed->energy();
332  int ieta = ebseedid.ieta();
333  int iphi = ebseedid.iphi();
334  eval[27] = ieta;
335  eval[28] = iphi;
336  int signieta = ieta > 0 ? +1 : -1;
337  eval[29] = (ieta - signieta) % 5;
338  eval[30] = (iphi - 1) % 2;
339  eval[31] = (abs(ieta) <= 25) * ((ieta - signieta)) + (abs(ieta) > 25) * ((ieta - 26 * signieta) % 20);
340  eval[32] = (iphi - 1) % 20;
341  eval[33] = ieta;
342  eval[34] = iphi;
343  } else {
344  EEDetId eeseedid(theseed->seed());
345  eval[26] = superClus->preshowerEnergy() / rawEnergy;
346  eval[27] = superClus->preshowerEnergyPlane1() / rawEnergy;
347  eval[28] = superClus->preshowerEnergyPlane2() / rawEnergy;
348  eval[29] = eeseedid.ix();
349  eval[30] = eeseedid.iy();
350  }
351 
352  //magic numbers for MINUIT-like transformation of BDT output onto limited range
353  //(These should be stored inside the conditions object in the future as well)
354  const double meanlimlow = 0.2;
355  const double meanlimhigh = 2.0;
356  const double meanoffset = meanlimlow + 0.5 * (meanlimhigh - meanlimlow);
357  const double meanscale = 0.5 * (meanlimhigh - meanlimlow);
358 
359  const double sigmalimlow = 0.0002;
360  const double sigmalimhigh = 0.5;
361  const double sigmaoffset = sigmalimlow + 0.5 * (sigmalimhigh - sigmalimlow);
362  const double sigmascale = 0.5 * (sigmalimhigh - sigmalimlow);
363 
364  const int coridx = isEB ? 0 : 1;
365 
366  //these are the actual BDT responses
367  const double rawmean = phoForestsMean_[coridx]->GetResponse(eval.data());
368  const double rawsigma = phoForestsSigma_[coridx]->GetResponse(eval.data());
369  //apply transformation to limited output range (matching the training)
370  const double mean = meanoffset + meanscale * vdt::fast_sin(rawmean);
371  const double sigma = sigmaoffset + sigmascale * vdt::fast_sin(rawsigma);
372 
373  //regression target is ln(Etrue/Eraw)
374  //so corrected energy is ecor=exp(mean)*e, uncertainty is exp(mean)*eraw*sigma=ecor*sigma
375  const double ecor = isEB ? mean * eval[0] : mean * (eval[0] + superClus->preshowerEnergy());
376 
377  const double sigmacor = sigma * ecor;
378  pho.setCorrectedEnergy(reco::Photon::P4type::regression2, ecor, sigmacor, true);
379 }
reco::CaloCluster::phi
double phi() const
azimuthal angle of cluster centroid
Definition: CaloCluster.h:184
reco::Photon::hadronicOverEm
float hadronicOverEm(int depth=0) const
Definition: Photon.h:233
EgHLTOffHistBins_cfi.eOverP
eOverP
Definition: EgHLTOffHistBins_cfi.py:37
EGRegressionModifierV1::autoDetectBunchSpacing_
const bool autoDetectBunchSpacing_
Definition: EGRegressionModifierV1.cc:41
EGRegressionModifierV1
Definition: EGRegressionModifierV1.cc:18
EGRegressionModifierV1::phoCond50nsTokens_
EGRegressionModifierCondTokens phoCond50nsTokens_
Definition: EGRegressionModifierV1.cc:34
EGRegressionModifierV1::bunchSpacingToken_
edm::EDGetTokenT< unsigned int > bunchSpacingToken_
Definition: EGRegressionModifierV1.cc:43
edm::ESInputTag
Definition: ESInputTag.h:87
EGRegressionModifierHelpers.h
reco::GsfElectron::isEB
bool isEB() const
Definition: GsfElectron.h:336
reco::Photon::showerShapeVariables
const ShowerShape & showerShapeVariables() const
Definition: Photon.h:209
reco::Photon::superCluster
reco::SuperClusterRef superCluster() const override
Ref to SuperCluster.
SiStripPI::mean
Definition: SiStripPayloadInspectorHelper.h:169
EGRegressionModifierV1::modifyObject
void modifyObject(reco::GsfElectron &) const final
Definition: EGRegressionModifierV1.cc:103
EGRegressionModifierV1::caloGeomToken_
edm::ESGetToken< CaloGeometry, CaloGeometryRecord > caloGeomToken_
Definition: EGRegressionModifierV1.cc:48
reco::GsfElectron::correctMomentum
void correctMomentum(const LorentzVector &p4, float trackMomentumError, float p4Error)
Definition: GsfElectron.h:892
edm::Ref::isAvailable
bool isAvailable() const
Definition: Ref.h:537
reco::deltaPhi
constexpr double deltaPhi(double phi1, double phi2)
Definition: deltaPhi.h:26
edm::EDGetTokenT< unsigned int >
EBDetId
Definition: EBDetId.h:17
GBRWrapperRcd.h
mps_merge.weight
weight
Definition: mps_merge.py:88
GBRForest
Definition: GBRForest.h:24
pat::Photon
Analysis-level Photon class.
Definition: Photon.h:46
EBDetId.h
EEDetId.h
egammaTools::localEcalClusterCoordsEE
void localEcalClusterCoordsEE(const reco::CaloCluster &bclus, const CaloGeometry &geom, float &xcry, float &ycry, int &ix, int &iy, float &thetatilt, float &phitilt)
Definition: EcalClusterLocal.cc:87
EGRegressionModifierV1::setEvent
void setEvent(const edm::Event &) final
Definition: EGRegressionModifierV1.cc:81
EGRegressionModifierV1::vtxToken_
edm::EDGetTokenT< reco::VertexCollection > vtxToken_
Definition: EGRegressionModifierV1.cc:47
GBRForest.h
reco::GsfElectron::showerShape
const ShowerShape & showerShape() const
Definition: GsfElectron.h:475
edm::Ref< SuperClusterCollection >
EGRegressionModifierV1::nVtx_
int nVtx_
Definition: EGRegressionModifierV1.cc:46
LEDCalibrationChannels.iphi
iphi
Definition: LEDCalibrationChannels.py:64
retrieveGBRForests
std::vector< const GBRForestD * > retrieveGBRForests(edm::EventSetup const &evs, std::vector< edm::ESGetToken< GBRForestD, GBRDWrapperRcd >> const &tokens)
Definition: EGRegressionModifierHelpers.cc:15
GBRForest::GetResponse
double GetResponse(const float *vector) const
Definition: GBRForest.h:48
Photon.h
CaloGeometry
Definition: CaloGeometry.h:21
EGRegressionModifierCondTokens::mean
std::vector< edm::ESGetToken< GBRForestD, GBRDWrapperRcd > > mean
Definition: EGRegressionModifierHelpers.h:18
ModifyObjectValueBase.h
EGRegressionModifierV1::modifyObject
void modifyObject(pat::Electron &ele) const final
Definition: EGRegressionModifierV1.cc:29
EGRegressionModifierV1::phoForestsMean_
std::vector< const GBRForestD * > phoForestsMean_
Definition: EGRegressionModifierV1.cc:52
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:19
EGRegressionModifierV1::EGRegressionModifierV1
EGRegressionModifierV1(const edm::ParameterSet &conf, edm::ConsumesCollector &cc)
Definition: EGRegressionModifierV1.cc:61
reco::GsfElectron::trackMomentumAtVtx
math::XYZVectorF trackMomentumAtVtx() const
Definition: GsfElectron.h:276
EGRegressionModifierV1::phoCond25nsTokens_
EGRegressionModifierCondTokens phoCond25nsTokens_
Definition: EGRegressionModifierV1.cc:36
reco::GsfElectron
Definition: GsfElectron.h:35
GsfElectron.h
DEFINE_EDM_PLUGIN
#define DEFINE_EDM_PLUGIN(factory, type, name)
Definition: PluginFactory.h:124
EDGetToken.h
EGRegressionModifierV1::eleForestsSigma_
std::vector< const GBRForestD * > eleForestsSigma_
Definition: EGRegressionModifierV1.cc:55
EEDetId
Definition: EEDetId.h:14
CaloGeometryRecord.h
EGRegressionModifierV1::epForest_
const GBRForest * epForest_
Definition: EGRegressionModifierV1.cc:56
EGRegressionModifierV1::bunchspacing_
int bunchspacing_
Definition: EGRegressionModifierV1.cc:42
egammaTools::localEcalClusterCoordsEB
void localEcalClusterCoordsEB(const reco::CaloCluster &bclus, const CaloGeometry &geom, float &etacry, float &phicry, int &ieta, int &iphi, float &thetatilt, float &phitilt)
Definition: EcalClusterLocal.cc:13
LEDCalibrationChannels.ieta
ieta
Definition: LEDCalibrationChannels.py:63
reco::GsfElectron::setCorrectedEcalEnergyError
void setCorrectedEcalEnergyError(float newEnergyError)
Definition: GsfElectron.cc:170
edm::ParameterSet
Definition: ParameterSet.h:47
reco::Photon::r9
float r9() const
Definition: Photon.h:273
ModifyObjectValueBase
Definition: ModifyObjectValueBase.h:18
reco::CaloCluster::eta
double eta() const
pseudorapidity of cluster centroid
Definition: CaloCluster.h:181
SiStripPI::max
Definition: SiStripPayloadInspectorHelper.h:169
EGRegressionModifierV1::caloGeom_
CaloGeometry const * caloGeom_
Definition: EGRegressionModifierV1.cc:49
edmplugin::PluginFactory
Definition: PluginFactory.h:34
reco::Photon::isEB
bool isEB() const
Definition: Photon.h:120
createfilelist.int
int
Definition: createfilelist.py:10
reco::GsfElectron::ecalDriven
bool ecalDriven() const
Definition: GsfElectron.cc:168
EGRegressionModifierCondTokens
Definition: EGRegressionModifierHelpers.h:13
edm::EventSetup
Definition: EventSetup.h:58
EGRegressionModifierV1::rhoToken_
edm::EDGetTokenT< double > rhoToken_
Definition: EGRegressionModifierV1.cc:45
reco::Photon::setCorrectedEnergy
void setCorrectedEnergy(P4type type, float E, float dE, bool toCand=true)
EGRegressionModifierV1::condNamesWeight25nsToken_
edm::ESGetToken< GBRForest, GBRWrapperRcd > condNamesWeight25nsToken_
Definition: EGRegressionModifierV1.cc:39
cc
reco::Photon::e5x5
float e5x5() const
Definition: Photon.h:267
reco::GsfElectron::p4
const LorentzVector & p4(P4Kind kind) const
Definition: GsfElectron.cc:217
reco::CaloCluster::seed
DetId seed() const
return DetId of seed
Definition: CaloCluster.h:219
edm::ESGetToken< GBRForest, GBRWrapperRcd >
EGRegressionModifierV1::setEventContent
void setEventContent(const edm::EventSetup &) final
Definition: EGRegressionModifierV1.cc:89
InputTag.h
edm::Ptr
Definition: AssociationVector.h:31
edm::EventSetup::getData
bool getData(T &iHolder) const
Definition: EventSetup.h:127
ValueMap.h
EcalClusterLocal.h
reco::Photon
Definition: Photon.h:22
EGRegressionModifierCondTokens::sigma
std::vector< edm::ESGetToken< GBRForestD, GBRDWrapperRcd > > sigma
Definition: EGRegressionModifierHelpers.h:19
pwdgSkimBPark_cfi.electrons
electrons
Definition: pwdgSkimBPark_cfi.py:6
math::XYZTLorentzVector
XYZTLorentzVectorD XYZTLorentzVector
Lorentz vector with cylindrical internal representation using pseudorapidity.
Definition: LorentzVector.h:29
EGRegressionModifierV1::condNamesWeight50nsToken_
edm::ESGetToken< GBRForest, GBRWrapperRcd > condNamesWeight50nsToken_
Definition: EGRegressionModifierV1.cc:38
EGRegressionModifierV1::modifyObject
void modifyObject(pat::Photon &pho) const final
Definition: EGRegressionModifierV1.cc:30
CaloGeometry.h
reco::GsfElectron::classification
Classification classification() const
Definition: GsfElectron.h:792
reco::deltaR
constexpr auto deltaR(const T1 &t1, const T2 &t2) -> decltype(t1.eta())
Definition: deltaR.h:30
reco::GsfElectron::trackerDrivenSeed
bool trackerDrivenSeed() const
Definition: GsfElectron.h:167
reco::GsfElectron::superCluster
SuperClusterRef superCluster() const override
reference to a SuperCluster
Definition: GsfElectron.h:163
EGRegressionModifierV1::rhoValue_
float rhoValue_
Definition: EGRegressionModifierV1.cc:44
EGRegressionModifierV1::eleForestsMean_
std::vector< const GBRForestD * > eleForestsMean_
Definition: EGRegressionModifierV1.cc:54
EGRegressionModifierV1::eleCond25nsTokens_
EGRegressionModifierCondTokens eleCond25nsTokens_
Definition: EGRegressionModifierV1.cc:35
L1TowerCalibrationProducer_cfi.iEta
iEta
Definition: L1TowerCalibrationProducer_cfi.py:60
reco::GsfElectron::trackMomentumError
float trackMomentumError() const
Definition: GsfElectron.h:871
pat::Electron
Analysis-level electron class.
Definition: Electron.h:51
funct::abs
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
jets_cff.maxDR
maxDR
Definition: jets_cff.py:130
dummy
Definition: DummySelector.h:38
reco::GsfElectron::setCorrectedEcalEnergy
void setCorrectedEcalEnergy(float newEnergy)
Definition: GsfElectron.cc:174
EGRegressionModifierV1::phoForestsSigma_
std::vector< const GBRForestD * > phoForestsSigma_
Definition: EGRegressionModifierV1.cc:53
edm::Event
Definition: Event.h:73
EGRegressionModifierV1::applyExtraHighEnergyProtection_
const bool applyExtraHighEnergyProtection_
Definition: EGRegressionModifierV1.cc:50
EGRegressionModifierV1::eleCond50nsTokens_
EGRegressionModifierCondTokens eleCond50nsTokens_
Definition: EGRegressionModifierV1.cc:33
reco::CaloCluster::energy
double energy() const
cluster energy
Definition: CaloCluster.h:149
edm::Event::get
bool get(ProductID const &oid, Handle< PROD > &result) const
Definition: Event.h:346
SiStripBadComponentsDQMServiceTemplate_cfg.ep
ep
Definition: SiStripBadComponentsDQMServiceTemplate_cfg.py:86
edm::InputTag
Definition: InputTag.h:15
edm::ConsumesCollector
Definition: ConsumesCollector.h:45
weight
Definition: weight.py:1
edm::ParameterSet::getParameterSet
ParameterSet const & getParameterSet(std::string const &) const
Definition: ParameterSet.cc:2128
Generated for CMSSW Reference Manual by doxygen 1.8.16