test
CMS 3D CMS Logo

 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
L1MuGlobalMuonTrigger.cc
Go to the documentation of this file.
1 //---------------------------------------------
2 //
3 // \class L1MuGlobalMuonTrigger
4 //
5 // Description: L1 Global Muon Trigger
6 //
7 //
8 //
9 // Author :
10 // Norbert Neumeister CERN EP
11 // Hannes Sakulin HEPHY Vienna
12 // Ivan Mikulec HEPHY Vienna
13 //
14 //--------------------------------------------------
15 
16 //-----------------------
17 // This Class's Header --
18 //-----------------------
19 
21 
22 //---------------
23 // C++ Headers --
24 //---------------
25 
26 #include <iostream>
27 #include <sys/stat.h>
28 #include <sys/types.h>
29 
30 //-------------------------------
31 // Collaborating Class Headers --
32 //-------------------------------
33 
45 
47 
49 
62 
65 
66 //----------------
67 // Constructors --
68 //----------------
70  produces<std::vector<L1MuGMTCand> >();
71  produces<L1MuGMTReadoutCollection>();
72  m_sendMipIso = ps.getUntrackedParameter<bool>("SendMipIso",false);
73  if( m_sendMipIso ) {
74  produces<std::vector<unsigned> >();
75  }
76 
83 
84  m_ExtendedCands.reserve(20);
85 
86  // set configuration parameters
87  if(!m_config) m_config = new L1MuGMTConfig(ps);
88  m_writeLUTsAndRegs = ps.getUntrackedParameter<bool>("WriteLUTsAndRegs",false);
89 
90  // build GMT
91  if ( L1MuGMTConfig::Debug(1) ) edm::LogVerbatim("GMT_info");
92  if ( L1MuGMTConfig::Debug(1) ) edm::LogVerbatim("GMT_info") << "**** L1GlobalMuonTrigger building ****";
93  if ( L1MuGMTConfig::Debug(1) ) edm::LogVerbatim("GMT_info");
94 
95  // create new PSB
96  if ( L1MuGMTConfig::Debug(2) ) edm::LogVerbatim("GMT_info") << "creating GMT PSB";
97  m_PSB = new L1MuGMTPSB(*this,consumesCollector());
98 
99  // create new matcher
100  if ( L1MuGMTConfig::Debug(2) ) edm::LogVerbatim("GMT_info") << "creating GMT Matcher (0,1)";
101  m_Matcher[0] = new L1MuGMTMatcher(*this,0); // barrel
102  m_Matcher[1] = new L1MuGMTMatcher(*this,1); // endcap
103 
104  // create new cancel-out units
105  if ( L1MuGMTConfig::Debug(2) ) edm::LogVerbatim("GMT_info") << "creating GMT Cancel Out Unit (0,1,2,3)";
106  m_CancelOutUnit[0] = new L1MuGMTCancelOutUnit(*this,0); // barrel
107  m_CancelOutUnit[1] = new L1MuGMTCancelOutUnit(*this,1); // endcap
108  m_CancelOutUnit[2] = new L1MuGMTCancelOutUnit(*this,2); // CSC/bRPC
109  m_CancelOutUnit[3] = new L1MuGMTCancelOutUnit(*this,3); // DT/fRPC
110 
111  // create new MIP & ISO bit assignment units
112  if ( L1MuGMTConfig::Debug(2) ) edm::LogVerbatim("GMT_info") << "creating GMT MIP & ISO bit Assigment Unit (0,1)";
113  m_MipIsoAU[0] = new L1MuGMTMipIsoAU(*this,0); // barrel
114  m_MipIsoAU[1] = new L1MuGMTMipIsoAU(*this,1); // endcap
115 
116  // create new Merger
117  if ( L1MuGMTConfig::Debug(2) ) edm::LogVerbatim("GMT_info") << "creating GMT Merger (0,1)";
118  m_Merger[0] = new L1MuGMTMerger(*this,0); // barrel
119  m_Merger[1] = new L1MuGMTMerger(*this,1); // endcap
120 
121  // create new sorter
122  if ( L1MuGMTConfig::Debug(2) ) edm::LogVerbatim("GMT_info") << "creating GMT Sorter";
123  m_Sorter = new L1MuGMTSorter(*this); // barrel
124 
126  usesResource("L1MuGlobalMuonTrigger");
127 }
128 
129 //--------------
130 // Destructor --
131 //--------------
133 
134  if(m_db) delete m_db;
135  m_db = 0;
136 
137  delete m_Sorter;
138  delete m_Merger[1]; // endcap Merger
139  delete m_Merger[0]; // barrel Merger
140  delete m_MipIsoAU[1]; // barrel MIP & ISO bit assignment unit
141  delete m_MipIsoAU[0]; // barrel MIP & ISO bit assignment unit
142  delete m_CancelOutUnit[3]; // DT/fRPC cancel-out unit (in endcap chip)
143  delete m_CancelOutUnit[2]; // CSC/bRPC cancel-out unit (in barrel chip)
144  delete m_CancelOutUnit[1]; // endcap DT/CSC cancel out unit
145  delete m_CancelOutUnit[0]; // barrel DT/CSC cancel out unit
146  delete m_Matcher[1]; // endcap matcher
147  delete m_Matcher[0]; // barrel matcher
148  delete m_PSB;
149 
150  if(m_config) delete m_config;
151  m_config = 0;
152 
153  // copied from produce() by Jim B, 7 Aug 2007
154  std::vector<L1MuGMTReadoutRecord*>::iterator irr = m_ReadoutRingbuffer.begin();
155  for ( ;irr!=m_ReadoutRingbuffer.end(); irr++) delete (*irr);
156  m_ReadoutRingbuffer.clear();
157  // end Jim B edit
158 
159 }
160 
161 //--------------
162 // Operations --
163 //--------------
164 
166 
167 }
168 
170 
171  // configure from the event setup
172 
173  unsigned long long L1MuGMTScalesCacheID = es.get< L1MuGMTScalesRcd >().cacheIdentifier();
174  if(L1MuGMTScalesCacheID != m_L1MuGMTScalesCacheID) {
175  edm::ESHandle< L1MuGMTScales > gmtscales_h;
176  es.get< L1MuGMTScalesRcd >().get( gmtscales_h );
177  m_config->setGMTScales( gmtscales_h.product() );
178  }
179 
180  unsigned long long L1MuTriggerScalesCacheID = es.get< L1MuTriggerScalesRcd >().cacheIdentifier();
181  if(L1MuTriggerScalesCacheID != m_L1MuTriggerScalesCacheID) {
183  es.get< L1MuTriggerScalesRcd >().get( trigscales_h );
184  m_config->setTriggerScales( trigscales_h.product() );
185  }
186 
187  unsigned long long L1MuTriggerPtScaleCacheID = es.get< L1MuTriggerPtScaleRcd >().cacheIdentifier();
188  if(L1MuTriggerPtScaleCacheID != m_L1MuTriggerPtScaleCacheID) {
190  es.get< L1MuTriggerPtScaleRcd >().get( trigptscale_h );
191  m_config->setTriggerPtScale( trigptscale_h.product() );
192  }
193 
194  unsigned long long L1MuGMTParametersCacheID = es.get< L1MuGMTParametersRcd >().cacheIdentifier();
195  if(L1MuGMTParametersCacheID != m_L1MuGMTParametersCacheID) {
197  es.get< L1MuGMTParametersRcd >().get( gmtparams_h );
198  m_config->setGMTParams( gmtparams_h.product() );
200  }
201 
202  unsigned long long L1MuGMTChannelMaskCacheID = es.get< L1MuGMTChannelMaskRcd >().cacheIdentifier();
203  if(L1MuGMTChannelMaskCacheID != m_L1MuGMTChannelMaskCacheID) {
205  es.get< L1MuGMTChannelMaskRcd >().get( gmtchanmask_h );
206  m_config->setGMTChanMask( gmtchanmask_h.product() );
207  if ( L1MuGMTConfig::Debug(1) ) {
208  std::string onoff;
209  const L1MuGMTChannelMask* theChannelMask = L1MuGMTConfig::getGMTChanMask();
210  unsigned mask = theChannelMask->getSubsystemMask();
211 
212  edm::LogVerbatim("GMT_info");
213  edm::LogVerbatim("GMT_info") << " GMT input Channel Mask:" << std::hex << mask << std::dec;
214  onoff = mask&1 ? "OFF" : "ON";
215  edm::LogVerbatim("GMT_info") << " DT input " << onoff;
216  onoff = mask&2 ? "OFF" : "ON";
217  edm::LogVerbatim("GMT_info") << " RPCb input " << onoff;
218  onoff = mask&4 ? "OFF" : "ON";
219  edm::LogVerbatim("GMT_info") << " CSC input " << onoff;
220  onoff = mask&8 ? "OFF" : "ON";
221  edm::LogVerbatim("GMT_info") << " RPCf input " << onoff;
222  edm::LogVerbatim("GMT_info");
223  }
224  }
225 
226  unsigned long long L1CaloGeometryCacheID = es.get< L1CaloGeometryRecord >().cacheIdentifier();
227  if(L1CaloGeometryCacheID != m_L1CaloGeometryCacheID) {
229  es.get< L1CaloGeometryRecord >().get( caloGeom_h ) ;
230  m_config->setCaloGeom( caloGeom_h.product() ) ;
231  }
232 
234 
235  // write LUTs and Regs if required
236 
237  if(m_writeLUTsAndRegs) {
238  std::string dir = "gmtconfig";
239 
240  mkdir(dir.c_str(), S_ISUID|S_ISGID|S_ISVTX|S_IRUSR|S_IWUSR|S_IXUSR);
241 
242  m_config->dumpLUTs(dir);
243  m_config->dumpRegs(dir);
244  }
245 
246  // process the event
247  if ( L1MuGMTConfig::Debug(2) ) edm::LogVerbatim("GMT_info");
248  if ( L1MuGMTConfig::Debug(2) ) edm::LogVerbatim("GMT_info") << "**** L1GlobalMuonTrigger processing ****";
249  if ( L1MuGMTConfig::Debug(2) ) edm::LogVerbatim("GMT_info");
250 
251  int bx_min = L1MuGMTConfig::getBxMin();
252  int bx_max = L1MuGMTConfig::getBxMax();
253 
254  m_ExtendedCands.clear();
255 
256  // clear readout ring buffer
257  std::vector<L1MuGMTReadoutRecord*>::iterator irr = m_ReadoutRingbuffer.begin();
258  for ( ;irr!=m_ReadoutRingbuffer.end(); irr++) delete (*irr);
259  m_ReadoutRingbuffer.clear();
260 
261  if(m_db) m_db->reset(); // reset debug block
262 
263  for ( int bx = bx_min; bx <= bx_max; bx++ ) {
264  m_db->SetBX(bx);
265 
266  // create new element in readout ring buffer
267  m_ReadoutRingbuffer.push_back( new L1MuGMTReadoutRecord(bx) );
268 
269  if ( L1MuGMTConfig::Debug(2) ) edm::LogVerbatim("GMT_info") << "L1GlobalMuonTrigger processing bunch-crossing : " << bx;
270 
271  // get data into the data buffer
272  if ( L1MuGMTConfig::Debug(2) ) edm::LogVerbatim("GMT_info") << "running GMT PSB";
273  if ( m_PSB ) {
274  m_PSB->receiveData(e,bx);
275  if ( L1MuGMTConfig::Debug(4) ) m_PSB->print();
276  }
277 
278  if ( m_PSB && !m_PSB->empty() ) {
279 
280  // run matcher
281  if ( L1MuGMTConfig::Debug(2) ) edm::LogVerbatim("GMT_info") << "running GMT barrel Matcher";
282  if ( m_Matcher[0] ) m_Matcher[0]->run();
283  if ( L1MuGMTConfig::Debug(3) && m_Matcher[0] ) m_Matcher[0]->print();
284  if ( L1MuGMTConfig::Debug(2) ) edm::LogVerbatim("GMT_info") << "running GMT endcap Matcher";
285  if ( m_Matcher[1] ) m_Matcher[1]->run();
286  if ( L1MuGMTConfig::Debug(3) && m_Matcher[1] ) m_Matcher[1]->print();
287 
288  // run cancel-out units
289  if ( L1MuGMTConfig::Debug(2) ) edm::LogVerbatim("GMT_info") << "running GMT barrel Cancel Out Unit";
290  if ( m_CancelOutUnit[0] ) m_CancelOutUnit[0]->run();
292 
293  if ( L1MuGMTConfig::Debug(2) ) edm::LogVerbatim("GMT_info") << "running GMT endcap Cancel Out Unit";
294  if ( m_CancelOutUnit[1] ) m_CancelOutUnit[1]->run();
296 
297  if ( L1MuGMTConfig::Debug(2) ) edm::LogVerbatim("GMT_info") << "running GMT CSC/fRPC Cancel Out Unit";
298  if ( m_CancelOutUnit[2] ) m_CancelOutUnit[2]->run();
300 
301  if ( L1MuGMTConfig::Debug(2) ) edm::LogVerbatim("GMT_info") << "running GMT DT/bRPC Cancel Out Unit";
302  if ( m_CancelOutUnit[3] ) m_CancelOutUnit[3]->run();
304 
305  // run MIP & ISO bit assignment units
306  if ( L1MuGMTConfig::Debug(2) ) edm::LogVerbatim("GMT_info") << "running GMT barrel MIP & ISO bit Assignment Unit";
307  if ( m_MipIsoAU[0] ) m_MipIsoAU[0]->run();
308  if ( L1MuGMTConfig::Debug(3) && m_MipIsoAU[0] ) m_MipIsoAU[0]->print();
309  if ( L1MuGMTConfig::Debug(2) ) edm::LogVerbatim("GMT_info") << "running GMT endcap MIP & ISO bit Assignment Unit";
310  if ( m_MipIsoAU[1] ) m_MipIsoAU[1]->run();
311  if ( L1MuGMTConfig::Debug(3) && m_MipIsoAU[1] ) m_MipIsoAU[1]->print();
312 
313  // run Merger
314  if ( L1MuGMTConfig::Debug(2) ) edm::LogVerbatim("GMT_info") << "running GMT barrel Merger";
315  if ( m_Merger[0] ) m_Merger[0]->run();
316  if ( L1MuGMTConfig::Debug(3) && m_Merger[0] ) m_Merger[0]->print();
317  if ( L1MuGMTConfig::Debug(2) ) edm::LogVerbatim("GMT_info") << "running GMT endcap Merger";
318  if ( m_Merger[1] ) m_Merger[1]->run();
319  if ( L1MuGMTConfig::Debug(3) && m_Merger[1] ) m_Merger[1]->print();
320 
321  // run sorter
322  if ( L1MuGMTConfig::Debug(2) ) edm::LogVerbatim("GMT_info") << "running GMT Sorter";
323  if ( m_Sorter ) m_Sorter->run();
325 
326  // store found track candidates in a container
327  if ( m_Sorter->numberOfCands() > 0 ) {
328  const std::vector<const L1MuGMTExtendedCand*>& gmt_cont = m_Sorter->Cands();
329  std::vector<const L1MuGMTExtendedCand*>::const_iterator iexc;
330  for ( iexc = gmt_cont.begin(); iexc != gmt_cont.end(); iexc++ ) {
331  if ( *iexc ) m_ExtendedCands.push_back( **iexc );
332  }
333  }
334 
335  // reset GMT
336  reset();
337 
338  }
339  }
340 
341  // produce the output
342  std::auto_ptr<std::vector<L1MuGMTCand> > GMTCands(new std::vector<L1MuGMTCand>);
343  std::vector<L1MuGMTExtendedCand>::const_iterator iexc;
344  for(iexc=m_ExtendedCands.begin(); iexc!=m_ExtendedCands.end(); iexc++) {
345  GMTCands->push_back(*iexc);
346  }
347  e.put(GMTCands);
348 
349  std::auto_ptr<L1MuGMTReadoutCollection> GMTRRC(getReadoutCollection());
350  e.put(GMTRRC);
351 
352  if( m_sendMipIso ) {
353  std::auto_ptr<std::vector<unsigned> > mipiso(new std::vector<unsigned>);
354  for(int i=0; i<32; i++) {
355  mipiso->push_back(m_db->IsMIPISO(0,i));
356  }
357  e.put(mipiso);
358  }
359 
360 // delete registers and LUTs
362 }
363 
364 //
365 // reset GMT
366 //
368 
369  if ( m_PSB ) m_PSB->reset();
370  if ( m_Matcher[0] ) m_Matcher[0]->reset();
371  if ( m_Matcher[1] ) m_Matcher[1]->reset();
372  if ( m_CancelOutUnit[0] ) m_CancelOutUnit[0]->reset();
373  if ( m_CancelOutUnit[1] ) m_CancelOutUnit[1]->reset();
374  if ( m_CancelOutUnit[2] ) m_CancelOutUnit[2]->reset();
375  if ( m_CancelOutUnit[3] ) m_CancelOutUnit[3]->reset();
376  if ( m_MipIsoAU[0] ) m_MipIsoAU[0]->reset();
377  if ( m_MipIsoAU[1] ) m_MipIsoAU[1]->reset();
378  if ( m_Merger[0] ) m_Merger[0]->reset();
379  if ( m_Merger[1] ) m_Merger[1]->reset();
380  if ( m_Sorter) m_Sorter->reset();
381 
382 }
383 
384 // get the GMT readout data for the triggered bx
385 std::auto_ptr<L1MuGMTReadoutCollection> L1MuGlobalMuonTrigger::getReadoutCollection() {
386 
387  int bx_min_ro = L1MuGMTConfig::getBxMinRo();
388  int bx_max_ro = L1MuGMTConfig::getBxMaxRo();
389  int bx_size = bx_max_ro - bx_min_ro + 1;
390 
391  std::auto_ptr<L1MuGMTReadoutCollection> rrc(new L1MuGMTReadoutCollection(bx_size));
392 
393  for (int bx = bx_min_ro; bx <= bx_max_ro; bx++) {
394  std::vector<L1MuGMTReadoutRecord*>::const_iterator iter = m_ReadoutRingbuffer.begin();
395 
396  for ( ;iter!=m_ReadoutRingbuffer.end(); iter++) {
397 
398  if ( (*iter)->getBxInEvent() == bx ) {
399  rrc->addRecord(**iter);
400  break;
401  }
402  }
403 
404  }
405 
406  return rrc;
407 }
408 
409 // static data members
410 
void setGMTChanMask(const L1MuGMTChannelMask *gmtchanmask)
T getUntrackedParameter(std::string const &, T const &) const
int i
Definition: DBlmapReader.cc:9
unsigned long long m_L1MuTriggerScalesCacheID
std::vector< L1MuGMTReadoutRecord * > m_ReadoutRingbuffer
unsigned long long m_L1MuGMTChannelMaskCacheID
void print() const
print PSB
Definition: L1MuGMTPSB.cc:219
void SetBX(int bx)
Set the current bunch crossing.
void setTriggerPtScale(const L1MuTriggerPtScale *trigptscale)
void dumpRegs(std::string dir)
std::vector< L1MuGMTExtendedCand > m_ExtendedCands
void reset()
clear cancel-out unit
void print() const
print results after MIP &amp; ISO bit assignment
void print()
print cancel-out bits
bool empty() const
are there any data in the PSB
Definition: L1MuGMTPSB.cc:311
virtual void reset()
reset Sorter
unsigned getSubsystemMask() const
void reset()
clear MIP &amp; ISO bit assignment unit
L1MuGlobalMuonTrigger(const edm::ParameterSet &)
void reset()
clear PSB
Definition: L1MuGMTPSB.cc:198
void run()
run Matcher
void print()
print matching results
const std::vector< const L1MuGMTExtendedCand * > & Cands() const
return std::vector with all muon candidates
Definition: L1MuGMTSorter.h:64
unsigned long long m_L1CaloGeometryCacheID
virtual void run()
run Sorter
void run()
run GMT MIP &amp; ISO bit assignment unit
std::auto_ptr< L1MuGMTReadoutCollection > getReadoutCollection()
void setGMTScales(const L1MuGMTScales *gmtscales)
unsigned long long m_L1MuTriggerPtScaleCacheID
static bool Debug()
L1MuGMTMipIsoAU * m_MipIsoAU[2]
OrphanHandle< PROD > put(std::auto_ptr< PROD > product)
Put a new product.
Definition: Event.h:121
ConsumesCollector consumesCollector()
Use a ConsumesCollector to gather consumes information from helper functions.
static const L1MuGMTChannelMask * getGMTChanMask()
unsigned IsMIPISO(int bx, int idx)
Get stored MIP/ISO select bits.
static int getBxMaxRo()
L1MuGMTMatcher * m_Matcher[2]
void dumpLUTs(std::string dir)
static int getBxMinRo()
L1MuGMTMerger * m_Merger[2]
void setCaloGeom(const L1CaloGeometry *caloGeom)
void print() const
print results after selection
void reset()
clear Merger
unsigned long long m_L1MuGMTScalesCacheID
static int getBxMin()
void print()
print results after sorting
L1MuGMTCancelOutUnit * m_CancelOutUnit[4]
const T & get() const
Definition: EventSetup.h:56
void reset()
clear Matcher
T const * product() const
Definition: ESHandle.h:86
void run()
run GMT Merger
void setGMTParams(const L1MuGMTParameters *gmtparams)
void clearLUTsRegs()
def mkdir
Definition: eostools.py:250
void receiveData(edm::Event &e, int bx)
receive muon candidates
Definition: L1MuGMTPSB.cc:90
void reset()
Reset the debug block.
static L1MuGMTConfig * m_config
unsigned long long m_L1MuGMTParametersCacheID
void createLUTsRegs()
static int getBxMax()
static L1MuGMTDebugBlock * m_db
dbl *** dir
Definition: mlp_gen.cc:35
virtual void produce(edm::Event &, const edm::EventSetup &)
void setTriggerScales(const L1MuTriggerScales *trigscales)
void run()
run cancel-out unit
int numberOfCands() const
return number of found muon candidates after sorter
Definition: L1MuGMTSorter.h:61