CMS 3D CMS Logo
L1GTProducer.cc
Go to the documentation of this file.
1 #include "FWCore/ParameterSet/interface/ParameterSetDescription.h"
2 #include "DataFormats/Common/interface/Handle.h"
3 #include "FWCore/Utilities/interface/EDGetToken.h"
4 
5 #include "FWCore/Framework/interface/MakerMacros.h"
6 #include "DataFormats/Common/interface/View.h"
7 
8 #include "FWCore/Framework/interface/Frameworkfwd.h"
9 #include "FWCore/Framework/interface/global/EDProducer.h"
10 
11 #include "FWCore/Framework/interface/Event.h"
12 #include "FWCore/ParameterSet/interface/ParameterSet.h"
13 #include "DataFormats/L1Trigger/interface/P2GTCandidate.h"
14 #include "FWCore/ParameterSet/interface/ConfigurationDescriptions.h"
15 #include "L1Trigger/Phase2L1GT/interface/L1GTScales.h"
16 
17 #include "DataFormats/L1Trigger/interface/TkJetWord.h"
18 #include "DataFormats/L1Trigger/interface/VertexWord.h"
19 
20 #include "DataFormats/L1TMuonPhase2/interface/SAMuon.h"
21 #include "DataFormats/L1TMuonPhase2/interface/TrackerMuon.h"
22 
23 #include "DataFormats/L1TParticleFlow/interface/PFJet.h"
24 #include "DataFormats/L1TCorrelator/interface/TkEmFwd.h"
25 #include "DataFormats/L1TCorrelator/interface/TkEm.h"
26 #include "DataFormats/L1TCorrelator/interface/TkElectronFwd.h"
27 #include "DataFormats/L1TCorrelator/interface/TkElectron.h"
28 #include "DataFormats/L1TParticleFlow/interface/PFTau.h"
29 #include "DataFormats/L1TParticleFlow/interface/gt_datatypes.h"
30 
31 #include "DataFormats/L1Trigger/interface/EtSum.h"
32 
33 #include <vector>
34 #include <array>
35 #include <string>
36 #include <type_traits>
37 
38 namespace l1t {
39 
40  class L1GTProducer : public edm::global::EDProducer<> {
41  public:
42  explicit L1GTProducer(const edm::ParameterSet &);
43  ~L1GTProducer() override = default;
44 
45  static void fillDescriptions(edm::ConfigurationDescriptions &);
46 
47  private:
48  void produceGTTPromptJets(edm::Event &event) const;
49  void produceGTTDisplacedJets(edm::Event &event) const;
50  void produceGTTPrimaryVert(edm::Event &event) const;
51 
52  void produceGMTSaPromptMuons(edm::Event &event) const;
53  void produceGMTSaDisplacedMuons(edm::Event &event) const;
54  void produceGMTTkMuons(edm::Event &event) const;
55 
56  void produceCL2Jets(edm::Event &event) const;
57  void produceCL2Photons(edm::Event &event) const;
58  void produceCL2Electrons(edm::Event &event) const;
59  void produceCL2Taus(edm::Event &event) const;
60  void produceCL2EtSum(edm::Event &event) const;
61  void produceCl2HtSum(edm::Event &event) const;
62 
63  void produce(edm::StreamID, edm::Event &, const edm::EventSetup &) const override;
64 
65  const L1GTScales scales_;
66 
67  const edm::EDGetTokenT<TkJetWordCollection> gttPromptJetToken_;
68  const edm::EDGetTokenT<TkJetWordCollection> gttDisplacedJetToken_;
69  const edm::EDGetTokenT<VertexWordCollection> gttPrimaryVertexToken_;
70 
71  const edm::EDGetTokenT<SAMuonCollection> gmtSaPromptMuonToken_;
72  const edm::EDGetTokenT<SAMuonCollection> gmtSaDisplacedMuonToken_;
73  const edm::EDGetTokenT<TrackerMuonCollection> gmtTkMuonToken_;
74 
75  const edm::EDGetTokenT<PFJetCollection> cl2JetToken_;
76  const edm::EDGetTokenT<TkEmCollection> cl2PhotonToken_;
77  const edm::EDGetTokenT<TkElectronCollection> cl2ElectronToken_;
78  const edm::EDGetTokenT<PFTauCollection> cl2TauToken_;
79  const edm::EDGetTokenT<std::vector<l1t::EtSum>> cl2EtSumToken_;
80  const edm::EDGetTokenT<std::vector<l1t::EtSum>> cl2HtSumToken_;
81  };
82 
83  L1GTProducer::L1GTProducer(const edm::ParameterSet &config)
84  : scales_(config.getParameter<edm::ParameterSet>("scales")),
85  gttPromptJetToken_(consumes<TkJetWordCollection>(config.getParameter<edm::InputTag>("GTTPromptJets"))),
86  gttDisplacedJetToken_(consumes<TkJetWordCollection>(config.getParameter<edm::InputTag>("GTTDisplacedJets"))),
87  gttPrimaryVertexToken_(consumes<VertexWordCollection>(config.getParameter<edm::InputTag>("GTTPrimaryVert"))),
88  gmtSaPromptMuonToken_(consumes<SAMuonCollection>(config.getParameter<edm::InputTag>("GMTSaPromptMuons"))),
89  gmtSaDisplacedMuonToken_(consumes<SAMuonCollection>(config.getParameter<edm::InputTag>("GMTSaDisplacedMuons"))),
90  gmtTkMuonToken_(consumes<TrackerMuonCollection>(config.getParameter<edm::InputTag>("GMTTkMuons"))),
91  cl2JetToken_(consumes<PFJetCollection>(config.getParameter<edm::InputTag>("CL2Jets"))),
92  cl2PhotonToken_(consumes<TkEmCollection>(config.getParameter<edm::InputTag>("CL2Photons"))),
93  cl2ElectronToken_(consumes<TkElectronCollection>(config.getParameter<edm::InputTag>("CL2Electrons"))),
94  cl2TauToken_(consumes<PFTauCollection>(config.getParameter<edm::InputTag>("CL2Taus"))),
95  cl2EtSumToken_(consumes<std::vector<l1t::EtSum>>(config.getParameter<edm::InputTag>("CL2EtSum"))),
96  cl2HtSumToken_(consumes<std::vector<l1t::EtSum>>(config.getParameter<edm::InputTag>("CL2HtSum"))) {
97  produces<P2GTCandidateCollection>("GTTPromptJets");
98  produces<P2GTCandidateCollection>("GTTDisplacedJets");
99  produces<P2GTCandidateCollection>("GTTPrimaryVert");
100 
101  produces<P2GTCandidateCollection>("GMTSaPromptMuons");
102  produces<P2GTCandidateCollection>("GMTSaDisplacedMuons");
103  produces<P2GTCandidateCollection>("GMTTkMuons");
104 
105  produces<P2GTCandidateCollection>("CL2Jets");
106  produces<P2GTCandidateCollection>("CL2Photons");
107  produces<P2GTCandidateCollection>("CL2Electrons");
108  produces<P2GTCandidateCollection>("CL2Taus");
109  produces<P2GTCandidateCollection>("CL2EtSum");
110  produces<P2GTCandidateCollection>("CL2HtSum");
111  }
112 
113  void L1GTProducer::fillDescriptions(edm::ConfigurationDescriptions &description) {
114  edm::ParameterSetDescription desc;
115 
116  edm::ParameterSetDescription scalesDesc;
117  L1GTScales::fillPSetDescription(scalesDesc);
118  desc.add<edm::ParameterSetDescription>("scales", scalesDesc);
119 
120  desc.add<edm::InputTag>("GTTPromptJets");
121  desc.add<edm::InputTag>("GTTDisplacedJets");
122  desc.add<edm::InputTag>("GTTPrimaryVert");
123 
124  desc.add<edm::InputTag>("GMTSaPromptMuons");
125  desc.add<edm::InputTag>("GMTSaDisplacedMuons");
126  desc.add<edm::InputTag>("GMTTkMuons");
127 
128  desc.add<edm::InputTag>("CL2Jets");
129  desc.add<edm::InputTag>("CL2Photons");
130  desc.add<edm::InputTag>("CL2Electrons");
131  desc.add<edm::InputTag>("CL2Taus");
132  desc.add<edm::InputTag>("CL2EtSum");
133  desc.add<edm::InputTag>("CL2HtSum");
134 
135  description.addWithDefaultLabel(desc);
136  }
137 
138  void L1GTProducer::produceGTTPrimaryVert(edm::Event &event) const {
139  std::unique_ptr<P2GTCandidateCollection> outputCollection = std::make_unique<P2GTCandidateCollection>();
140  const VertexWordCollection &collection = event.get(gttPrimaryVertexToken_);
141  for (std::size_t i = 0; i < collection.size() && i < 10; i++) {
142  const VertexWord &obj = collection[i];
143  int hwZ0 = obj.z0Word().V.to_int() * 5;
144  P2GTCandidate gtObj(
145  0, reco::ParticleState::PolarLorentzVector(), reco::ParticleState::Point(0, 0, scales_.to_z0(hwZ0)));
146  gtObj.hwZ0_ = hwZ0;
147  gtObj.hwQual_ = obj.qualityWord().V.to_int();
148  gtObj.hwSum_pT_pv_ = obj.multiplicityWord().V.to_int();
149  gtObj.hwNumber_of_tracks_in_pv_ = obj.multiplicityWord().V.to_int();
150  gtObj.hwNumber_of_tracks_not_in_pv_ = obj.inverseMultiplicityWord().V.to_int();
151  gtObj.objectType_ = P2GTCandidate::GTTPrimaryVert;
152 
153  outputCollection->push_back(gtObj);
154  }
155  event.put(std::move(outputCollection), "GTTPrimaryVert");
156  }
157 
158  void L1GTProducer::produceGTTPromptJets(edm::Event &event) const {
159  std::unique_ptr<P2GTCandidateCollection> outputCollection = std::make_unique<P2GTCandidateCollection>();
160  const TkJetWordCollection &collection = event.get(gttPromptJetToken_);
161  for (std::size_t i = 0; i < collection.size() && i < 12; i++) {
162  const TkJetWord &obj = collection[i];
163  int hwZ0 = obj.z0Word().V.to_int() << 7;
164  P2GTCandidate gtObj(0,
165  reco::ParticleState::PolarLorentzVector(scales_.to_pT(obj.ptWord().V.to_int()),
166  scales_.to_eta(obj.glbEtaWord().V.to_int()),
167  scales_.to_phi(obj.glbPhiWord().V.to_int()),
168  0),
169  reco::ParticleState::Point(0, 0, scales_.to_z0(hwZ0)));
170  gtObj.hwPT_ = obj.ptWord().V.to_int();
171  gtObj.hwPhi_ = obj.glbPhiWord().V.to_int();
172  gtObj.hwEta_ = obj.glbEtaWord().V.to_int();
173  gtObj.hwZ0_ = hwZ0;
174  gtObj.hwNumber_of_tracks_ = obj.ntWord().V.to_int();
175  gtObj.objectType_ = P2GTCandidate::GTTPromptJets;
176 
177  outputCollection->push_back(gtObj);
178  }
179  event.put(std::move(outputCollection), "GTTPromptJets");
180  }
181 
182  void L1GTProducer::produceGTTDisplacedJets(edm::Event &event) const {
183  std::unique_ptr<P2GTCandidateCollection> outputCollection = std::make_unique<P2GTCandidateCollection>();
184  const TkJetWordCollection &collection = event.get(gttDisplacedJetToken_);
185  for (std::size_t i = 0; i < collection.size() && i < 12; i++) {
186  const TkJetWord &obj = collection[i];
187  int hwZ0 = obj.z0Word().V.to_int() << 7;
188  P2GTCandidate gtObj(0,
189  reco::ParticleState::PolarLorentzVector(scales_.to_pT(obj.ptWord().V.to_int()),
190  scales_.to_eta(obj.glbEtaWord().V.to_int()),
191  scales_.to_phi(obj.glbPhiWord().V.to_int()),
192  0),
193  reco::ParticleState::Point(0, 0, scales_.to_z0(hwZ0)));
194  gtObj.hwPT_ = obj.ptWord().V.to_int();
195  gtObj.hwPhi_ = obj.glbPhiWord().V.to_int();
196  gtObj.hwEta_ = obj.glbEtaWord().V.to_int();
197  gtObj.hwZ0_ = hwZ0;
198  gtObj.hwNumber_of_tracks_ = obj.ntWord().V.to_int();
199  gtObj.objectType_ = P2GTCandidate::GTTDisplacedJets;
200 
201  outputCollection->push_back(gtObj);
202  }
203  event.put(std::move(outputCollection), "GTTDisplacedJets");
204  }
205 
206  void L1GTProducer::produceGMTSaPromptMuons(edm::Event &event) const {
207  std::unique_ptr<P2GTCandidateCollection> outputCollection = std::make_unique<P2GTCandidateCollection>();
208  const SAMuonCollection &collection = event.get(gmtSaPromptMuonToken_);
209  for (size_t i = 0; i < collection.size() && i < 12; i++) {
210  const SAMuon &obj = collection[i];
211  int hwZ0 = obj.apZ0().to_int() << 12;
212  P2GTCandidate gtObj(scales_.to_chg(obj.apCharge().to_int()),
213  reco::ParticleState::PolarLorentzVector(scales_.to_pT(obj.apPt().to_int()),
214  scales_.to_eta(obj.apEta().to_int()),
215  scales_.to_phi(obj.apPhi().to_int()),
216  0),
217  reco::ParticleState::Point(0, 0, scales_.to_z0(hwZ0)));
218  gtObj.hwPT_ = obj.apPt().to_int();
219  gtObj.hwPhi_ = obj.apPhi().to_int();
220  gtObj.hwEta_ = obj.apEta().to_int();
221  gtObj.hwZ0_ = hwZ0;
222  gtObj.hwQual_ = obj.apQual().to_int();
223  gtObj.hwCharge_ = obj.apCharge().to_int();
224  gtObj.hwD0_ = obj.apD0().to_int();
225  gtObj.objectType_ = P2GTCandidate::GMTSaPromptMuons;
226 
227  outputCollection->push_back(gtObj);
228  }
229  event.put(std::move(outputCollection), "GMTSaPromptMuons");
230  }
231 
232  void L1GTProducer::produceGMTSaDisplacedMuons(edm::Event &event) const {
233  std::unique_ptr<P2GTCandidateCollection> outputCollection = std::make_unique<P2GTCandidateCollection>();
234  const SAMuonCollection &collection = event.get(gmtSaDisplacedMuonToken_);
235  for (size_t i = 0; i < collection.size() && i < 12; i++) {
236  const SAMuon &obj = collection[i];
237  int hwZ0 = obj.apZ0().to_int() << 12;
238  P2GTCandidate gtObj(scales_.to_chg(obj.apCharge().to_int()),
239  reco::ParticleState::PolarLorentzVector(scales_.to_pT(obj.apPt().to_int()),
240  scales_.to_eta(obj.apEta().to_int()),
241  scales_.to_phi(obj.apPhi().to_int()),
242  0),
243  reco::ParticleState::Point(0, 0, scales_.to_z0(hwZ0)));
244  gtObj.hwPT_ = obj.apPt().to_int();
245  gtObj.hwPhi_ = obj.apPhi().to_int();
246  gtObj.hwEta_ = obj.apEta().to_int();
247  gtObj.hwZ0_ = hwZ0;
248  gtObj.hwQual_ = obj.apQual().to_int();
249  gtObj.hwCharge_ = obj.apCharge().to_int();
250  gtObj.hwD0_ = obj.apD0().to_int();
251  gtObj.objectType_ = P2GTCandidate::GMTSaDisplacedMuons;
252 
253  outputCollection->push_back(gtObj);
254  }
255  event.put(std::move(outputCollection), "GMTSaDisplacedMuons");
256  }
257 
258  void L1GTProducer::produceGMTTkMuons(edm::Event &event) const {
259  std::unique_ptr<P2GTCandidateCollection> outputCollection = std::make_unique<P2GTCandidateCollection>();
260  const TrackerMuonCollection &collection = event.get(gmtTkMuonToken_);
261  for (size_t i = 0; i < collection.size() && i < 12; i++) {
262  const TrackerMuon &obj = collection[i];
263  int hwZ0 = obj.apZ0().to_int() << 7;
264  P2GTCandidate gtObj(scales_.to_chg(obj.apCharge().to_int()),
265  reco::ParticleState::PolarLorentzVector(scales_.to_pT(obj.apPt().to_int()),
266  scales_.to_eta(obj.apEta().to_int()),
267  scales_.to_phi(obj.apPhi().to_int()),
268  0),
269  reco::ParticleState::Point(0, 0, scales_.to_z0(hwZ0)));
270  gtObj.hwPT_ = obj.apPt().to_int();
271  gtObj.hwPhi_ = obj.apPhi().to_int();
272  gtObj.hwEta_ = obj.apEta().to_int();
273  gtObj.hwZ0_ = hwZ0;
274  gtObj.hwIso_ = obj.apIso().to_int();
275  gtObj.hwQual_ = obj.apQual().to_int();
276  gtObj.hwCharge_ = obj.apCharge().to_int();
277  gtObj.hwD0_ = obj.apD0().to_int();
278  gtObj.hwBeta_ = obj.apBeta().to_int();
279  gtObj.objectType_ = P2GTCandidate::GMTTkMuons;
280 
281  outputCollection->push_back(gtObj);
282  }
283  event.put(std::move(outputCollection), "GMTTkMuons");
284  }
285 
286  void L1GTProducer::produceCL2Jets(edm::Event &event) const {
287  std::unique_ptr<P2GTCandidateCollection> outputCollection = std::make_unique<P2GTCandidateCollection>();
288  const PFJetCollection &collection = event.get(cl2JetToken_);
289  for (size_t i = 0; i < collection.size() && i < 12; i++) {
290  l1gt::Jet gtJet = l1gt::Jet::unpack(collection[i].getHWJetGT());
291  int hwZ0 = gtJet.z0.V.to_int() << 7;
292  P2GTCandidate gtObj(0,
293  reco::ParticleState::PolarLorentzVector(scales_.to_pT(gtJet.v3.pt.V.to_int()),
294  scales_.to_eta(gtJet.v3.eta.V.to_int()),
295  scales_.to_phi(gtJet.v3.phi.V.to_int()),
296  0),
297  reco::ParticleState::Point(0, 0, scales_.to_z0(hwZ0)));
298  gtObj.hwPT_ = gtJet.v3.pt.V.to_int();
299  gtObj.hwPhi_ = gtJet.v3.phi.V.to_int();
300  gtObj.hwEta_ = gtJet.v3.eta.V.to_int();
301  gtObj.hwZ0_ = hwZ0;
302  gtObj.objectType_ = P2GTCandidate::CL2Jets;
303 
304  outputCollection->push_back(gtObj);
305  }
306  event.put(std::move(outputCollection), "CL2Jets");
307  }
308 
309  void L1GTProducer::produceCL2Photons(edm::Event &event) const {
310  std::unique_ptr<P2GTCandidateCollection> outputCollection = std::make_unique<P2GTCandidateCollection>();
311  const TkEmCollection &collection = event.get(cl2PhotonToken_);
312  for (size_t i = 0; i < collection.size() && i < 12; i++) {
313  l1gt::Photon gtPhoton = collection[i].hwObj();
314  P2GTCandidate gtObj(0,
315  reco::ParticleState::PolarLorentzVector(scales_.to_pT(gtPhoton.v3.pt.V.to_int()),
316  scales_.to_eta(gtPhoton.v3.eta.V.to_int()),
317  scales_.to_phi(gtPhoton.v3.phi.V.to_int()),
318  0));
319  gtObj.hwPT_ = gtPhoton.v3.pt.V.to_int();
320  gtObj.hwPhi_ = gtPhoton.v3.phi.V.to_int();
321  gtObj.hwEta_ = gtPhoton.v3.eta.V.to_int();
322  gtObj.hwIso_ = gtPhoton.isolationPT.V.to_int();
323  gtObj.hwQual_ = gtPhoton.qualityFlags.V.to_int();
324  gtObj.objectType_ = P2GTCandidate::CL2Photons;
325 
326  outputCollection->push_back(gtObj);
327  }
328  event.put(std::move(outputCollection), "CL2Photons");
329  }
330 
331  void L1GTProducer::produceCL2Electrons(edm::Event &event) const {
332  std::unique_ptr<P2GTCandidateCollection> outputCollection = std::make_unique<P2GTCandidateCollection>();
333  const TkElectronCollection &collection = event.get(cl2ElectronToken_);
334  for (size_t i = 0; i < collection.size() && i < 12; i++) {
335  l1gt::Electron gtElectron = collection[i].hwObj();
336  int hwZ0 = gtElectron.z0.V.to_int() << 7;
337  P2GTCandidate gtObj(scales_.to_chg(gtElectron.charge.V.to_int()),
338  reco::ParticleState::PolarLorentzVector(scales_.to_pT(gtElectron.v3.pt.V.to_int()),
339  scales_.to_eta(gtElectron.v3.eta.V.to_int()),
340  scales_.to_phi(gtElectron.v3.phi.V.to_int()),
341  0),
342  reco::ParticleState::Point(0, 0, scales_.to_z0(hwZ0)));
343  gtObj.hwPT_ = gtElectron.v3.pt.V.to_int();
344  gtObj.hwPhi_ = gtElectron.v3.phi.V.to_int();
345  gtObj.hwEta_ = gtElectron.v3.eta.V.to_int();
346  gtObj.hwZ0_ = hwZ0;
347  gtObj.hwIso_ = gtElectron.isolationPT.V.to_int();
348  gtObj.hwQual_ = gtElectron.qualityFlags.V.to_int();
349  gtObj.hwCharge_ = gtElectron.charge.V.to_int();
350  gtObj.objectType_ = P2GTCandidate::CL2Electrons;
351 
352  outputCollection->push_back(gtObj);
353  }
354  event.put(std::move(outputCollection), "CL2Electrons");
355  }
356 
357  void L1GTProducer::produceCL2Taus(edm::Event &event) const {
358  std::unique_ptr<P2GTCandidateCollection> outputCollection = std::make_unique<P2GTCandidateCollection>();
359  const PFTauCollection &collection = event.get(cl2TauToken_);
360  for (size_t i = 0; i < collection.size() && i < 12; i++) {
361  l1gt::Tau gtTau = collection[i].getHWTauGT();
362  P2GTCandidate gtObj(scales_.to_chg(gtTau.charge.V.to_int()),
363  reco::ParticleState::PolarLorentzVector(scales_.to_pT(gtTau.v3.pt.V.to_int()),
364  scales_.to_eta(gtTau.v3.eta.V.to_int()),
365  scales_.to_phi(gtTau.v3.phi.V.to_int()),
366  0));
367  gtObj.hwPT_ = gtTau.v3.pt.V.to_int();
368  gtObj.hwPhi_ = gtTau.v3.phi.V.to_int();
369  gtObj.hwEta_ = gtTau.v3.eta.V.to_int();
370  gtObj.hwSeed_pT_ = gtTau.seed_pt.V.to_int();
371  gtObj.hwSeed_z0_ = gtTau.seed_z0.V.to_int();
372  gtObj.hwCharge_ = gtTau.charge.V.to_int();
373  gtObj.hwType_ = gtTau.type.V.to_int();
374  gtObj.hwIso_ = gtTau.quality.V.to_int();
375  gtObj.objectType_ = P2GTCandidate::CL2Taus;
376 
377  outputCollection->push_back(gtObj);
378  }
379  event.put(std::move(outputCollection), "CL2Taus");
380  }
381 
382  void L1GTProducer::produceCL2EtSum(edm::Event &event) const {
383  std::unique_ptr<P2GTCandidateCollection> outputCollection = std::make_unique<P2GTCandidateCollection>();
384  const std::vector<EtSum> &collection = event.get(cl2EtSumToken_);
385  const EtSum &met = collection[0];
386 
387  l1gt::Sum sum{true /* valid */, met.pt(), met.phi() / l1gt::Scales::ETAPHI_LSB, 0 /* scalar sum */};
388 
389  P2GTCandidate gtObj(0,
390  reco::ParticleState::PolarLorentzVector(
391  scales_.to_pT(sum.vector_pt.V.to_int()), 0, scales_.to_phi(sum.vector_phi.V.to_int()), 0));
392  gtObj.hwPT_ = sum.vector_pt.V.to_int();
393  gtObj.hwPhi_ = sum.vector_phi.V.to_int();
394  gtObj.hwSca_sum_ = sum.scalar_pt.V.to_int();
395  gtObj.objectType_ = P2GTCandidate::CL2EtSum;
396 
397  outputCollection->push_back(gtObj);
398  event.put(std::move(outputCollection), "CL2EtSum");
399  }
400 
401  void L1GTProducer::produceCl2HtSum(edm::Event &event) const {
402  std::unique_ptr<P2GTCandidateCollection> outputCollection = std::make_unique<P2GTCandidateCollection>();
403  const std::vector<EtSum> &collection = event.get(cl2HtSumToken_);
404  const EtSum &ht = collection[0];
405  const EtSum &mht = collection[1];
406 
407  P2GTCandidate gtObj(
408  0, reco::ParticleState::PolarLorentzVector(scales_.to_pT(mht.hwPt()), 0, scales_.to_phi(mht.hwPhi()), 0));
409  gtObj.hwPT_ = mht.hwPt();
410  gtObj.hwPhi_ = mht.hwPhi();
411  gtObj.hwSca_sum_ = ht.hwPt();
412  gtObj.objectType_ = P2GTCandidate::CL2HtSum;
413 
414  outputCollection->push_back(gtObj);
415  event.put(std::move(outputCollection), "CL2HtSum");
416  }
417 
418  void L1GTProducer::produce(edm::StreamID, edm::Event &event, const edm::EventSetup &setup) const {
419  produceGTTPromptJets(event);
420  produceGTTDisplacedJets(event);
421  produceGTTPrimaryVert(event);
422 
423  produceGMTSaPromptMuons(event);
424  produceGMTSaDisplacedMuons(event);
425  produceGMTTkMuons(event);
426 
427  produceCL2Jets(event);
428  produceCL2Photons(event);
429  produceCL2Electrons(event);
430  produceCL2Taus(event);
431  produceCL2EtSum(event);
432  produceCl2HtSum(event);
433  }
434 } // namespace l1t
435 
436 using namespace l1t;
437 
438 DEFINE_FWK_MODULE(L1GTProducer);
l1t::L1GTProducer::~L1GTProducer
~L1GTProducer() override=default
l1gt::Photon::qualityFlags
egquality_t qualityFlags
Definition: gt_datatypes.h:316
l1t::P2GTCandidate::hwIso_
Optional< int > hwIso_
Definition: P2GTCandidate.h:307
ProducerED_cfi.InputTag
InputTag
Definition: ProducerED_cfi.py:5
l1t::L1GTProducer::gttPrimaryVertexToken_
const edm::EDGetTokenT< VertexWordCollection > gttPrimaryVertexToken_
Definition: L1GTProducer.cc:69
l1t::TkJetWordCollection
std::vector< l1t::TkJetWord > TkJetWordCollection
Definition: TkJetWord.h:191
l1t::TkElectronCollection
std::vector< TkElectron > TkElectronCollection
Definition: TkElectronFwd.h:16
l1t::PFTauCollection
std::vector< l1t::PFTau > PFTauCollection
Definition: PFTau.h:90
mps_fire.i
i
Definition: mps_fire.py:429
gt_datatypes.h
l1t::L1GTProducer::cl2ElectronToken_
const edm::EDGetTokenT< TkElectronCollection > cl2ElectronToken_
Definition: L1GTProducer.cc:77
l1t::L1GTProducer::L1GTProducer
L1GTProducer(const edm::ParameterSet &)
Definition: L1GTProducer.cc:83
l1t::P2GTCandidate::hwNumber_of_tracks_in_pv_
Optional< int > hwNumber_of_tracks_in_pv_
Definition: P2GTCandidate.h:322
l1t::L1GTProducer::scales_
const L1GTScales scales_
Definition: L1GTProducer.cc:65
reco::ParticleState::PolarLorentzVector
math::PtEtaPhiMLorentzVector PolarLorentzVector
Lorentz vector.
Definition: ParticleState.h:29
l1t::L1GTProducer::gttPromptJetToken_
const edm::EDGetTokenT< TkJetWordCollection > gttPromptJetToken_
Definition: L1GTProducer.cc:67
l1t::L1GTScales::to_phi
double to_phi(int value) const
Definition: L1GTScales.h:56
l1t::P2GTCandidate::objectType_
ObjectType objectType_
Definition: P2GTCandidate.h:325
l1t::L1GTScales::to_chg
int to_chg(int value) const
Definition: L1GTScales.h:60
l1t::L1GTScales::fillPSetDescription
static void fillPSetDescription(edm::ParameterSetDescription &)
Definition: L1GTScales.cc:55
l1t
delete x;
Definition: CaloConfig.h:22
config
Definition: config.py:1
reco::ParticleState::Point
math::XYZPoint Point
point in the space
Definition: ParticleState.h:31
l1t::P2GTCandidate::hwD0_
Optional< int > hwD0_
Definition: P2GTCandidate.h:310
l1t::L1GTProducer::produceGMTSaPromptMuons
void produceGMTSaPromptMuons(edm::Event &event) const
Definition: L1GTProducer.cc:206
l1t::L1GTProducer::cl2HtSumToken_
const edm::EDGetTokenT< std::vector< l1t::EtSum > > cl2HtSumToken_
Definition: L1GTProducer.cc:80
l1gt::Photon::isolationPT
iso_t isolationPT
Definition: gt_datatypes.h:317
l1t::L1GTProducer::cl2TauToken_
const edm::EDGetTokenT< PFTauCollection > cl2TauToken_
Definition: L1GTProducer.cc:78
l1t::L1GTScales::to_z0
double to_z0(int value) const
Definition: L1GTScales.h:58
l1t::L1GTProducer::produceCL2Taus
void produceCL2Taus(edm::Event &event) const
Definition: L1GTProducer.cc:357
l1t::L1GTProducer::produceCL2Photons
void produceCL2Photons(edm::Event &event) const
Definition: L1GTProducer.cc:309
l1t::P2GTCandidate::hwCharge_
Optional< int > hwCharge_
Definition: P2GTCandidate.h:309
l1gt::Jet::unpack
static Jet unpack(const std::array< uint64_t, 2 > &src)
Definition: gt_datatypes.h:129
l1t::PFJetCollection
std::vector< l1t::PFJet > PFJetCollection
Definition: PFJet.h:61
l1t::TkEmCollection
std::vector< TkEm > TkEmCollection
Definition: TkEmFwd.h:17
l1t::L1GTProducer::cl2JetToken_
const edm::EDGetTokenT< PFJetCollection > cl2JetToken_
Definition: L1GTProducer.cc:75
l1t::L1GTScales::to_eta
double to_eta(int value) const
Definition: L1GTScales.h:57
l1t::L1GTProducer::produceCL2Electrons
void produceCL2Electrons(edm::Event &event) const
Definition: L1GTProducer.cc:331
l1t::L1GTProducer::produceGMTTkMuons
void produceGMTTkMuons(edm::Event &event) const
Definition: L1GTProducer.cc:258
l1t::P2GTCandidate::hwPT_
Optional< int > hwPT_
Definition: P2GTCandidate.h:303
l1t::L1GTProducer::produceGTTPromptJets
void produceGTTPromptJets(edm::Event &event) const
Definition: L1GTProducer.cc:158
l1gt::Electron::qualityFlags
egquality_t qualityFlags
Definition: gt_datatypes.h:264
DEFINE_FWK_MODULE
#define DEFINE_FWK_MODULE(type)
Definition: MakerMacros.h:16
l1t::L1GTProducer::gmtTkMuonToken_
const edm::EDGetTokenT< TrackerMuonCollection > gmtTkMuonToken_
Definition: L1GTProducer.cc:73
l1t::VertexWordCollection
std::vector< VertexWord > VertexWordCollection
Definition: VertexWord.h:197
l1t::L1GTProducer::cl2EtSumToken_
const edm::EDGetTokenT< std::vector< l1t::EtSum > > cl2EtSumToken_
Definition: L1GTProducer.cc:79
singleTopDQM_cfi.setup
setup
SETUP
Definition: singleTopDQM_cfi.py:37
l1gt::Electron::v3
ThreeVector v3
Definition: gt_datatypes.h:263
l1t::L1GTProducer::produceCL2EtSum
void produceCL2EtSum(edm::Event &event) const
Definition: L1GTProducer.cc:382
l1t::L1GTProducer::fillDescriptions
static void fillDescriptions(edm::ConfigurationDescriptions &)
Definition: L1GTProducer.cc:113
l1t::L1GTProducer::gmtSaDisplacedMuonToken_
const edm::EDGetTokenT< SAMuonCollection > gmtSaDisplacedMuonToken_
Definition: L1GTProducer.cc:72
l1t::L1GTProducer::produceGTTDisplacedJets
void produceGTTDisplacedJets(edm::Event &event) const
Definition: L1GTProducer.cc:182
l1t::TrackerMuonCollection
std::vector< TrackerMuon > TrackerMuonCollection
Definition: TrackerMuon.h:16
l1t::SAMuonCollection
std::vector< SAMuon > SAMuonCollection
Definition: SAMuon.h:14
l1t::L1GTProducer::produce
void produce(edm::StreamID, edm::Event &, const edm::EventSetup &) const override
Definition: L1GTProducer.cc:418
l1t::P2GTCandidate::hwNumber_of_tracks_
Optional< int > hwNumber_of_tracks_
Definition: P2GTCandidate.h:317
l1t::L1GTProducer::gmtSaPromptMuonToken_
const edm::EDGetTokenT< SAMuonCollection > gmtSaPromptMuonToken_
Definition: L1GTProducer.cc:71
l1t::L1GTScales::to_pT
double to_pT(int value) const
Definition: L1GTScales.h:55
l1gt::Scales::ETAPHI_LSB
constexpr float ETAPHI_LSB
Definition: gt_datatypes.h:45
BTaggingMonitor_cfi.met
met
Definition: BTaggingMonitor_cfi.py:74
l1gt::Photon::v3
ThreeVector v3
Definition: gt_datatypes.h:315
edm
HLT enums.
Definition: AlignableModifier.h:19
l1t::L1GTProducer::produceCl2HtSum
void produceCl2HtSum(edm::Event &event) const
Definition: L1GTProducer.cc:401
l1t::P2GTCandidate::hwSum_pT_pv_
Optional< int > hwSum_pT_pv_
Definition: P2GTCandidate.h:320
l1gt::Electron::charge
ap_uint< 1 > charge
Definition: gt_datatypes.h:265
l1t::P2GTCandidate::hwNumber_of_tracks_not_in_pv_
Optional< int > hwNumber_of_tracks_not_in_pv_
Definition: P2GTCandidate.h:323
l1t::L1GTProducer::produceGTTPrimaryVert
void produceGTTPrimaryVert(edm::Event &event) const
Definition: L1GTProducer.cc:138
l1t::P2GTCandidate::hwQual_
Optional< int > hwQual_
Definition: P2GTCandidate.h:308
l1t::P2GTCandidate::hwPhi_
Optional< int > hwPhi_
Definition: P2GTCandidate.h:304
l1t::L1GTProducer::produceCL2Jets
void produceCL2Jets(edm::Event &event) const
Definition: L1GTProducer.cc:286
l1t::L1GTProducer::cl2PhotonToken_
const edm::EDGetTokenT< TkEmCollection > cl2PhotonToken_
Definition: L1GTProducer.cc:76
l1t::P2GTCandidate::hwEta_
Optional< int > hwEta_
Definition: P2GTCandidate.h:305
eostools.move
def move(src, dest)
Definition: eostools.py:511
l1t::P2GTCandidate::hwBeta_
Optional< int > hwBeta_
Definition: P2GTCandidate.h:311
l1t::L1GTProducer::produceGMTSaDisplacedMuons
void produceGMTSaDisplacedMuons(edm::Event &event) const
Definition: L1GTProducer.cc:232
event
Definition: event.py:1
l1t::L1GTProducer::gttDisplacedJetToken_
const edm::EDGetTokenT< TkJetWordCollection > gttDisplacedJetToken_
Definition: L1GTProducer.cc:68
l1t::P2GTCandidate::hwZ0_
Optional< int > hwZ0_
Definition: P2GTCandidate.h:306
Generated for CMSSW Reference Manual by doxygen 1.8.14