CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
L1TCSCTF.cc
Go to the documentation of this file.
1 /*
2  * \file L1TCSCTF.cc
3  *
4  * \author J. Berryhill
5  *
6  */
7 
8 #include "DQM/L1TMonitor/interface/L1TCSCTF.h"
9 #include "DQMServices/Core/interface/DQMStore.h"
10 
11 // includes to fetch all reguired data products from the edm::Event
12 #include "DataFormats/CSCDigi/interface/CSCCorrelatedLCTDigiCollection.h"
13 #include "DataFormats/L1CSCTrackFinder/interface/L1CSCTrackCollection.h"
14 #include "DataFormats/L1CSCTrackFinder/interface/L1CSCStatusDigiCollection.h"
15 #include "L1Trigger/CSCCommonTrigger/interface/CSCTriggerGeometry.h"
16 #include "Geometry/Records/interface/MuonGeometryRecord.h"
17 
18 #include "DataFormats/L1CSCTrackFinder/interface/CSCTriggerContainer.h"
19 #include "DataFormats/L1CSCTrackFinder/interface/TrackStub.h"
20 
21 
22 using namespace std;
23 using namespace edm;
24 
25 L1TCSCTF::L1TCSCTF(const ParameterSet& ps)
26 // if some piece of data is absent - configure corresponding source with 'null:'
27 // : csctfSource_( ps.getParameter< InputTag >("csctfSource") )
28  : gmtProducer( ps.getParameter< InputTag >("gmtProducer") ),
29  lctProducer( ps.getParameter< InputTag >("lctProducer") ),
30  trackProducer( ps.getParameter< InputTag >("trackProducer") ),
31  statusProducer( ps.getParameter< InputTag >("statusProducer") ),
32  mbProducer( ps.getParameter< InputTag >("mbProducer") )
33 {
34 
35  // verbosity switch
36  verbose_ = ps.getUntrackedParameter<bool>("verbose", false);
37 
38  if(verbose_) edm::LogInfo("DataNotFound") << "L1TCSCTF: constructor...." << endl;
39 
40 
41  dbe = NULL;
42  if ( ps.getUntrackedParameter<bool>("DQMStore", false) )
43  {
44  dbe = Service<DQMStore>().operator->();
45  dbe->setVerbose(0);
46  }
47 
48  outputFile_ = ps.getUntrackedParameter<string>("outputFile", "");
49  if ( outputFile_.size() != 0 )
50  {
51  edm::LogInfo("DataNotFound") << "L1T Monitoring histograms will be saved to " << outputFile_.c_str() << endl;
52  }
53 
54  bool disable = ps.getUntrackedParameter<bool>("disableROOToutput", false);
55  if(disable){
56  outputFile_="";
57  }
58 
59 
60  if ( dbe !=NULL )
61  {
62  dbe->setCurrentFolder("L1T/L1TCSCTF");
63  }
64 
65  // instantiate standard on-fly SR LUTs from CSC TF emulator package
66  bzero(srLUTs_,sizeof(srLUTs_));
67  int endcap=1, sector=1; // assume SR LUTs are all same for every sector in either of endcaps
68  bool TMB07=true; // specific TMB firmware
69  // Create a dummy pset for SR LUTs
70  edm::ParameterSet srLUTset;
71  srLUTset.addUntrackedParameter<bool>("ReadLUTs", false);
72  srLUTset.addUntrackedParameter<bool>("Binary", false);
73  srLUTset.addUntrackedParameter<std::string>("LUTPath", "./");
74  for(int station=1,fpga=0; station<=4 && fpga<5; station++)
75  {
76  if(station==1)
77  for(int subSector=0; subSector<2 && fpga<5; subSector++)
78  srLUTs_[fpga++] = new CSCSectorReceiverLUT(endcap, sector, subSector+1, station, srLUTset, TMB07);
79  else
80  srLUTs_[fpga++] = new CSCSectorReceiverLUT(endcap, sector, 0, station, srLUTset, TMB07);
81  }
82 
83  //set Token(-s)
84  edm::InputTag statusTag_(statusProducer.label(),statusProducer.instance());
85  edm::InputTag corrlctsTag_(lctProducer.label(),lctProducer.instance());
86  edm::InputTag tracksTag_(trackProducer.label(),trackProducer.instance());
87  edm::InputTag dtStubsTag_(mbProducer.label(), mbProducer.instance());
88  edm::InputTag mbtracksTag_(trackProducer.label(),trackProducer.instance());
89 
90  gmtProducerToken_ = consumes<L1MuGMTReadoutCollection>(ps.getParameter< InputTag >("gmtProducer"));
91  statusToken_ = consumes<L1CSCStatusDigiCollection>(statusTag_);
92  corrlctsToken_ = consumes<CSCCorrelatedLCTDigiCollection>(corrlctsTag_);
93  tracksToken_ = consumes<L1CSCTrackCollection>(tracksTag_);
94  dtStubsToken_ = consumes<CSCTriggerContainer<csctf::TrackStub> >(dtStubsTag_);
95  mbtracksToken_ = consumes<L1CSCTrackCollection>(mbtracksTag_);
96 }
97 
98 L1TCSCTF::~L1TCSCTF()
99 {
100 
101  for(int i=0; i<5; i++)
102  delete srLUTs_[i]; //free the array of pointers
103 }
104 
105 void L1TCSCTF::beginJob(void)
106 {
107  m_scalesCacheID = -999;
108  m_ptScaleCacheID = -999;
109 
110  nev_ = 0;
111 
112  // get hold of back-end interface
113  DQMStore* dbe = 0;
114  dbe = Service<DQMStore>().operator->();
115 
116  if( dbe )
117  {
118  dbe->setCurrentFolder("L1T/L1TCSCTF");
119  dbe->rmdir("L1T/L1TCSCTF");
120  }
121 
122 
123  if( dbe )
124  {
125  dbe->setCurrentFolder("L1T/L1TCSCTF");
126 
127  // Error counting histogram:
128  // 1) checks TF data integrity (error rate - first bin),
129  // 2) monitors sychronization on input links (4 errors types: SE/SM/BX/AF; ORed for all time bins, links, and SPs),
130  // 3) reports FMM status (if in any SP FMM status != "Ready" - fill the last bin)
131  csctferrors = dbe->book1D("CSCTF_errors","CSCTF Errors",6,0,6);
132  csctferrors->setAxisTitle("Error type",1);
133  csctferrors->setAxisTitle("Number of Errors",2);
134  csctferrors->setBinLabel(1,"Corruptions",1);
135  csctferrors->setBinLabel(2,"Synch. Err.",1);
136  csctferrors->setBinLabel(3,"Synch. Mod.",1);
137  csctferrors->setBinLabel(4,"BX mismatch",1);
138  csctferrors->setBinLabel(5,"Time misalign.",1);
139  csctferrors->setBinLabel(6,"FMM != Ready",1);
140 
141  // Occupancy histogram Eta x Y, where Y:
142  // 1) Phi_packed of input LCTs from 1st, 2nd, 3rd, and 4th stations
143  // 2) Phi_packed of output tracks
144  // (all 12 SPs - 360 degree coveradge)
145  csctfoccupancies = dbe->book2D("CSCTF_occupancies", "CSCTF Occupancies", 64,-32,31,32,0,6.2);
146  csctfoccupancies->setAxisTitle("#eta",1);
147  csctfoccupancies->setAxisTitle("#phi",2);
148  csctfoccupancies->setBinLabel( 1,"-2.5", 1);
149  csctfoccupancies->setBinLabel( 8,"-2.1", 1);
150  csctfoccupancies->setBinLabel(18,"-1.6", 1);
151  csctfoccupancies->setBinLabel(26,"-1.2", 1);
152  csctfoccupancies->setBinLabel(32,"-0.9", 1);
153  csctfoccupancies->setBinLabel(33, "0.9", 1);
154  csctfoccupancies->setBinLabel(39, "1.2", 1);
155  csctfoccupancies->setBinLabel(47, "1.6", 1);
156  csctfoccupancies->setBinLabel(57, "2.1", 1);
157  csctfoccupancies->setBinLabel(64, "2.5", 1);
158 
159  // ... and for halo muons only
160  csctfoccupancies_H = dbe->book2D("CSCTF_occupancies_H", "CSCTF Halo Occupancies", 64,-32,31,32,0,6.2);
161  csctfoccupancies_H->setAxisTitle("#eta",1);
162  csctfoccupancies_H->setAxisTitle("#phi",2);
163  csctfoccupancies_H->setBinLabel( 1,"-2.5", 1);
164  csctfoccupancies_H->setBinLabel( 8,"-2.1", 1);
165  csctfoccupancies_H->setBinLabel(18,"-1.6", 1);
166  csctfoccupancies_H->setBinLabel(26,"-1.2", 1);
167  csctfoccupancies_H->setBinLabel(32,"-0.9", 1);
168  csctfoccupancies_H->setBinLabel(33, "0.9", 1);
169  csctfoccupancies_H->setBinLabel(39, "1.2", 1);
170  csctfoccupancies_H->setBinLabel(47, "1.6", 1);
171  csctfoccupancies_H->setBinLabel(57, "2.1", 1);
172  csctfoccupancies_H->setBinLabel(64, "2.5", 1);
173 
174  //haloDelEta12 = dbe->book1D("CSCTF_Halo_Eta12", "#Delta #eta_{12} for Halo Muons", 40, -0.20,0.30);
175  //haloDelEta112 = dbe->book1D("CSCTF_Halo_Eta112","#Delta #eta_{112} for Halo Muons", 40, -0.20,0.30);
176  //haloDelEta13 = dbe->book1D("CSCTF_Halo_Eta13", "#Delta #eta_{13} for Halo Muons", 40, -0.20,0.30);
177  //haloDelEta113 = dbe->book1D("CSCTF_Halo_Eta113","#Delta #eta_{113} for Halo Muons", 40, -0.20,0.30);
178 
179  // Quality VS Mode
180  trackModeVsQ = dbe->book2D("CSCTF_Track_ModeVsQual","CSC Track Mode Vs Quality", 19, -0.5, 18.5, 4, 0, 4);
181  trackModeVsQ->setAxisTitle("Track Type", 1);
182  trackModeVsQ->setBinLabel(1,"No Track",1);
183  trackModeVsQ->setBinLabel(2,"Bad Phi/Single",1);
184  trackModeVsQ->setBinLabel(3,"ME1-2-3",1);
185  trackModeVsQ->setBinLabel(4,"ME1-2-4",1);
186  trackModeVsQ->setBinLabel(5,"ME1-3-4",1);
187  trackModeVsQ->setBinLabel(6,"ME2-3-4",1);
188  trackModeVsQ->setBinLabel(7,"ME1-2",1);
189  trackModeVsQ->setBinLabel(8,"ME1-3",1);
190  trackModeVsQ->setBinLabel(9,"ME2-3",1);
191  trackModeVsQ->setBinLabel(10,"ME2-4",1);
192  trackModeVsQ->setBinLabel(11,"ME3-4",1);
193  trackModeVsQ->setBinLabel(12,"MB1-ME3",1);
194  trackModeVsQ->setBinLabel(13,"MB1-ME2",1);
195  trackModeVsQ->setBinLabel(14,"ME1-4",1);
196  trackModeVsQ->setBinLabel(15,"MB1-ME1",1);
197  trackModeVsQ->setBinLabel(16,"Halo Trigger",1);
198  trackModeVsQ->setBinLabel(17,"MB1-ME1-2",1);
199  trackModeVsQ->setBinLabel(18,"MB1-ME1-3",1);
200  trackModeVsQ->setBinLabel(19,"MB1-ME2-3",1);
201 
202  trackModeVsQ->setAxisTitle("Quality",2);
203  trackModeVsQ->setBinLabel(1,"0",2);
204  trackModeVsQ->setBinLabel(2,"1",2);
205  trackModeVsQ->setBinLabel(3,"2",2);
206  trackModeVsQ->setBinLabel(4,"3",2);
207 
208  // Mode
209  csctfTrackM = dbe->book1D("CSCTF_Track_Mode","CSC Track Mode", 19, -0.5, 18.5);
210  csctfTrackM->setAxisTitle("Track Type", 1);
211  csctfTrackM->setBinLabel(1,"No Track",1);
212  csctfTrackM->setBinLabel(2,"Bad Phi/Single",1);
213  csctfTrackM->setBinLabel(3,"ME1-2-3",1);
214  csctfTrackM->setBinLabel(4,"ME1-2-4",1);
215  csctfTrackM->setBinLabel(5,"ME1-3-4",1);
216  csctfTrackM->setBinLabel(6,"ME2-3-4",1);
217  csctfTrackM->setBinLabel(7,"ME1-2",1);
218  csctfTrackM->setBinLabel(8,"ME1-3",1);
219  csctfTrackM->setBinLabel(9,"ME2-3",1);
220  csctfTrackM->setBinLabel(10,"ME2-4",1);
221  csctfTrackM->setBinLabel(11,"ME3-4",1);
222  csctfTrackM->setBinLabel(12,"MB1-ME3",1);
223  csctfTrackM->setBinLabel(13,"MB1-ME2",1);
224  csctfTrackM->setBinLabel(14,"ME1-4",1);
225  csctfTrackM->setBinLabel(15,"MB1-ME1",1);
226  csctfTrackM->setBinLabel(16,"Halo Trigger",1);
227  csctfTrackM->setBinLabel(17,"MB1-ME1-2",1);
228  csctfTrackM->setBinLabel(18,"MB1-ME1-3",1);
229  csctfTrackM->setBinLabel(19,"MB1-ME2-3",1);
230 
231  // Chamber Occupancy
232  csctfChamberOccupancies = dbe->book2D("CSCTF_Chamber_Occupancies","CSCTF Chamber Occupancies", 54, -0.05, 5.35, 10, -5.5, 4.5);
233  csctfChamberOccupancies->setAxisTitle("Sector, (chambers 1-9 not labeled)",1);
234  csctfChamberOccupancies->setBinLabel(1,"ME-4",2);
235  csctfChamberOccupancies->setBinLabel(2,"ME-3",2);
236  csctfChamberOccupancies->setBinLabel(3,"ME-2",2);
237  csctfChamberOccupancies->setBinLabel(4,"ME-1b",2);
238  csctfChamberOccupancies->setBinLabel(5,"ME-1a",2);
239  csctfChamberOccupancies->setBinLabel(6,"ME+1a",2);
240  csctfChamberOccupancies->setBinLabel(7,"ME+1b",2);
241  csctfChamberOccupancies->setBinLabel(8,"ME+2",2);
242  csctfChamberOccupancies->setBinLabel(9,"ME+3",2);
243  csctfChamberOccupancies->setBinLabel(10,"ME+4",2);
244  csctfChamberOccupancies->setBinLabel(1, "1",1);
245  csctfChamberOccupancies->setBinLabel(10,"2",1);
246  csctfChamberOccupancies->setBinLabel(19,"3",1);
247  csctfChamberOccupancies->setBinLabel(28,"4",1);
248  csctfChamberOccupancies->setBinLabel(37,"5",1);
249  csctfChamberOccupancies->setBinLabel(46,"6",1);
250 
251  // Track Phi
252  csctfTrackPhi = dbe->book1D("CSCTF_Track_Phi", "CSCTF Track #phi",144,0,2*M_PI);
253  csctfTrackPhi->setAxisTitle("Track #phi", 1);
254 
255  // Track Eta
256  csctfTrackEta = dbe->book1D("CSCTF_Track_Eta", "CSCTF Track #eta",64,-32,32);
257  csctfTrackEta->setAxisTitle("Track #eta", 1);
258  csctfTrackEta->setBinLabel( 1,"-2.5", 1);
259  csctfTrackEta->setBinLabel( 8,"-2.1", 1);
260  csctfTrackEta->setBinLabel(18,"-1.6", 1);
261  csctfTrackEta->setBinLabel(26,"-1.2", 1);
262  csctfTrackEta->setBinLabel(32,"-0.9", 1);
263  csctfTrackEta->setBinLabel(33, "0.9", 1);
264  csctfTrackEta->setBinLabel(39, "1.2", 1);
265  csctfTrackEta->setBinLabel(47, "1.6", 1);
266  csctfTrackEta->setBinLabel(57, "2.1", 1);
267  csctfTrackEta->setBinLabel(64, "2.5", 1);
268 
269  // Track Eta Low Quality
270  csctfTrackEtaLowQ = dbe->book1D("CSCTF_Track_Eta_LowQ", "CSCTF Track #eta LQ",64,-32,32);
271  csctfTrackEtaLowQ->setAxisTitle("Track #eta", 1);
272  csctfTrackEtaLowQ->setBinLabel( 1,"-2.5", 1);
273  csctfTrackEtaLowQ->setBinLabel( 8,"-2.1", 1);
274  csctfTrackEtaLowQ->setBinLabel(18,"-1.6", 1);
275  csctfTrackEtaLowQ->setBinLabel(26,"-1.2", 1);
276  csctfTrackEtaLowQ->setBinLabel(32,"-0.9", 1);
277  csctfTrackEtaLowQ->setBinLabel(33, "0.9", 1);
278  csctfTrackEtaLowQ->setBinLabel(39, "1.2", 1);
279  csctfTrackEtaLowQ->setBinLabel(47, "1.6", 1);
280  csctfTrackEtaLowQ->setBinLabel(57, "2.1", 1);
281  csctfTrackEtaLowQ->setBinLabel(64, "2.5", 1);
282 
283 
284  // Track Eta High Quality
285  csctfTrackEtaHighQ = dbe->book1D("CSCTF_Track_Eta_HighQ", "CSCTF Track #eta HQ",64,-32,32);
286  csctfTrackEtaHighQ->setAxisTitle("Track #eta", 1);
287  csctfTrackEtaHighQ->setBinLabel( 1,"-2.5", 1);
288  csctfTrackEtaHighQ->setBinLabel( 8,"-2.1", 1);
289  csctfTrackEtaHighQ->setBinLabel(18,"-1.6", 1);
290  csctfTrackEtaHighQ->setBinLabel(26,"-1.2", 1);
291  csctfTrackEtaHighQ->setBinLabel(32,"-0.9", 1);
292  csctfTrackEtaHighQ->setBinLabel(33, "0.9", 1);
293  csctfTrackEtaHighQ->setBinLabel(39, "1.2", 1);
294  csctfTrackEtaHighQ->setBinLabel(47, "1.6", 1);
295  csctfTrackEtaHighQ->setBinLabel(57, "2.1", 1);
296  csctfTrackEtaHighQ->setBinLabel(64, "2.5", 1);
297 
298 
299  // Halo Phi
300  csctfTrackPhi_H = dbe->book1D("CSCTF_Track_Phi_H", "CSCTF Halo #phi",144,0,2*M_PI);
301  csctfTrackPhi_H->setAxisTitle("Track #phi", 1);
302 
303  // Halo Eta
304  csctfTrackEta_H = dbe->book1D("CSCTF_Track_Eta_H", "CSCTF Halo #eta",64,-32,32);
305  csctfTrackEta_H->setAxisTitle("Track #eta", 1);
306  csctfTrackEta_H->setBinLabel( 1,"-2.5", 1);
307  csctfTrackEta_H->setBinLabel( 8,"-2.1", 1);
308  csctfTrackEta_H->setBinLabel(18,"-1.6", 1);
309  csctfTrackEta_H->setBinLabel(26,"-1.2", 1);
310  csctfTrackEta_H->setBinLabel(32,"-0.9", 1);
311  csctfTrackEta_H->setBinLabel(33, "0.9", 1);
312  csctfTrackEta_H->setBinLabel(39, "1.2", 1);
313  csctfTrackEta_H->setBinLabel(47, "1.6", 1);
314  csctfTrackEta_H->setBinLabel(57, "2.1", 1);
315  csctfTrackEta_H->setBinLabel(64, "2.5", 1);
316 
317  // Track Timing
318  csctfbx = dbe->book2D("CSCTF_bx","CSCTF BX", 12,1,13, 7,-3,3) ;
319  csctfbx->setAxisTitle("Sector (Endcap)", 1);
320  csctfbx->setBinLabel( 1," 1 (+)",1);
321  csctfbx->setBinLabel( 2," 2 (+)",1);
322  csctfbx->setBinLabel( 3," 3 (+)",1);
323  csctfbx->setBinLabel( 4," 4 (+)",1);
324  csctfbx->setBinLabel( 5," 5 (+)",1);
325  csctfbx->setBinLabel( 6," 6 (+)",1);
326  csctfbx->setBinLabel( 7," 7 (-)",1);
327  csctfbx->setBinLabel( 8," 8 (-)",1);
328  csctfbx->setBinLabel( 9," 9 (-)",1);
329  csctfbx->setBinLabel(10,"10 (-)",1);
330  csctfbx->setBinLabel(11,"11 (-)",1);
331  csctfbx->setBinLabel(12,"12 (-)",1);
332 
333  csctfbx->setAxisTitle("CSCTF BX", 2);
334  csctfbx->setBinLabel( 1, "-3", 2);
335  csctfbx->setBinLabel( 2, "-2", 2);
336  csctfbx->setBinLabel( 3, "-1", 2);
337  csctfbx->setBinLabel( 4, "-0", 2);
338  csctfbx->setBinLabel( 5, " 1", 2);
339  csctfbx->setBinLabel( 6, " 2", 2);
340  csctfbx->setBinLabel( 7, " 3", 2);
341 
342  // Halo Timing
343  csctfbx_H = dbe->book2D("CSCTF_bx_H","CSCTF HALO BX", 12,1,13, 7,-3,3) ;
344  csctfbx_H->setAxisTitle("Sector (Endcap)", 1);
345  csctfbx_H->setBinLabel( 1," 1 (+)",1);
346  csctfbx_H->setBinLabel( 2," 2 (+)",1);
347  csctfbx_H->setBinLabel( 3," 3 (+)",1);
348  csctfbx_H->setBinLabel( 4," 4 (+)",1);
349  csctfbx_H->setBinLabel( 5," 5 (+)",1);
350  csctfbx_H->setBinLabel( 6," 6 (+)",1);
351  csctfbx_H->setBinLabel( 7," 7 (-)",1);
352  csctfbx_H->setBinLabel( 8," 8 (-)",1);
353  csctfbx_H->setBinLabel( 9," 9 (-)",1);
354  csctfbx_H->setBinLabel(10,"10 (-)",1);
355  csctfbx_H->setBinLabel(11,"11 (-)",1);
356  csctfbx_H->setBinLabel(12,"12 (-)",1);
357 
358  csctfbx_H->setAxisTitle("CSCTF BX", 2);
359  csctfbx_H->setBinLabel( 1, "-3", 2);
360  csctfbx_H->setBinLabel( 2, "-2", 2);
361  csctfbx_H->setBinLabel( 3, "-1", 2);
362  csctfbx_H->setBinLabel( 4, "-0", 2);
363  csctfbx_H->setBinLabel( 5, " 1", 2);
364  csctfbx_H->setBinLabel( 6, " 2", 2);
365  csctfbx_H->setBinLabel( 7, " 3", 2);
366 
367  // Number of Tracks Stubs
368  cscTrackStubNumbers = dbe->book1D("CSCTF_TrackStubs", "Number of Stubs in CSCTF Tracks", 5, 0, 5);
369  cscTrackStubNumbers->setBinLabel( 1, "0", 1);
370  cscTrackStubNumbers->setBinLabel( 2, "1", 1);
371  cscTrackStubNumbers->setBinLabel( 3, "2", 1);
372  cscTrackStubNumbers->setBinLabel( 4, "3", 1);
373  cscTrackStubNumbers->setBinLabel( 5, "4", 1);
374 
375  // Number of Tracks
376  csctfntrack = dbe->book1D("CSCTF_ntrack","Number of CSCTracks found per event", 5, 0, 5 ) ;
377  csctfntrack->setBinLabel( 1, "0", 1);
378  csctfntrack->setBinLabel( 2, "1", 1);
379  csctfntrack->setBinLabel( 3, "2", 1);
380  csctfntrack->setBinLabel( 4, "3", 1);
381  csctfntrack->setBinLabel( 5, "4", 1);
382  }
383 
384  char hname [200];
385  char htitle[200];
386 
387  for(int i=0; i<12; i++) {
388 
389  sprintf(hname ,"DTstubsTimeTrackMenTimeArrival_%d",i+1);
390  sprintf(htitle,"T_{track} - T_{DT stub} sector %d",i+1);
391 
392  DTstubsTimeTrackMenTimeArrival[i] = dbe->book2D(hname,htitle, 7,-3,3, 2,1,3);
393  DTstubsTimeTrackMenTimeArrival[i]->getTH2F()->SetMinimum(0);
394 
395  // axis makeup
396  DTstubsTimeTrackMenTimeArrival[i]->setAxisTitle("bx_{CSC track} - bx_{DT stub}",1);
397  DTstubsTimeTrackMenTimeArrival[i]->setAxisTitle("subsector",2);
398 
399  DTstubsTimeTrackMenTimeArrival[i]->setBinLabel(1,"-3",1);
400  DTstubsTimeTrackMenTimeArrival[i]->setBinLabel(2,"-2",1);
401  DTstubsTimeTrackMenTimeArrival[i]->setBinLabel(3,"-1",1);
402  DTstubsTimeTrackMenTimeArrival[i]->setBinLabel(4, "0",1);
403  DTstubsTimeTrackMenTimeArrival[i]->setBinLabel(5,"+1",1);
404  DTstubsTimeTrackMenTimeArrival[i]->setBinLabel(6,"+2",1);
405  DTstubsTimeTrackMenTimeArrival[i]->setBinLabel(7,"+3",1);
406 
407  DTstubsTimeTrackMenTimeArrival[i]->setBinLabel(1,"sub1",2);
408  DTstubsTimeTrackMenTimeArrival[i]->setBinLabel(2,"sub2",2);
409 
410  }
411 
412 
413  // NEW: CSC EVENT LCT PLOTS, Renjie Wang
414  csctflcts = dbe->book2D("CSCTF_LCT", "CSCTF LCTs", 12,1,13, 18,0,18);
415  csctflcts->setAxisTitle("CSCTF LCT BX",1);
416  csctflcts->setBinLabel(1,"1",1);
417  csctflcts->setBinLabel(2,"2",1);
418  csctflcts->setBinLabel(3,"3",1);
419  csctflcts->setBinLabel(4,"4",1);
420  csctflcts->setBinLabel(5,"5",1);
421  csctflcts->setBinLabel(6,"6",1);
422  csctflcts->setBinLabel(7,"7",1);
423  csctflcts->setBinLabel(8,"8",1);
424  csctflcts->setBinLabel(9,"9",1);
425  csctflcts->setBinLabel(10,"10",1);
426  csctflcts->setBinLabel(11,"11",1);
427  csctflcts->setBinLabel(12,"12",1);
428 
429  int ihist = 0;
430  for (int iEndcap = 0; iEndcap < 2; iEndcap++) {
431  for (int iStation = 1; iStation < 5; iStation++) {
432  for (int iRing = 1; iRing < 4; iRing++) {
433  if (iStation != 1 && iRing > 2) continue;
434  TString signEndcap="+";
435  if(iEndcap==0) signEndcap="-";
436 
437  char lcttitle[200];
438  snprintf(lcttitle,200,"ME%s%d/%d", signEndcap.Data(), iStation, iRing);
439  if(ihist<=8){
440  csctflcts -> setBinLabel(9-ihist,lcttitle,2);
441  }
442  else csctflcts -> setBinLabel(ihist+1,lcttitle,2);
443 
444  ihist++;
445  }
446  }
447  }
448 
449 
450  // plots for ME1/1 chambers
451  me11_lctStrip = dbe->book1D("CSC_ME11_LCT_Strip", "CSC_ME11_LCT_Strip", 223, 0, 223);
452  me11_lctStrip->setAxisTitle("Cathode HalfStrip, ME1/1", 1);
453 
454  me11_lctWire = dbe->book1D("CSC_ME11_LCT_Wire", "CSC_ME11_LCT_Wire", 112, 0, 112);
455  me11_lctWire->setAxisTitle("Anode Wiregroup, ME1/1", 1);
456 
457  me11_lctLocalPhi = dbe->book1D("CSC_ME11_LCT_LocalPhi", "CSC_ME11_LCT_LocalPhi", 200,0,1024);
458  me11_lctLocalPhi ->setAxisTitle("LCT Local #it{#phi}, ME1/1", 1);
459 
460  me11_lctPackedPhi = dbe->book1D("CSC_ME11_LCT_PackedPhi", "CSC_ME11_LCT_PackedPhi", 200,0,4096);
461  me11_lctPackedPhi ->setAxisTitle("LCT Packed #it{#phi}, ME1/1",1);
462 
463  me11_lctGblPhi = dbe->book1D("CSC_ME11_LCT_GblPhi", "CSC_ME11_LCT_GblPhi", 200, 0, 2*M_PI);
464  me11_lctGblPhi ->setAxisTitle("LCT Global #it{#phi}, ME1/1", 1);
465 
466  me11_lctGblEta = dbe->book1D("CSC_ME11_LCT_GblEta", "CSC_ME11_LCT_GblEta", 50, 0.9, 2.5);
467  me11_lctGblEta ->setAxisTitle("LCT Global #eta, ME1/1", 1);
468 
469 
470  // plots for ME4/2 chambers
471  me42_lctGblPhi = dbe->book1D("CSC_ME42_LCT_GblPhi", "CSC_ME42_LCT_GblPhi", 200, 0, 2*M_PI);
472  me42_lctGblPhi ->setAxisTitle("LCT Global #it{#phi}, ME4/2", 1);
473 
474  me42_lctGblEta = dbe->book1D("CSC_ME42_LCT_GblEta", "CSC_ME42_LCT_GblEta", 50, 0.9, 2.5);
475  me42_lctGblEta ->setAxisTitle("LCT Global #eta, ME4/2", 1);
476 
477 
478 
479 }
480 
481 
482 void L1TCSCTF::endJob(void)
483 {
484 
485  if(verbose_) edm::LogInfo("DataNotFound") << "L1TCSCTF: end job...." << endl;
486  LogInfo("EndJob") << "analyzed " << nev_ << " events";
487 
488  if ( outputFile_.size() != 0 && dbe ) dbe->save(outputFile_);
489 
490  return;
491 }
492 
493 void L1TCSCTF::analyze(const Event& e, const EventSetup& c)
494 {
495 
496  if( c.get< L1MuTriggerScalesRcd > ().cacheIdentifier() != m_scalesCacheID ||
497  c.get< L1MuTriggerPtScaleRcd >().cacheIdentifier() != m_ptScaleCacheID ){
498 
499  ESHandle< L1MuTriggerScales > scales;
500  c.get< L1MuTriggerScalesRcd >().get(scales);
501  ts = scales.product();
502  ESHandle< L1MuTriggerPtScale > ptscales;
503  c.get< L1MuTriggerPtScaleRcd >().get(ptscales);
504  tpts = ptscales.product();
505  m_scalesCacheID = c.get< L1MuTriggerScalesRcd >().cacheIdentifier();
506  m_ptScaleCacheID = c.get< L1MuTriggerPtScaleRcd >().cacheIdentifier();
507 
508  edm::LogInfo("L1TCSCTF") << "Changing triggerscales and triggerptscales...";
509  }
510 
511  int NumCSCTfTracksRep = 0;
512  nev_++;
513  if(verbose_) edm::LogInfo("DataNotFound") << "L1TCSCTF: analyze...." << endl;
514 
515  edm::Handle<L1MuGMTReadoutCollection> pCollection;
516  if( gmtProducer.label() != "null" )
517  { // GMT block
518  e.getByToken(gmtProducerToken_, pCollection);
519  if (!pCollection.isValid())
520  {
521  edm::LogInfo("DataNotFound") << "can't find L1MuGMTReadoutCollection with label "; // << csctfSource_.label() ;
522  return;
523  }
524 
525  L1MuGMTReadoutCollection const* gmtrc = pCollection.product();
526  vector<L1MuGMTReadoutRecord> gmt_records = gmtrc->getRecords();
527  vector<L1MuGMTReadoutRecord>::const_iterator RRItr;
528 
529  // Look if the readout window contains one (and only one CSC cands)
530  // to make it simpler I reject events with more than a CSC cand in the
531  // same readout window
532 
533  // count non-empty candidates in this bx
534  int bxWindow = 0;
535  int nCands = 0;
536 
537  for( RRItr = gmt_records.begin(); RRItr != gmt_records.end(); RRItr++ ) {
538  bxWindow++;
539 
540  // get the csc candidates
541  vector<L1MuRegionalCand> INPCands = RRItr->getCSCCands();
542  vector<L1MuRegionalCand>::const_iterator INPItr;
543 
544  BxInEvent_ = 0;
545  isCSCcand_ = false;
546  int nCandsBx = 0;
547 
548  for( INPItr = INPCands.begin(); INPItr != INPCands.end(); ++INPItr ) {
549  if(!INPItr->empty())
550  {
551  nCandsBx++;
552  nCands++;
553  BxInEvent_ = RRItr->getBxInEvent();
554  if (verbose_) edm::LogInfo("DataNotFound") << "cand " << nCandsBx << " -> assigned CSCTF bx: " << INPItr->bx() << endl;
555  }
556  }
557  if (verbose_)
558  if(nCandsBx) edm::LogInfo("DataNotFound") << nCandsBx << " cands in bx: " << BxInEvent_ << endl;
559  }
560 
561  if (nCands != 1) return;
562  else isCSCcand_ = true;
563  if (verbose_) edm::LogInfo("DataNotFound") << "bxWindow: " << bxWindow << endl;
564 
565  int ncsctftrack = 0;
566  if (verbose_)
567  {
568  edm::LogInfo("DataNotFound") << "\tCSCTFCand ntrack " << ncsctftrack << endl;
569  }
570  } // end of GMT block
571 
572  L1ABXN = -999;
573  if( statusProducer.label() != "null" )
574  {
575  edm::Handle<L1CSCStatusDigiCollection> status;
576  e.getByToken(statusToken_, status);
577  bool integrity=status->first, se=false, sm=false, bx=false, af=false, fmm=false;
578  int nStat = 0;
579 
580  for(std::vector<L1CSCSPStatusDigi>::const_iterator stat=status->second.begin(); stat!=status->second.end(); stat++)
581  {
582  se |= stat->SEs()&0xFFF;
583  sm |= stat->SMs()&0xFFF;
584  bx |= stat->BXs()&0xFFF;
585  af |= stat->AFs()&0xFFF;
586  fmm|= stat->FMM()!=8;
587 
588  if(stat->VPs() != 0)
589  {
590  L1ABXN += stat->BXN();
591  nStat++;
592  }
593  }
594  // compute the average
595  if(nStat!=0) L1ABXN /= nStat;
596  if(integrity) csctferrors->Fill(0.5);
597  if(se) csctferrors->Fill(1.5);
598  if(sm) csctferrors->Fill(2.5);
599  if(bx) csctferrors->Fill(3.5);
600  if(af) csctferrors->Fill(4.5);
601  if(fmm) csctferrors->Fill(5.5);
602  }
603 
604  if( lctProducer.label() != "null" )
605  {
606  edm::ESHandle<CSCGeometry> pDD;
607  c.get<MuonGeometryRecord>().get( pDD );
608  CSCTriggerGeometry::setGeometry(pDD);
609 
610  edm::Handle<CSCCorrelatedLCTDigiCollection> corrlcts;
611  e.getByToken(corrlctsToken_, corrlcts);
612 
613  for(CSCCorrelatedLCTDigiCollection::DigiRangeIterator csc=corrlcts.product()->begin(); csc!=corrlcts.product()->end(); csc++)
614  {
615  CSCCorrelatedLCTDigiCollection::Range range1 = corrlcts.product()->get((*csc).first);
616  for(CSCCorrelatedLCTDigiCollection::const_iterator lct=range1.first; lct!=range1.second; lct++)
617  {
618  int endcap = (*csc).first.endcap()-1;
619  int station = (*csc).first.station()-1;
620  int sector = (*csc).first.triggerSector()-1;
621  int subSector = CSCTriggerNumbering::triggerSubSectorFromLabels((*csc).first);
622  int ring = (*csc).first.ring();
623  int cscId = (*csc).first.triggerCscId()-1;
624  int fpga = ( subSector ? subSector-1 : station+1 );
625  int strip = lct -> getStrip();
626  int keyWire = lct -> getKeyWG();
627  int bx = lct -> getBX();
628 
629  int endcapAssignment = 1;
630  int shift = 1;
631  float sectorArg = sector;
632  //float sectorArg = j;
633 
634  if( endcap == 1 ){
635  endcapAssignment = -1;
636  shift = 2;
637  //sectorArg = sector - 6;
638  }
639 
640  int signedStation = (station + shift)* endcapAssignment;
641  if( (station == 0) && (endcap == 0)) signedStation = subSector - 1;
642  if( (station == 0) && (endcap == 1)) signedStation = (-1)*subSector;
643 
644  float chamberArg1 = cscId * 0.1 + sectorArg;
645  //float chamberArg1 = i*0.1 + sectorArg;
646  //std::cout << "First" << i << " " << sectorArg << " " << chamberArg1 << std::endl;
647 
648  float chamberArg11 = chamberArg1;
649  if(sectorArg == 1) chamberArg1 = chamberArg11 - 0.1;
650  if(sectorArg == 2) chamberArg1 = chamberArg11 - 0.2;
651  if(sectorArg == 3) chamberArg1 = chamberArg11 - 0.3;
652  if(sectorArg == 4) chamberArg1 = chamberArg11 - 0.4;
653  if(sectorArg == 5) chamberArg1 = chamberArg11 - 0.5;
654 
655  //std::cout << "cscId, station, sector, endcap, sectorArg, chamber Arg: " << cscId << ", " << station << ", " <<sector << ", " << endcap << ", " << chamberArg1 << ", " << signedStation << std::endl;
656 
657  csctfChamberOccupancies->Fill(chamberArg1, signedStation);
658  //int bunchX = ( (lct->getBX()) - 6 );
659 
660  //int timingSectorArg = 3*(sector) + (lct->getMPCLink());
661  //if( endcap == 1) timingSectorArg = 3*(sector + 6) + (lct->getMPCLink());
662  //std::cout << "Sector, MPCLink, TSA, endcap: " << sector << ", " << lct->getMPCLink() << ", " << timingSectorArg << ", " << endcap << std::endl;
663 
664  //csctfbx->Fill(timingSectorArg, bunchX );
665  //std::cout << "LCT'S, encap: " << endcap << ", station: " << station << ", sector: " << sector << ", subSector: " << subSector << ", cscId: " << cscId << std:: endl;
666  //End JAG
667 
668  // Check if Det Id is within pysical range:
669  if( endcap<0||endcap>1 || sector<0||sector>6 || station<0||station>3 || cscId<0||cscId>8 || fpga<0||fpga>4)
670  {
671  edm::LogError("L1CSCTF: CSC TP are out of range: ")<<" endcap: "<<(endcap+1)<<" station: "<<(station+1) <<" sector: "<<(sector+1)<<" subSector: "<<subSector<<" fpga: "<<fpga<<" cscId: "<<(cscId+1);
672  continue;
673  }
674  lclphidat lclPhi;
675  try {
676  lclPhi = srLUTs_[fpga]->localPhi(lct->getStrip(), lct->getPattern(), lct->getQuality(), lct->getBend());
677  } catch(cms::Exception &) {
678  bzero(&lclPhi,sizeof(lclPhi));
679  }
680 
681  gblphidat gblPhi;
682  try {
683  gblPhi = srLUTs_[fpga]->globalPhiME(lclPhi.phi_local, lct->getKeyWG(), cscId+1);
684  } catch(cms::Exception &) {
685  bzero(&gblPhi,sizeof(gblPhi));
686  }
687 
688  gbletadat gblEta;
689  try {
690  gblEta = srLUTs_[fpga]->globalEtaME(lclPhi.phi_bend_local, lclPhi.phi_local, lct->getKeyWG(), cscId+1);
691  } catch(cms::Exception &) {
692  bzero(&gblEta,sizeof(gblEta));
693  }
694 
695 
696  //TrackStub
697  csctf::TrackStub theStub((*lct), (*csc).first);
698  theStub.setPhiPacked(gblPhi.global_phi);
699  theStub.setEtaPacked(gblEta.global_eta);
700 
701  float etaG = theStub.etaValue();
702  float phiG = fmod( theStub.phiValue()+15.0*M_PI/180+(sector)*60.0*M_PI/180, 2.*M_PI );
703 
704  // BX plots
705  // minus side
706  if (endcap == 0 && station == 0 && ring == 1) { csctflcts -> Fill(bx, 8.5); }
707  if (endcap == 0 && station == 0 && ring == 2) { csctflcts -> Fill(bx, 7.5); }
708  if (endcap == 0 && station == 0 && ring == 3) { csctflcts -> Fill(bx, 6.5); }
709  if (endcap == 0 && station == 1 && ring == 1) { csctflcts -> Fill(bx, 5.5); }
710  if (endcap == 0 && station == 1 && ring == 2) { csctflcts -> Fill(bx, 4.5); }
711  if (endcap == 0 && station == 2 && ring == 1) { csctflcts -> Fill(bx, 3.5); }
712  if (endcap == 0 && station == 2 && ring == 2) { csctflcts -> Fill(bx, 2.5); }
713  if (endcap == 0 && station == 3 && ring == 1) { csctflcts -> Fill(bx, 1.5); }
714  if (endcap == 0 && station == 3 && ring == 2) { csctflcts -> Fill(bx, 0.5); }
715 
716  // plus side
717  if (endcap == 1 && station == 0 && ring == 1) { csctflcts -> Fill(bx, 9.5); }
718  if (endcap == 1 && station == 0 && ring == 2) { csctflcts -> Fill(bx, 10.5); }
719  if (endcap == 1 && station == 0 && ring == 3) { csctflcts -> Fill(bx, 11.5); }
720  if (endcap == 1 && station == 1 && ring == 1) { csctflcts -> Fill(bx, 12.5); }
721  if (endcap == 1 && station == 1 && ring == 2) { csctflcts -> Fill(bx, 13.5); }
722  if (endcap == 1 && station == 2 && ring == 1) { csctflcts -> Fill(bx, 14.5); }
723  if (endcap == 1 && station == 2 && ring == 2) { csctflcts -> Fill(bx, 15.5); }
724  if (endcap == 1 && station == 3 && ring == 1) { csctflcts -> Fill(bx, 16.5); }
725  if (endcap == 1 && station == 3 && ring == 2) { csctflcts -> Fill(bx, 17.5); }
726 
727  // only for ME1/1
728  if(station == 0 && ring == 1){
729  me11_lctStrip -> Fill(strip);
730  me11_lctWire -> Fill(keyWire);
731  me11_lctLocalPhi -> Fill(lclPhi.phi_local);
732  me11_lctPackedPhi-> Fill(theStub.phiPacked());
733  me11_lctGblPhi -> Fill(phiG);
734  me11_lctGblEta -> Fill(etaG);
735  }
736 
737  // only for ME4/2
738  if(station == 3 && ring == 2){
739  me42_lctGblPhi -> Fill(phiG);
740  me42_lctGblEta -> Fill(etaG);
741  }
742 
743 
744 
745 
746 
747  // SR LUT gives packed eta and phi values -> normilize them to 1 by scale them to 'max' and shift by 'min'
748  //float etaP = gblEta.global_eta/127*1.5 + 0.9;
749  //float phiP = (gblPhi.global_phi);// + ( sector )*4096 + station*4096*12) * 1./(4*4096*12);
750  //std::cout << "LCT Eta & Phi Coordinates: " << etaP << ", " << phiP << "." << std::endl;
751  //csctfoccupancies->Fill( gblEta.global_eta/127. * 1.5 + 0.9, (gblPhi.global_phi + ( sector + (endcap?0:6) )*4096 + station*4096*12) * 1./(4*4096*12) );
752  }//lct != range1.scond
753  }//csc!=corrlcts.product()->end()
754  }// lctProducer.label() != "null"
755 
756 
757 
758  if( trackProducer.label() != "null" )
759  {
760  edm::Handle<L1CSCTrackCollection> tracks;
761  e.getByToken(tracksToken_, tracks);
762  for(L1CSCTrackCollection::const_iterator trk=tracks->begin(); trk<tracks->end(); trk++)
763  {
764 
765  NumCSCTfTracksRep++;
766  long LUTAdd = trk->first.ptLUTAddress();
767  int trigMode = ( (LUTAdd)&0xf0000 ) >> 16;
768  int trEta = (trk->first.eta_packed() );
769 
770 
771  // trk->first.endcap() = 2 for - endcap
772  // = 1 for + endcap
773  //int trEndcap = (trk->first.endcap()==2 ? trk->first.endcap()-3 : trk->first.endcap());
774  if( trk->first.endcap() != 1)
775  {
776  int holder = trEta;
777  trEta = -1*holder;
778  trEta -= 1;
779  }
780 
781  int trSector = 6*(trk->first.endcap()-1)+trk->first.sector();
782  int trBX = trk->first.BX();
783 
784  //Here is what is done with output phi value:
785  //output_phi = (phi / 32) * 3 /16
786  //where:
787  //phi is 12-bit phi, 4096 bins covering 62 degrees
788  //output_phi is 5-bit value
789 
790  //Easy to see that output_phi can have values from 0 to 23, or 24 total combinations.
791  //This gives per-bin phi value of 62/24 = 2.583333 degrees.
792 
793  // Sector 1 nominally starts at 15 degrees but there 1 degree overlap between sectors so 14 degrees effectively
794  //double trPhi = trk->first.localPhi() * 62. / 24.;
795  double trPhi = ts->getPhiScale()->getLowEdge(trk->first.localPhi());
796  double trPhi02PI = fmod(trPhi +
797  ((trSector-1)*M_PI/3) +
798  (M_PI*14/180.), 2*M_PI);
799 
800  if (trigMode == 15) {
801  csctfTrackPhi_H -> Fill( trPhi02PI );
802  csctfTrackEta_H -> Fill( trEta );
803  csctfoccupancies_H -> Fill( trEta, trPhi02PI );
804  csctfbx_H -> Fill( trSector, trBX );
805  }
806  else{
807  csctfTrackPhi -> Fill( trPhi02PI );
808  csctfTrackEta -> Fill( trEta );
809  csctfoccupancies -> Fill( trEta, trPhi02PI );
810  csctfbx -> Fill( trSector, trBX );
811 
812  // Low Quality / High Quality Eta Distributions
813  //|eta| < 2.1
814  if (abs(trEta) < 24) {
815  if (trigMode == 2 ||
816  trigMode == 3 ||
817  trigMode == 4 ||
818  trigMode == 5 ||
819  trigMode == 6 ||
820  trigMode == 7 ||
821  trigMode == 11 ||
822  trigMode == 12 ||
823  trigMode == 13 ||
824  trigMode == 14 ) csctfTrackEtaHighQ -> Fill (trEta);
825 
826  if (trigMode == 8 ||
827  trigMode == 9 ||
828  trigMode == 10 ) csctfTrackEtaLowQ -> Fill (trEta);
829  }
830  else {//|eta| > 2.1
831  if (trigMode == 2 ||
832  trigMode == 3 ||
833  trigMode == 4 ||
834  trigMode == 5 ) csctfTrackEtaHighQ -> Fill (trEta);
835  else
836  csctfTrackEtaLowQ -> Fill (trEta);
837  }
838  }
839 
840  csctfTrackM->Fill( trk->first.modeExtended() );
841 
842  // we monitor the track quality only on the first link
843  // so let's make sure to fill the plot if there is something that
844  // is read from the hardware
845  int trRank = trk->first.rank();
846  if (trRank) {
847  int trQuality = ((trRank>>5)&0x3);
848  trackModeVsQ->Fill( trk->first.modeExtended(), trQuality );
849  }
850 
851  /*
852  OLD METHOD FOR FILLING HALO PLOTS, IMPROVED METHOD USING ASSOCIATED TRACK STUBS
853  BELOW ~LINE 605
854  if( trigMode == 15 )
855  {
856 
857  double haloVals[4][4];
858  for( int i = 0; i < 4; i++)
859  {
860  haloVals[i][0] = 0;
861  }
862 
863  edm::Handle<CSCCorrelatedLCTDigiCollection> corrlcts;
864  for(CSCCorrelatedLCTDigiCollection::DigiRangeIterator csc=corrlcts.product()->begin(); csc!=corrlcts.product()->end(); csc++)
865  {
866  CSCCorrelatedLCTDigiCollection::Range range1 = corrlcts.product()->get((*csc).first);
867  for(CSCCorrelatedLCTDigiCollection::const_iterator lct=range1.first; lct!=range1.second; lct++)
868  {
869  int endcap = (*csc).first.endcap()-1;
870  int station = (*csc).first.station()-1;
871  int sector = (*csc).first.triggerSector()-1;
872  int cscId = (*csc).first.triggerCscId()-1;
873  int subSector = CSCTriggerNumbering::triggerSubSectorFromLabels((*csc).first);
874  int fpga = ( subSector ? subSector-1 : station+1 );
875 
876  if(station != 4)
877  {
878  int modEnd = 1;
879  if( endcap == 0 ) modEnd = -1;
880  int indexHalo = modEnd + station;
881  if(haloVals[indexHalo][0] == 1.0) haloVals[indexHalo][3] = 1.0;
882  if(haloVals[indexHalo][0] == 0) haloVals[indexHalo][0] = 1.0;
883  haloVals[indexHalo][1] = sector*1.0;
884 
885  lclphidat lclPhi;
886  lclPhi = srLUTs_[fpga]->localPhi(lct->getStrip(), lct->getPattern(), lct->getQuality(), lct->getBend());
887  gblphidat gblPhi;
888  gblPhi = srLUTs_[fpga]->globalPhiME(lclPhi.phi_local, lct->getKeyWG(), cscId+1);
889  gbletadat gblEta;
890  gblEta = srLUTs_[fpga]->globalEtaME(lclPhi.phi_bend_local, lclPhi.phi_local, lct->getKeyWG(), cscId+1);
891 
892  haloVals[indexHalo][2] = gblEta.global_eta/127. * 1.5 + 0.9;
893  } //station1 or 2
894  } //lct first to second
895  } //corrlcts
896 
897  if( (haloVals[0][0] == 1.) && (haloVals[1][0] == 1.) && (haloVals[0][3] != 1.) && (haloVals[1][3] != 1.) )
898  {
899  if( haloVals[0][1] == haloVals[1][1] ){
900  double delEta23 = haloVals[1][2] - haloVals[0][2];
901  haloDelEta23->Fill( delEta23 );
902  }
903  }
904 
905  if( (haloVals[2][0] == 1.) && (haloVals[3][0] == 1.) && (haloVals[2][3] != 1.) && (haloVals[3][3] != 1.) )
906  {
907  if( haloVals[2][1] == haloVals[3][1] ){
908  double delEta23 = haloVals[3][2] - haloVals[2][2];
909  haloDelEta23->Fill( delEta23 );
910  }
911  }
912  } //halo trigger
913  */
914 
915  int cscTrackStub = 0;
916  //float haloEta[3];
917  //for(int i=0; i<3; i++) haloEta[i]=-1.0;
918  //bool haloME11 = false;
919  CSCCorrelatedLCTDigiCollection lctsOfTracks=trk->second;
920  for(CSCCorrelatedLCTDigiCollection::DigiRangeIterator trackStub=lctsOfTracks.begin(); trackStub!=lctsOfTracks.end(); trackStub++)
921  {
922  CSCCorrelatedLCTDigiCollection::Range range2 = lctsOfTracks.get((*trackStub).first);
923  for(CSCCorrelatedLCTDigiCollection::const_iterator lct=range2.first; lct!=range2.second; lct++)
924  {
925 // int station = (*trackStub).first.station()-1;
926 // if(station != 4)
927 // {
928 // // int endcap = (*trackStub).first.endcap()-1;
929 // // int sector = (*trackStub).first.triggerSector()-1;
930 // int cscId = (*trackStub).first.triggerCscId()-1;
931 // int subSector = CSCTriggerNumbering::triggerSubSectorFromLabels((*trackStub).first);
932 // int fpga = ( subSector ? subSector-1 : station+1 );
933 
934 // lclphidat lclPhi;
935 // lclPhi = srLUTs_[fpga]->localPhi(lct->getStrip(), lct->getPattern(), lct->getQuality(), lct->getBend());
936 // gblphidat gblPhi;
937 // gblPhi = srLUTs_[fpga]->globalPhiME(lclPhi.phi_local, lct->getKeyWG(), cscId+1);
938 // gbletadat gblEta;
939 // gblEta = srLUTs_[fpga]->globalEtaME(lclPhi.phi_bend_local, lclPhi.phi_local, lct->getKeyWG(), cscId+1);
940 // haloEta[station-1] = gblEta.global_eta/127. * 1.5 + 0.9;
941 // if(station==1 && cscId<2) haloME11 = true;
942 // }
943  cscTrackStub++;
944  }
945  }
946  cscTrackStubNumbers->Fill(cscTrackStub);
947 
948 // if(trigMode == 15)
949 // {
950 // float dEta13 = haloEta[2]-haloEta[0];
951 // float dEta12 = haloEta[1]-haloEta[0];
952 // if(haloME11)
953 // {
954 // if(haloEta[1]!=-1.0) haloDelEta112->Fill(dEta12);
955 // if(haloEta[2]!=-1.0) haloDelEta113->Fill(dEta13);
956 // } else {
957 // if(haloEta[1]!=-1.0) haloDelEta12->Fill(dEta12);
958 // if(haloEta[2]!=-1.0) haloDelEta13->Fill(dEta13);
959 // }
960 // }
961  //
962 
963 
964 
965  }
966  }
967  csctfntrack->Fill(NumCSCTfTracksRep);
968 
969 
970  if( mbProducer.label() != "null" )
971  {
972  // handle to needed collections
973  edm::Handle<CSCTriggerContainer<csctf::TrackStub> > dtStubs;
974  e.getByToken(dtStubsToken_, dtStubs);
975  edm::Handle<L1CSCTrackCollection> tracks;
976  e.getByToken(mbtracksToken_, tracks);
977 
978  // loop on the DT stubs
979  std::vector<csctf::TrackStub> vstubs = dtStubs->get();
980  for(std::vector<csctf::TrackStub>::const_iterator stub=vstubs.begin();
981  stub!=vstubs.end(); stub++)
982  {
983  if (verbose_)
984  {
985  edm::LogInfo("DataNotFound") << "\n mbEndcap: " << stub->endcap();
986  edm::LogInfo("DataNotFound") << "\n stub->getStrip()[FLAG]: " << stub->getStrip();
987  edm::LogInfo("DataNotFound") << "\n stub->getKeyWG()[CAL]: " << stub->getKeyWG();
988  edm::LogInfo("DataNotFound") << "\n stub->BX(): " << stub->BX();
989  edm::LogInfo("DataNotFound") << "\n stub->sector(): " << stub->sector();
990  edm::LogInfo("DataNotFound") << "\n stub->subsector(): " << stub->subsector();
991  edm::LogInfo("DataNotFound") << "\n stub->station(): " << stub->station();
992  edm::LogInfo("DataNotFound") << "\n stub->phiPacked(): " << stub->phiPacked();
993  edm::LogInfo("DataNotFound") << "\n stub->getBend(): " << stub->getBend();
994  edm::LogInfo("DataNotFound") << "\n stub->getQuality(): " << stub->getQuality();
995  edm::LogInfo("DataNotFound") << "\n stub->cscid(): " << stub->cscid() << endl;
996  }
997  // define the sector ID
998  int mbId = (stub->endcap()==2) ? 6 : 0;
999  mbId += stub->sector();
1000  // *** do not fill if CalMB variable is set ***
1001  // horrible! They used the same class to write up the LCT and MB info,
1002  // but given the MB does not have strip and WG they replaced this two
1003  // with the flag and cal bits... :S
1004  if (stub->getKeyWG() == 0)
1005  {
1006  // if FLAG =1, muon belong to previous BX
1007  int bxDT = stub->BX()-stub->getStrip(); // correct by the FLAG
1008  int subDT = stub->subsector();
1009 
1010  // Fill the event only if CSC had or would have triggered
1011  if (isCSCcand_)
1012  {
1013  //look for tracks in the event and compare the matching DT stubs
1014  int trkBX = 0;
1015  for(L1CSCTrackCollection::const_iterator trk=tracks->begin(); trk<tracks->end(); trk++)
1016  {
1017  trkBX = trk->first.BX();
1018  int trkId = (trk->first.endcap()==2) ? 6 : 0;
1019  trkId += trk->first.sector();
1020  if (verbose_){
1021  edm::LogInfo("DataNotFound") << "\n trk BX: " << trkBX
1022  << " Sector: " << trkId
1023  << " SubSector: " << trk->first.subsector()
1024  << " Endcap: " << trk->first.endcap();
1025 
1026  edm::LogInfo("DataNotFound") << "\n DT BX: " << stub->BX()
1027  << " Sector: " << mbId
1028  << " SubSector: " << stub->subsector()
1029  << " Endcap: " << stub->endcap() << endl;
1030  }
1031 
1032  if (mbId == trkId)
1033  {
1034  if (verbose_) {
1035  edm::LogInfo("DataNotFound") << " --> MATCH" << endl;
1036  edm::LogInfo("DataNotFound") << "Fill :" << trkBX+6-bxDT << " -- " << subDT << " -- cands" << endl;
1037  }
1038  // DT bx ranges from 3 to 9
1039  // trk bx ranges from -3 to 3
1040  DTstubsTimeTrackMenTimeArrival[mbId-1]->Fill(bxDT-trkBX-6,subDT);//subsec
1041  }
1042  }// loop on the tracks
1043  }//if (isCSCcand_){
1044  }//if (stub->getKeyWG() == 0) {
1045  }
1046  }
1047 }
L1TCSCTF::DTstubsTimeTrackMenTimeArrival
MonitorElement * DTstubsTimeTrackMenTimeArrival[12]
Definition: L1TCSCTF.h:124
L1TCSCTF::csctfoccupancies
MonitorElement * csctfoccupancies
Definition: L1TCSCTF.h:85
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
edm::eventsetup::EventSetupRecord::cacheIdentifier
unsigned long long cacheIdentifier() const
Definition: EventSetupRecord.h:148
L1TCSCTF::csctfTrackEtaHighQ
MonitorElement * csctfTrackEtaHighQ
Definition: L1TCSCTF.h:97
edm::ParameterSet::getUntrackedParameter
T getUntrackedParameter(std::string const &, T const &) const
i
int i
Definition: DBlmapReader.cc:9
L1TCSCTF::csctflcts
MonitorElement * csctflcts
Definition: L1TCSCTF.h:106
L1TCSCTF::dbe
DQMStore * dbe
Definition: L1TCSCTF.h:78
L1TCSCTF::csctfTrackPhi
MonitorElement * csctfTrackPhi
Definition: L1TCSCTF.h:94
L1TCSCTF::verbose_
bool verbose_
Definition: L1TCSCTF.h:132
L1TCSCTF::m_scalesCacheID
unsigned long long m_scalesCacheID
Definition: L1TCSCTF.h:141
L1TCSCTF::csctfTrackPhi_H
MonitorElement * csctfTrackPhi_H
Definition: L1TCSCTF.h:98
DQMStore::book1D
MonitorElement * book1D(const char *name, const char *title, int nchX, double lowX, double highX)
Book 1D histogram.
Definition: DQMStore.cc:954
DQMStore::rmdir
void rmdir(const std::string &fullpath)
Definition: DQMStore.cc:3101
L1TCSCTF::beginJob
void beginJob(void)
Definition: L1TCSCTF.cc:105
edm::Event::getByToken
bool getByToken(EDGetToken token, Handle< PROD > &result) const
Definition: Event.h:434
Reference_intrackfit_cff.endcap
endcap
Definition: Reference_intrackfit_cff.py:82
L1MuScale::getLowEdge
virtual float getLowEdge(unsigned packed) const =0
get the low edge of bin represented by packed
L1TCSCTF::me11_lctGblPhi
MonitorElement * me11_lctGblPhi
Definition: L1TCSCTF.h:113
MonitorElement::setBinLabel
void setBinLabel(int bin, const std::string &label, int axis=1)
set bin label for x, y or z axis (axis=1, 2, 3 respectively)
Definition: MonitorElement.cc:942
NULL
#define NULL
Definition: scimark2.h:8
L1TCSCTF::csctfbx
MonitorElement * csctfbx
Definition: L1TCSCTF.h:81
L1TCSCTF::csctfTrackEtaLowQ
MonitorElement * csctfTrackEtaLowQ
Definition: L1TCSCTF.h:96
L1TCSCTF::mbProducer
edm::InputTag mbProducer
Definition: L1TCSCTF.h:135
CSCSectorReceiverLUT::globalEtaME
gbletadat globalEtaME(int phi_bend, int phi_local, int wire_group, int cscid) const
Definition: CSCSectorReceiverLUT.cc:772
CSCTriggerGeometry::setGeometry
static void setGeometry(const edm::ESHandle< CSCGeometry > &thegeom)
Definition: CSCTriggerGeometry.h:19
L1TCSCTF::csctfTrackEta
MonitorElement * csctfTrackEta
Definition: L1TCSCTF.h:95
L1TCSCTF::analyze
void analyze(const edm::Event &e, const edm::EventSetup &c)
Definition: L1TCSCTF.cc:493
L1TCSCTF::me11_lctWire
MonitorElement * me11_lctWire
Definition: L1TCSCTF.h:110
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:255
L1TCSCTF::outputFile_
std::string outputFile_
Definition: L1TCSCTF.h:131
csctf::TrackStub::phiPacked
unsigned phiPacked() const
Return the binned phi for this stub.
Definition: TrackStub.h:44
MonitorElement::Fill
void Fill(long long x)
Definition: MonitorElement.h:166
L1TCSCTF::statusProducer
edm::InputTag statusProducer
Definition: L1TCSCTF.h:135
L1TCSCTF::statusToken_
edm::EDGetTokenT< L1CSCStatusDigiCollection > statusToken_
Definition: L1TCSCTF.h:146
L1TCSCTF::csctfTrackEta_H
MonitorElement * csctfTrackEta_H
Definition: L1TCSCTF.h:99
L1TCSCTF::me42_lctGblPhi
MonitorElement * me42_lctGblPhi
Definition: L1TCSCTF.h:117
L1TCSCTF::ts
const L1MuTriggerScales * ts
Definition: L1TCSCTF.h:139
L1TCSCTF::cscTrackStubNumbers
MonitorElement * cscTrackStubNumbers
Definition: L1TCSCTF.h:100
CSCSectorReceiverLUT::localPhi
lclphidat localPhi(int strip, int pattern, int quality, int lr) const
Geometry Lookup Tables.
Definition: CSCSectorReceiverLUT.cc:205
L1TCSCTF::L1ABXN
int L1ABXN
Definition: L1TCSCTF.h:128
L1TCSCTF::me11_lctStrip
MonitorElement * me11_lctStrip
Definition: L1TCSCTF.h:109
HcalObjRepresent::Fill
void Fill(HcalDetId &id, double val, std::vector< TH2F > &depth)
Definition: HcalObjRepresent.h:611
L1TCSCTF::BxInEvent_
int BxInEvent_
Definition: L1TCSCTF.h:125
L1TCSCTF::~L1TCSCTF
virtual ~L1TCSCTF()
Definition: L1TCSCTF.cc:98
csctf::TrackStub::etaValue
double etaValue() const
return the Eta Value of this stub&#39;s position.
Definition: TrackStub.h:36
L1TCSCTF::endJob
void endJob(void)
Definition: L1TCSCTF.cc:482
L1MuTriggerScales::getPhiScale
const L1MuScale * getPhiScale() const
get the phi scale
Definition: L1MuTriggerScales.h:201
funct::abs
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
L1TCSCTF::nev_
int nev_
Definition: L1TCSCTF.h:130
L1TCSCTF::L1TCSCTF
L1TCSCTF(const edm::ParameterSet &ps)
Definition: L1TCSCTF.cc:25
DQMStore::setVerbose
void setVerbose(unsigned level)
Definition: DQMStore.cc:631
L1TCSCTF::csctfChamberOccupancies
MonitorElement * csctfChamberOccupancies
Definition: L1TCSCTF.h:93
edm::HandleBase::isValid
bool isValid() const
Definition: HandleBase.h:76
L1TCSCTF::lctProducer
edm::InputTag lctProducer
Definition: L1TCSCTF.h:135
M_PI
#define M_PI
Definition: readTestVector.cc:25
L1TCSCTF::tracksToken_
edm::EDGetTokenT< L1CSCTrackCollection > tracksToken_
Definition: L1TCSCTF.h:148
L1TCSCTF::me42_lctGblEta
MonitorElement * me42_lctGblEta
Definition: L1TCSCTF.h:118
L1TCSCTF::srLUTs_
CSCSectorReceiverLUT * srLUTs_[5]
Definition: L1TCSCTF.h:137
L1TCSCTF::corrlctsToken_
edm::EDGetTokenT< CSCCorrelatedLCTDigiCollection > corrlctsToken_
Definition: L1TCSCTF.h:147
csctf::TrackStub::phiValue
double phiValue() const
return the Phi Value of this stub&#39;s position in local coordinates.
Definition: TrackStub.h:38
L1TCSCTF::csctferrors
MonitorElement * csctferrors
Definition: L1TCSCTF.h:84
testEve_cfg.tracks
tuple tracks
Definition: testEve_cfg.py:39
edm::ParameterSet::addUntrackedParameter
void addUntrackedParameter(std::string const &name, T const &value)
Definition: ParameterSet.h:206
L1TCSCTF::isCSCcand_
bool isCSCcand_
Definition: L1TCSCTF.h:126
gblphidat
class global_phi_data gblphidat
edm::EventSetup::get
const T & get() const
Definition: EventSetup.h:55
L1TCSCTF::csctfTrackM
MonitorElement * csctfTrackM
Definition: L1TCSCTF.h:101
L1TCSCTF::mbtracksToken_
edm::EDGetTokenT< L1CSCTrackCollection > mbtracksToken_
Definition: L1TCSCTF.h:150
edm::ESHandle::product
T const * product() const
Definition: ESHandle.h:62
L1TCSCTF::csctfbx_H
MonitorElement * csctfbx_H
Definition: L1TCSCTF.h:82
MuonDigiCollection< CSCDetId, CSCCorrelatedLCTDigi >::const_iterator
std::vector< CSCCorrelatedLCTDigi >::const_iterator const_iterator
Definition: MuonDigiCollection.h:110
lclphidat
class local_phi_data lclphidat
Data Types.
CSCTriggerNumbering::triggerSubSectorFromLabels
static int triggerSubSectorFromLabels(int station, int chamber)
Definition: CSCTriggerNumbering.cc:118
edm::Handle::product
T const * product() const
Definition: Handle.h:81
L1TCSCTF::csctfntrack
MonitorElement * csctfntrack
Definition: L1TCSCTF.h:80
L1TCSCTF::trackProducer
edm::InputTag trackProducer
Definition: L1TCSCTF.h:135
L1TCSCTF::gmtProducerToken_
edm::EDGetTokenT< L1MuGMTReadoutCollection > gmtProducerToken_
Definition: L1TCSCTF.h:145
edm::InputTag::label
std::string const & label() const
Definition: InputTag.h:42
L1TCSCTF::trackModeVsQ
MonitorElement * trackModeVsQ
Definition: L1TCSCTF.h:102
L1TCSCTF::me11_lctGblEta
MonitorElement * me11_lctGblEta
Definition: L1TCSCTF.h:114
DQMStore::save
void save(const std::string &filename, const std::string &path="", const std::string &pattern="", const std::string &rewrite="", const uint32_t run=0, const uint32_t lumi=0, SaveReferenceTag ref=SaveWithReference, int minStatus=dqm::qstatus::STATUS_OK, const std::string &fileupdate="RECREATE", const bool resetMEsAfterWriting=false)
Definition: DQMStore.cc:2540
L1MuGMTReadoutCollection::getRecords
std::vector< L1MuGMTReadoutRecord > const & getRecords() const
Definition: L1MuGMTReadoutCollection.h:51
L1TCSCTF::csctfoccupancies_H
MonitorElement * csctfoccupancies_H
Definition: L1TCSCTF.h:86
edm::shift
static unsigned int const shift
Definition: LuminosityBlockID.cc:9
csctf::TrackStub::setPhiPacked
void setPhiPacked(const unsigned &phi_)
Definition: TrackStub.h:33
MuonDigiCollection< CSCDetId, CSCCorrelatedLCTDigi >::Range
std::pair< const_iterator, const_iterator > Range
Definition: MuonDigiCollection.h:111
L1TCSCTF::tpts
const L1MuTriggerPtScale * tpts
Definition: L1TCSCTF.h:140
L1TCSCTF::me11_lctLocalPhi
MonitorElement * me11_lctLocalPhi
Definition: L1TCSCTF.h:111
L1TCSCTF::m_ptScaleCacheID
unsigned long long m_ptScaleCacheID
Definition: L1TCSCTF.h:142
ntuplemaker.status
tuple status
Definition: ntuplemaker.py:245
MonitorElement::getTH2F
TH2F * getTH2F(void) const
Definition: MonitorElement.cc:1563
DQMStore::book2D
MonitorElement * book2D(const char *name, const char *title, int nchX, double lowX, double highX, int nchY, double lowY, double highY)
Book 2D histogram.
Definition: DQMStore.cc:1082
gbletadat
class global_eta_data gbletadat
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:969
edm::InputTag::instance
std::string const & instance() const
Definition: InputTag.h:43
L1TCSCTF::dtStubsToken_
edm::EDGetTokenT< CSCTriggerContainer< csctf::TrackStub > > dtStubsToken_
Definition: L1TCSCTF.h:149
L1TCSCTF::gmtProducer
edm::InputTag gmtProducer
Definition: L1TCSCTF.h:135
csctf::TrackStub::setEtaPacked
void setEtaPacked(const unsigned &eta_)
set Eta and Phi from integer values.
Definition: TrackStub.h:32
DQMStore::setCurrentFolder
void setCurrentFolder(const std::string &fullpath)
Definition: DQMStore.cc:667
L1TCSCTF::me11_lctPackedPhi
MonitorElement * me11_lctPackedPhi
Definition: L1TCSCTF.h:112
CSCSectorReceiverLUT::globalPhiME
gblphidat globalPhiME(int phi_local, int wire_group, int cscid) const
Definition: CSCSectorReceiverLUT.cc:502
Generated for CMSSW Reference Manual by doxygen 1.8.5