CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
HLTDoubletDZ.cc
Go to the documentation of this file.
1 #include "FWCore/ParameterSet/interface/ParameterSet.h"
2 #include "HLTrigger/HLTfilters/interface/HLTDoubletDZ.h"
3 
4 #include "DataFormats/Common/interface/Handle.h"
5 
6 #include "DataFormats/Common/interface/Ref.h"
7 #include "DataFormats/HLTReco/interface/TriggerFilterObjectWithRefs.h"
8 
9 //#include "DataFormats/Candidate/interface/Particle.h"
10 #include "DataFormats/RecoCandidate/interface/RecoEcalCandidate.h"
11 #include "DataFormats/RecoCandidate/interface/RecoEcalCandidateFwd.h"
12 #include "DataFormats/EgammaCandidates/interface/Electron.h"
13 #include "DataFormats/EgammaCandidates/interface/ElectronFwd.h"
14 #include "DataFormats/RecoCandidate/interface/RecoChargedCandidate.h"
15 #include "DataFormats/RecoCandidate/interface/RecoChargedCandidateFwd.h"
16 #include "FWCore/MessageLogger/interface/MessageLogger.h"
17 
18 #include "FWCore/Framework/interface/MakerMacros.h"
19 
20 #include "DataFormats/Math/interface/deltaR.h"
21 
22 #include "HLTrigger/HLTcore/interface/defaultModuleLabel.h"
23 #include <cmath>
24 
25 //
26 // constructors and destructor
27 //
28 template<typename T1, typename T2>
29 HLTDoubletDZ<T1,T2>::HLTDoubletDZ(const edm::ParameterSet& iConfig) :
30  HLTFilter(iConfig),
31  originTag1_(iConfig.template getParameter<std::vector<edm::InputTag> >("originTag1")),
32  originTag2_(iConfig.template getParameter<std::vector<edm::InputTag> >("originTag2")),
33  inputTag1_(iConfig.template getParameter<edm::InputTag>("inputTag1")),
34  inputTag2_(iConfig.template getParameter<edm::InputTag>("inputTag2")),
35  inputToken1_(consumes<trigger::TriggerFilterObjectWithRefs>(inputTag1_)),
36  inputToken2_(consumes<trigger::TriggerFilterObjectWithRefs>(inputTag2_)),
37  //electronToken_ (consumes<reco::ElectronCollection>(iConfig.template getParameter<edm::InputTag>("electronTag"))),
38  triggerType1_(iConfig.template getParameter<int>("triggerType1")),
39  triggerType2_(iConfig.template getParameter<int>("triggerType2")),
40  minDR_ (iConfig.template getParameter<double>("MinDR")),
41  maxDZ_ (iConfig.template getParameter<double>("MaxDZ")),
42  minPixHitsForDZ_ (iConfig.template getParameter<int>("MinPixHitsForDZ")),
43  min_N_ (iConfig.template getParameter<int>("MinN")),
44  checkSC_ (iConfig.template getParameter<bool>("checkSC")),
45  same_ (inputTag1_.encode()==inputTag2_.encode()) // same collections to be compared?
46 {}
47 
48 template<>
49 HLTDoubletDZ<reco::RecoEcalCandidate, reco::RecoEcalCandidate>::HLTDoubletDZ(const edm::ParameterSet& iConfig) :
50  HLTFilter(iConfig),
51  originTag1_(iConfig.template getParameter<std::vector<edm::InputTag> >("originTag1")),
52  originTag2_(iConfig.template getParameter<std::vector<edm::InputTag> >("originTag2")),
53  inputTag1_(iConfig.template getParameter<edm::InputTag>("inputTag1")),
54  inputTag2_(iConfig.template getParameter<edm::InputTag>("inputTag2")),
55  inputToken1_(consumes<trigger::TriggerFilterObjectWithRefs>(inputTag1_)),
56  inputToken2_(consumes<trigger::TriggerFilterObjectWithRefs>(inputTag2_)),
57  electronToken_ (consumes<reco::ElectronCollection>(iConfig.template getParameter<edm::InputTag>("electronTag"))),
58  triggerType1_(iConfig.template getParameter<int>("triggerType1")),
59  triggerType2_(iConfig.template getParameter<int>("triggerType2")),
60  minDR_ (iConfig.template getParameter<double>("MinDR")),
61  maxDZ_ (iConfig.template getParameter<double>("MaxDZ")),
62  minPixHitsForDZ_ (iConfig.template getParameter<int>("MinPixHitsForDZ")),
63  min_N_ (iConfig.template getParameter<int>("MinN")),
64  checkSC_ (iConfig.template getParameter<bool>("checkSC")),
65  same_ (inputTag1_.encode()==inputTag2_.encode()) // same collections to be compared?
66 {}
67 
68 template<>
69 HLTDoubletDZ<reco::RecoChargedCandidate, reco::RecoEcalCandidate>::HLTDoubletDZ(const edm::ParameterSet& iConfig) :
70  HLTFilter(iConfig),
71  originTag1_(iConfig.template getParameter<std::vector<edm::InputTag> >("originTag1")),
72  originTag2_(iConfig.template getParameter<std::vector<edm::InputTag> >("originTag2")),
73  inputTag1_(iConfig.template getParameter<edm::InputTag>("inputTag1")),
74  inputTag2_(iConfig.template getParameter<edm::InputTag>("inputTag2")),
75  inputToken1_(consumes<trigger::TriggerFilterObjectWithRefs>(inputTag1_)),
76  inputToken2_(consumes<trigger::TriggerFilterObjectWithRefs>(inputTag2_)),
77  electronToken_ (consumes<reco::ElectronCollection>(iConfig.template getParameter<edm::InputTag>("electronTag"))),
78  triggerType1_(iConfig.template getParameter<int>("triggerType1")),
79  triggerType2_(iConfig.template getParameter<int>("triggerType2")),
80  minDR_ (iConfig.template getParameter<double>("MinDR")),
81  maxDZ_ (iConfig.template getParameter<double>("MaxDZ")),
82  minPixHitsForDZ_ (iConfig.template getParameter<int>("MinPixHitsForDZ")),
83  min_N_ (iConfig.template getParameter<int>("MinN")),
84  checkSC_ (iConfig.template getParameter<bool>("checkSC")),
85  same_ (inputTag1_.encode()==inputTag2_.encode()) // same collections to be compared?
86 {}
87 
88 template<>
89 HLTDoubletDZ<reco::RecoEcalCandidate, reco::RecoChargedCandidate>::HLTDoubletDZ(const edm::ParameterSet& iConfig) :
90  HLTFilter(iConfig),
91  originTag1_(iConfig.template getParameter<std::vector<edm::InputTag> >("originTag1")),
92  originTag2_(iConfig.template getParameter<std::vector<edm::InputTag> >("originTag2")),
93  inputTag1_(iConfig.template getParameter<edm::InputTag>("inputTag1")),
94  inputTag2_(iConfig.template getParameter<edm::InputTag>("inputTag2")),
95  inputToken1_(consumes<trigger::TriggerFilterObjectWithRefs>(inputTag1_)),
96  inputToken2_(consumes<trigger::TriggerFilterObjectWithRefs>(inputTag2_)),
97  electronToken_ (consumes<reco::ElectronCollection>(iConfig.template getParameter<edm::InputTag>("electronTag"))),
98  triggerType1_(iConfig.template getParameter<int>("triggerType1")),
99  triggerType2_(iConfig.template getParameter<int>("triggerType2")),
100  minDR_ (iConfig.template getParameter<double>("MinDR")),
101  maxDZ_ (iConfig.template getParameter<double>("MaxDZ")),
102  minPixHitsForDZ_ (iConfig.template getParameter<int>("MinPixHitsForDZ")),
103  min_N_ (iConfig.template getParameter<int>("MinN")),
104  checkSC_ (iConfig.template getParameter<bool>("checkSC")),
105  same_ (inputTag1_.encode()==inputTag2_.encode()) // same collections to be compared?
106 {}
107 
108 template<>
109 HLTDoubletDZ<reco::RecoChargedCandidate, reco::RecoChargedCandidate>::HLTDoubletDZ(const edm::ParameterSet& iConfig) :
110  HLTFilter(iConfig),
111  originTag1_(iConfig.template getParameter<std::vector<edm::InputTag> >("originTag1")),
112  originTag2_(iConfig.template getParameter<std::vector<edm::InputTag> >("originTag2")),
113  inputTag1_(iConfig.template getParameter<edm::InputTag>("inputTag1")),
114  inputTag2_(iConfig.template getParameter<edm::InputTag>("inputTag2")),
115  inputToken1_(consumes<trigger::TriggerFilterObjectWithRefs>(inputTag1_)),
116  inputToken2_(consumes<trigger::TriggerFilterObjectWithRefs>(inputTag2_)),
117  //electronToken_ (consumes<reco::ElectronCollection>(iConfig.template getParameter<edm::InputTag>("electronTag"))),
118  triggerType1_(iConfig.template getParameter<int>("triggerType1")),
119  triggerType2_(iConfig.template getParameter<int>("triggerType2")),
120  minDR_ (iConfig.template getParameter<double>("MinDR")),
121  maxDZ_ (iConfig.template getParameter<double>("MaxDZ")),
122  minPixHitsForDZ_ (iConfig.template getParameter<int>("MinPixHitsForDZ")),
123  min_N_ (iConfig.template getParameter<int>("MinN")),
124  checkSC_ (iConfig.template getParameter<bool>("checkSC")),
125  same_ (inputTag1_.encode()==inputTag2_.encode()) // same collections to be compared?
126 {}
127 
128 template<typename T1, typename T2>
129 HLTDoubletDZ<T1,T2>::~HLTDoubletDZ()
130 {}
131 
132 template<typename T1, typename T2>
133 void
134 HLTDoubletDZ<T1,T2>::fillDescriptions(edm::ConfigurationDescriptions& descriptions) {
135  edm::ParameterSetDescription desc;
136  makeHLTFilterDescription(desc);
137  std::vector<edm::InputTag> originTag1(1,edm::InputTag("hltOriginal1"));
138  std::vector<edm::InputTag> originTag2(1,edm::InputTag("hltOriginal2"));
139  desc.add<std::vector<edm::InputTag> >("originTag1",originTag1);
140  desc.add<std::vector<edm::InputTag> >("originTag2",originTag2);
141  desc.add<edm::InputTag>("inputTag1",edm::InputTag("hltFiltered1"));
142  desc.add<edm::InputTag>("inputTag2",edm::InputTag("hltFiltered2"));
143  desc.add<int>("triggerType1",0);
144  desc.add<int>("triggerType2",0);
145  desc.add<double>("MinDR",-1.0);
146  desc.add<double>("MaxDZ",0.2);
147  desc.add<int>("MinPixHitsForDZ",0);
148  desc.add<bool>("checkSC",false);
149  desc.add<int>("MinN",1);
150  descriptions.add(defaultModuleLabel<HLTDoubletDZ<T1,T2>>(), desc);
151 }
152 
153 template<>
154 void
155 HLTDoubletDZ<reco::RecoEcalCandidate, reco::RecoEcalCandidate>::fillDescriptions(edm::ConfigurationDescriptions& descriptions) {
156  edm::ParameterSetDescription desc;
157  makeHLTFilterDescription(desc);
158  std::vector<edm::InputTag> originTag1(1,edm::InputTag("hltOriginal1"));
159  std::vector<edm::InputTag> originTag2(1,edm::InputTag("hltOriginal2"));
160  desc.add<std::vector<edm::InputTag> >("originTag1",originTag1);
161  desc.add<std::vector<edm::InputTag> >("originTag2",originTag2);
162  desc.add<edm::InputTag>("inputTag1",edm::InputTag("hltFiltered1"));
163  desc.add<edm::InputTag>("inputTag2",edm::InputTag("hltFiltered2"));
164  desc.add<edm::InputTag>("electronTag",edm::InputTag("electronTag"));
165  desc.add<int>("triggerType1",0);
166  desc.add<int>("triggerType2",0);
167  desc.add<double>("MinDR",-1.0);
168  desc.add<double>("MaxDZ",0.2);
169  desc.add<int>("MinPixHitsForDZ",0);
170  desc.add<bool>("checkSC",false);
171  desc.add<int>("MinN",1);
172  descriptions.add(defaultModuleLabel<HLTDoubletDZ<reco::RecoEcalCandidate, reco::RecoEcalCandidate>>(), desc);
173 }
174 
175 template<>
176 void
177 HLTDoubletDZ<reco::RecoChargedCandidate, reco::RecoEcalCandidate>::fillDescriptions(edm::ConfigurationDescriptions& descriptions) {
178  edm::ParameterSetDescription desc;
179  makeHLTFilterDescription(desc);
180  std::vector<edm::InputTag> originTag1(1,edm::InputTag("hltOriginal1"));
181  std::vector<edm::InputTag> originTag2(1,edm::InputTag("hltOriginal2"));
182  desc.add<std::vector<edm::InputTag> >("originTag1",originTag1);
183  desc.add<std::vector<edm::InputTag> >("originTag2",originTag2);
184  desc.add<edm::InputTag>("inputTag1",edm::InputTag("hltFiltered1"));
185  desc.add<edm::InputTag>("inputTag2",edm::InputTag("hltFiltered2"));
186  desc.add<edm::InputTag>("electronTag",edm::InputTag("electronTag"));
187  desc.add<int>("triggerType1",0);
188  desc.add<int>("triggerType2",0);
189  desc.add<double>("MinDR",-1.0);
190  desc.add<double>("MaxDZ",0.2);
191  desc.add<int>("MinPixHitsForDZ",0);
192  desc.add<bool>("checkSC",false);
193  desc.add<int>("MinN",1);
194  descriptions.add(defaultModuleLabel<HLTDoubletDZ<reco::RecoChargedCandidate, reco::RecoEcalCandidate>>(), desc);
195 }
196 
197 template<>
198 void
199 HLTDoubletDZ<reco::RecoEcalCandidate, reco::RecoChargedCandidate>::fillDescriptions(edm::ConfigurationDescriptions& descriptions) {
200  edm::ParameterSetDescription desc;
201  makeHLTFilterDescription(desc);
202  std::vector<edm::InputTag> originTag1(1,edm::InputTag("hltOriginal1"));
203  std::vector<edm::InputTag> originTag2(1,edm::InputTag("hltOriginal2"));
204  desc.add<std::vector<edm::InputTag> >("originTag1",originTag1);
205  desc.add<std::vector<edm::InputTag> >("originTag2",originTag2);
206  desc.add<edm::InputTag>("inputTag1",edm::InputTag("hltFiltered1"));
207  desc.add<edm::InputTag>("inputTag2",edm::InputTag("hltFiltered2"));
208  desc.add<edm::InputTag>("electronTag",edm::InputTag("electronTag"));
209  desc.add<int>("triggerType1",0);
210  desc.add<int>("triggerType2",0);
211  desc.add<double>("MinDR",-1.0);
212  desc.add<double>("MaxDZ",0.2);
213  desc.add<int>("MinPixHitsForDZ",0);
214  desc.add<bool>("checkSC",false);
215  desc.add<int>("MinN",1);
216  descriptions.add(defaultModuleLabel<HLTDoubletDZ<reco::RecoEcalCandidate, reco::RecoChargedCandidate>>(), desc);
217 }
218 
219 template<>
220 void
221 HLTDoubletDZ<reco::RecoChargedCandidate, reco::RecoChargedCandidate>::fillDescriptions(edm::ConfigurationDescriptions& descriptions) {
222  edm::ParameterSetDescription desc;
223  makeHLTFilterDescription(desc);
224  std::vector<edm::InputTag> originTag1(1,edm::InputTag("hltOriginal1"));
225  std::vector<edm::InputTag> originTag2(1,edm::InputTag("hltOriginal2"));
226  desc.add<std::vector<edm::InputTag> >("originTag1",originTag1);
227  desc.add<std::vector<edm::InputTag> >("originTag2",originTag2);
228  desc.add<edm::InputTag>("inputTag1",edm::InputTag("hltFiltered1"));
229  desc.add<edm::InputTag>("inputTag2",edm::InputTag("hltFiltered2"));
230  desc.add<int>("triggerType1",0);
231  desc.add<int>("triggerType2",0);
232  desc.add<double>("MinDR",-1.0);
233  desc.add<double>("MaxDZ",0.2);
234  desc.add<int>("MinPixHitsForDZ",0);
235  desc.add<bool>("checkSC",false);
236  desc.add<int>("MinN",1);
237  descriptions.add(defaultModuleLabel<HLTDoubletDZ<reco::RecoChargedCandidate, reco::RecoChargedCandidate>>(), desc);
238 }
239 
240 template<typename T1, typename T2>
241 bool
242 //HLTDoubletDZ<T1, T2>::getCollections(edm::Event& iEvent, std::vector<T1Ref>& coll1, std::vector<T2Ref>& coll2) const {
243 HLTDoubletDZ<T1, T2>::getCollections(edm::Event& iEvent, std::vector<T1Ref>& coll1, std::vector<T2Ref>& coll2, trigger::TriggerFilterObjectWithRefs & filterproduct) const {
244 
245  edm::Handle<trigger::TriggerFilterObjectWithRefs> handle1,handle2;
246  if (iEvent.getByToken(inputToken1_, handle1) and iEvent.getByToken(inputToken2_, handle2)) {
247 
248  // get hold of pre-filtered object collections
249  handle1->getObjects(triggerType1_, coll1);
250  handle2->getObjects(triggerType2_, coll2);
251  const trigger::size_type n1(coll1.size());
252  const trigger::size_type n2(coll2.size());
253 
254  if (saveTags()) {
255  edm::InputTag tagOld;
256  for (unsigned int i=0; i<originTag1_.size(); ++i) {
257  filterproduct.addCollectionTag(originTag1_[i]);
258  }
259  tagOld=edm::InputTag();
260  for (trigger::size_type i1=0; i1!=n1; ++i1) {
261  const edm::ProductID pid(coll1[i1].id());
262  const std::string& label(iEvent.getProvenance(pid).moduleLabel());
263  const std::string& instance(iEvent.getProvenance(pid).productInstanceName());
264  const std::string& process(iEvent.getProvenance(pid).processName());
265  edm::InputTag tagNew(edm::InputTag(label,instance,process));
266  if (tagOld.encode()!=tagNew.encode()) {
267  filterproduct.addCollectionTag(tagNew);
268  tagOld=tagNew;
269  }
270  }
271  for (unsigned int i=0; i<originTag2_.size(); ++i) {
272  filterproduct.addCollectionTag(originTag2_[i]);
273  }
274  tagOld=edm::InputTag();
275  for (trigger::size_type i2=0; i2!=n2; ++i2) {
276  const edm::ProductID pid(coll2[i2].id());
277  const std::string& label(iEvent.getProvenance(pid).moduleLabel());
278  const std::string& instance(iEvent.getProvenance(pid).productInstanceName());
279  const std::string& process(iEvent.getProvenance(pid).processName());
280  edm::InputTag tagNew(edm::InputTag(label,instance,process));
281  if (tagOld.encode()!=tagNew.encode()) {
282  filterproduct.addCollectionTag(tagNew);
283  tagOld=tagNew;
284  }
285  }
286  }
287 
288  return true;
289  } else
290  return false;
291 }
292 
293 template<typename T1, typename T2>
294 bool
295 HLTDoubletDZ<T1,T2>::computeDZ(edm::Event& iEvent, T1Ref& r1, T2Ref& r2) const
296 {
297 
298  const reco::Candidate& candidate1(*r1);
299  const reco::Candidate& candidate2(*r2);
300  if ( reco::deltaR(candidate1, candidate2) < minDR_ )
301  return false;
302  if ( std::abs(candidate1.vz()-candidate2.vz()) > maxDZ_ )
303  return false;
304 
305  return true;
306 }
307 
308 template<>
309 bool
310 HLTDoubletDZ<reco::RecoEcalCandidate, reco::RecoChargedCandidate>::computeDZ(edm::Event& iEvent, T1Ref& r1, T2Ref& r2) const
311 {
312  edm::Handle<reco::ElectronCollection> electronHandle_;
313  iEvent.getByToken(electronToken_, electronHandle_);
314  if (!electronHandle_.isValid())
315  edm::LogError("HLTDoubletDZ") << "HLTDoubletDZ: Electron Handle not valid.";
316 
317  if ( reco::deltaR(*r1, *r2) < minDR_ )
318  return false;
319  reco::Electron e1;
320  for(reco::ElectronCollection::const_iterator eleIt = electronHandle_->begin(); eleIt != electronHandle_->end(); eleIt++) {
321  if (eleIt->superCluster() == r1->superCluster())
322  e1 = *(eleIt);
323  }
324 
325  const reco::RecoChargedCandidate& candidate2(*r2);
326  bool skipDZ = false;
327  if ( minPixHitsForDZ_ > 0 &&
328  (e1.gsfTrack()->hitPattern().numberOfValidPixelHits() < minPixHitsForDZ_
329  || candidate2.track()->hitPattern().numberOfValidPixelHits() < minPixHitsForDZ_) )
330  skipDZ = true;
331  if ( !skipDZ && std::abs(e1.vz()-candidate2.vz()) > maxDZ_ )
332  return false;
333 
334  return true;
335 }
336 
337 template<>
338 bool
339 HLTDoubletDZ<reco::RecoChargedCandidate, reco::RecoEcalCandidate>::computeDZ(edm::Event& iEvent, T1Ref& r1, T2Ref& r2) const
340 {
341  edm::Handle<reco::ElectronCollection> electronHandle_;
342  iEvent.getByToken(electronToken_, electronHandle_);
343  if (!electronHandle_.isValid())
344  edm::LogError("HLTDoubletDZ") << "HLTDoubletDZ: Electron Handle not valid.";
345 
346  if ( reco::deltaR(*r1, *r2) < minDR_ )
347  return false;
348  reco::Electron e2;
349  for(reco::ElectronCollection::const_iterator eleIt = electronHandle_->begin(); eleIt != electronHandle_->end(); eleIt++) {
350  if (eleIt->superCluster() == r2->superCluster())
351  e2 = *(eleIt);
352  }
353 
354  const reco::RecoChargedCandidate& candidate1(*r1);
355  bool skipDZ = false;
356  if ( minPixHitsForDZ_ > 0 &&
357  (candidate1.track()->hitPattern().numberOfValidPixelHits() < minPixHitsForDZ_
358  || e2.gsfTrack()->hitPattern().numberOfValidPixelHits() < minPixHitsForDZ_) )
359  skipDZ = true;
360  if ( !skipDZ && std::abs(e2.vz()-candidate1.vz()) > maxDZ_ )
361  return false;
362 
363  return true;
364 }
365 
366 template<>
367 bool
368 HLTDoubletDZ<reco::RecoEcalCandidate, reco::RecoEcalCandidate>::computeDZ(edm::Event& iEvent, T1Ref& r1, T2Ref& r2) const
369 {
370  edm::Handle<reco::ElectronCollection> electronHandle_;
371  iEvent.getByToken(electronToken_, electronHandle_);
372  if (!electronHandle_.isValid())
373  edm::LogError("HLTDoubletDZ") << "HLTDoubletDZ: Electron Handle not valid.";
374 
375  if ( reco::deltaR(*r1, *r2) < minDR_ )
376  return false;
377  reco::Electron e1, e2;
378  for(reco::ElectronCollection::const_iterator eleIt = electronHandle_->begin(); eleIt != electronHandle_->end(); eleIt++) {
379  if (eleIt->superCluster() == r2->superCluster())
380  e2 = *(eleIt);
381  if (eleIt->superCluster() == r1->superCluster())
382  e1 = *(eleIt);
383  }
384 
385  bool skipDZ = false;
386  if ( minPixHitsForDZ_ > 0 &&
387  (e1.gsfTrack()->hitPattern().numberOfValidPixelHits() < minPixHitsForDZ_
388  || e2.gsfTrack()->hitPattern().numberOfValidPixelHits() < minPixHitsForDZ_) )
389  skipDZ = true;
390  if ( !skipDZ && std::abs(e2.vz()-e1.vz()) > maxDZ_ )
391  return false;
392 
393  return true;
394 }
395 
396 template<>
397 bool
398 HLTDoubletDZ<reco::RecoChargedCandidate, reco::RecoChargedCandidate>::computeDZ(edm::Event& iEvent, T1Ref& r1, T2Ref& r2) const
399 {
400 
401  const reco::RecoChargedCandidate& candidate1(*r1);
402  const reco::RecoChargedCandidate& candidate2(*r2);
403  if ( reco::deltaR(candidate1, candidate2) < minDR_ )
404  return false;
405  bool skipDZ = false;
406  if ( minPixHitsForDZ_ > 0 &&
407  (candidate1.track()->hitPattern().numberOfValidPixelHits() < minPixHitsForDZ_
408  || candidate2.track()->hitPattern().numberOfValidPixelHits() < minPixHitsForDZ_) )
409  skipDZ = true;
410  if ( !skipDZ && std::abs(candidate1.vz()-candidate2.vz()) > maxDZ_ )
411  return false;
412 
413  return true;
414 }
415 
416 // ------------ method called to produce the data ------------
417 template<typename T1, typename T2>
418 bool
419 HLTDoubletDZ<T1,T2>::hltFilter(edm::Event& iEvent, const edm::EventSetup& iSetup, trigger::TriggerFilterObjectWithRefs & filterproduct) const
420 {
421 
422  // All HLT filters must create and fill an HLT filter object,
423  // recording any reconstructed physics objects satisfying (or not)
424  // this HLT filter, and place it in the Event.
425  bool accept(false);
426 
427  std::vector<T1Ref> coll1;
428  std::vector<T2Ref> coll2;
429 
430  if (getCollections(iEvent, coll1, coll2, filterproduct)) {
431  int n(0);
432  T1Ref r1;
433  T2Ref r2;
434 
435  for (unsigned int i1=0; i1!=coll1.size(); i1++) {
436  r1=coll1[i1];
437  //const reco::Candidate& candidate1(*r1);
438  unsigned int I(0);
439  if (same_) {I=i1+1;}
440  for (unsigned int i2=I; i2!=coll2.size(); i2++) {
441  r2=coll2[i2];
442  if (checkSC_) {
443  if (r1->superCluster().isNonnull() && r2->superCluster().isNonnull()) {
444  if (r1->superCluster() == r2->superCluster()) continue;
445  }
446  }
447 
448  if (!computeDZ(iEvent, r1, r2))
449  continue;
450 
451  n++;
452  filterproduct.addObject(triggerType1_,r1);
453  filterproduct.addObject(triggerType2_,r2);
454  }
455  }
456 
457  accept = accept || (n>=min_N_);
458  }
459 
460  return accept;
461 }
462 
463 typedef HLTDoubletDZ<reco::Electron , reco::Electron> HLT2ElectronElectronDZ;
464 typedef HLTDoubletDZ<reco::RecoChargedCandidate, reco::RecoChargedCandidate> HLT2MuonMuonDZ;
465 typedef HLTDoubletDZ<reco::Electron , reco::RecoChargedCandidate> HLT2ElectronMuonDZ;
466 typedef HLTDoubletDZ<reco::RecoEcalCandidate , reco::RecoEcalCandidate> HLT2PhotonPhotonDZ;
467 typedef HLTDoubletDZ<reco::RecoChargedCandidate, reco::RecoEcalCandidate> HLT2MuonPhotonDZ;
468 typedef HLTDoubletDZ<reco::RecoEcalCandidate , reco::RecoChargedCandidate> HLT2PhotonMuonDZ;
469 
470 DEFINE_FWK_MODULE(HLT2ElectronElectronDZ);
471 DEFINE_FWK_MODULE(HLT2MuonMuonDZ);
472 DEFINE_FWK_MODULE(HLT2ElectronMuonDZ);
473 DEFINE_FWK_MODULE(HLT2PhotonPhotonDZ);
474 DEFINE_FWK_MODULE(HLT2PhotonMuonDZ);
475 DEFINE_FWK_MODULE(HLT2MuonPhotonDZ);
defaultModuleLabel
std::string defaultModuleLabel()
Definition: defaultModuleLabel.h:16
i
int i
Definition: DBlmapReader.cc:9
HLTDoubletDZ::HLTDoubletDZ
HLTDoubletDZ(const edm::ParameterSet &)
Definition: HLTDoubletDZ.cc:29
CommPDSkim_cfg.saveTags
tuple saveTags
Definition: CommPDSkim_cfg.py:92
edm::Ref::isNonnull
bool isNonnull() const
Checks for non-null.
Definition: Ref.h:252
HLTDoubletDZ::fillDescriptions
static void fillDescriptions(edm::ConfigurationDescriptions &descriptions)
Definition: HLTDoubletDZ.cc:134
instance
static PFTauRenderPlugin instance
Definition: PFTauRenderPlugin.cc:70
edm::Event::getByToken
bool getByToken(EDGetToken token, Handle< PROD > &result) const
Definition: Event.h:462
DEFINE_FWK_MODULE
#define DEFINE_FWK_MODULE(type)
Definition: MakerMacros.h:17
reco::Electron::gsfTrack
virtual reco::GsfTrackRef gsfTrack() const
reference to a GsfTrack
HLT2PhotonMuonDZ
HLTDoubletDZ< reco::RecoEcalCandidate, reco::RecoChargedCandidate > HLT2PhotonMuonDZ
Definition: HLTDoubletDZ.cc:468
HLT2ElectronElectronDZ
HLTDoubletDZ< reco::Electron, reco::Electron > HLT2ElectronElectronDZ
Definition: HLTDoubletDZ.cc:463
dt_dqm_sourceclient_common_cff.reco
tuple reco
Definition: dt_dqm_sourceclient_common_cff.py:107
edm::Provenance::processName
std::string const & processName() const
Definition: Provenance.h:52
accept
bool accept(const edm::Event &event, const edm::TriggerResults &triggerTable, const std::string &triggerPath)
Definition: TopDQMHelpers.h:30
trigger::size_type
uint16_t size_type
Definition: TriggerTypeDefs.h:19
HLT_25ns10e33_v2_cff.originTag1
tuple originTag1
Definition: HLT_25ns10e33_v2_cff.py:21182
edm::InputTag::encode
std::string encode() const
Definition: InputTag.cc:164
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:255
HLTDoubletDZ::getCollections
bool getCollections(edm::Event &iEvent, std::vector< T1Ref > &coll1, std::vector< T2Ref > &coll2, trigger::TriggerFilterObjectWithRefs &filterproduct) const
Definition: HLTDoubletDZ.cc:243
trigger::TriggerRefsCollections::addObject
void addObject(int id, const reco::RecoEcalCandidateRef &ref)
setters for L3 collections: (id=physics type, and Ref&lt;C&gt;)
Definition: TriggerRefsCollections.h:217
HLT_25ns10e33_v2_cff.InputTag
tuple InputTag
Definition: HLT_25ns10e33_v2_cff.py:58467
iEvent
int iEvent
Definition: GenABIO.cc:230
HLTDoubletDZ::computeDZ
bool computeDZ(edm::Event &iEvent, T1Ref &c1, T2Ref &c2) const
Definition: HLTDoubletDZ.cc:295
funct::abs
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
reco::deltaR
auto deltaR(const T1 &t1, const T2 &t2) -> decltype(t1.eta())
Definition: deltaR.h:28
Exhume::I
const std::complex< double > I
Definition: I.h:8
edm::ParameterSetDescription::add
ParameterDescriptionBase * add(U const &iLabel, T const &value)
Definition: ParameterSetDescription.h:89
edm::HandleBase::isValid
bool isValid() const
Definition: HandleBase.h:75
reco::LeafCandidate::vz
virtual double vz() const
z coordinate of vertex position
Definition: LeafCandidate.h:156
reco::ElectronCollection
std::vector< Electron > ElectronCollection
collectin of Electron objects
Definition: ElectronFwd.h:9
HLTDoubletDZ::hltFilter
virtual bool hltFilter(edm::Event &, const edm::EventSetup &, trigger::TriggerFilterObjectWithRefs &filterproduct) const override
Definition: HLTDoubletDZ.cc:419
HLT2MuonMuonDZ
HLTDoubletDZ< reco::RecoChargedCandidate, reco::RecoChargedCandidate > HLT2MuonMuonDZ
Definition: HLTDoubletDZ.cc:464
HLT_25ns10e33_v2_cff.originTag2
tuple originTag2
Definition: HLT_25ns10e33_v2_cff.py:21178
HLT2ElectronMuonDZ
HLTDoubletDZ< reco::Electron, reco::RecoChargedCandidate > HLT2ElectronMuonDZ
Definition: HLTDoubletDZ.cc:465
reco::e1
Float e1
Definition: deltaR.h:20
HLT2PhotonPhotonDZ
HLTDoubletDZ< reco::RecoEcalCandidate, reco::RecoEcalCandidate > HLT2PhotonPhotonDZ
Definition: HLTDoubletDZ.cc:466
sysUtil.pid
tuple pid
Definition: sysUtil.py:22
trigger::TriggerFilterObjectWithRefs::addCollectionTag
void addCollectionTag(const edm::InputTag &collectionTag)
collectionTags
Definition: TriggerFilterObjectWithRefs.h:64
edm::ConfigurationDescriptions::add
void add(std::string const &label, ParameterSetDescription const &psetDescription)
Definition: ConfigurationDescriptions.cc:65
edm::Provenance::moduleLabel
std::string const & moduleLabel() const
Definition: Provenance.h:51
reco::RecoChargedCandidate::track
virtual reco::TrackRef track() const
reference to a track
Definition: RecoChargedCandidate.cc:11
HLT2MuonPhotonDZ
HLTDoubletDZ< reco::RecoChargedCandidate, reco::RecoEcalCandidate > HLT2MuonPhotonDZ
Definition: HLTDoubletDZ.cc:467
reco::e2
Float e2
Definition: deltaR.h:21
reco::Candidate::vz
virtual double vz() const =0
z coordinate of vertex position
edm::Provenance::productInstanceName
std::string const & productInstanceName() const
Definition: Provenance.h:53
LaserDQM_cfg.process
tuple process
Definition: LaserDQM_cfg.py:3
edm::Event::getProvenance
Provenance getProvenance(BranchID const &theID) const
Definition: Event.cc:80
svgfig.template
def template
Definition: svgfig.py:520
Generated for CMSSW Reference Manual by doxygen 1.8.5