CMS 3D CMS Logo

 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
L1TGCT.cc
Go to the documentation of this file.
4 
5 // Trigger Headers
6 
7 // GCT and RCT data formats
11 
12 using namespace edm;
13 
14 // Define statics for bins etc.
15 const unsigned int JETETABINS = 22;
16 const float JETETAMIN = -0.5;
17 const float JETETAMAX = 21.5;
18 
19 const unsigned int EMETABINS = 22;
20 const float EMETAMIN = -0.5;
21 const float EMETAMAX = 21.5;
22 
23 const unsigned int METPHIBINS = 72;
24 const float METPHIMIN = -0.5;
25 const float METPHIMAX = 71.5;
26 
27 const unsigned int MHTPHIBINS = 18;
28 const float MHTPHIMIN = -0.5;
29 const float MHTPHIMAX = 17.5;
30 
31 const unsigned int PHIBINS = 18;
32 const float PHIMIN = -0.5;
33 const float PHIMAX = 17.5;
34 
35 const unsigned int OFBINS = 2;
36 const float OFMIN = -0.5;
37 const float OFMAX = 1.5;
38 
39 const unsigned int BXBINS = 5;
40 const float BXMIN = -2.5;
41 const float BXMAX = 2.5;
42 
43 // Bins for 3, 5, 6, 7, 10 and 12 bits
44 const unsigned int R3BINS = 8;
45 const float R3MIN = -0.5;
46 const float R3MAX = 7.5;
47 const unsigned int R5BINS = 32;
48 const float R5MIN = -0.5;
49 const float R5MAX = 31.5;
50 const unsigned int R6BINS = 64;
51 const float R6MIN = -0.5;
52 const float R6MAX = 63.5;
53 const unsigned int R7BINS = 128;
54 const float R7MIN = -0.5;
55 const float R7MAX = 127.5;
56 const unsigned int R12BINS = 4096;
57 const float R12MIN = -0.5;
58 const float R12MAX = 4095.5;
59 
61  monitorDir_(ps.getUntrackedParameter<std::string>("monitorDir","")),
62  gctCenJetsSource_(ps.getParameter<edm::InputTag>("gctCentralJetsSource")),
63  gctForJetsSource_(ps.getParameter<edm::InputTag>("gctForwardJetsSource")),
64  gctTauJetsSource_(ps.getParameter<edm::InputTag>("gctTauJetsSource")),
65  gctIsoTauJetsSource_(ps.getParameter<edm::InputTag>("gctIsoTauJetsSource")),
66  gctEnergySumsSource_(ps.getParameter<edm::InputTag>("gctEnergySumsSource")),
67  gctIsoEmSource_(ps.getParameter<edm::InputTag>("gctIsoEmSource")),
68  gctNonIsoEmSource_(ps.getParameter<edm::InputTag>("gctNonIsoEmSource")),
69  m_stage1_layer2_(ps.getParameter<bool>("stage1_layer2_")),
70  filterTriggerType_ (ps.getParameter< int >("filterTriggerType"))
71 {
72 
73  // verbosity switch
74  verbose_ = ps.getUntrackedParameter < bool > ("verbose", false);
75 
76  if (verbose_)
77  edm::LogInfo("L1TGCT") << "L1TGCT: constructor...." << std::endl;
78 
79  outputFile_ = ps.getUntrackedParameter < std::string > ("outputFile", "");
80  if (outputFile_.size() != 0) {
81  edm::LogInfo("L1TGCT") << "L1T Monitoring histograms will be saved to "
82  << outputFile_ << std::endl;
83  }
84 
85  bool disable = ps.getUntrackedParameter<bool>("disableROOToutput", false);
86  if(disable){
87  outputFile_="";
88  }
89 
90  //set Token(-s)
91  gctIsoEmSourceToken_ = consumes<L1GctEmCandCollection>(ps.getParameter<edm::InputTag>("gctIsoEmSource"));
92  gctNonIsoEmSourceToken_ = consumes<L1GctEmCandCollection>(ps.getParameter<edm::InputTag>("gctNonIsoEmSource"));
93  gctCenJetsSourceToken_ = consumes<L1GctJetCandCollection>(ps.getParameter<edm::InputTag>("gctCentralJetsSource"));
94  gctForJetsSourceToken_ = consumes<L1GctJetCandCollection>(ps.getParameter<edm::InputTag>("gctForwardJetsSource"));
95  gctTauJetsSourceToken_ = consumes<L1GctJetCandCollection>(ps.getParameter<edm::InputTag>("gctTauJetsSource"));
96  if(m_stage1_layer2_ == true){
97  gctIsoTauJetsSourceToken_=consumes<L1GctJetCandCollection>(ps.getParameter<edm::InputTag>("gctIsoTauJetsSource"));
98  }
99  gctEnergySumsSourceToken_ = consumes<L1GctHFRingEtSumsCollection>(ps.getParameter<edm::InputTag>("gctEnergySumsSource"));
100  l1HFCountsToken_ = consumes<L1GctHFBitCountsCollection>(ps.getParameter<edm::InputTag>("gctEnergySumsSource"));
101  l1EtMissToken_ = consumes<L1GctEtMissCollection>(ps.getParameter<edm::InputTag>("gctEnergySumsSource"));
102  l1HtMissToken_ = consumes<L1GctHtMissCollection>(ps.getParameter<edm::InputTag>("gctEnergySumsSource"));
103  l1EtHadToken_ = consumes<L1GctEtHadCollection>(ps.getParameter<edm::InputTag>("gctEnergySumsSource"));
104  l1EtTotalToken_ = consumes<L1GctEtTotalCollection>(ps.getParameter<edm::InputTag>("gctEnergySumsSource"));
105 }
106 
108 {
109 }
110 
112 {
113 
114  nev_ = 0;
115 
116  ibooker.setCurrentFolder(monitorDir_);
117 
118  runId_ = ibooker.bookInt("iRun");
119  runId_->Fill(-1);
120  lumisecId_ = ibooker.bookInt("iLumiSection");
121  lumisecId_->Fill(-1);
122 
123  triggerType_ =ibooker.book1D("TriggerType", "TriggerType", 17, -0.5, 16.5);
124 
125  l1GctAllJetsEtEtaPhi_ = ibooker.book2D("AllJetsEtEtaPhi", "CENTRAL AND FORWARD JET E_{T}",JETETABINS, JETETAMIN, JETETAMAX,PHIBINS, PHIMIN, PHIMAX);
126  l1GctCenJetsEtEtaPhi_ = ibooker.book2D("CenJetsEtEtaPhi", "CENTRAL JET E_{T}",JETETABINS, JETETAMIN, JETETAMAX, PHIBINS, PHIMIN, PHIMAX);
127  l1GctForJetsEtEtaPhi_ = ibooker.book2D("ForJetsEtEtaPhi", "FORWARD JET E_{T}", JETETABINS, JETETAMIN, JETETAMAX, PHIBINS, PHIMIN, PHIMAX);
128  l1GctTauJetsEtEtaPhi_ = ibooker.book2D("TauJetsEtEtaPhi", "TAU JET E_{T}", EMETABINS, EMETAMIN, EMETAMAX, PHIBINS, PHIMIN, PHIMAX);
129  if (m_stage1_layer2_ == true){
130  l1GctIsoTauJetsEtEtaPhi_ = ibooker.book2D("IsoTauJetsEtEtaPhi", "ISOTAU JET E_{T}", EMETABINS, EMETAMIN, EMETAMAX, PHIBINS, PHIMIN, PHIMAX);
131  }
132  l1GctIsoEmRankEtaPhi_ = ibooker.book2D("IsoEmRankEtaPhi", "ISO EM E_{T}", EMETABINS, EMETAMIN, EMETAMAX, PHIBINS, PHIMIN, PHIMAX);
133  l1GctNonIsoEmRankEtaPhi_ = ibooker.book2D("NonIsoEmRankEtaPhi", "NON-ISO EM E_{T}", EMETABINS, EMETAMIN, EMETAMAX,PHIBINS, PHIMIN, PHIMAX);
134  l1GctAllJetsOccEtaPhi_ = ibooker.book2D("AllJetsOccEtaPhi", "CENTRAL AND FORWARD JET OCCUPANCY", JETETABINS, JETETAMIN, JETETAMAX, PHIBINS, PHIMIN, PHIMAX);
135  l1GctCenJetsOccEtaPhi_ = ibooker.book2D("CenJetsOccEtaPhi", "CENTRAL JET OCCUPANCY", JETETABINS, JETETAMIN, JETETAMAX, PHIBINS, PHIMIN, PHIMAX);
136  l1GctForJetsOccEtaPhi_ = ibooker.book2D("ForJetsOccEtaPhi", "FORWARD JET OCCUPANCY",JETETABINS, JETETAMIN, JETETAMAX, PHIBINS, PHIMIN, PHIMAX);
137  l1GctTauJetsOccEtaPhi_ = ibooker.book2D("TauJetsOccEtaPhi", "TAU JET OCCUPANCY", EMETABINS, EMETAMIN, EMETAMAX, PHIBINS, PHIMIN, PHIMAX);
138  if (m_stage1_layer2_ == true){
139  l1GctIsoTauJetsOccEtaPhi_ = ibooker.book2D("IsoTauJetsOccEtaPhi", "ISOTAU JET OCCUPANCY", EMETABINS, EMETAMIN, EMETAMAX, PHIBINS, PHIMIN, PHIMAX);
140  }
141  l1GctIsoEmOccEtaPhi_ = ibooker.book2D("IsoEmOccEtaPhi", "ISO EM OCCUPANCY", EMETABINS, EMETAMIN, EMETAMAX, PHIBINS, PHIMIN, PHIMAX);
142  l1GctNonIsoEmOccEtaPhi_ = ibooker.book2D("NonIsoEmOccEtaPhi", "NON-ISO EM OCCUPANCY", EMETABINS, EMETAMIN, EMETAMAX, PHIBINS, PHIMIN, PHIMAX);
143 
144 
145  //HF Ring stuff
146 
147  l1GctHFRing1TowerCountPosEtaNegEta_ = ibooker.book2D("HFRing1TowerCountCorr", "HF RING1 TOWER COUNT CORRELATION +/- #eta", R3BINS, R3MIN, R3MAX, R3BINS, R3MIN, R3MAX);
148  l1GctHFRing2TowerCountPosEtaNegEta_ = ibooker.book2D("HFRing2TowerCountCorr", "HF RING2 TOWER COUNT CORRELATION +/- #eta", R3BINS, R3MIN, R3MAX, R3BINS, R3MIN, R3MAX);
149 
150  l1GctHFRing1TowerCountPosEta_ = ibooker.book1D("HFRing1TowerCountPosEta", "HF RING1 TOWER COUNT #eta +", R3BINS, R3MIN, R3MAX);
151  l1GctHFRing1TowerCountNegEta_ = ibooker.book1D("HFRing1TowerCountNegEta", "HF RING1 TOWER COUNT #eta -", R3BINS, R3MIN, R3MAX);
152  l1GctHFRing2TowerCountPosEta_ = ibooker.book1D("HFRing2TowerCountPosEta", "HF RING2 TOWER COUNT #eta +", R3BINS, R3MIN, R3MAX);
153  l1GctHFRing2TowerCountNegEta_ = ibooker.book1D("HFRing2TowerCountNegEta", "HF RING2 TOWER COUNT #eta -", R3BINS, R3MIN, R3MAX);
154 
155  l1GctHFRingTowerCountOccBx_ = ibooker.book2D("HFRingTowerCountOccBx", "HF RING TOWER COUNT PER BX",BXBINS, BXMIN, BXMAX, R3BINS, R3MIN, R3MAX);
156 
157  if (m_stage1_layer2_ == false){
158  l1GctHFRing1PosEtaNegEta_ = ibooker.book2D("HFRing1Corr", "HF RING1 E_{T} CORRELATION +/- #eta", R3BINS, R3MIN, R3MAX, R3BINS, R3MIN, R3MAX);
159  l1GctHFRing2PosEtaNegEta_ = ibooker.book2D("HFRing2Corr", "HF RING2 E_{T} CORRELATION +/- #eta", R3BINS, R3MIN, R3MAX, R3BINS, R3MIN, R3MAX);
160  l1GctHFRing1ETSumPosEta_ = ibooker.book1D("HFRing1ETSumPosEta", "HF RING1 E_{T} #eta +", R3BINS, R3MIN, R3MAX);
161  l1GctHFRing1ETSumNegEta_ = ibooker.book1D("HFRing1ETSumNegEta", "HF RING1 E_{T} #eta -", R3BINS, R3MIN, R3MAX);
162  l1GctHFRing2ETSumPosEta_ = ibooker.book1D("HFRing2ETSumPosEta", "HF RING2 E_{T} #eta +", R3BINS, R3MIN, R3MAX);
163  l1GctHFRing2ETSumNegEta_ = ibooker.book1D("HFRing2ETSumNegEta", "HF RING2 E_{T} #eta -", R3BINS, R3MIN, R3MAX);
164  l1GctHFRingETSumOccBx_ = ibooker.book2D("HFRingETSumOccBx", "HF RING E_{T} PER BX",BXBINS, BXMIN, BXMAX, R3BINS, R3MIN, R3MAX);
165  l1GctHFRingRatioPosEta_ = ibooker.book1D("HFRingRatioPosEta", "HF RING E_{T} RATIO #eta +", R5BINS, R5MIN, R5MAX);
166  l1GctHFRingRatioNegEta_ = ibooker.book1D("HFRingRatioNegEta", "HF RING E_{T} RATIO #eta -", R5BINS, R5MIN, R5MAX);
167  }
168 
169  if (m_stage1_layer2_ == true){
170  l1GctHFRing1PosEtaNegEta_ = ibooker.book2D("IsoTau 1 2 Corr", "IsoTau 1 IsoTau 2 E_{T} CORRELATION", R3BINS, R3MIN, R3MAX, R3BINS, R3MIN, R3MAX);
171  l1GctHFRing2PosEtaNegEta_ = ibooker.book2D("IsoTau 3 4 Corr", "IsoTau 3 IsoTau 4 CORRELATION", R3BINS, R3MIN, R3MAX, R3BINS, R3MIN, R3MAX);
172  l1GctHFRing1ETSumPosEta_ = ibooker.book1D("Iso Tau 1 Et", "Isolated Tau1 E_{T}", 9, -0.5, 8.5);
173  l1GctHFRing1ETSumNegEta_ = ibooker.book1D("Iso Tau 2 Et", "Isolated Tau2 E_{T}", 9, -0.5, 8.5);
174  l1GctHFRing2ETSumPosEta_ = ibooker.book1D("Iso Tau 3 Et", "Isolated Tau3 E_{T}", 9, -0.5, 8.5);
175  l1GctHFRing2ETSumNegEta_ = ibooker.book1D("Iso Tau 4 Et", "Isolated Tau4 E_{T}", 9, -0.5, 8.5);
176  l1GctHFRingETSumOccBx_ = ibooker.book2D("IsoTau HFRingSum OccBx", "Iso Tau PER BX",BXBINS, BXMIN, BXMAX, R3BINS, R3MIN, R3MAX);
177  l1GctHFRingRatioPosEta_ = ibooker.book1D("IsoTau Ratio 1 2", "IsoTau E_{T} RATIO", 9, -0.5, 8.5);
178  l1GctHFRingRatioNegEta_ = ibooker.book1D("IsoTau Ratio 1 2", "IsoTau E_{T} RATIO", 9, -0.5, 8.5);
179  }
180 
181  // Rank histograms
182  l1GctCenJetsRank_ = ibooker.book1D("CenJetsRank", "CENTRAL JET E_{T}", R6BINS, R6MIN, R6MAX);
183  l1GctForJetsRank_ = ibooker.book1D("ForJetsRank", "FORWARD JET E_{T}", R6BINS, R6MIN, R6MAX);
184  l1GctTauJetsRank_ = ibooker.book1D("TauJetsRank", "TAU JET E_{T}", R6BINS, R6MIN, R6MAX);
185  if (m_stage1_layer2_ == true){
186  l1GctIsoTauJetsRank_ = ibooker.book1D("IsoTauJetsRank", "ISOTAU JET E_{T}", R6BINS, R6MIN, R6MAX);
187  }
188  l1GctIsoEmRank_ = ibooker.book1D("IsoEmRank", "ISO EM E_{T}", R6BINS, R6MIN, R6MAX);
189  l1GctNonIsoEmRank_ = ibooker.book1D("NonIsoEmRank", "NON-ISO EM E_{T}", R6BINS, R6MIN, R6MAX);
190 
191  l1GctAllJetsOccRankBx_ = ibooker.book2D("AllJetsOccRankBx","ALL JETS E_{T} PER BX",BXBINS,BXMIN,BXMAX,R6BINS,R6MIN,R6MAX);
192  l1GctAllEmOccRankBx_ = ibooker.book2D("AllEmOccRankBx","ALL EM E_{T} PER BX",BXBINS,BXMIN,BXMAX,R6BINS,R6MIN,R6MAX);
193 
194  // Energy sums
195  l1GctEtMiss_ = ibooker.book1D("EtMiss", "MET", R12BINS, R12MIN, R12MAX);
196  l1GctEtMissPhi_ = ibooker.book1D("EtMissPhi", "MET #phi", METPHIBINS, METPHIMIN, METPHIMAX);
197  l1GctEtMissOf_ = ibooker.book1D("EtMissOf", "MET OVERFLOW", OFBINS, OFMIN, OFMAX);
198  l1GctEtMissOccBx_ = ibooker.book2D("EtMissOccBx","MET PER BX",BXBINS,BXMIN,BXMAX,R12BINS,R12MIN,R12MAX);
199  if (m_stage1_layer2_ == false) {
200  l1GctHtMiss_ = ibooker.book1D("HtMiss", "MHT", R7BINS, R7MIN, R7MAX);
201  l1GctHtMissPhi_ = ibooker.book1D("HtMissPhi", "MHT #phi", MHTPHIBINS, MHTPHIMIN, MHTPHIMAX);
202  l1GctHtMissOf_ = ibooker.book1D("HtMissOf", "MHT OVERFLOW", OFBINS, OFMIN, OFMAX);
203  l1GctHtMissOccBx_ = ibooker.book2D("HtMissOccBx","MHT PER BX",BXBINS,BXMIN,BXMAX,R7BINS,R7MIN,R7MAX);
204  }
205  if (m_stage1_layer2_ == true) {
206  l1GctHtMiss_ = ibooker.book1D("HtMissHtTotal", "MHTHT", R7BINS, R7MIN, R7MAX);
207  l1GctHtMissPhi_ = ibooker.book1D("HtMissHtTotal Phi", "MHTHT #phi", MHTPHIBINS, MHTPHIMIN, MHTPHIMAX);
208  l1GctHtMissOf_ = ibooker.book1D("HtMissHtTotal Of", "MHTHT OVERFLOW", OFBINS, OFMIN, OFMAX);
209  l1GctHtMissOccBx_ = ibooker.book2D("HtMissHtTotal OccBx","MHTHT PER BX",BXBINS,BXMIN,BXMAX,R7BINS,R7MIN,R7MAX);
210  }
211  l1GctEtMissHtMissCorr_ = ibooker.book2D("EtMissHtMissCorr", "MET MHT CORRELATION", R6BINS, R12MIN, R12MAX, R6BINS, R7MIN, R7MAX);
212  l1GctEtMissHtMissCorrPhi_ = ibooker.book2D("EtMissHtMissPhiCorr", "MET MHT #phi CORRELATION", METPHIBINS, METPHIMIN, METPHIMAX, MHTPHIBINS, MHTPHIMIN, MHTPHIMAX);
213  l1GctEtTotal_ = ibooker.book1D("EtTotal", "SUM E_{T}", R12BINS, R12MIN, R12MAX);
214  l1GctEtTotalOf_ = ibooker.book1D("EtTotalOf", "SUM E_{T} OVERFLOW", OFBINS, OFMIN, OFMAX);
215  l1GctEtTotalOccBx_ = ibooker.book2D("EtTotalOccBx","SUM E_{T} PER BX",BXBINS,BXMIN,BXMAX,R12BINS,R12MIN,R12MAX);
216  l1GctEtHad_ = ibooker.book1D("EtHad", "H_{T}", R12BINS, R12MIN, R12MAX);
217  l1GctEtHadOf_ = ibooker.book1D("EtHadOf", "H_{T} OVERFLOW", OFBINS, OFMIN, OFMAX);
218  l1GctEtHadOccBx_ = ibooker.book2D("EtHadOccBx","H_{T} PER BX",BXBINS,BXMIN,BXMAX,R12BINS,R12MIN,R12MAX);
219  l1GctEtTotalEtHadCorr_ = ibooker.book2D("EtTotalEtHadCorr", "Sum E_{T} H_{T} CORRELATION", R6BINS, R12MIN, R12MAX, R6BINS, R12MIN, R12MAX);
220  //}
221 }
222 
223 
224 void L1TGCT::dqmBeginRun(edm::Run const& iRrun, edm::EventSetup const& evSetup) {
225  //runId_->Fill(iRrun.id().run());
226 }
227 
229  //lumisecId_->Fill(iLumi.id().luminosityBlock());
230 }
231 
233 {
234  nev_++;
235  if (verbose_) {
236  edm::LogInfo("L1TGCT") << "L1TGCT: analyze...." << std::endl;
237  }
238 
239 
240  // filter according trigger type
241  // enum ExperimentType {
242  // Undefined = 0,
243  // PhysicsTrigger = 1,
244  // CalibrationTrigger = 2,
245  // RandomTrigger = 3,
246  // Reserved = 4,
247  // TracedEvent = 5,
248  // TestTrigger = 6,
249  // ErrorTrigger = 15
250 
251  // fill a histogram with the trigger type, for normalization fill also last bin
252  // ErrorTrigger + 1
253  double triggerType = static_cast<double> (e.experimentType()) + 0.001;
254  double triggerTypeLast = static_cast<double> (edm::EventAuxiliary::ExperimentType::ErrorTrigger)
255  + 0.001;
256  triggerType_->Fill(triggerType);
257  triggerType_->Fill(triggerTypeLast + 1);
258 
259  // filter only if trigger type is greater than 0, negative values disable filtering
260  if (filterTriggerType_ >= 0) {
261 
262  // now filter, for real data only
263  if (e.isRealData()) {
264  if (!(e.experimentType() == filterTriggerType_)) {
265 
266  edm::LogInfo("L1TGCT") << "\n Event of TriggerType "
267  << e.experimentType() << " rejected" << std::endl;
268  return;
269 
270  }
271  }
272 
273  }
274 
275  // Get all the collections
281  if(m_stage1_layer2_ == true) {
283  e.getByToken(gctIsoTauJetsSourceToken_, l1IsoTauJets);
284  }
291 
292  e.getByToken(gctIsoEmSourceToken_, l1IsoEm);
293  e.getByToken(gctNonIsoEmSourceToken_, l1NonIsoEm);
294  e.getByToken(gctCenJetsSourceToken_, l1CenJets);
295  e.getByToken(gctForJetsSourceToken_, l1ForJets);
296  e.getByToken(gctTauJetsSourceToken_, l1TauJets);
298  e.getByToken(l1HFCountsToken_, l1HFCounts);
299  e.getByToken(l1EtMissToken_, l1EtMiss);
300  e.getByToken(l1HtMissToken_, l1HtMiss);
301  e.getByToken(l1EtHadToken_, l1EtHad);
302  e.getByToken(l1EtTotalToken_, l1EtTotal);
303 
304  // Fill histograms
305 
306  // Central jets
307  if (l1CenJets.isValid()) {
308  for (L1GctJetCandCollection::const_iterator cj = l1CenJets->begin();cj != l1CenJets->end(); cj++) {
309  // only plot central BX
310  if (cj->bx()==0) {
311  l1GctCenJetsRank_->Fill(cj->rank());
312  // only plot eta and phi maps for non-zero candidates
313  if (cj->rank()) {
314  l1GctAllJetsEtEtaPhi_->Fill(cj->regionId().ieta(),cj->regionId().iphi(),cj->rank());
315  l1GctAllJetsOccEtaPhi_->Fill(cj->regionId().ieta(),cj->regionId().iphi());
316  l1GctCenJetsEtEtaPhi_->Fill(cj->regionId().ieta(),cj->regionId().iphi(),cj->rank());
317  l1GctCenJetsOccEtaPhi_->Fill(cj->regionId().ieta(),cj->regionId().iphi());
318  }
319  }
320  if (cj->rank()) l1GctAllJetsOccRankBx_->Fill(cj->bx(),cj->rank()); // for all BX
321  }
322  } else {
323  edm::LogWarning("DataNotFound") << " Could not find l1CenJets label was " << gctCenJetsSource_ ;
324  }
325 
326  // Forward jets
327  if (l1ForJets.isValid()) {
328  for (L1GctJetCandCollection::const_iterator fj = l1ForJets->begin(); fj != l1ForJets->end(); fj++) {
329  // only plot central BX
330  if (fj->bx()==0) {
331  l1GctForJetsRank_->Fill(fj->rank());
332  // only plot eta and phi maps for non-zero candidates
333  if (fj->rank()) {
334  l1GctAllJetsEtEtaPhi_->Fill(fj->regionId().ieta(),fj->regionId().iphi(),fj->rank());
335  l1GctAllJetsOccEtaPhi_->Fill(fj->regionId().ieta(),fj->regionId().iphi());
336  l1GctForJetsEtEtaPhi_->Fill(fj->regionId().ieta(),fj->regionId().iphi(),fj->rank());
337  l1GctForJetsOccEtaPhi_->Fill(fj->regionId().ieta(),fj->regionId().iphi());
338  }
339  }
340  if (fj->rank()) l1GctAllJetsOccRankBx_->Fill(fj->bx(),fj->rank()); // for all BX
341  }
342  } else {
343  edm::LogWarning("DataNotFound") << " Could not find l1ForJets label was " << gctForJetsSource_ ;
344  }
345 
346  // Tau jets
347  if (l1TauJets.isValid()) {
348  for (L1GctJetCandCollection::const_iterator tj = l1TauJets->begin(); tj != l1TauJets->end(); tj++) {
349  // only plot central BX
350  if (tj->bx()==0) {
351  l1GctTauJetsRank_->Fill(tj->rank());
352  // only plot eta and phi maps for non-zero candidates
353  if (tj->rank()) {
354  l1GctTauJetsEtEtaPhi_->Fill(tj->regionId().ieta(),tj->regionId().iphi(),tj->rank());
355  l1GctTauJetsOccEtaPhi_->Fill(tj->regionId().ieta(),tj->regionId().iphi());
356  }
357  }
358  if (tj->rank()) l1GctAllJetsOccRankBx_->Fill(tj->bx(),tj->rank()); // for all BX
359  }
360  } else {
361  edm::LogWarning("DataNotFound") << " Could not find l1TauJets label was " << gctTauJetsSource_ ;
362  }
363 
364  // IsoTau jets
365  if (m_stage1_layer2_ == true) {
367  e.getByToken(gctIsoTauJetsSourceToken_, l1IsoTauJets);
368  if (l1IsoTauJets.isValid()) {
369  for (L1GctJetCandCollection::const_iterator itj = l1IsoTauJets->begin(); itj != l1IsoTauJets->end(); itj++) {
370  // only plot central BX
371  if (itj->bx()==0) {
372  l1GctIsoTauJetsRank_->Fill(itj->rank());
373  // only plot eta and phi maps for non-zero candidates
374  if (itj->rank()) {
375  l1GctIsoTauJetsEtEtaPhi_->Fill(itj->regionId().ieta(),itj->regionId().iphi(),itj->rank());
376  l1GctIsoTauJetsOccEtaPhi_->Fill(itj->regionId().ieta(),itj->regionId().iphi());
377  }
378  }
379  if (itj->rank()) l1GctAllJetsOccRankBx_->Fill(itj->bx(),itj->rank()); // for all BX
380  }
381  } else {
382  edm::LogWarning("DataNotFound") << " Could not find l1IsoTauJets label was " << gctIsoTauJetsSource_ ;
383  }
384  }
385 
386 
387  // Missing ET
388  if (l1EtMiss.isValid()) {
389  for (L1GctEtMissCollection::const_iterator met = l1EtMiss->begin(); met != l1EtMiss->end(); met++) {
390  // only plot central BX
391  if (met->bx()==0) {
392  if (met->overFlow() == 0 && met->et() > 0) {
393  //Avoid problems with met=0 candidates affecting MET_PHI plots
394  l1GctEtMiss_->Fill(met->et());
395  l1GctEtMissPhi_->Fill(met->phi());
396  }
397  l1GctEtMissOf_->Fill(met->overFlow());
398  }
399  if (met->overFlow() == 0 && met->et() > 0) l1GctEtMissOccBx_->Fill(met->bx(),met->et()); // for all BX
400  }
401  } else {
402  edm::LogWarning("DataNotFound") << " Could not find l1EtMiss label was " << gctEnergySumsSource_ ;
403  }
404 
405  // Missing HT
406  if (l1HtMiss.isValid()) {
407  for (L1GctHtMissCollection::const_iterator mht = l1HtMiss->begin(); mht != l1HtMiss->end(); mht++) {
408  // only plot central BX
409  if (mht->bx()==0) {
410  if (mht->overFlow() == 0 && mht->et() > 0) {
411  //Avoid problems with mht=0 candidates affecting MHT_PHI plots
412  l1GctHtMiss_->Fill(mht->et());
413  l1GctHtMissPhi_->Fill(mht->phi());
414  }
415  l1GctHtMissOf_->Fill(mht->overFlow());
416  }
417  if (mht->overFlow() == 0 && mht->et() > 0) l1GctHtMissOccBx_->Fill(mht->bx(),mht->et()); // for all BX
418  }
419  } else {
420  edm::LogWarning("DataNotFound") << " Could not find l1HtMiss label was " << gctEnergySumsSource_ ;
421  }
422 
423  // Missing ET HT correlations
424  if (l1HtMiss.isValid() && l1EtMiss.isValid()) {
425  if (l1HtMiss->size() == l1EtMiss->size()) {
426  for (unsigned i=0; i<l1HtMiss->size(); i++) {
427  if (l1HtMiss->at(i).overFlow() == 0 && l1EtMiss->at(i).overFlow() == 0 &&
428  l1HtMiss->at(i).bx() == 0 && l1EtMiss->at(i).bx() == 0) {
429  // Avoid problems overflows and only plot central BX
430  l1GctEtMissHtMissCorr_->Fill(l1EtMiss->at(i).et(),l1HtMiss->at(i).et());
431  if (l1HtMiss->at(i).et() && l1EtMiss->at(i).et()){ // Don't plot phi if one or other is zero
432  l1GctEtMissHtMissCorrPhi_->Fill(l1EtMiss->at(i).phi(),l1HtMiss->at(i).phi());
433  }
434  }
435  }
436  }
437  } else {
438  edm::LogWarning("DataNotFound") << " Could not find l1EtMiss or l1HtMiss label was " << gctEnergySumsSource_ ;
439  }
440 
441  // HT
442  if (l1EtHad.isValid()) {
443  for (L1GctEtHadCollection::const_iterator ht = l1EtHad->begin(); ht != l1EtHad->end(); ht++) {
444  // only plot central BX
445  if (ht->bx()==0) {
446  l1GctEtHad_->Fill(ht->et());
447  l1GctEtHadOf_->Fill(ht->overFlow());
448  }
449  l1GctEtHadOccBx_->Fill(ht->bx(),ht->et()); // for all BX
450  }
451  } else {
452  edm::LogWarning("DataNotFound") << " Could not find l1EtHad label was " << gctEnergySumsSource_ ;
453  }
454 
455  // Total ET
456  if (l1EtTotal.isValid()) {
457  for (L1GctEtTotalCollection::const_iterator et = l1EtTotal->begin(); et != l1EtTotal->end(); et++) {
458  // only plot central BX
459  if (et->bx()==0) {
460  l1GctEtTotal_->Fill(et->et());
461  l1GctEtTotalOf_->Fill(et->overFlow());
462  }
463  l1GctEtTotalOccBx_->Fill(et->bx(),et->et()); // for all BX
464  }
465  } else {
466  edm::LogWarning("DataNotFound") << " Could not find l1EtTotal label was " << gctEnergySumsSource_ ;
467  }
468 
469  // Total ET HT correlations
470  if (l1EtTotal.isValid() && l1EtHad.isValid()) {
471  if (l1EtTotal->size() == l1EtHad->size()) {
472  for (unsigned i=0; i<l1EtHad->size(); i++) {
473  if (l1EtHad->at(i).overFlow() == 0 && l1EtTotal->at(i).overFlow() == 0 &&
474  l1EtHad->at(i).bx() == 0 && l1EtTotal->at(i).bx() == 0) {
475  // Avoid problems overflows and only plot central BX
476  l1GctEtTotalEtHadCorr_->Fill(l1EtTotal->at(i).et(),l1EtHad->at(i).et());
477  }
478  }
479  }
480  } else {
481  edm::LogWarning("DataNotFound") << " Could not find l1EtTotal or l1EtHad label was " << gctEnergySumsSource_ ;
482  }
483 
484  //HF Ring Et Sums
485  if (l1HFSums.isValid()) {
486  for (L1GctHFRingEtSumsCollection::const_iterator hfs=l1HFSums->begin(); hfs!=l1HFSums->end(); hfs++){
487  // only plot central BX
488  if (hfs->bx()==0) {
489  // Individual ring Et sums
490  l1GctHFRing1ETSumPosEta_->Fill(hfs->etSum(0));
491  l1GctHFRing1ETSumNegEta_->Fill(hfs->etSum(1));
492  l1GctHFRing2ETSumPosEta_->Fill(hfs->etSum(2));
493  l1GctHFRing2ETSumNegEta_->Fill(hfs->etSum(3));
494  // Ratio of ring Et sums
495  if (hfs->etSum(2)!=0) l1GctHFRingRatioPosEta_->Fill((hfs->etSum(0))/(hfs->etSum(2)));
496  if (hfs->etSum(3)!=0) l1GctHFRingRatioNegEta_->Fill((hfs->etSum(1))/(hfs->etSum(3)));
497  // Correlate positive and neagative eta
498  l1GctHFRing1PosEtaNegEta_->Fill(hfs->etSum(0),hfs->etSum(1));
499  l1GctHFRing2PosEtaNegEta_->Fill(hfs->etSum(2),hfs->etSum(3));
500  }
501  // Occupancy vs BX
502  for (unsigned i=0; i<4; i++){
503  l1GctHFRingETSumOccBx_->Fill(hfs->bx(),hfs->etSum(i));
504  }
505  }
506  } else {
507  edm::LogWarning("DataNotFound") << " Could not find l1HFSums label was " << gctEnergySumsSource_ ;
508  }
509 
510  // HF Ring Counts
511  if (l1HFCounts.isValid()) {
512  for (L1GctHFBitCountsCollection::const_iterator hfc=l1HFCounts->begin(); hfc!=l1HFCounts->end(); hfc++){
513  // only plot central BX
514  if (hfc->bx()==0) {
515  // Individual ring counts
516  l1GctHFRing1TowerCountPosEta_->Fill(hfc->bitCount(0));
517  l1GctHFRing1TowerCountNegEta_->Fill(hfc->bitCount(1));
518  l1GctHFRing2TowerCountPosEta_->Fill(hfc->bitCount(2));
519  l1GctHFRing2TowerCountNegEta_->Fill(hfc->bitCount(3));
520  // Correlate positive and negative eta
521  l1GctHFRing1TowerCountPosEtaNegEta_->Fill(hfc->bitCount(0),hfc->bitCount(1));
522  l1GctHFRing2TowerCountPosEtaNegEta_->Fill(hfc->bitCount(2),hfc->bitCount(3));
523  }
524  // Occupancy vs BX
525  for (unsigned i=0; i<4; i++){
526  l1GctHFRingTowerCountOccBx_->Fill(hfc->bx(),hfc->bitCount(i));
527  }
528  }
529  } else {
530  edm::LogWarning("DataNotFound") << " Could not find l1HFCounts label was " << gctEnergySumsSource_ ;
531  }
532 
533  // Isolated EM
534  if (l1IsoEm.isValid()) {
535  for (L1GctEmCandCollection::const_iterator ie=l1IsoEm->begin(); ie!=l1IsoEm->end(); ie++) {
536  // only plot central BX
537  if (ie->bx()==0) {
538  l1GctIsoEmRank_->Fill(ie->rank());
539  // only plot eta and phi maps for non-zero candidates
540  if (ie->rank()){
541  l1GctIsoEmRankEtaPhi_->Fill(ie->regionId().ieta(),ie->regionId().iphi(),ie->rank());
542  l1GctIsoEmOccEtaPhi_->Fill(ie->regionId().ieta(),ie->regionId().iphi());
543  }
544  }
545  if (ie->rank()) l1GctAllEmOccRankBx_->Fill(ie->bx(),ie->rank()); // for all BX
546  }
547  } else {
548  edm::LogWarning("DataNotFound") << " Could not find l1IsoEm label was " << gctIsoEmSource_ ;
549  }
550 
551  // Non-isolated EM
552  if (l1NonIsoEm.isValid()) {
553  for (L1GctEmCandCollection::const_iterator ne=l1NonIsoEm->begin(); ne!=l1NonIsoEm->end(); ne++) {
554  // only plot central BX
555  if (ne->bx()==0) {
556  l1GctNonIsoEmRank_->Fill(ne->rank());
557  // only plot eta and phi maps for non-zero candidates
558  if (ne->rank()){
559  l1GctNonIsoEmRankEtaPhi_->Fill(ne->regionId().ieta(),ne->regionId().iphi(),ne->rank());
560  l1GctNonIsoEmOccEtaPhi_->Fill(ne->regionId().ieta(),ne->regionId().iphi());
561  }
562  }
563  if (ne->rank()) l1GctAllEmOccRankBx_->Fill(ne->bx(),ne->rank()); // for all BX
564  }
565  } else {
566  edm::LogWarning("DataNotFound") << " Could not find l1NonIsoEm label was " << gctNonIsoEmSource_ ;
567  }
568  edm::LogInfo("L1TGCT") << "L1TGCT: end job...." << std::endl;
569  edm::LogInfo("EndJob") << "analyzed " << nev_ << " events";
570 }
571 
572 
const float METPHIMAX
Definition: L1TGCT.cc:25
T getParameter(std::string const &) const
const unsigned int R6BINS
Definition: L1TCompare.cc:71
const unsigned int BXBINS
Definition: L1TGCT.cc:39
T getUntrackedParameter(std::string const &, T const &) const
#define PHIMAX
int i
Definition: DBlmapReader.cc:9
MonitorElement * l1GctHtMiss_
Definition: L1TGCT.h:182
MonitorElement * l1GctEtHad_
Definition: L1TGCT.h:191
MonitorElement * l1GctHtMissPhi_
Definition: L1TGCT.h:183
const unsigned int OFBINS
Definition: L1TGCT.cc:35
edm::InputTag gctEnergySumsSource_
Definition: L1TGCT.h:226
MonitorElement * l1GctHFRing2TowerCountNegEta_
Definition: L1TGCT.h:204
MonitorElement * l1GctEtMissOccBx_
Definition: L1TGCT.h:181
MonitorElement * l1GctCenJetsRank_
Definition: L1TGCT.h:167
MonitorElement * l1GctHFRing1ETSumNegEta_
Definition: L1TGCT.h:206
MonitorElement * l1GctEtTotalEtHadCorr_
Definition: L1TGCT.h:194
MonitorElement * l1GctHFRing1TowerCountPosEta_
Definition: L1TGCT.h:201
const unsigned int EMETABINS
Definition: L1TGCT.cc:19
MonitorElement * l1GctTauJetsEtEtaPhi_
Definition: L1TGCT.h:154
const float BXMIN
Definition: L1TGCT.cc:40
MonitorElement * l1GctHFRingETSumOccBx_
Definition: L1TGCT.h:211
MonitorElement * l1GctNonIsoEmRank_
Definition: L1TGCT.h:172
bool getByToken(EDGetToken token, Handle< PROD > &result) const
Definition: Event.h:464
const unsigned int MHTPHIBINS
Definition: L1TGCT.cc:27
const float OFMAX
Definition: L1TGCT.cc:37
MonitorElement * l1GctEtMiss_
Definition: L1TGCT.h:178
MonitorElement * bookInt(Args &&...args)
Definition: DQMStore.h:103
std::string monitorDir_
Definition: L1TGCT.h:145
MonitorElement * l1GctHFRing2ETSumNegEta_
Definition: L1TGCT.h:208
MonitorElement * l1GctAllEmOccRankBx_
Definition: L1TGCT.h:175
edm::EDGetTokenT< L1GctJetCandCollection > gctCenJetsSourceToken_
Definition: L1TGCT.h:237
const float BXMAX
Definition: L1TGCT.cc:41
MonitorElement * l1GctEtTotal_
Definition: L1TGCT.h:188
MonitorElement * l1GctHtMissOf_
Definition: L1TGCT.h:184
edm::EDGetTokenT< L1GctHtMissCollection > l1HtMissToken_
Definition: L1TGCT.h:244
MonitorElement * l1GctIsoTauJetsEtEtaPhi_
Definition: L1TGCT.h:155
MonitorElement * l1GctEtHadOccBx_
Definition: L1TGCT.h:193
MonitorElement * l1GctEtTotalOccBx_
Definition: L1TGCT.h:190
edm::InputTag gctIsoTauJetsSource_
Definition: L1TGCT.h:225
bool isRealData() const
Definition: EventBase.h:64
edm::EDGetTokenT< L1GctEtTotalCollection > l1EtTotalToken_
Definition: L1TGCT.h:246
MonitorElement * l1GctHFRingTowerCountOccBx_
Definition: L1TGCT.h:212
MonitorElement * l1GctForJetsRank_
Definition: L1TGCT.h:168
MonitorElement * l1GctTauJetsRank_
Definition: L1TGCT.h:169
MonitorElement * l1GctHFRing2ETSumPosEta_
Definition: L1TGCT.h:207
virtual void beginLuminosityBlock(const edm::LuminosityBlock &, const edm::EventSetup &)
Definition: L1TGCT.cc:228
void Fill(long long x)
edm::InputTag gctCenJetsSource_
Definition: L1TGCT.h:222
L1TGCT(const edm::ParameterSet &ps)
Definition: L1TGCT.cc:60
MonitorElement * l1GctEtMissPhi_
Definition: L1TGCT.h:179
MonitorElement * l1GctCenJetsEtEtaPhi_
Definition: L1TGCT.h:152
virtual void dqmBeginRun(const edm::Run &, const edm::EventSetup &)
Definition: L1TGCT.cc:224
edm::InputTag gctTauJetsSource_
Definition: L1TGCT.h:224
edm::EDGetTokenT< L1GctEmCandCollection > gctIsoEmSourceToken_
Definition: L1TGCT.h:235
MonitorElement * l1GctAllJetsEtEtaPhi_
Definition: L1TGCT.h:151
const float JETETAMAX
Definition: L1TGCT.cc:17
MonitorElement * l1GctHFRing1ETSumPosEta_
Definition: L1TGCT.h:205
const float R6MAX
Definition: L1TCompare.cc:73
edm::EDGetTokenT< L1GctJetCandCollection > gctTauJetsSourceToken_
Definition: L1TGCT.h:239
const unsigned int R3BINS
Definition: L1TGCT.cc:44
edm::EDGetTokenT< L1GctHFBitCountsCollection > l1HFCountsToken_
Definition: L1TGCT.h:242
MonitorElement * l1GctHFRingRatioNegEta_
Definition: L1TGCT.h:210
const float EMETAMIN
Definition: L1TGCT.cc:20
MonitorElement * l1GctCenJetsOccEtaPhi_
Definition: L1TGCT.h:160
MonitorElement * l1GctIsoTauJetsRank_
Definition: L1TGCT.h:170
MonitorElement * l1GctEtMissOf_
Definition: L1TGCT.h:180
MonitorElement * l1GctForJetsEtEtaPhi_
Definition: L1TGCT.h:153
const unsigned int JETETABINS
Definition: L1TGCT.cc:15
MonitorElement * l1GctIsoEmRank_
Definition: L1TGCT.h:171
const float R12MAX
Definition: L1TGCT.cc:58
MonitorElement * book1D(Args &&...args)
Definition: DQMStore.h:115
MonitorElement * l1GctHtMissOccBx_
Definition: L1TGCT.h:185
void analyze(const edm::Event &e, const edm::EventSetup &c)
Definition: L1TGCT.cc:232
const float R7MAX
Definition: L1TGCT.cc:55
edm::EDGetTokenT< L1GctEmCandCollection > gctNonIsoEmSourceToken_
Definition: L1TGCT.h:236
MonitorElement * runId_
Definition: L1TGCT.h:213
const float R5MAX
Definition: L1TGCT.cc:49
MonitorElement * l1GctHFRing2TowerCountPosEta_
Definition: L1TGCT.h:203
edm::EDGetTokenT< L1GctEtHadCollection > l1EtHadToken_
Definition: L1TGCT.h:245
const float R12MIN
Definition: L1TGCT.cc:57
MonitorElement * l1GctHFRing2TowerCountPosEtaNegEta_
Definition: L1TGCT.h:200
bool isValid() const
Definition: HandleBase.h:75
edm::InputTag gctIsoEmSource_
Definition: L1TGCT.h:227
#define PHIMIN
MonitorElement * triggerType_
Definition: L1TGCT.h:148
const float OFMIN
Definition: L1TGCT.cc:36
MonitorElement * l1GctEtMissHtMissCorrPhi_
Definition: L1TGCT.h:187
const float MHTPHIMAX
Definition: L1TGCT.cc:29
const float R5MIN
Definition: L1TGCT.cc:48
MonitorElement * l1GctNonIsoEmOccEtaPhi_
Definition: L1TGCT.h:165
const float MHTPHIMIN
Definition: L1TGCT.cc:28
int filterTriggerType_
Definition: L1TGCT.h:232
const unsigned int R7BINS
Definition: L1TGCT.cc:53
MonitorElement * lumisecId_
Definition: L1TGCT.h:214
const float R6MIN
Definition: L1TCompare.cc:72
void setCurrentFolder(const std::string &fullpath)
Definition: DQMStore.cc:274
edm::InputTag gctForJetsSource_
Definition: L1TGCT.h:223
MonitorElement * l1GctNonIsoEmRankEtaPhi_
Definition: L1TGCT.h:157
int nev_
Definition: L1TGCT.h:216
MonitorElement * l1GctIsoTauJetsOccEtaPhi_
Definition: L1TGCT.h:163
MonitorElement * book2D(Args &&...args)
Definition: DQMStore.h:133
bool m_stage1_layer2_
filter TriggerType
Definition: L1TGCT.h:231
MonitorElement * l1GctHFRing1PosEtaNegEta_
Definition: L1TGCT.h:197
const unsigned int R12BINS
Definition: L1TGCT.cc:56
const unsigned int R5BINS
Definition: L1TGCT.cc:47
virtual ~L1TGCT()
Definition: L1TGCT.cc:107
edm::EDGetTokenT< L1GctHFRingEtSumsCollection > gctEnergySumsSourceToken_
Definition: L1TGCT.h:241
MonitorElement * l1GctAllJetsOccEtaPhi_
Definition: L1TGCT.h:159
edm::EDGetTokenT< L1GctJetCandCollection > gctIsoTauJetsSourceToken_
Definition: L1TGCT.h:240
MonitorElement * l1GctHFRingRatioPosEta_
Definition: L1TGCT.h:209
virtual void bookHistograms(DQMStore::IBooker &ibooker, edm::Run const &, edm::EventSetup const &) override
Definition: L1TGCT.cc:111
MonitorElement * l1GctIsoEmRankEtaPhi_
Definition: L1TGCT.h:156
MonitorElement * l1GctHFRing1TowerCountPosEtaNegEta_
Definition: L1TGCT.h:199
const float METPHIMIN
Definition: L1TGCT.cc:24
edm::EventAuxiliary::ExperimentType experimentType() const
Definition: EventBase.h:65
edm::EDGetTokenT< L1GctJetCandCollection > gctForJetsSourceToken_
Definition: L1TGCT.h:238
#define PHIBINS
MonitorElement * l1GctEtHadOf_
Definition: L1TGCT.h:192
std::string outputFile_
Definition: L1TGCT.h:217
MonitorElement * l1GctForJetsOccEtaPhi_
Definition: L1TGCT.h:161
const float EMETAMAX
Definition: L1TGCT.cc:21
MonitorElement * l1GctHFRing1TowerCountNegEta_
Definition: L1TGCT.h:202
edm::EDGetTokenT< L1GctEtMissCollection > l1EtMissToken_
Definition: L1TGCT.h:243
bool verbose_
Definition: L1TGCT.h:218
MonitorElement * l1GctEtTotalOf_
Definition: L1TGCT.h:189
const float R3MAX
Definition: L1TGCT.cc:46
MonitorElement * l1GctAllJetsOccRankBx_
Definition: L1TGCT.h:174
MonitorElement * l1GctEtMissHtMissCorr_
Definition: L1TGCT.h:186
MonitorElement * l1GctHFRing2PosEtaNegEta_
Definition: L1TGCT.h:198
const unsigned int METPHIBINS
Definition: L1TGCT.cc:23
edm::InputTag gctNonIsoEmSource_
Definition: L1TGCT.h:228
MonitorElement * l1GctTauJetsOccEtaPhi_
Definition: L1TGCT.h:162
Definition: Run.h:43
MonitorElement * l1GctIsoEmOccEtaPhi_
Definition: L1TGCT.h:164
const float R3MIN
Definition: L1TGCT.cc:45
const float R7MIN
Definition: L1TGCT.cc:54
const float JETETAMIN
Definition: L1TGCT.cc:16