CMS 3D CMS Logo
L1TObjectsTiming.cc
Go to the documentation of this file.
1 #include "DQM/L1TMonitor/interface/L1TObjectsTiming.h"
2 
3 
4 L1TObjectsTiming::L1TObjectsTiming(const edm::ParameterSet& ps)
5  : ugmtMuonToken_(consumes<l1t::MuonBxCollection>(ps.getParameter<edm::InputTag>("muonProducer"))),
6  stage2CaloLayer2JetToken_(consumes<l1t::JetBxCollection>(ps.getParameter<edm::InputTag>("stage2CaloLayer2JetProducer"))),
7  stage2CaloLayer2EGammaToken_(consumes<l1t::EGammaBxCollection>(ps.getParameter<edm::InputTag>("stage2CaloLayer2EGammaProducer"))),
8  stage2CaloLayer2TauToken_(consumes<l1t::TauBxCollection>(ps.getParameter<edm::InputTag>("stage2CaloLayer2TauProducer"))),
9  stage2CaloLayer2EtSumToken_(consumes<l1t::EtSumBxCollection>(ps.getParameter<edm::InputTag>("stage2CaloLayer2EtSumProducer"))),
10  l1tStage2uGtProducer_(consumes<GlobalAlgBlkBxCollection>(ps.getParameter<edm::InputTag>("ugtProducer"))),
11  monitorDir_(ps.getUntrackedParameter<std::string>("monitorDir")),
12  verbose_(ps.getUntrackedParameter<bool>("verbose")),
13  gtUtil_(new l1t::L1TGlobalUtil(ps, consumesCollector(), *this, ps.getParameter<edm::InputTag>("ugtProducer"), ps.getParameter<edm::InputTag>("ugtProducer"))),
14  algoBitFirstBxInTrain_(-1),
15  algoBitLastBxInTrain_(-1),
16  algoNameFirstBxInTrain_(ps.getUntrackedParameter<std::string>("firstBXInTrainAlgo","")),
17  algoNameLastBxInTrain_(ps.getUntrackedParameter<std::string>("lastBXInTrainAlgo","")),
18  bxrange_(5)
19 
20 {
21 
22  muons_eta_phi.reserve(bxrange_);
23  jet_eta_phi.reserve(bxrange_);
24  egamma_eta_phi.reserve(bxrange_);
25  tau_eta_phi.reserve(bxrange_);
26  etsum_eta_phi_MET.reserve(bxrange_);
27  etsum_eta_phi_METHF.reserve(bxrange_);
28  etsum_eta_phi_MHT.reserve(bxrange_);
29  etsum_eta_phi_MHTHF.reserve(bxrange_);
30 
31  muons_eta_phi_isolated.reserve(bxrange_);
32  jet_eta_phi_isolated.reserve(bxrange_);
33  egamma_eta_phi_isolated.reserve(bxrange_);
34  tau_eta_phi_isolated.reserve(bxrange_);
35  etsum_eta_phi_MET_isolated.reserve(bxrange_);
36  etsum_eta_phi_METHF_isolated.reserve(bxrange_);
37  etsum_eta_phi_MHT_isolated.reserve(bxrange_);
38  etsum_eta_phi_MHTHF_isolated.reserve(bxrange_);
39 
40  muons_eta_phi_firstbunch.reserve(bxrange_);
41  jet_eta_phi_firstbunch.reserve(bxrange_);
42  egamma_eta_phi_firstbunch.reserve(bxrange_);
43  tau_eta_phi_firstbunch.reserve(bxrange_);
44  etsum_eta_phi_MET_firstbunch.reserve(bxrange_);
45  etsum_eta_phi_METHF_firstbunch.reserve(bxrange_);
46  etsum_eta_phi_MHT_firstbunch.reserve(bxrange_);
47  etsum_eta_phi_MHTHF_firstbunch.reserve(bxrange_);
48 
49  muons_eta_phi_lastbunch.reserve(bxrange_);
50  jet_eta_phi_lastbunch.reserve(bxrange_);
51  egamma_eta_phi_lastbunch.reserve(bxrange_);
52  tau_eta_phi_lastbunch.reserve(bxrange_);
53  etsum_eta_phi_MET_lastbunch.reserve(bxrange_);
54  etsum_eta_phi_METHF_lastbunch.reserve(bxrange_);
55  etsum_eta_phi_MHT_lastbunch.reserve(bxrange_);
56  etsum_eta_phi_MHTHF_lastbunch.reserve(bxrange_);
57 
58 }
59 
60 L1TObjectsTiming::~L1TObjectsTiming() {}
61 
62 void L1TObjectsTiming::fillDescriptions(edm::ConfigurationDescriptions& descriptions) {
63  edm::ParameterSetDescription desc;
64  desc.add<edm::InputTag>("muonProducer")->setComment("L1T muons");;
65  desc.add<edm::InputTag>("stage2CaloLayer2JetProducer")->setComment("L1T jets");
66  desc.add<edm::InputTag>("stage2CaloLayer2EGammaProducer")->setComment("L1T egamma");
67  desc.add<edm::InputTag>("stage2CaloLayer2TauProducer")->setComment("L1T taus");
68  desc.add<edm::InputTag>("stage2CaloLayer2EtSumProducer")->setComment("L1T etsums");
69  desc.add<edm::InputTag>("ugtProducer")->setComment("uGT output");
70  desc.addUntracked<std::string>("monitorDir", "")->setComment("Target directory in the DQM file. Will be created if not existing.");
71  desc.addUntracked<bool>("verbose", false);
72  desc.addUntracked<std::string>("firstBXInTrainAlgo","")->setComment("Pick the right algo name for L1 First Collision In Train");
73  desc.addUntracked<std::string>("lastBXInTrainAlgo","")->setComment("Pick the right algo name for L1 Last Collision In Train");
74  descriptions.add("l1tObjectsTiming", desc);
75 }
76 
77 void L1TObjectsTiming::dqmBeginRun(const edm::Run& r, const edm::EventSetup& c) {
78  // Get the trigger menu information
79  gtUtil_->retrieveL1Setup(c);
80  // Get the algo bits needed for the timing histograms
81  if (!gtUtil_->getAlgBitFromName(algoNameFirstBxInTrain_, algoBitFirstBxInTrain_)) {
82  edm::LogWarning("L1TObjectsTiming") << "Algo \"" << algoNameFirstBxInTrain_ << "\" not found in the trigger menu " << gtUtil_->gtTriggerMenuName() << ". Could not retrieve algo bit number.";
83  }
84  if (!gtUtil_->getAlgBitFromName(algoNameLastBxInTrain_, algoBitLastBxInTrain_)) {
85  edm::LogWarning("L1TObjectsTiming") << "Algo \"" << algoNameLastBxInTrain_ << "\" not found in the trigger menu " << gtUtil_->gtTriggerMenuName() << ". Could not retrieve algo bit number.";
86  }
87 
88 
89 }
90 
91 void L1TObjectsTiming::beginLuminosityBlock(const edm::LuminosityBlock&, const edm::EventSetup&) {}
92 
93 void L1TObjectsTiming::bookHistograms(DQMStore::IBooker& ibooker, const edm::Run&, const edm::EventSetup&) {
94 
95 
96  // Subsystem Monitoring and Muon Output
97  std::string bx_obj[5]={"minus2","minus1","0","plus1","plus2"};
98 
99  ibooker.setCurrentFolder(monitorDir_+"/L1TMuon"+"/timing");
100  denominator_muons = ibooker.book2D("denominator_muons","Denominator for L1TMuon",25, -2.5, 2.5, 25, -3.2, 3.2);
101  denominator_muons->setAxisTitle("#eta", 1);
102  denominator_muons->setAxisTitle("#phi", 2);
103  for(unsigned int i=0; i<bxrange_; ++i) {
104  muons_eta_phi.push_back(ibooker.book2D("muons_eta_phi_bx_"+bx_obj[i],"L1T Muon #eta vs #phi BX="+bx_obj[i],25, -2.5, 2.5, 25, -3.2, 3.2));
105  muons_eta_phi.at(i)->setAxisTitle("#eta", 1);
106  muons_eta_phi.at(i)->setAxisTitle("#phi", 2);
107  }
108 
109  ibooker.setCurrentFolder(monitorDir_+"/L1TJet"+"/timing");
110  denominator_jet = ibooker.book2D("denominator_jet","Denominator for L1TJet",25, -2.5, 2.5, 25, -3.2, 3.2);
111  denominator_jet->setAxisTitle("#eta", 1);
112  denominator_jet->setAxisTitle("#phi", 2);
113  for(unsigned int i=0; i<bxrange_; ++i) {
114  jet_eta_phi.push_back(ibooker.book2D("jet_eta_phi_bx_"+bx_obj[i],"L1T Jet #eta vs #phi BX="+bx_obj[i],25, -2.5, 2.5, 25, -3.2, 3.2));
115  jet_eta_phi.at(i)->setAxisTitle("#eta", 1);
116  jet_eta_phi.at(i)->setAxisTitle("#phi", 2);
117  }
118 
119  ibooker.setCurrentFolder(monitorDir_+"/L1TEGamma"+"/timing");
120  denominator_egamma = ibooker.book2D("denominator_egamma","Denominator for L1TEGamma",25, -2.5, 2.5, 25, -3.2, 3.2);
121  denominator_egamma->setAxisTitle("#eta", 1);
122  denominator_egamma->setAxisTitle("#phi", 2);
123  for(unsigned int i=0; i<bxrange_; ++i) {
124  egamma_eta_phi.push_back(ibooker.book2D("egamma_eta_phi_bx_"+bx_obj[i],"L1T EGamma #eta vs #phi BX="+bx_obj[i],25, -2.5, 2.5, 25, -3.2, 3.2));
125  egamma_eta_phi.at(i)->setAxisTitle("#eta", 1);
126  egamma_eta_phi.at(i)->setAxisTitle("#phi", 2);
127  }
128 
129  ibooker.setCurrentFolder(monitorDir_+"/L1TTau"+"/timing");
130  denominator_tau = ibooker.book2D("denominator_tau","Denominator for L1TTau",25, -2.5, 2.5, 25, -3.2, 3.2);
131  denominator_tau->setAxisTitle("#eta", 1);
132  denominator_tau->setAxisTitle("#phi", 2);
133  for(unsigned int i=0; i<bxrange_; ++i) {
134  tau_eta_phi.push_back(ibooker.book2D("tau_eta_phi_bx_"+bx_obj[i],"L1T Tau #eta vs #phi BX="+bx_obj[i],25, -2.5, 2.5, 25, -3.2, 3.2));
135  tau_eta_phi.at(i)->setAxisTitle("#eta", 1);
136  tau_eta_phi.at(i)->setAxisTitle("#phi", 2);
137  }
138 
139  ibooker.setCurrentFolder(monitorDir_+"/L1TEtSum"+"/timing");
140  denominator_etsum_MET = ibooker.book1D("denominator_etsum_MET","Denominator for L1TEtSum MET", 25, -3.2, 3.2);
141  denominator_etsum_MET->setAxisTitle("#phi", 1);
142  denominator_etsum_METHF = ibooker.book1D("denominator_etsum_METHF","Denominator for L1TEtSum METHF", 25, -3.2, 3.2);
143  denominator_etsum_METHF->setAxisTitle("#phi", 1);
144  denominator_etsum_MHT = ibooker.book1D("denominator_etsum_MHT","Denominator for L1TEtSum MHT", 25, -3.2, 3.2);
145  denominator_etsum_MHT->setAxisTitle("#phi", 1);
146  denominator_etsum_MHTHF = ibooker.book1D("denominator_etsum_MHTHF","Denominator for L1TEtSum MHTHF", 25, -3.2, 3.2);
147  denominator_etsum_MHTHF->setAxisTitle("#phi", 1);
148  for(unsigned int i=0; i<bxrange_; ++i) {
149  etsum_eta_phi_MET.push_back(ibooker.book1D("etsum_phi_bx_MET_"+bx_obj[i],"L1T MET #phi BX="+bx_obj[i],25, -3.2, 3.2));
150  etsum_eta_phi_MET.at(i)->setAxisTitle("#phi", 1);
151  etsum_eta_phi_METHF.push_back(ibooker.book1D("etsum_phi_bx_METHF_"+bx_obj[i],"L1T METHF #phi BX="+bx_obj[i],25, -3.2, 3.2));
152  etsum_eta_phi_METHF.at(i)->setAxisTitle("#phi", 1);
153  etsum_eta_phi_MHT.push_back(ibooker.book1D("etsum_phi_bx_MHT_"+bx_obj[i],"L1T MHT #phi BX="+bx_obj[i],25, -3.2, 3.2));
154  etsum_eta_phi_MHT.at(i)->setAxisTitle("#phi", 1);
155  etsum_eta_phi_MHTHF.push_back(ibooker.book1D("etsum_phi_bx_MHTHF_"+bx_obj[i],"L1T MHTHF #phi BX="+bx_obj[i],25, -3.2, 3.2));
156  etsum_eta_phi_MHTHF.at(i)->setAxisTitle("#phi", 1);
157  }
158 
159 if(algoBitFirstBxInTrain_ > -1 && algoBitLastBxInTrain_ > -1) {
160  ibooker.setCurrentFolder(monitorDir_+"/L1TMuon"+"/timing"+"/Isolated_bunch");
161  denominator_muons_isolated = ibooker.book2D("denominator_muons_isolated","Denominator for Isolated Bunch for L1TMuon",25, -2.5, 2.5, 25, -3.2, 3.2);
162  denominator_muons_isolated->setAxisTitle("#eta", 1);
163  denominator_muons_isolated->setAxisTitle("#phi", 2);
164  ibooker.setCurrentFolder(monitorDir_+"/L1TJet"+"/timing"+"/Isolated_bunch");
165  denominator_jet_isolated = ibooker.book2D("denominator_jet_isolated","Denominator for Isolated Bunch for L1TJet",25, -2.5, 2.5, 25, -3.2, 3.2);
166  denominator_jet_isolated->setAxisTitle("#eta", 1);
167  denominator_jet_isolated->setAxisTitle("#phi", 2);
168  ibooker.setCurrentFolder(monitorDir_+"/L1TEGamma"+"/timing"+"/Isolated_bunch");
169  denominator_egamma_isolated = ibooker.book2D("denominator_egamma_isolated","Denominator for Isolated Bunch for L1EGamma",25, -2.5, 2.5, 25, -3.2, 3.2);
170  denominator_egamma_isolated->setAxisTitle("#eta", 1);
171  denominator_egamma_isolated->setAxisTitle("#phi", 2);
172  ibooker.setCurrentFolder(monitorDir_+"/L1TTau"+"/timing"+"/Isolated_bunch");
173  denominator_tau_isolated = ibooker.book2D("denominator_tau_isolated","Denominator for Isolated Bunch for L1TTau",25, -2.5, 2.5, 25, -3.2, 3.2);
174  denominator_tau_isolated->setAxisTitle("#eta", 1);
175  denominator_tau_isolated->setAxisTitle("#phi", 2);
176  ibooker.setCurrentFolder(monitorDir_+"/L1TEtSum"+"/timing"+"/Isolated_bunch");
177  denominator_etsum_isolated_MET = ibooker.book1D("denominator_etsum_isolated_MET","Denominator for Isolated Bunch for L1TEtSum MET", 25, -3.2, 3.2);
178  denominator_etsum_isolated_MET->setAxisTitle("#phi", 1);
179  denominator_etsum_isolated_METHF = ibooker.book1D("denominator_etsum_isolated_METHF","Denominator for Isolated Bunch for L1TEtSum METHF", 25, -3.2, 3.2);
180  denominator_etsum_isolated_METHF->setAxisTitle("#phi", 1);
181  denominator_etsum_isolated_MHT = ibooker.book1D("denominator_etsum_isolated_MHT","Denominator for Isolated Bunch for L1TEtSum MHT", 25, -3.2, 3.2);
182  denominator_etsum_isolated_MHT->setAxisTitle("#phi", 1);
183  denominator_etsum_isolated_MHTHF = ibooker.book1D("denominator_etsum_isolated_MHTHF","Denominator for Isolated Bunch for L1TEtSum MHTHF", 25, -3.2, 3.2);
184  denominator_etsum_isolated_MHTHF->setAxisTitle("#phi", 1);
185 
186  ibooker.setCurrentFolder(monitorDir_+"/L1TMuon"+"/timing"+"/Isolated_bunch");
187  for(unsigned int i=0; i<bxrange_; ++i) {
188  muons_eta_phi_isolated.push_back(ibooker.book2D("muons_eta_phi_bx_isolated_"+bx_obj[i],"L1T Muon #eta vs #phi for isolated bunch BX="+bx_obj[i],25, -2.5, 2.5, 25, -3.2, 3.2));
189  muons_eta_phi_isolated.at(i)->setAxisTitle("#eta", 1);
190  muons_eta_phi_isolated.at(i)->setAxisTitle("#phi", 2);
191  }
192 
193  ibooker.setCurrentFolder(monitorDir_+"/L1TJet"+"/timing"+"/Isolated_bunch");
194  for(unsigned int i=0; i<bxrange_; ++i) {
195  jet_eta_phi_isolated.push_back(ibooker.book2D("jet_eta_phi_bx_isolated_"+bx_obj[i],"L1T Jet #eta vs #phi for isolated bunch BX="+bx_obj[i],25, -2.5, 2.5, 25, -3.2, 3.2));
196  jet_eta_phi_isolated.at(i)->setAxisTitle("#eta", 1);
197  jet_eta_phi_isolated.at(i)->setAxisTitle("#phi", 2);
198  }
199 
200  ibooker.setCurrentFolder(monitorDir_+"/L1TEGamma"+"/timing"+"/Isolated_bunch");
201  for(unsigned int i=0; i<bxrange_; ++i) {
202  egamma_eta_phi_isolated.push_back(ibooker.book2D("egamma_eta_phi_bx_isolated_"+bx_obj[i],"L1T EGamma #eta vs #phi for isolated bunch BX="+bx_obj[i],25, -2.5, 2.5, 25, -3.2, 3.2));
203  egamma_eta_phi_isolated.at(i)->setAxisTitle("#eta", 1);
204  egamma_eta_phi_isolated.at(i)->setAxisTitle("#phi", 2);
205  }
206 
207  ibooker.setCurrentFolder(monitorDir_+"/L1TTau"+"/timing"+"/Isolated_bunch");
208  for(unsigned int i=0; i<bxrange_; ++i) {
209  tau_eta_phi_isolated.push_back(ibooker.book2D("tau_eta_phi_bx_isolated_"+bx_obj[i],"L1T Tau #eta vs #phi for isolated bunch BX="+bx_obj[i],25, -2.5, 2.5, 25, -3.2, 3.2));
210  tau_eta_phi_isolated.at(i)->setAxisTitle("#eta", 1);
211  tau_eta_phi_isolated.at(i)->setAxisTitle("#phi", 2);
212  }
213 
214  ibooker.setCurrentFolder(monitorDir_+"/L1TEtSum"+"/timing"+"/Isolated_bunch");
215  for(unsigned int i=0; i<bxrange_; ++i) {
216  etsum_eta_phi_MET_isolated.push_back(ibooker.book1D("etsum_phi_bx_MET_isolated_"+bx_obj[i],"L1T MET #phi for isolated bunch BX="+bx_obj[i],25, -3.2, 3.2));
217  etsum_eta_phi_MET_isolated.at(i)->setAxisTitle("#phi", 1);
218  etsum_eta_phi_METHF_isolated.push_back(ibooker.book1D("etsum_phi_bx_METHF_isolated_"+bx_obj[i],"L1T METHF #phi for isolated bunch BX="+bx_obj[i],25, -3.2, 3.2));
219  etsum_eta_phi_METHF_isolated.at(i)->setAxisTitle("#phi", 1);
220  etsum_eta_phi_MHT_isolated.push_back(ibooker.book1D("etsum_phi_bx_MHT_isolated_"+bx_obj[i],"L1T MHT #phi for isolated bunch BX="+bx_obj[i],25, -3.2, 3.2));
221  etsum_eta_phi_MHT_isolated.at(i)->setAxisTitle("#phi", 1);
222  etsum_eta_phi_MHTHF_isolated.push_back(ibooker.book1D("etsum_phi_bx_MHTHF_isolated_"+bx_obj[i],"L1T MHTHF #phi for isolated bunch BX="+bx_obj[i],25, -3.2, 3.2));
223  etsum_eta_phi_MHTHF_isolated.at(i)->setAxisTitle("#phi", 1);
224  }
225 }
226 
227 if(algoBitFirstBxInTrain_ > -1) {
228  ibooker.setCurrentFolder(monitorDir_+"/L1TMuon"+"/timing"+"/First_bunch");
229  denominator_muons_firstbunch = ibooker.book2D("denominator_muons_firstbunch","Denominator for First Bunch for L1TMuon",25, -2.5, 2.5, 25, -3.2, 3.2);
230  denominator_muons_firstbunch->setAxisTitle("#eta", 1);
231  denominator_muons_firstbunch->setAxisTitle("#phi", 2);
232  ibooker.setCurrentFolder(monitorDir_+"/L1TJet"+"/timing"+"/First_bunch");
233  denominator_jet_firstbunch = ibooker.book2D("denominator_jet_firstbunch","Denominator for First Bunch for L1TJet",25, -2.5, 2.5, 25, -3.2, 3.2);
234  denominator_jet_firstbunch->setAxisTitle("#eta", 1);
235  denominator_jet_firstbunch->setAxisTitle("#phi", 2);
236  ibooker.setCurrentFolder(monitorDir_+"/L1TEGamma"+"/timing"+"/First_bunch");
237  denominator_egamma_firstbunch = ibooker.book2D("denominator_egamma_firstbunch","Denominator for First Bunch for L1TEGamma",25, -2.5, 2.5, 25, -3.2, 3.2);
238  denominator_egamma_firstbunch->setAxisTitle("#eta", 1);
239  denominator_egamma_firstbunch->setAxisTitle("#phi", 2);
240  ibooker.setCurrentFolder(monitorDir_+"/L1TTau"+"/timing"+"/First_bunch");
241  denominator_tau_firstbunch = ibooker.book2D("denominator_tau_firstbunch","Denominator for First Bunch for L1TTau",25, -2.5, 2.5, 25, -3.2, 3.2);
242  denominator_tau_firstbunch->setAxisTitle("#eta", 1);
243  denominator_tau_firstbunch->setAxisTitle("#phi", 2);
244  ibooker.setCurrentFolder(monitorDir_+"/L1TEtSum"+"/timing"+"/First_bunch");
245  denominator_etsum_firstbunch_MET = ibooker.book1D("denominator_etsum_firstbunch_MET","Denominator for First Bunch for L1TEtSum MET", 25, -3.2, 3.2);
246  denominator_etsum_firstbunch_MET->setAxisTitle("#phi", 1);
247  denominator_etsum_firstbunch_METHF = ibooker.book1D("denominator_etsum_firstbunch_METHF","Denominator for First Bunch for L1TEtSum METHF", 25, -3.2, 3.2);
248  denominator_etsum_firstbunch_METHF->setAxisTitle("#phi", 1);
249  denominator_etsum_firstbunch_MHT = ibooker.book1D("denominator_etsum_firstbunch_MHT","Denominator for First Bunch for L1TEtSum MHT", 25, -3.2, 3.2);
250  denominator_etsum_firstbunch_MHT->setAxisTitle("#phi", 1);
251  denominator_etsum_firstbunch_MHTHF = ibooker.book1D("denominator_etsum_firstbunch_MHTHF","Denominator for First Bunch for L1TEtSum MHTHF", 25, -3.2, 3.2);
252  denominator_etsum_firstbunch_MHTHF->setAxisTitle("#phi", 1);
253 
254  ibooker.setCurrentFolder(monitorDir_+"/L1TMuon"+"/timing"+"/First_bunch");
255  for(unsigned int i=0; i<bxrange_-2; ++i) {
256  muons_eta_phi_firstbunch.push_back(ibooker.book2D("muons_eta_phi_bx_firstbunch_"+bx_obj[i],"L1T Muon #eta vs #phi for first bunch BX="+bx_obj[i],25, -2.5, 2.5, 25, -3.2, 3.2));
257  muons_eta_phi_firstbunch.at(i)->setAxisTitle("#eta", 1);
258  muons_eta_phi_firstbunch.at(i)->setAxisTitle("#phi", 2);
259  }
260 
261  ibooker.setCurrentFolder(monitorDir_+"/L1TJet"+"/timing"+"/First_bunch");
262  for(unsigned int i=0; i<bxrange_-2; ++i) {
263  jet_eta_phi_firstbunch.push_back(ibooker.book2D("jet_eta_phi_bx_firstbunch_"+bx_obj[i],"L1T Jet #eta vs #phi for first bunch BX="+bx_obj[i],25, -2.5, 2.5, 25, -3.2, 3.2));
264  jet_eta_phi_firstbunch.at(i)->setAxisTitle("#eta", 1);
265  jet_eta_phi_firstbunch.at(i)->setAxisTitle("#phi", 2);
266  }
267 
268  ibooker.setCurrentFolder(monitorDir_+"/L1TEGamma"+"/timing"+"/First_bunch");
269  for(unsigned int i=0; i<bxrange_-2; ++i) {
270  egamma_eta_phi_firstbunch.push_back(ibooker.book2D("egamma_eta_phi_bx_firstbunch_"+bx_obj[i],"L1T EGamma #eta vs #phi for first bunch BX="+bx_obj[i],25, -2.5, 2.5, 25, -3.2, 3.2));
271  egamma_eta_phi_firstbunch.at(i)->setAxisTitle("#eta", 1);
272  egamma_eta_phi_firstbunch.at(i)->setAxisTitle("#phi", 2);
273  }
274 
275  ibooker.setCurrentFolder(monitorDir_+"/L1TTau"+"/timing"+"/First_bunch");
276  for(unsigned int i=0; i<bxrange_-2; ++i) {
277  tau_eta_phi_firstbunch.push_back(ibooker.book2D("tau_eta_phi_bx_firstbunch_"+bx_obj[i],"L1T Tau #eta vs #phi for first bunch BX="+bx_obj[i],25, -2.5, 2.5, 25, -3.2, 3.2));
278  tau_eta_phi_firstbunch.at(i)->setAxisTitle("#eta", 1);
279  tau_eta_phi_firstbunch.at(i)->setAxisTitle("#phi", 2);
280  }
281 
282  ibooker.setCurrentFolder(monitorDir_+"/L1TEtSum"+"/timing"+"/First_bunch");
283  for(unsigned int i=0; i<bxrange_-2; ++i) {
284  etsum_eta_phi_MET_firstbunch.push_back(ibooker.book1D("etsum_phi_bx_MET_firstbunch_"+bx_obj[i],"L1T MET #phi for firstbunch bunch BX="+bx_obj[i],25, -3.2, 3.2));
285  etsum_eta_phi_MET_firstbunch.at(i)->setAxisTitle("#phi", 1);
286  etsum_eta_phi_METHF_firstbunch.push_back(ibooker.book1D("etsum_phi_bx_METHF_firstbunch_"+bx_obj[i],"L1T METHF #phi for firstbunch bunch BX="+bx_obj[i],25, -3.2, 3.2));
287  etsum_eta_phi_METHF_firstbunch.at(i)->setAxisTitle("#phi", 1);
288  etsum_eta_phi_MHT_firstbunch.push_back(ibooker.book1D("etsum_phi_bx_MHT_firstbunch_"+bx_obj[i],"L1T MHT #phi for firstbunch bunch BX="+bx_obj[i],25, -3.2, 3.2));
289  etsum_eta_phi_MHT_firstbunch.at(i)->setAxisTitle("#phi", 1);
290  etsum_eta_phi_MHTHF_firstbunch.push_back(ibooker.book1D("etsum_phi_bx_MHTHF_firstbunch_"+bx_obj[i],"L1T MHTHF #phi for firstbunch bunch BX="+bx_obj[i],25, -3.2, 3.2));
291  etsum_eta_phi_MHTHF_firstbunch.at(i)->setAxisTitle("#phi", 1);
292  }
293 }
294 
295 if(algoBitLastBxInTrain_ > -1) {
296  ibooker.setCurrentFolder(monitorDir_+"/L1TMuon"+"/timing"+"/Last_bunch");
297  denominator_muons_lastbunch = ibooker.book2D("denominator_muons_lastbunch","Denominator for Last Bunch for L1TMuon",25, -2.5, 2.5, 25, -3.2, 3.2);
298  denominator_muons_lastbunch->setAxisTitle("#eta", 1);
299  denominator_muons_lastbunch->setAxisTitle("#phi", 2);
300  ibooker.setCurrentFolder(monitorDir_+"/L1TJet"+"/timing"+"/Last_bunch");
301  denominator_jet_lastbunch = ibooker.book2D("denominator_jet_lastbunch","Denominator for Last Bunch for L1TJet",25, -2.5, 2.5, 25, -3.2, 3.2);
302  denominator_jet_lastbunch->setAxisTitle("#eta", 1);
303  denominator_jet_lastbunch->setAxisTitle("#phi", 2);
304  ibooker.setCurrentFolder(monitorDir_+"/L1TEGamma"+"/timing"+"/Last_bunch");
305  denominator_egamma_lastbunch = ibooker.book2D("denominator_egamma_lastbunch","Denominator for Last Bunch for L1TEGamma",25, -2.5, 2.5, 25, -3.2, 3.2);
306  denominator_egamma_lastbunch->setAxisTitle("#eta", 1);
307  denominator_egamma_lastbunch->setAxisTitle("#phi", 2);
308  ibooker.setCurrentFolder(monitorDir_+"/L1TTau"+"/timing"+"/Last_bunch");
309  denominator_tau_lastbunch = ibooker.book2D("denominator_tau_lastbunch","Denominator for Last Bunch for L1TTau",25, -2.5, 2.5, 25, -3.2, 3.2);
310  denominator_tau_lastbunch->setAxisTitle("#eta", 1);
311  denominator_tau_lastbunch->setAxisTitle("#phi", 2);
312  ibooker.setCurrentFolder(monitorDir_+"/L1TEtSum"+"/timing"+"/Last_bunch");
313  denominator_etsum_lastbunch_MET = ibooker.book1D("denominator_etsum_lastbunch_MET","Denominator for Last Bunch for L1TEtSum MET", 25, -3.2, 3.2);
314  denominator_etsum_lastbunch_MET->setAxisTitle("#phi", 1);
315  denominator_etsum_lastbunch_METHF = ibooker.book1D("denominator_etsum_lastbunch_METHF","Denominator for Last Bunch for L1TEtSum METHF", 25, -3.2, 3.2);
316  denominator_etsum_lastbunch_METHF->setAxisTitle("#phi", 1);
317  denominator_etsum_lastbunch_MHT = ibooker.book1D("denominator_etsum_lastbunch_MHT","Denominator for Last Bunch for L1TEtSum MHT", 25, -3.2, 3.2);
318  denominator_etsum_lastbunch_MHT->setAxisTitle("#phi", 1);
319  denominator_etsum_lastbunch_MHTHF = ibooker.book1D("denominator_etsum_lastbunch_MHTHF","Denominator for Last Bunch for L1TEtSum MHTHF", 25, -3.2, 3.2);
320  denominator_etsum_lastbunch_MHTHF->setAxisTitle("#phi", 1);
321 
322  ibooker.setCurrentFolder(monitorDir_+"/L1TMuon"+"/timing"+"/Last_bunch");
323  for(unsigned int i=0; i<bxrange_-2; ++i) {
324  muons_eta_phi_lastbunch.push_back(ibooker.book2D("muons_eta_phi_bx_lastbunch_"+bx_obj[i+2],"L1T Muon #eta vs #phi for last bunch BX="+bx_obj[i+2],25, -2.5, 2.5, 25, -3.2, 3.2));
325  muons_eta_phi_lastbunch.at(i)->setAxisTitle("#eta", 1);
326  muons_eta_phi_lastbunch.at(i)->setAxisTitle("#phi", 2);
327  }
328 
329  ibooker.setCurrentFolder(monitorDir_+"/L1TJet"+"/timing"+"/Last_bunch");
330  for(unsigned int i=0; i<bxrange_-2; ++i) {
331  jet_eta_phi_lastbunch.push_back(ibooker.book2D("jet_eta_phi_bx_lastbunch_"+bx_obj[i+2],"L1T Jet #eta vs #phi for last bunch BX="+bx_obj[i+2],25, -2.5, 2.5, 25, -3.2, 3.2));
332  jet_eta_phi_lastbunch.at(i)->setAxisTitle("#eta", 1);
333  jet_eta_phi_lastbunch.at(i)->setAxisTitle("#phi", 2);
334  }
335 
336  ibooker.setCurrentFolder(monitorDir_+"/L1TEGamma"+"/timing"+"/Last_bunch");
337  for(unsigned int i=0; i<bxrange_-2; ++i) {
338  egamma_eta_phi_lastbunch.push_back(ibooker.book2D("egamma_eta_phi_bx_lastbunch_"+bx_obj[i+2],"L1T EGamma #eta vs #phi for last bunch BX="+bx_obj[i+2],25, -2.5, 2.5, 25, -3.2, 3.2));
339  egamma_eta_phi_lastbunch.at(i)->setAxisTitle("#eta", 1);
340  egamma_eta_phi_lastbunch.at(i)->setAxisTitle("#phi", 2);
341  }
342 
343  ibooker.setCurrentFolder(monitorDir_+"/L1TTau"+"/timing"+"/Last_bunch");
344  for(unsigned int i=0; i<bxrange_-2; ++i) {
345  tau_eta_phi_lastbunch.push_back(ibooker.book2D("tau_eta_phi_bx_lastbunch_"+bx_obj[i+2],"L1T Tau #eta vs #phi for last bunch BX="+bx_obj[i+2],25, -2.5, 2.5, 25, -3.2, 3.2));
346  tau_eta_phi_lastbunch.at(i)->setAxisTitle("#eta", 1);
347  tau_eta_phi_lastbunch.at(i)->setAxisTitle("#phi", 2);
348  }
349 
350  ibooker.setCurrentFolder(monitorDir_+"/L1TEtSum"+"/timing"+"/Last_bunch");
351  for(unsigned int i=0; i<bxrange_-2; ++i) {
352  etsum_eta_phi_MET_lastbunch.push_back(ibooker.book1D("etsum_phi_bx_MET_lastbunch_"+bx_obj[i+2],"L1T MET #phi for lastbunch bunch BX="+bx_obj[i+2],25, -3.2, 3.2));
353  etsum_eta_phi_MET_lastbunch.at(i)->setAxisTitle("#phi", 1);
354  etsum_eta_phi_METHF_lastbunch.push_back(ibooker.book1D("etsum_phi_bx_METHF_lastbunch_"+bx_obj[i+2],"L1T METHF #phi for lastbunch bunch BX="+bx_obj[i+2],25, -3.2, 3.2));
355  etsum_eta_phi_METHF_lastbunch.at(i)->setAxisTitle("#phi", 1);
356  etsum_eta_phi_MHT_lastbunch.push_back(ibooker.book1D("etsum_phi_bx_MHT_lastbunch_"+bx_obj[i+2],"L1T MHT #phi for lastbunch bunch BX="+bx_obj[i+2],25, -3.2, 3.2));
357  etsum_eta_phi_MHT_lastbunch.at(i)->setAxisTitle("#phi", 1);
358  etsum_eta_phi_MHTHF_lastbunch.push_back(ibooker.book1D("etsum_phi_bx_MHTHF_lastbunch_"+bx_obj[i+2],"L1T MHTHF #phi for lastbunch bunch BX="+bx_obj[i+2],25, -3.2, 3.2));
359  etsum_eta_phi_MHTHF_lastbunch.at(i)->setAxisTitle("#phi", 1);
360  }
361 }
362 
363 }
364 
365 
366 
367 void L1TObjectsTiming::analyze(const edm::Event& e, const edm::EventSetup& c) {
368 
369  if (verbose_) edm::LogInfo("L1TObjectsTiming") << "L1TObjectsTiming: analyze..." << std::endl;
370 
371 
372  // Muon Collection
373  edm::Handle<l1t::MuonBxCollection> MuonBxCollection;
374  e.getByToken(ugmtMuonToken_, MuonBxCollection);
375  // Jet Collection
376  edm::Handle<l1t::JetBxCollection> JetBxCollection;
377  e.getByToken(stage2CaloLayer2JetToken_, JetBxCollection);
378  // EGamma Collection
379  edm::Handle<l1t::EGammaBxCollection> EGammaBxCollection;
380  e.getByToken(stage2CaloLayer2EGammaToken_, EGammaBxCollection);
381  // Tau Collection
382  edm::Handle<l1t::TauBxCollection> TauBxCollection;
383  e.getByToken(stage2CaloLayer2TauToken_, TauBxCollection);
384  // EtSum Collection
385  edm::Handle<l1t::EtSumBxCollection> EtSumBxCollection;
386  e.getByToken(stage2CaloLayer2EtSumToken_, EtSumBxCollection);
387 
388  // Open uGT readout record
389  edm::Handle<GlobalAlgBlkBxCollection> uGtAlgs;
390  e.getByToken(l1tStage2uGtProducer_, uGtAlgs);
391 
392 
393  // Filling eta-phi map for muons for BX=-2,-1,0,+1,+2
394  for (int itBX = MuonBxCollection->getFirstBX(); itBX <= MuonBxCollection->getLastBX(); ++itBX) {
395  for (l1t::MuonBxCollection::const_iterator Muon = MuonBxCollection->begin(itBX); Muon != MuonBxCollection->end(itBX); ++Muon) {
396  denominator_muons->Fill(Muon->eta(), Muon->phi());
397  int index = (int)itBX - std::min(0, 1 - (int)bxrange_%2 - (int)std::floor(bxrange_/2.)); // the correlation from itBX to respective index of the vector
398  muons_eta_phi.at(index)->Fill(Muon->eta(), Muon->phi());
399  }
400  }
401 
402  // Filling eta-phi map for jets for BX=-2,-1,0,+1,+2
403  for (int itBX = JetBxCollection->getFirstBX(); itBX <= JetBxCollection->getLastBX(); ++itBX) {
404  for (l1t::JetBxCollection::const_iterator jet = JetBxCollection->begin(itBX); jet != JetBxCollection->end(itBX); ++jet) {
405  denominator_jet->Fill(jet->eta(), jet->phi());
406  int index = itBX - std::min(0, 1 - (int)bxrange_%2 - (int)std::floor(bxrange_/2.)); // the correlation from itBX to respective index of the vector
407  jet_eta_phi.at(index)->Fill(jet->eta(), jet->phi());
408  }
409  }
410 
411  // Filling eta-phi map for egamma for BX=-2,-1,0,+1,+2
412  for (int itBX = EGammaBxCollection->getFirstBX(); itBX <= EGammaBxCollection->getLastBX(); ++itBX) {
413  for (l1t::EGammaBxCollection::const_iterator egamma = EGammaBxCollection->begin(itBX); egamma != EGammaBxCollection->end(itBX); ++egamma) {
414  denominator_egamma->Fill(egamma->eta(), egamma->phi());
415  int index = itBX - std::min(0, 1 - (int)bxrange_%2 - (int)std::floor(bxrange_/2.)); // the correlation from itBX to respective index of the vector
416  egamma_eta_phi.at(index)->Fill(egamma->eta(), egamma->phi());
417  }
418  }
419 
420  // Filling eta-phi map for tau for BX=-2,-1,0,+1,+2
421  for (int itBX = TauBxCollection->getFirstBX(); itBX <= TauBxCollection->getLastBX(); ++itBX) {
422  for (l1t::TauBxCollection::const_iterator tau = TauBxCollection->begin(itBX); tau != TauBxCollection->end(itBX); ++tau) {
423  denominator_tau->Fill(tau->eta(), tau->phi());
424  int index = itBX - std::min(0, 1 - (int)bxrange_%2 - (int)std::floor(bxrange_/2.)); // the correlation from itBX to respective index of the vector
425  tau_eta_phi.at(index)->Fill(tau->eta(), tau->phi());
426  }
427  }
428 
429  // Filling eta-phi map for etsum for BX=-2,-1,0,+1,+2
430  for (int itBX = EtSumBxCollection->getFirstBX(); itBX <= EtSumBxCollection->getLastBX(); ++itBX) {
431  for (l1t::EtSumBxCollection::const_iterator EtSum = EtSumBxCollection->begin(itBX); EtSum != EtSumBxCollection->end(itBX); ++EtSum) {
432  int index = itBX - std::min(0, 1 - (int)bxrange_%2 - (int)std::floor(bxrange_/2.)); // the correlation from itBX to respective index of the vector
433  if (l1t::EtSum::EtSumType::kMissingEt == EtSum->getType()) {
434  etsum_eta_phi_MET.at(index)->Fill(EtSum->phi());
435  denominator_etsum_MET->Fill(EtSum->phi());
436  }
437  else if (l1t::EtSum::EtSumType::kMissingEtHF == EtSum->getType()) {
438  etsum_eta_phi_METHF.at(index)->Fill(EtSum->phi());
439  denominator_etsum_METHF->Fill(EtSum->phi());
440 
441  }
442  else if(l1t::EtSum::EtSumType::kMissingHt == EtSum->getType()) {
443  etsum_eta_phi_MHT.at(index)->Fill(EtSum->phi());
444  denominator_etsum_MHT->Fill(EtSum->phi());
445 
446  }
447  else if(l1t::EtSum::EtSumType::kMissingHtHF == EtSum->getType()) {
448  etsum_eta_phi_MHTHF.at(index)->Fill(EtSum->phi());
449  denominator_etsum_MHTHF->Fill(EtSum->phi());
450  }
451  }
452  }
453 
454  for(GlobalAlgBlkBxCollection::const_iterator itr = uGtAlgs->begin(0); itr != uGtAlgs->end(0); ++itr) {
455  if((algoBitFirstBxInTrain_ != -1 && algoBitLastBxInTrain_ != -1) && (itr->getAlgoDecisionInitial(algoBitFirstBxInTrain_) && itr->getAlgoDecisionInitial(algoBitLastBxInTrain_))) {
456 // for(GlobalAlgBlkBxCollection::const_iterator itr = uGtAlgs->begin(0); itr != uGtAlgs->end(0); ++itr)
457  for (int itBX = MuonBxCollection->getFirstBX(); itBX <= MuonBxCollection->getLastBX(); ++itBX) {
458  for (l1t::MuonBxCollection::const_iterator muon = MuonBxCollection->begin(itBX); muon != MuonBxCollection->end(itBX); ++muon) { // Starting with Muons
459  denominator_muons_isolated->Fill(muon->eta(),muon->phi());
460  int index = itBX - std::min(0, 1 - (int)bxrange_%2 - (int)std::floor(bxrange_/2.)); // the correlation from itBX to respective index of the vector
461  muons_eta_phi_isolated.at(index)->Fill(muon->eta(),muon->phi());
462  }
463  }
464  for (int itBX = JetBxCollection->getFirstBX(); itBX <= JetBxCollection->getLastBX(); ++itBX) {
465  for (l1t::JetBxCollection::const_iterator jet = JetBxCollection->begin(itBX); jet != JetBxCollection->end(itBX); ++jet) {
466  denominator_jet_isolated->Fill(jet->eta(), jet->phi());
467  int index = itBX - std::min(0, 1 - (int)bxrange_%2 - (int)std::floor(bxrange_/2.)); // the correlation from itBX to respective index of the vector
468  jet_eta_phi_isolated.at(index)->Fill(jet->eta(), jet->phi());
469  }
470  }
471  for (int itBX = EGammaBxCollection->getFirstBX(); itBX <= EGammaBxCollection->getLastBX(); ++itBX) {
472  for (l1t::EGammaBxCollection::const_iterator egamma = EGammaBxCollection->begin(itBX); egamma != EGammaBxCollection->end(itBX); ++egamma) {
473  denominator_egamma_isolated->Fill(egamma->eta(), egamma->phi());
474  int index = itBX - std::min(0, 1 - (int)bxrange_%2 - (int)std::floor(bxrange_/2.)); // the correlation from itBX to respective index of the vector
475  egamma_eta_phi_isolated.at(index)->Fill(egamma->eta(), egamma->phi());
476  }
477  }
478  for (int itBX = TauBxCollection->getFirstBX(); itBX <= TauBxCollection->getLastBX(); ++itBX) {
479  for (l1t::TauBxCollection::const_iterator tau = TauBxCollection->begin(itBX); tau != TauBxCollection->end(itBX); ++tau) {
480  denominator_tau_isolated->Fill(tau->eta(), tau->phi());
481  int index = itBX - std::min(0, 1 - (int)bxrange_%2 - (int)std::floor(bxrange_/2.)); // the correlation from itBX to respective index of the vector
482  tau_eta_phi_isolated.at(index)->Fill(tau->eta(), tau->phi());
483  }
484  }
485  for (int itBX = EtSumBxCollection->getFirstBX(); itBX <= EtSumBxCollection->getLastBX(); ++itBX) {
486  for (l1t::EtSumBxCollection::const_iterator EtSum = EtSumBxCollection->begin(itBX); EtSum != EtSumBxCollection->end(itBX); ++EtSum) {
487  int index = itBX - std::min(0, 1 - (int)bxrange_%2 - (int)std::floor(bxrange_/2.)); // the correlation from itBX to respective index of the vector
488  if (l1t::EtSum::EtSumType::kMissingEt == EtSum->getType()) {
489  etsum_eta_phi_MET_isolated.at(index)->Fill(EtSum->phi());
490  denominator_etsum_isolated_MET->Fill(EtSum->phi());
491 
492  }
493  else if (l1t::EtSum::EtSumType::kMissingEtHF == EtSum->getType()) {
494  etsum_eta_phi_METHF_isolated.at(index)->Fill(EtSum->phi());
495  denominator_etsum_isolated_METHF->Fill(EtSum->phi());
496 
497  }
498  else if(l1t::EtSum::EtSumType::kMissingHt == EtSum->getType()) {
499  etsum_eta_phi_MHT_isolated.at(index)->Fill(EtSum->phi());
500  denominator_etsum_isolated_MHT->Fill(EtSum->phi());
501  }
502  else if(l1t::EtSum::EtSumType::kMissingHtHF == EtSum->getType()) {
503  etsum_eta_phi_MHTHF_isolated.at(index)->Fill(EtSum->phi());
504  denominator_etsum_isolated_MHTHF->Fill(EtSum->phi());
505 
506  }
507  }
508  }
509  }
510 
511  //if(itr->getAlgoDecisionInitial(algoBitFirstBxInTrain_))
512  if(algoBitFirstBxInTrain_ != -1 && itr->getAlgoDecisionInitial(algoBitFirstBxInTrain_)) { // Filling eta-phi map for all objects for first bunch each BX
513  for (int itBX = MuonBxCollection->getFirstBX(); itBX <= 0 ; ++itBX) {
514  for (l1t::MuonBxCollection::const_iterator muon = MuonBxCollection->begin(itBX); muon != MuonBxCollection->end(itBX); ++muon) { // Starting with Muons
515  denominator_muons_firstbunch->Fill(muon->eta(),muon->phi());
516  int index = itBX - std::min(0, 1 - (int)bxrange_%2 - (int)std::floor(bxrange_/2.)); // the correlation from itBX to respective index of the vector
517  muons_eta_phi_firstbunch.at(index)->Fill(muon->eta(),muon->phi());
518  }
519  }
520  for (int itBX = JetBxCollection->getFirstBX(); itBX <= 0; ++itBX) {
521  for (l1t::JetBxCollection::const_iterator jet = JetBxCollection->begin(itBX); jet != JetBxCollection->end(itBX); ++jet) {
522  denominator_jet_firstbunch->Fill(jet->eta(), jet->phi());
523  int index = itBX - std::min(0, 1 - (int)bxrange_%2 - (int)std::floor(bxrange_/2.)); // the correlation from itBX to respective index of the vector
524  jet_eta_phi_firstbunch.at(index)->Fill(jet->eta(), jet->phi());
525  }
526  }
527  for (int itBX = EGammaBxCollection->getFirstBX(); itBX <= 0; ++itBX) {
528  for (l1t::EGammaBxCollection::const_iterator egamma = EGammaBxCollection->begin(itBX); egamma != EGammaBxCollection->end(itBX); ++egamma) {
529  denominator_egamma_firstbunch->Fill(egamma->eta(), egamma->phi());
530  int index = itBX - std::min(0, 1 - (int)bxrange_%2 - (int)std::floor(bxrange_/2.)); // the correlation from itBX to respective index of the vector
531  egamma_eta_phi_firstbunch.at(index)->Fill(egamma->eta(), egamma->phi());
532  }
533  }
534  for (int itBX = TauBxCollection->getFirstBX(); itBX <= 0; ++itBX) {
535  for (l1t::TauBxCollection::const_iterator tau = TauBxCollection->begin(itBX); tau != TauBxCollection->end(itBX); ++tau) {
536  denominator_tau_firstbunch->Fill(tau->eta(), tau->phi());
537  int index = itBX - std::min(0, 1 - (int)bxrange_%2 - (int)std::floor(bxrange_/2.)); // the correlation from itBX to respective index of the vector
538  tau_eta_phi_firstbunch.at(index)->Fill(tau->eta(), tau->phi());
539  }
540  }
541  for (int itBX = EtSumBxCollection->getFirstBX(); itBX <= 0; ++itBX) {
542  for (l1t::EtSumBxCollection::const_iterator EtSum = EtSumBxCollection->begin(itBX); EtSum != EtSumBxCollection->end(itBX); ++EtSum) {
543  int index = itBX - std::min(0, 1 - (int)bxrange_%2 - (int)std::floor(bxrange_/2.)); // the correlation from itBX to respective index of the vector
544  if (l1t::EtSum::EtSumType::kMissingEt == EtSum->getType()) {
545  etsum_eta_phi_MET_firstbunch.at(index)->Fill(EtSum->phi());
546  denominator_etsum_firstbunch_MET->Fill(EtSum->phi());
547 
548  }
549  else if (l1t::EtSum::EtSumType::kMissingEtHF == EtSum->getType()) {
550  etsum_eta_phi_METHF_firstbunch.at(index)->Fill(EtSum->phi());
551  denominator_etsum_firstbunch_METHF->Fill(EtSum->phi());
552 
553  }
554  else if(l1t::EtSum::EtSumType::kMissingHt == EtSum->getType()) {
555  etsum_eta_phi_MHT_firstbunch.at(index)->Fill(EtSum->phi());
556  denominator_etsum_firstbunch_MHT->Fill(EtSum->phi());
557 
558  }
559  else if(l1t::EtSum::EtSumType::kMissingHtHF == EtSum->getType()) {
560  etsum_eta_phi_MHTHF_firstbunch.at(index)->Fill(EtSum->phi());
561  denominator_etsum_firstbunch_MHTHF->Fill(EtSum->phi());
562 
563  }
564  }
565  }
566  }
567  //if(itr->getAlgoDecisionInitial(algoBitLastBxInTrain_))
568  if(algoBitLastBxInTrain_ != -1 && itr->getAlgoDecisionInitial(algoBitLastBxInTrain_)) { // Filling eta-phi map for all objects for last bunch each BX
569  for (int itBX = 0; itBX <= MuonBxCollection->getLastBX(); ++itBX) {
570  for (l1t::MuonBxCollection::const_iterator muon = MuonBxCollection->begin(itBX); muon != MuonBxCollection->end(itBX); ++muon) { // Starting with Muons
571  denominator_muons_lastbunch->Fill(muon->eta(),muon->phi());
572  int index = itBX - std::min(0, 1 - (int)bxrange_%2 + (int)std::floor(bxrange_/2.)); // the correlation from itBX to respective index of the vector
573  muons_eta_phi_lastbunch.at(index)->Fill(muon->eta(),muon->phi());
574  }
575  }
576 
577  for (int itBX = 0; itBX <= JetBxCollection->getLastBX(); ++itBX) {
578  for (l1t::JetBxCollection::const_iterator jet = JetBxCollection->begin(itBX); jet != JetBxCollection->end(itBX); ++jet) {
579  denominator_jet_lastbunch->Fill(jet->eta(), jet->phi());
580  int index = itBX - std::min(0, 1 - (int)bxrange_%2 + (int)std::floor(bxrange_/2.)); // the correlation from itBX to respective index of the vector
581  jet_eta_phi_lastbunch.at(index)->Fill(jet->eta(), jet->phi());
582  }
583  }
584  for (int itBX = 0; itBX <= EGammaBxCollection->getLastBX(); ++itBX) {
585  for (l1t::EGammaBxCollection::const_iterator egamma = EGammaBxCollection->begin(itBX); egamma != EGammaBxCollection->end(itBX); ++egamma) {
586  denominator_egamma_lastbunch->Fill(egamma->eta(), egamma->phi());
587  int index = itBX - std::min(0, 1 - (int)bxrange_%2 + (int)std::floor(bxrange_/2.)); // the correlation from itBX to respective index of the vector
588  egamma_eta_phi_lastbunch.at(index)->Fill(egamma->eta(), egamma->phi());
589  }
590  }
591  for (int itBX = 0; itBX <= TauBxCollection->getLastBX(); ++itBX) {
592  for (l1t::TauBxCollection::const_iterator tau = TauBxCollection->begin(itBX); tau != TauBxCollection->end(itBX); ++tau) {
593  denominator_tau_lastbunch->Fill(tau->eta(), tau->phi());
594  int index = itBX - std::min(0, 1 - (int)bxrange_%2 + (int)std::floor(bxrange_/2.)); // the correlation from itBX to respective index of the vector
595  tau_eta_phi_lastbunch.at(index)->Fill(tau->eta(), tau->phi());
596  }
597  }
598  for (int itBX = 0; itBX <= EtSumBxCollection->getLastBX(); ++itBX) {
599  for (l1t::EtSumBxCollection::const_iterator EtSum = EtSumBxCollection->begin(itBX); EtSum != EtSumBxCollection->end(itBX); ++EtSum) {
600  int index = itBX - std::min(0, 1 - (int)bxrange_%2 + (int)std::floor(bxrange_/2.)); // the correlation from itBX to respective index of the vector
601  if (l1t::EtSum::EtSumType::kMissingEt == EtSum->getType()) {
602  etsum_eta_phi_MET_lastbunch.at(index)->Fill(EtSum->phi());
603  denominator_etsum_lastbunch_MET->Fill(EtSum->phi());
604 
605  }
606  else if (l1t::EtSum::EtSumType::kMissingEtHF == EtSum->getType()) {
607  etsum_eta_phi_METHF_lastbunch.at(index)->Fill(EtSum->phi());
608  denominator_etsum_lastbunch_METHF->Fill(EtSum->phi());
609 
610  }
611  else if(l1t::EtSum::EtSumType::kMissingHt == EtSum->getType()) {
612  etsum_eta_phi_MHT_lastbunch.at(index)->Fill(EtSum->phi());
613  denominator_etsum_lastbunch_MHT->Fill(EtSum->phi());
614 
615  }
616  else if(l1t::EtSum::EtSumType::kMissingHtHF == EtSum->getType()) {
617  etsum_eta_phi_MHTHF_lastbunch.at(index)->Fill(EtSum->phi());
618  denominator_etsum_lastbunch_MHTHF->Fill(EtSum->phi());
619 
620  }
621  }
622  }
623  }
624  }
625 }
626 
L1TObjectsTiming::egamma_eta_phi_isolated
std::vector< MonitorElement * > egamma_eta_phi_isolated
Definition: L1TObjectsTiming.h:81
L1TObjectsTiming::denominator_etsum_METHF
MonitorElement * denominator_etsum_METHF
Definition: L1TObjectsTiming.h:127
l1t::EGammaBxCollection
BXVector< EGamma > EGammaBxCollection
Definition: EGamma.h:11
BXVector::end
const_iterator end(int bx) const
L1TObjectsTiming::jet_eta_phi_firstbunch
std::vector< MonitorElement * > jet_eta_phi_firstbunch
Definition: L1TObjectsTiming.h:89
L1TObjectsTiming::denominator_egamma_firstbunch
MonitorElement * denominator_egamma_firstbunch
Definition: L1TObjectsTiming.h:117
L1TObjectsTiming::denominator_etsum_lastbunch_MHT
MonitorElement * denominator_etsum_lastbunch_MHT
Definition: L1TObjectsTiming.h:134
L1TObjectsTiming::stage2CaloLayer2EGammaToken_
edm::EDGetTokenT< l1t::EGammaBxCollection > stage2CaloLayer2EGammaToken_
Definition: L1TObjectsTiming.h:49
L1TObjectsTiming::denominator_muons_isolated
MonitorElement * denominator_muons_isolated
Definition: L1TObjectsTiming.h:108
L1TObjectsTiming::beginLuminosityBlock
void beginLuminosityBlock(const edm::LuminosityBlock &, const edm::EventSetup &) override
Definition: L1TObjectsTiming.cc:91
L1TObjectsTiming::l1tStage2uGtProducer_
edm::EDGetTokenT< GlobalAlgBlkBxCollection > l1tStage2uGtProducer_
Definition: L1TObjectsTiming.h:53
mps_fire.i
i
Definition: mps_fire.py:269
edm::ParameterSetDescription::addUntracked
ParameterDescriptionBase * addUntracked(U const &iLabel, T const &value)
Definition: ParameterSetDescription.h:94
L1TObjectsTiming::etsum_eta_phi_MET_isolated
std::vector< MonitorElement * > etsum_eta_phi_MET_isolated
Definition: L1TObjectsTiming.h:83
L1TObjectsTiming::etsum_eta_phi_METHF_isolated
std::vector< MonitorElement * > etsum_eta_phi_METHF_isolated
Definition: L1TObjectsTiming.h:84
L1TObjectsTiming::muons_eta_phi_isolated
std::vector< MonitorElement * > muons_eta_phi_isolated
Definition: L1TObjectsTiming.h:79
L1TObjectsTiming::denominator_etsum_MET
MonitorElement * denominator_etsum_MET
Definition: L1TObjectsTiming.h:123
L1TObjectsTiming::egamma_eta_phi
std::vector< MonitorElement * > egamma_eta_phi
Definition: L1TObjectsTiming.h:72
L1TObjectsTiming::denominator_etsum_isolated_MHTHF
MonitorElement * denominator_etsum_isolated_MHTHF
Definition: L1TObjectsTiming.h:136
edm::Event::getByToken
bool getByToken(EDGetToken token, Handle< PROD > &result) const
Definition: Event.h:519
L1TObjectsTiming::denominator_muons_firstbunch
MonitorElement * denominator_muons_firstbunch
Definition: L1TObjectsTiming.h:109
L1TObjectsTiming::jet_eta_phi_lastbunch
std::vector< MonitorElement * > jet_eta_phi_lastbunch
Definition: L1TObjectsTiming.h:99
AlCaHLTBitMon_QueryRunRegistry.string
string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:255
L1TObjectsTiming::dqmBeginRun
void dqmBeginRun(const edm::Run &, const edm::EventSetup &) override
Definition: L1TObjectsTiming.cc:77
L1TObjectsTiming::tau_eta_phi_lastbunch
std::vector< MonitorElement * > tau_eta_phi_lastbunch
Definition: L1TObjectsTiming.h:101
L1TObjectsTiming::analyze
void analyze(const edm::Event &, const edm::EventSetup &) override
Definition: L1TObjectsTiming.cc:367
L1TObjectsTiming::denominator_etsum_MHT
MonitorElement * denominator_etsum_MHT
Definition: L1TObjectsTiming.h:131
L1TObjectsTiming::bxrange_
const unsigned int bxrange_
Definition: L1TObjectsTiming.h:66
L1TObjectsTiming::denominator_egamma_isolated
MonitorElement * denominator_egamma_isolated
Definition: L1TObjectsTiming.h:116
l1t
delete x;
Definition: CaloConfig.h:22
MillePedeFileConverter_cfg.e
e
Definition: MillePedeFileConverter_cfg.py:37
l1t::TauBxCollection
BXVector< Tau > TauBxCollection
Definition: Tau.h:11
L1TObjectsTiming::etsum_eta_phi_MET_firstbunch
std::vector< MonitorElement * > etsum_eta_phi_MET_firstbunch
Definition: L1TObjectsTiming.h:93
L1TObjectsTiming::denominator_etsum_isolated_METHF
MonitorElement * denominator_etsum_isolated_METHF
Definition: L1TObjectsTiming.h:128
L1TObjectsTiming::denominator_tau_isolated
MonitorElement * denominator_tau_isolated
Definition: L1TObjectsTiming.h:120
L1TObjectsTiming::etsum_eta_phi_MHT_lastbunch
std::vector< MonitorElement * > etsum_eta_phi_MHT_lastbunch
Definition: L1TObjectsTiming.h:104
L1TObjectsTiming::denominator_etsum_firstbunch_MET
MonitorElement * denominator_etsum_firstbunch_MET
Definition: L1TObjectsTiming.h:125
MonitorElement::Fill
void Fill(long long x)
Definition: MonitorElement.h:166
L1TObjectsTiming::denominator_etsum_firstbunch_MHT
MonitorElement * denominator_etsum_firstbunch_MHT
Definition: L1TObjectsTiming.h:133
L1TObjectsTiming::muons_eta_phi_lastbunch
std::vector< MonitorElement * > muons_eta_phi_lastbunch
Definition: L1TObjectsTiming.h:98
L1TObjectsTiming::denominator_etsum_MHTHF
MonitorElement * denominator_etsum_MHTHF
Definition: L1TObjectsTiming.h:135
L1TObjectsTiming::etsum_eta_phi_MHT
std::vector< MonitorElement * > etsum_eta_phi_MHT
Definition: L1TObjectsTiming.h:76
L1TObjectsTiming::bookHistograms
void bookHistograms(DQMStore::IBooker &, const edm::Run &, const edm::EventSetup &) override
Definition: L1TObjectsTiming.cc:93
l1t::EtSumBxCollection
BXVector< EtSum > EtSumBxCollection
Definition: EtSum.h:11
L1TObjectsTiming::denominator_egamma
MonitorElement * denominator_egamma
Definition: L1TObjectsTiming.h:115
L1TObjectsTiming::etsum_eta_phi_MET_lastbunch
std::vector< MonitorElement * > etsum_eta_phi_MET_lastbunch
Definition: L1TObjectsTiming.h:102
Muon
Definition: Muon.py:1
L1TObjectsTiming::fillDescriptions
static void fillDescriptions(edm::ConfigurationDescriptions &descriptions)
Definition: L1TObjectsTiming.cc:62
L1TObjectsTiming::denominator_tau_lastbunch
MonitorElement * denominator_tau_lastbunch
Definition: L1TObjectsTiming.h:122
L1TObjectsTiming::jet_eta_phi
std::vector< MonitorElement * > jet_eta_phi
Definition: L1TObjectsTiming.h:71
L1TObjectsTiming::etsum_eta_phi_METHF_firstbunch
std::vector< MonitorElement * > etsum_eta_phi_METHF_firstbunch
Definition: L1TObjectsTiming.h:94
L1TObjectsTiming::etsum_eta_phi_MHTHF_isolated
std::vector< MonitorElement * > etsum_eta_phi_MHTHF_isolated
Definition: L1TObjectsTiming.h:86
DQMStore::IBooker::book1D
MonitorElement * book1D(Args &&...args)
Definition: DQMStore.h:118
L1TObjectsTiming::L1TObjectsTiming
L1TObjectsTiming(const edm::ParameterSet &ps)
Definition: L1TObjectsTiming.cc:4
L1TObjectsTiming::ugmtMuonToken_
edm::EDGetTokenT< l1t::MuonBxCollection > ugmtMuonToken_
Definition: L1TObjectsTiming.h:47
L1TObjectsTiming::denominator_tau_firstbunch
MonitorElement * denominator_tau_firstbunch
Definition: L1TObjectsTiming.h:121
L1TObjectsTiming::egamma_eta_phi_lastbunch
std::vector< MonitorElement * > egamma_eta_phi_lastbunch
Definition: L1TObjectsTiming.h:100
L1TObjectsTiming::denominator_etsum_lastbunch_MHTHF
MonitorElement * denominator_etsum_lastbunch_MHTHF
Definition: L1TObjectsTiming.h:138
L1TObjectsTiming::denominator_muons
MonitorElement * denominator_muons
Definition: L1TObjectsTiming.h:107
L1TObjectsTiming::muons_eta_phi_firstbunch
std::vector< MonitorElement * > muons_eta_phi_firstbunch
Definition: L1TObjectsTiming.h:88
min
T min(T a, T b)
Definition: MathUtil.h:58
edm::ParameterSetDescription::add
ParameterDescriptionBase * add(U const &iLabel, T const &value)
Definition: ParameterSetDescription.h:89
nanoDQM_cfi.Muon
Muon
Definition: nanoDQM_cfi.py:268
L1TObjectsTiming::etsum_eta_phi_MHTHF_lastbunch
std::vector< MonitorElement * > etsum_eta_phi_MHTHF_lastbunch
Definition: L1TObjectsTiming.h:105
L1TObjectsTiming::stage2CaloLayer2EtSumToken_
edm::EDGetTokenT< l1t::EtSumBxCollection > stage2CaloLayer2EtSumToken_
Definition: L1TObjectsTiming.h:51
L1TObjectsTiming::stage2CaloLayer2JetToken_
edm::EDGetTokenT< l1t::JetBxCollection > stage2CaloLayer2JetToken_
Definition: L1TObjectsTiming.h:48
L1TObjectsTiming::stage2CaloLayer2TauToken_
edm::EDGetTokenT< l1t::TauBxCollection > stage2CaloLayer2TauToken_
Definition: L1TObjectsTiming.h:50
electrons_cff.bool
bool
Definition: electrons_cff.py:231
L1TObjectsTiming::tau_eta_phi_isolated
std::vector< MonitorElement * > tau_eta_phi_isolated
Definition: L1TObjectsTiming.h:82
L1TObjectsTiming::etsum_eta_phi_METHF
std::vector< MonitorElement * > etsum_eta_phi_METHF
Definition: L1TObjectsTiming.h:75
L1TObjectsTiming::etsum_eta_phi_MHTHF
std::vector< MonitorElement * > etsum_eta_phi_MHTHF
Definition: L1TObjectsTiming.h:77
l1t::MuonBxCollection
BXVector< Muon > MuonBxCollection
Definition: Muon.h:11
L1TObjectsTiming::etsum_eta_phi_MHTHF_firstbunch
std::vector< MonitorElement * > etsum_eta_phi_MHTHF_firstbunch
Definition: L1TObjectsTiming.h:96
L1TObjectsTiming::tau_eta_phi
std::vector< MonitorElement * > tau_eta_phi
Definition: L1TObjectsTiming.h:73
BXVector::getFirstBX
int getFirstBX() const
DQMStore::IBooker::setCurrentFolder
void setCurrentFolder(const std::string &fullpath)
Definition: DQMStore.cc:279
l1t::JetBxCollection
BXVector< Jet > JetBxCollection
Definition: Jet.h:11
L1TObjectsTiming::egamma_eta_phi_firstbunch
std::vector< MonitorElement * > egamma_eta_phi_firstbunch
Definition: L1TObjectsTiming.h:90
L1TObjectsTiming::denominator_etsum_firstbunch_METHF
MonitorElement * denominator_etsum_firstbunch_METHF
Definition: L1TObjectsTiming.h:129
L1TObjectsTiming::denominator_jet_lastbunch
MonitorElement * denominator_jet_lastbunch
Definition: L1TObjectsTiming.h:114
DQMStore::IBooker::book2D
MonitorElement * book2D(Args &&...args)
Definition: DQMStore.h:136
L1TObjectsTiming::monitorDir_
std::string monitorDir_
Definition: L1TObjectsTiming.h:55
L1TObjectsTiming::denominator_jet
MonitorElement * denominator_jet
Definition: L1TObjectsTiming.h:111
L1TObjectsTiming::etsum_eta_phi_MHT_isolated
std::vector< MonitorElement * > etsum_eta_phi_MHT_isolated
Definition: L1TObjectsTiming.h:85
edm::ConfigurationDescriptions::add
void add(std::string const &label, ParameterSetDescription const &psetDescription)
Definition: ConfigurationDescriptions.cc:67
L1TObjectsTiming::etsum_eta_phi_METHF_lastbunch
std::vector< MonitorElement * > etsum_eta_phi_METHF_lastbunch
Definition: L1TObjectsTiming.h:103
L1TObjectsTiming::denominator_muons_lastbunch
MonitorElement * denominator_muons_lastbunch
Definition: L1TObjectsTiming.h:110
L1TObjectsTiming::algoNameLastBxInTrain_
const std::string algoNameLastBxInTrain_
Definition: L1TObjectsTiming.h:65
L1TObjectsTiming::etsum_eta_phi_MET
std::vector< MonitorElement * > etsum_eta_phi_MET
Definition: L1TObjectsTiming.h:74
L1TObjectsTiming::denominator_tau
MonitorElement * denominator_tau
Definition: L1TObjectsTiming.h:119
L1TObjectsTiming::~L1TObjectsTiming
~L1TObjectsTiming() override
Definition: L1TObjectsTiming.cc:60
L1TObjectsTiming::jet_eta_phi_isolated
std::vector< MonitorElement * > jet_eta_phi_isolated
Definition: L1TObjectsTiming.h:80
L1TObjectsTiming::denominator_egamma_lastbunch
MonitorElement * denominator_egamma_lastbunch
Definition: L1TObjectsTiming.h:118
edm
HLT enums.
Definition: AlignableModifier.h:17
L1TObjectsTiming::denominator_etsum_isolated_MET
MonitorElement * denominator_etsum_isolated_MET
Definition: L1TObjectsTiming.h:124
L1TObjectsTiming::denominator_etsum_firstbunch_MHTHF
MonitorElement * denominator_etsum_firstbunch_MHTHF
Definition: L1TObjectsTiming.h:137
BXVector::getLastBX
int getLastBX() const
L1TObjectsTiming::muons_eta_phi
std::vector< MonitorElement * > muons_eta_phi
Definition: L1TObjectsTiming.h:70
L1TObjectsTiming::denominator_etsum_lastbunch_METHF
MonitorElement * denominator_etsum_lastbunch_METHF
Definition: L1TObjectsTiming.h:130
L1TObjectsTiming::denominator_etsum_isolated_MHT
MonitorElement * denominator_etsum_isolated_MHT
Definition: L1TObjectsTiming.h:132
L1TObjectsTiming::algoNameFirstBxInTrain_
const std::string algoNameFirstBxInTrain_
Definition: L1TObjectsTiming.h:64
L1TObjectsTiming::etsum_eta_phi_MHT_firstbunch
std::vector< MonitorElement * > etsum_eta_phi_MHT_firstbunch
Definition: L1TObjectsTiming.h:95
BXVector::begin
const_iterator begin(int bx) const
MonitorElement::setAxisTitle
void setAxisTitle(const std::string &title, int axis=1)
set x-, y- or z-axis title (axis=1, 2, 3 respectively)
Definition: MonitorElement.cc:973
L1TObjectsTiming::gtUtil_
std::shared_ptr< l1t::L1TGlobalUtil > gtUtil_
Definition: L1TObjectsTiming.h:59
L1TObjectsTiming::denominator_jet_firstbunch
MonitorElement * denominator_jet_firstbunch
Definition: L1TObjectsTiming.h:113
L1TObjectsTiming::denominator_jet_isolated
MonitorElement * denominator_jet_isolated
Definition: L1TObjectsTiming.h:112
edm::Run
Definition: Run.h:43
L1TObjectsTiming::denominator_etsum_lastbunch_MET
MonitorElement * denominator_etsum_lastbunch_MET
Definition: L1TObjectsTiming.h:126
BXVector< Muon >::const_iterator
std::vector< Muon >::const_iterator const_iterator
Definition: BXVector.h:20
L1TObjectsTiming::tau_eta_phi_firstbunch
std::vector< MonitorElement * > tau_eta_phi_firstbunch
Definition: L1TObjectsTiming.h:91
Generated for CMSSW Reference Manual by doxygen 1.8.11