CMS 3D CMS Logo
DPFIsolation.cc
Go to the documentation of this file.
1 /*
2  * \class DPFIsolation
3  *
4  * Deep ParticleFlow tau isolation using Deep NN.
5  *
6  * \author Owen Colegrove, UCSB
7  */
8 
9 #include "RecoTauTag/RecoTau/interface/DeepTauBase.h"
10 
11 namespace {
12  inline int getPFCandidateIndex(const edm::Handle<pat::PackedCandidateCollection>& pfcands,
13  const reco::CandidatePtr& cptr) {
14  for (unsigned int i = 0; i < pfcands->size(); ++i) {
15  if (reco::CandidatePtr(pfcands, i) == cptr)
16  return i;
17  }
18  return -1;
19  }
20 } // anonymous namespace
21 
22 class DPFIsolation : public deep_tau::DeepTauBase {
23 public:
24  static const OutputCollection& GetOutputs() {
25  const size_t tau_index = 0;
26  static const OutputCollection outputs_ = {{"VSall", Output({tau_index}, {})}};
27  return outputs_;
28  };
29 
30  static unsigned getNumberOfParticles(unsigned graphVersion) {
31  static const std::map<unsigned, unsigned> nparticles{{0, 60}, {1, 36}};
32  return nparticles.at(graphVersion);
33  }
34 
35  static unsigned GetNumberOfFeatures(unsigned graphVersion) {
36  static const std::map<unsigned, unsigned> nfeatures{{0, 47}, {1, 51}};
37  return nfeatures.at(graphVersion);
38  }
39 
40  static void fillDescriptions(edm::ConfigurationDescriptions& descriptions) {
41  edm::ParameterSetDescription desc;
42  desc.add<edm::InputTag>("pfcands", edm::InputTag("packedPFCandidates"));
43  desc.add<edm::InputTag>("taus", edm::InputTag("slimmedTaus"));
44  desc.add<edm::InputTag>("vertices", edm::InputTag("offlineSlimmedPrimaryVertices"));
45  desc.add<std::string>("graph_file", "RecoTauTag/TrainingFiles/data/DPFTauId/DPFIsolation_2017v0_quantized.pb");
46  desc.add<unsigned>("version", 0);
47  desc.add<bool>("mem_mapped", false);
48 
49  edm::ParameterSetDescription descWP;
50  descWP.add<std::string>("VVVLoose", "0");
51  descWP.add<std::string>("VVLoose", "0");
52  descWP.add<std::string>("VLoose", "0");
53  descWP.add<std::string>("Loose", "0");
54  descWP.add<std::string>("Medium", "0");
55  descWP.add<std::string>("Tight", "0");
56  descWP.add<std::string>("VTight", "0");
57  descWP.add<std::string>("VVTight", "0");
58  descWP.add<std::string>("VVVTight", "0");
59  desc.add<edm::ParameterSetDescription>("VSallWP", descWP);
60  descriptions.add("DPFTau2016v0", desc);
61  }
62 
63  explicit DPFIsolation(const edm::ParameterSet& cfg, const deep_tau::DeepTauCache* cache)
64  : DeepTauBase(cfg, GetOutputs(), cache), graphVersion(cfg.getParameter<unsigned>("version")) {
65  const auto& shape = cache_->getGraph().node(0).attr().at("shape").shape();
66 
67  if (!(graphVersion == 1 || graphVersion == 0))
68  throw cms::Exception("DPFIsolation") << "unknown version of the graph file.";
69 
70  if (!(shape.dim(1).size() == getNumberOfParticles(graphVersion) &&
71  shape.dim(2).size() == GetNumberOfFeatures(graphVersion)))
72  throw cms::Exception("DPFIsolation")
73  << "number of inputs does not match the expected inputs for the given version";
74  }
75 
76 private:
77  tensorflow::Tensor getPredictions(edm::Event& event, edm::Handle<TauCollection> taus) override {
78  edm::Handle<pat::PackedCandidateCollection> pfcands;
79  event.getByToken(pfcandToken_, pfcands);
80 
81  edm::Handle<reco::VertexCollection> vertices;
82  event.getByToken(vtxToken_, vertices);
83 
84  tensorflow::Tensor tensor(
85  tensorflow::DT_FLOAT,
86  {1, static_cast<int>(getNumberOfParticles(graphVersion)), static_cast<int>(GetNumberOfFeatures(graphVersion))});
87 
88  tensorflow::Tensor predictions(tensorflow::DT_FLOAT, {static_cast<int>(taus->size()), 1});
89 
90  std::vector<tensorflow::Tensor> outputs_;
91 
92  float pfCandPt, pfCandPz, pfCandPtRel, pfCandPzRel, pfCandDr, pfCandDEta, pfCandDPhi, pfCandEta, pfCandDz,
93  pfCandDzErr, pfCandD0, pfCandD0D0, pfCandD0Dz, pfCandD0Dphi, pfCandPuppiWeight, pfCandPixHits, pfCandHits,
94  pfCandLostInnerHits, pfCandPdgID, pfCandCharge, pfCandFromPV, pfCandVtxQuality, pfCandHighPurityTrk,
95  pfCandTauIndMatch, pfCandDzSig, pfCandD0Sig, pfCandD0Err, pfCandPtRelPtRel, pfCandDzDz, pfCandDVx_1,
96  pfCandDVy_1, pfCandDVz_1, pfCandD_1;
97  float pvx = !vertices->empty() ? (*vertices)[0].x() : -1;
98  float pvy = !vertices->empty() ? (*vertices)[0].y() : -1;
99  float pvz = !vertices->empty() ? (*vertices)[0].z() : -1;
100 
101  bool pfCandIsBarrel;
102 
103  // These variables define ranges further used for standardization
104  static constexpr float pfCandPt_max = 500.f;
105  static constexpr float pfCandPz_max = 1000.f;
106  static constexpr float pfCandPtRel_max = 1.f;
107  static constexpr float pfCandPzRel_max = 100.f;
108  static constexpr float pfCandPtRelPtRel_max = 1.f;
109  static constexpr float pfCandD0_max = 5.f;
110  static constexpr float pfCandDz_max = 5.f;
111  static constexpr float pfCandDVx_y_z_1_max = 0.05f;
112  static constexpr float pfCandD_1_max = 0.1f;
113  static constexpr float pfCandD0_z_Err_max = 1.f;
114  static constexpr float pfCandDzSig_max = 3.f;
115  static constexpr float pfCandD0Sig_max = 1.f;
116  static constexpr float pfCandDr_max = 0.5f;
117  static constexpr float pfCandEta_max = 2.75f;
118  static constexpr float pfCandDEta_max = 0.5f;
119  static constexpr float pfCandDPhi_max = 0.5f;
120  static constexpr float pfCandPixHits_max = 7.f;
121  static constexpr float pfCandHits_max = 30.f;
122 
123  for (size_t tau_index = 0; tau_index < taus->size(); tau_index++) {
124  pat::Tau tau = taus->at(tau_index);
125  bool isGoodTau = false;
126  const float lepRecoPt = tau.pt();
127  const float lepRecoPz = std::abs(tau.pz());
128  const float lepRecoEta = tau.eta();
129  const float lepRecoPhi = tau.phi();
130 
131  if (lepRecoPt >= 30 && std::abs(lepRecoEta) < 2.3 && tau.isTauIDAvailable("againstMuonLoose3") &&
132  tau.isTauIDAvailable("againstElectronVLooseMVA6")) {
133  isGoodTau = (tau.tauID("againstElectronVLooseMVA6") && tau.tauID("againstMuonLoose3"));
134  }
135 
136  if (!isGoodTau) {
137  predictions.matrix<float>()(tau_index, 0) = -1;
138  continue;
139  }
140 
141  std::vector<unsigned int> signalCandidateInds;
142 
143  for (const auto& c : tau.signalCands())
144  signalCandidateInds.push_back(getPFCandidateIndex(pfcands, c));
145 
146  // Use of setZero results in warnings in eigen library during compilation.
147  //tensor.flat<float>().setZero();
148  const unsigned n_inputs = getNumberOfParticles(graphVersion) * GetNumberOfFeatures(graphVersion);
149  for (unsigned input_idx = 0; input_idx < n_inputs; ++input_idx)
150  tensor.flat<float>()(input_idx) = 0;
151 
152  unsigned int iPF = 0;
153  const unsigned max_iPF = getNumberOfParticles(graphVersion);
154 
155  std::vector<unsigned int> sorted_inds(pfcands->size());
156  std::size_t n = 0;
157  std::generate(std::begin(sorted_inds), std::end(sorted_inds), [&] { return n++; });
158 
159  std::sort(std::begin(sorted_inds), std::end(sorted_inds), [&](int i1, int i2) {
160  return pfcands->at(i1).pt() > pfcands->at(i2).pt();
161  });
162 
163  for (size_t pf_index = 0; pf_index < pfcands->size() && iPF < max_iPF; pf_index++) {
164  pat::PackedCandidate p = pfcands->at(sorted_inds.at(pf_index));
165  float deltaR_tau_p = deltaR(p.p4(), tau.p4());
166 
167  if (p.pt() < 0.5)
168  continue;
169  if (deltaR_tau_p > 0.5)
170  continue;
171  if (p.fromPV() < 1 && p.charge() != 0)
172  continue;
173  pfCandPt = p.pt();
174  pfCandPtRel = p.pt() / lepRecoPt;
175 
176  pfCandDr = deltaR_tau_p;
177  pfCandDEta = std::abs(lepRecoEta - p.eta());
178  pfCandDPhi = std::abs(deltaPhi(lepRecoPhi, p.phi()));
179  pfCandEta = p.eta();
180  pfCandIsBarrel = (std::abs(pfCandEta) < 1.4);
181  pfCandPz = std::abs(std::sinh(pfCandEta) * pfCandPt);
182  pfCandPzRel = pfCandPz / lepRecoPz;
183  pfCandPdgID = std::abs(p.pdgId());
184  pfCandCharge = p.charge();
185  pfCandDVx_1 = p.vx() - pvx;
186  pfCandDVy_1 = p.vy() - pvy;
187  pfCandDVz_1 = p.vz() - pvz;
188 
189  pfCandD_1 = std::sqrt(pfCandDVx_1 * pfCandDVx_1 + pfCandDVy_1 * pfCandDVy_1 + pfCandDVz_1 * pfCandDVz_1);
190 
191  if (pfCandCharge != 0 and p.hasTrackDetails()) {
192  pfCandDz = p.dz();
193  pfCandDzErr = p.dzError();
194  pfCandDzSig = (std::abs(p.dz()) + 0.000001) / (p.dzError() + 0.00001);
195  pfCandD0 = p.dxy();
196  pfCandD0Err = p.dxyError();
197  pfCandD0Sig = (std::abs(p.dxy()) + 0.000001) / (p.dxyError() + 0.00001);
198  pfCandPixHits = p.numberOfPixelHits();
199  pfCandHits = p.numberOfHits();
200  pfCandLostInnerHits = p.lostInnerHits();
201  } else {
202  float disp = 1;
203  int psudorand = p.pt() * 1000000;
204  if (psudorand % 2 == 0)
205  disp = -1;
206  pfCandDz = 5 * disp;
207  pfCandDzErr = 0;
208  pfCandDzSig = 0;
209  pfCandD0 = 5 * disp;
210  pfCandD0Err = 0;
211  pfCandD0Sig = 0;
212  pfCandPixHits = 0;
213  pfCandHits = 0;
214  pfCandLostInnerHits = 2.;
215  pfCandDVx_1 = 1;
216  pfCandDVy_1 = 1;
217  pfCandDVz_1 = 1;
218  pfCandD_1 = 1;
219  }
220 
221  pfCandPuppiWeight = p.puppiWeight();
222  pfCandFromPV = p.fromPV();
223  pfCandVtxQuality = p.pvAssociationQuality();
224  pfCandHighPurityTrk = p.trackHighPurity();
225  float pfCandTauIndMatch_temp = 0;
226 
227  for (auto i : signalCandidateInds) {
228  if (i == sorted_inds.at(pf_index))
229  pfCandTauIndMatch_temp = 1;
230  }
231 
232  pfCandTauIndMatch = pfCandTauIndMatch_temp;
233  pfCandPtRelPtRel = pfCandPtRel * pfCandPtRel;
234  pfCandPt = std::min(pfCandPt, pfCandPt_max);
235  pfCandPt = pfCandPt / pfCandPt_max;
236 
237  pfCandPz = std::min(pfCandPz, pfCandPz_max);
238  pfCandPz = pfCandPz / pfCandPz_max;
239 
240  pfCandPtRel = std::min(pfCandPtRel, pfCandPtRel_max);
241  pfCandPzRel = std::min(pfCandPzRel, pfCandPzRel_max);
242  pfCandPzRel = pfCandPzRel / pfCandPzRel_max;
243  pfCandDr = pfCandDr / pfCandDr_max;
244  pfCandEta = pfCandEta / pfCandEta_max;
245  pfCandDEta = pfCandDEta / pfCandDEta_max;
246  pfCandDPhi = pfCandDPhi / pfCandDPhi_max;
247  pfCandPixHits = pfCandPixHits / pfCandPixHits_max;
248  pfCandHits = pfCandHits / pfCandHits_max;
249 
250  pfCandPtRelPtRel = std::min(pfCandPtRelPtRel, pfCandPtRelPtRel_max);
251 
252  pfCandD0 = std::clamp(pfCandD0, -pfCandD0_max, pfCandD0_max);
253  pfCandD0 = pfCandD0 / pfCandD0_max;
254 
255  pfCandDz = std::clamp(pfCandDz, -pfCandDz_max, pfCandDz_max);
256  pfCandDz = pfCandDz / pfCandDz_max;
257 
258  pfCandD0Err = std::min(pfCandD0Err, pfCandD0_z_Err_max);
259  pfCandDzErr = std::min(pfCandDzErr, pfCandD0_z_Err_max);
260  pfCandDzSig = std::min(pfCandDzSig, pfCandDzSig_max);
261  pfCandDzSig = pfCandDzSig / pfCandDzSig_max;
262 
263  pfCandD0Sig = std::min(pfCandD0Sig, pfCandD0Sig_max);
264  pfCandD0D0 = pfCandD0 * pfCandD0;
265  pfCandDzDz = pfCandDz * pfCandDz;
266  pfCandD0Dz = pfCandD0 * pfCandDz;
267  pfCandD0Dphi = pfCandD0 * pfCandDPhi;
268 
269  pfCandDVx_1 = std::clamp(pfCandDVx_1, -pfCandDVx_y_z_1_max, pfCandDVx_y_z_1_max);
270  pfCandDVx_1 = pfCandDVx_1 / pfCandDVx_y_z_1_max;
271 
272  pfCandDVy_1 = std::clamp(pfCandDVy_1, -pfCandDVx_y_z_1_max, pfCandDVx_y_z_1_max);
273  pfCandDVy_1 = pfCandDVy_1 / pfCandDVx_y_z_1_max;
274 
275  pfCandDVz_1 = std::clamp(pfCandDVz_1, -pfCandDVx_y_z_1_max, pfCandDVx_y_z_1_max);
276  pfCandDVz_1 = pfCandDVz_1 / pfCandDVx_y_z_1_max;
277 
278  pfCandD_1 = std::clamp(pfCandD_1, -pfCandD_1_max, pfCandD_1_max);
279  pfCandD_1 = pfCandD_1 / pfCandD_1_max;
280 
281  if (graphVersion == 0) {
282  tensor.tensor<float, 3>()(0, 60 - 1 - iPF, 0) = pfCandPt;
283  tensor.tensor<float, 3>()(0, 60 - 1 - iPF, 1) = pfCandPz;
284  tensor.tensor<float, 3>()(0, 60 - 1 - iPF, 2) = pfCandPtRel;
285  tensor.tensor<float, 3>()(0, 60 - 1 - iPF, 3) = pfCandPzRel;
286  tensor.tensor<float, 3>()(0, 60 - 1 - iPF, 4) = pfCandDr;
287  tensor.tensor<float, 3>()(0, 60 - 1 - iPF, 5) = pfCandDEta;
288  tensor.tensor<float, 3>()(0, 60 - 1 - iPF, 6) = pfCandDPhi;
289  tensor.tensor<float, 3>()(0, 60 - 1 - iPF, 7) = pfCandEta;
290  tensor.tensor<float, 3>()(0, 60 - 1 - iPF, 8) = pfCandDz;
291  tensor.tensor<float, 3>()(0, 60 - 1 - iPF, 9) = pfCandDzSig;
292  tensor.tensor<float, 3>()(0, 60 - 1 - iPF, 10) = pfCandD0;
293  tensor.tensor<float, 3>()(0, 60 - 1 - iPF, 11) = pfCandD0Sig;
294  tensor.tensor<float, 3>()(0, 60 - 1 - iPF, 12) = pfCandDzErr;
295  tensor.tensor<float, 3>()(0, 60 - 1 - iPF, 13) = pfCandD0Err;
296  tensor.tensor<float, 3>()(0, 60 - 1 - iPF, 14) = pfCandD0D0;
297  tensor.tensor<float, 3>()(0, 60 - 1 - iPF, 15) = pfCandCharge == 0;
298  tensor.tensor<float, 3>()(0, 60 - 1 - iPF, 16) = pfCandCharge == 1;
299  tensor.tensor<float, 3>()(0, 60 - 1 - iPF, 17) = pfCandCharge == -1;
300  tensor.tensor<float, 3>()(0, 60 - 1 - iPF, 18) = pfCandPdgID > 22;
301  tensor.tensor<float, 3>()(0, 60 - 1 - iPF, 19) = pfCandPdgID == 22;
302  tensor.tensor<float, 3>()(0, 60 - 1 - iPF, 20) = pfCandDzDz;
303  tensor.tensor<float, 3>()(0, 60 - 1 - iPF, 21) = pfCandD0Dz;
304  tensor.tensor<float, 3>()(0, 60 - 1 - iPF, 22) = pfCandD0Dphi;
305  tensor.tensor<float, 3>()(0, 60 - 1 - iPF, 23) = pfCandPtRelPtRel;
306  tensor.tensor<float, 3>()(0, 60 - 1 - iPF, 24) = pfCandPixHits;
307  tensor.tensor<float, 3>()(0, 60 - 1 - iPF, 25) = pfCandHits;
308  tensor.tensor<float, 3>()(0, 60 - 1 - iPF, 26) = pfCandLostInnerHits == -1;
309  tensor.tensor<float, 3>()(0, 60 - 1 - iPF, 27) = pfCandLostInnerHits == 0;
310  tensor.tensor<float, 3>()(0, 60 - 1 - iPF, 28) = pfCandLostInnerHits == 1;
311  tensor.tensor<float, 3>()(0, 60 - 1 - iPF, 29) = pfCandLostInnerHits == 2;
312  tensor.tensor<float, 3>()(0, 60 - 1 - iPF, 30) = pfCandPuppiWeight;
313  tensor.tensor<float, 3>()(0, 60 - 1 - iPF, 31) = (pfCandVtxQuality == 1);
314  tensor.tensor<float, 3>()(0, 60 - 1 - iPF, 32) = (pfCandVtxQuality == 5);
315  tensor.tensor<float, 3>()(0, 60 - 1 - iPF, 33) = (pfCandVtxQuality == 6);
316  tensor.tensor<float, 3>()(0, 60 - 1 - iPF, 34) = (pfCandVtxQuality == 7);
317  tensor.tensor<float, 3>()(0, 60 - 1 - iPF, 35) = (pfCandFromPV == 1);
318  tensor.tensor<float, 3>()(0, 60 - 1 - iPF, 36) = (pfCandFromPV == 2);
319  tensor.tensor<float, 3>()(0, 60 - 1 - iPF, 37) = (pfCandFromPV == 3);
320  tensor.tensor<float, 3>()(0, 60 - 1 - iPF, 38) = pfCandIsBarrel;
321  tensor.tensor<float, 3>()(0, 60 - 1 - iPF, 39) = pfCandHighPurityTrk;
322  tensor.tensor<float, 3>()(0, 60 - 1 - iPF, 40) = pfCandPdgID == 1;
323  tensor.tensor<float, 3>()(0, 60 - 1 - iPF, 41) = pfCandPdgID == 2;
324  tensor.tensor<float, 3>()(0, 60 - 1 - iPF, 42) = pfCandPdgID == 11;
325  tensor.tensor<float, 3>()(0, 60 - 1 - iPF, 43) = pfCandPdgID == 13;
326  tensor.tensor<float, 3>()(0, 60 - 1 - iPF, 44) = pfCandPdgID == 130;
327  tensor.tensor<float, 3>()(0, 60 - 1 - iPF, 45) = pfCandPdgID == 211;
328  tensor.tensor<float, 3>()(0, 60 - 1 - iPF, 46) = pfCandTauIndMatch;
329  }
330 
331  if (graphVersion == 1) {
332  tensor.tensor<float, 3>()(0, 36 - 1 - iPF, 0) = pfCandPt;
333  tensor.tensor<float, 3>()(0, 36 - 1 - iPF, 1) = pfCandPz;
334  tensor.tensor<float, 3>()(0, 36 - 1 - iPF, 2) = pfCandPtRel;
335  tensor.tensor<float, 3>()(0, 36 - 1 - iPF, 3) = pfCandPzRel;
336  tensor.tensor<float, 3>()(0, 36 - 1 - iPF, 4) = pfCandDr;
337  tensor.tensor<float, 3>()(0, 36 - 1 - iPF, 5) = pfCandDEta;
338  tensor.tensor<float, 3>()(0, 36 - 1 - iPF, 6) = pfCandDPhi;
339  tensor.tensor<float, 3>()(0, 36 - 1 - iPF, 7) = pfCandEta;
340  tensor.tensor<float, 3>()(0, 36 - 1 - iPF, 8) = pfCandDz;
341  tensor.tensor<float, 3>()(0, 36 - 1 - iPF, 9) = pfCandDzSig;
342  tensor.tensor<float, 3>()(0, 36 - 1 - iPF, 10) = pfCandD0;
343  tensor.tensor<float, 3>()(0, 36 - 1 - iPF, 11) = pfCandD0Sig;
344  tensor.tensor<float, 3>()(0, 36 - 1 - iPF, 12) = pfCandDzErr;
345  tensor.tensor<float, 3>()(0, 36 - 1 - iPF, 13) = pfCandD0Err;
346  tensor.tensor<float, 3>()(0, 36 - 1 - iPF, 14) = pfCandD0D0;
347  tensor.tensor<float, 3>()(0, 36 - 1 - iPF, 15) = pfCandCharge == 0;
348  tensor.tensor<float, 3>()(0, 36 - 1 - iPF, 16) = pfCandCharge == 1;
349  tensor.tensor<float, 3>()(0, 36 - 1 - iPF, 17) = pfCandCharge == -1;
350  tensor.tensor<float, 3>()(0, 36 - 1 - iPF, 18) = pfCandPdgID > 22;
351  tensor.tensor<float, 3>()(0, 36 - 1 - iPF, 19) = pfCandPdgID == 22;
352  tensor.tensor<float, 3>()(0, 36 - 1 - iPF, 20) = pfCandDVx_1;
353  tensor.tensor<float, 3>()(0, 36 - 1 - iPF, 21) = pfCandDVy_1;
354  tensor.tensor<float, 3>()(0, 36 - 1 - iPF, 22) = pfCandDVz_1;
355  tensor.tensor<float, 3>()(0, 36 - 1 - iPF, 23) = pfCandD_1;
356  tensor.tensor<float, 3>()(0, 36 - 1 - iPF, 24) = pfCandDzDz;
357  tensor.tensor<float, 3>()(0, 36 - 1 - iPF, 25) = pfCandD0Dz;
358  tensor.tensor<float, 3>()(0, 36 - 1 - iPF, 26) = pfCandD0Dphi;
359  tensor.tensor<float, 3>()(0, 36 - 1 - iPF, 27) = pfCandPtRelPtRel;
360  tensor.tensor<float, 3>()(0, 36 - 1 - iPF, 28) = pfCandPixHits;
361  tensor.tensor<float, 3>()(0, 36 - 1 - iPF, 29) = pfCandHits;
362  tensor.tensor<float, 3>()(0, 36 - 1 - iPF, 30) = pfCandLostInnerHits == -1;
363  tensor.tensor<float, 3>()(0, 36 - 1 - iPF, 31) = pfCandLostInnerHits == 0;
364  tensor.tensor<float, 3>()(0, 36 - 1 - iPF, 32) = pfCandLostInnerHits == 1;
365  tensor.tensor<float, 3>()(0, 36 - 1 - iPF, 33) = pfCandLostInnerHits == 2;
366  tensor.tensor<float, 3>()(0, 36 - 1 - iPF, 34) = pfCandPuppiWeight;
367  tensor.tensor<float, 3>()(0, 36 - 1 - iPF, 35) = (pfCandVtxQuality == 1);
368  tensor.tensor<float, 3>()(0, 36 - 1 - iPF, 36) = (pfCandVtxQuality == 5);
369  tensor.tensor<float, 3>()(0, 36 - 1 - iPF, 37) = (pfCandVtxQuality == 6);
370  tensor.tensor<float, 3>()(0, 36 - 1 - iPF, 38) = (pfCandVtxQuality == 7);
371  tensor.tensor<float, 3>()(0, 36 - 1 - iPF, 39) = (pfCandFromPV == 1);
372  tensor.tensor<float, 3>()(0, 36 - 1 - iPF, 40) = (pfCandFromPV == 2);
373  tensor.tensor<float, 3>()(0, 36 - 1 - iPF, 41) = (pfCandFromPV == 3);
374  tensor.tensor<float, 3>()(0, 36 - 1 - iPF, 42) = pfCandIsBarrel;
375  tensor.tensor<float, 3>()(0, 36 - 1 - iPF, 43) = pfCandHighPurityTrk;
376  tensor.tensor<float, 3>()(0, 36 - 1 - iPF, 44) = pfCandPdgID == 1;
377  tensor.tensor<float, 3>()(0, 36 - 1 - iPF, 45) = pfCandPdgID == 2;
378  tensor.tensor<float, 3>()(0, 36 - 1 - iPF, 46) = pfCandPdgID == 11;
379  tensor.tensor<float, 3>()(0, 36 - 1 - iPF, 47) = pfCandPdgID == 13;
380  tensor.tensor<float, 3>()(0, 36 - 1 - iPF, 48) = pfCandPdgID == 130;
381  tensor.tensor<float, 3>()(0, 36 - 1 - iPF, 49) = pfCandPdgID == 211;
382  tensor.tensor<float, 3>()(0, 36 - 1 - iPF, 50) = pfCandTauIndMatch;
383  }
384  iPF++;
385  }
386  tensorflow::run(&(cache_->getSession()), {{"input_1", tensor}}, {"output_node0"}, &outputs_);
387  predictions.matrix<float>()(tau_index, 0) = outputs_[0].flat<float>()(0);
388  }
389  return predictions;
390  }
391 
392 private:
393  unsigned graphVersion;
394 };
395 
396 #include "FWCore/Framework/interface/MakerMacros.h"
397 DEFINE_FWK_MODULE(DPFIsolation);
deep_tau::DeepTauBase::cache_
const DeepTauCache * cache_
Definition: DeepTauBase.h:135
testProducerWithPsetDescEmpty_cfi.i2
i2
Definition: testProducerWithPsetDescEmpty_cfi.py:46
DeepTauBase.h
mps_fire.i
i
Definition: mps_fire.py:428
edm::ParameterSetDescription::add
ParameterDescriptionBase * add(U const &iLabel, T const &value)
Definition: ParameterSetDescription.h:95
dqmMemoryStats.float
float
Definition: dqmMemoryStats.py:127
deep_tau::DeepTauCache::getSession
tensorflow::Session & getSession(const std::string &name="") const
Definition: DeepTauBase.h:57
dqmiodumpmetadata.n
n
Definition: dqmiodumpmetadata.py:28
deep_tau::DeepTauBase::pfcandToken_
edm::EDGetTokenT< CandidateCollection > pfcandToken_
Definition: DeepTauBase.h:130
metsig::tau
Definition: SignAlgoResolutions.h:49
min
T min(T a, T b)
Definition: MathUtil.h:58
Tau3MuMonitor_cff.taus
taus
Definition: Tau3MuMonitor_cff.py:7
AlCaHLTBitMon_ParallelJobs.p
p
Definition: AlCaHLTBitMon_ParallelJobs.py:153
pat::Tau
Analysis-level tau class.
Definition: Tau.h:53
deep_tau::DeepTauBase::OutputCollection
std::map< std::string, Output > OutputCollection
Definition: DeepTauBase.h:91
testProducerWithPsetDescEmpty_cfi.i1
i1
Definition: testProducerWithPsetDescEmpty_cfi.py:45
HLT_FULL_cff.InputTag
InputTag
Definition: HLT_FULL_cff.py:89281
edm::ParameterSetDescription
Definition: ParameterSetDescription.h:52
DPFIsolation::DPFIsolation
DPFIsolation(const edm::ParameterSet &cfg, const deep_tau::DeepTauCache *cache)
Definition: DPFIsolation.cc:63
deep_tau::DeepTauBase::vtxToken_
edm::EDGetTokenT< reco::VertexCollection > vtxToken_
Definition: DeepTauBase.h:131
DPFIsolation::getPredictions
tensorflow::Tensor getPredictions(edm::Event &event, edm::Handle< TauCollection > taus) override
Definition: DPFIsolation.cc:77
edm::Handle< pat::PackedCandidateCollection >
GetRecoTauVFromDQM_MC_cff.Output
Output
Definition: GetRecoTauVFromDQM_MC_cff.py:12
DPFIsolation::GetOutputs
static const OutputCollection & GetOutputs()
Definition: DPFIsolation.cc:24
deep_tau::DeepTauBase
Definition: DeepTauBase.h:66
generate
void generate(uint32_t const nbins, float const *initValues, std::vector< float > &values)
Definition: generateQIEShapes.cc:48
deep_tau::DeepTauBase::DeepTauBase
DeepTauBase(const edm::ParameterSet &cfg, const OutputCollection &outputs, const DeepTauCache *cache)
Definition: DeepTauBase.cc:86
MakerMacros.h
DEFINE_FWK_MODULE
#define DEFINE_FWK_MODULE(type)
Definition: MakerMacros.h:16
edm::ConfigurationDescriptions::add
void add(std::string const &label, ParameterSetDescription const &psetDescription)
Definition: ConfigurationDescriptions.cc:57
SiPixelRawToDigiRegional_cfi.deltaPhi
deltaPhi
Definition: SiPixelRawToDigiRegional_cfi.py:9
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:19
mps_fire.end
end
Definition: mps_fire.py:242
DPFIsolation
Definition: DPFIsolation.cc:22
utilities.cache
def cache(function)
Definition: utilities.py:3
PbPb_ZMuSkimMuonDPG_cff.deltaR
deltaR
Definition: PbPb_ZMuSkimMuonDPG_cff.py:63
edm::ConfigurationDescriptions
Definition: ConfigurationDescriptions.h:28
AlCaHLTBitMon_QueryRunRegistry.string
string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
deep_tau::DeepTauBase::outputs_
OutputCollection outputs_
Definition: DeepTauBase.h:134
edm::ParameterSet
Definition: ParameterSet.h:47
pat::PackedCandidate
Definition: PackedCandidate.h:22
jetUpdater_cfi.sort
sort
Definition: jetUpdater_cfi.py:29
DPFIsolation::fillDescriptions
static void fillDescriptions(edm::ConfigurationDescriptions &descriptions)
Definition: DPFIsolation.cc:40
DPFIsolation::getNumberOfParticles
static unsigned getNumberOfParticles(unsigned graphVersion)
Definition: DPFIsolation.cc:30
edm::Ptr< Candidate >
looper.cfg
cfg
Definition: looper.py:297
submitPVResolutionJobs.desc
string desc
Definition: submitPVResolutionJobs.py:251
Exception
Definition: hltDiff.cc:245
tensorflow::run
void run(Session *session, const NamedTensorList &inputs, const std::vector< std::string > &outputNames, std::vector< Tensor > *outputs, const thread::ThreadPoolOptions &threadPoolOptions)
Definition: TensorFlow.cc:213
DPFIsolation::GetNumberOfFeatures
static unsigned GetNumberOfFeatures(unsigned graphVersion)
Definition: DPFIsolation.cc:35
funct::abs
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
deep_tau::DeepTauCache
Definition: DeepTauBase.h:48
c
auto & c
Definition: CAHitNtupletGeneratorKernelsImpl.h:46
event
Definition: event.py:1
edm::Event
Definition: Event.h:73
edm::InputTag
Definition: InputTag.h:15
deep_tau::DeepTauCache::getGraph
const tensorflow::GraphDef & getGraph(const std::string &name="") const
Definition: DeepTauBase.h:58
DPFIsolation::graphVersion
unsigned graphVersion
Definition: DPFIsolation.cc:393
pwdgSkimBPark_cfi.vertices
vertices
Definition: pwdgSkimBPark_cfi.py:7
Generated for CMSSW Reference Manual by doxygen 1.8.16