259 unsigned int nDataMuons = 0;
260 unsigned int nEmulMuons = 0;
261 int dataMuonArray[8][10], emuMuonArray[8][10];
264 for(
int par=0; par<3; par++)
266 dataMuonArray[
muon][par]=0;
267 emuMuonArray[
muon][par] =0;
269 emuMuonArray[
muon][3] =-1;
270 dataMuonArray[
muon][3]=-1;
272 emuMuonArray[
muon][4]=7;
273 dataMuonArray[
muon][4]=7;
275 for(
int par2=5; par2<10; par2++)
277 emuMuonArray[
muon][par2]= -1;
278 dataMuonArray[
muon][par2]= -1;
291 <<
"\n No valid [data tracks] product found: " 292 <<
" L1CSCTrackCollection" 298 for(L1CSCTrackCollection::const_iterator trk=tracks.
product()->begin(); trk!=tracks.
product()->end(); trk++)
302 if ( (trk->first.BX() <2) && (trk->first.BX() > -1) )
306 dataMuonArray[nDataMuons][0] = trk->first.ptLUTAddress();
307 dataMuonArray[nDataMuons][1] = trk->first.sector();
308 dataMuonArray[nDataMuons][2] = trk->first.endcap();
309 dataMuonArray[nDataMuons][8] = trk->first.outputLink();
310 dataMuonArray[nDataMuons][4] = trk->first.BX();
311 dataMuonArray[nDataMuons][5] = trk->first.rank();
312 dataMuonArray[nDataMuons][6] = trk->first.localPhi();
313 dataMuonArray[nDataMuons][7] = trk->first.eta_packed();
314 dataMuonArray[nDataMuons][9] = trk->first.modeExtended();
329 <<
"\n No valid [emulator tracks] product found: " 330 <<
" L1CSCTrackCollection" 335 for(L1CSCTrackCollection::const_iterator trk=tracks.
product()->begin(); trk!=tracks.
product()->end(); trk++)
339 if((trk->first.BX() <2) && (trk->first.BX() >-1))
343 emuMuonArray[nEmulMuons][0] = trk->first.ptLUTAddress();
344 emuMuonArray[nEmulMuons][1] = trk->first.sector();
345 emuMuonArray[nEmulMuons][2] = trk->first.endcap();
346 emuMuonArray[nEmulMuons][4] = trk->first.BX();
347 emuMuonArray[nEmulMuons][5] = trk->first.rank();
348 emuMuonArray[nEmulMuons][6] = trk->first.localPhi();
349 emuMuonArray[nEmulMuons][7] = trk->first.eta_packed();
350 emuMuonArray[nEmulMuons][9] = trk->first.modeExtended();
356 if( (nDataMuons!=0)||(nEmulMuons!=0) ) {
362 if(nDataMuons==nEmulMuons)
365 for(
unsigned int mu1=0; mu1<nDataMuons; mu1++)
367 for(
unsigned int mu2=0; mu2<nEmulMuons; mu2++)
368 if((emuMuonArray[mu2][1]==dataMuonArray[mu1][1])&&(emuMuonArray[mu2][2]==dataMuonArray[mu1][2]))
370 if(emuMuonArray[mu2][0]==dataMuonArray[mu1][0])
372 emuMuonArray[mu2][3]=mu1;
373 dataMuonArray[mu1][3]=1;
378 for(
unsigned int c2a=0; c2a<nEmulMuons; c2a++)
380 if(emuMuonArray[c2a][3]==-1)
382 for(
unsigned int cor_a=0; cor_a<nDataMuons; cor_a++)
384 if( (dataMuonArray[cor_a][1]==emuMuonArray[c2a][1]) && (dataMuonArray[cor_a][2]==emuMuonArray[c2a][2]))
386 emuMuonArray[c2a][3]=cor_a;
387 dataMuonArray[cor_a][3]=1;
393 bool multiMap =
false;
396 for(
unsigned int c1a=0; c1a<(nEmulMuons-1); c1a++)
398 for(
unsigned int c1b=(c1a+1); c1b<nEmulMuons; c1b++)
400 if(emuMuonArray[c1a][3]==emuMuonArray[c1b][3])
412 for(
unsigned int mu3=0; mu3<nEmulMuons; mu3++)
414 int mapping = emuMuonArray[mu3][3];
415 if((mapping!=-1)&&(multiMap==
false))
418 int emuPhi12 = (0x0000ff & emuMuonArray[mu3][0]);
419 int datPhi12 = (0x0000ff & dataMuonArray[
mapping][0]);
420 int emuPhi23 = (0x000f00 & emuMuonArray[mu3][0])>>8;
421 int datPhi23 = (0x000f00 & dataMuonArray[
mapping][0])>>8;
422 int emuEta = (0x00f000 & emuMuonArray[mu3][0])>>12;
423 int datEta = (0x00f000 & dataMuonArray[
mapping][0])>>12;
426 int emuFrSin = (0xf00000 & emuMuonArray[mu3][0])>>20;
427 int datFrSin = (0xf00000 & dataMuonArray[
mapping][0])>>20;
429 int emuQual = emuMuonArray[mu3][5]>>5;
430 int datQual = dataMuonArray[
mapping][5]>>5;
431 int emuPt = 0x1f & emuMuonArray[mu3][5];
432 int datPt = 0x1f & dataMuonArray[
mapping][5];
433 int emuModeExtended = emuMuonArray[mu3][9];
434 int datModeExtended = dataMuonArray[
mapping][9];
440 if(emuModeExtended==datModeExtended)
448 if(dataMuonArray[mapping][8]==1)
453 phiComp->
Fill(dataMuonArray[mapping][6],emuMuonArray[mu3][6]);
454 etaComp->
Fill(dataMuonArray[mapping][7],emuMuonArray[mu3][7]);
470 for(
int dJ=0; dJ<7; dJ++)
472 for(
int dK=0; dK<15; dK++)
474 eDtStub[dJ][dK] = -55;
475 dDtStub[dJ][dK] = -55;
476 dDtStub[7][dK] = -55;
488 <<
"\n No valid [Data Stubs] product found: " 489 <<
" L1CSCTrackCollection" 496 vector<csctf::TrackStub> stuList = stub_list.
get();
497 vector<csctf::TrackStub>::const_iterator stu= stuList.begin();
498 for(; stu!=stuList.end(); stu++)
502 if((stu->BX()>4) && (stu->BX()<9))
504 dDtStub[0][dDtCounter] = stu->phiPacked();
505 dDtStub[1][dDtCounter] = stu->getQuality();
506 dDtStub[2][dDtCounter] = stu->endcap();
507 dDtStub[3][dDtCounter] = stu->sector();
508 dDtStub[4][dDtCounter] = stu->subsector();
523 <<
"\n No valid [Data Stubs] product found: " 524 <<
" L1CSCTrackCollection" 529 vector<csctf::TrackStub> emuList = emulStub.
get();
530 vector<csctf::TrackStub>::const_iterator eStu=emuList.begin();
531 for(; eStu!=emuList.end(); eStu++)
535 if((eStu->BX()>4) && (eStu->BX()<9) )
537 eDtStub[0][eDtCounter] = eStu->phiPacked();
538 eDtStub[1][eDtCounter] = eStu->getQuality();
539 eDtStub[2][eDtCounter] = eStu->endcap();
540 eDtStub[3][eDtCounter] = eStu->sector();
541 eDtStub[4][eDtCounter] = eStu->subsector();
549 for(
int eS=0; eS<eDtCounter; eS++)
552 for(
int dS=0; dS<dDtCounter; dS++)
555 if(eDtStub[2][eS]==dDtStub[2][dS])
558 if(eDtStub[3][eS]==dDtStub[3][dS])
561 if(eDtStub[4][eS]==dDtStub[4][dS])
564 if( (eDtStub[0][eS]==dDtStub[0][dS]) && (eDtStub[1][eS]==dDtStub[1][dS]) && (eDtStub[6][eS]!=1) && (dDtStub[6][dS]!=1) )
579 for(
int eS2=0; eS2<eDtCounter; eS2++)
581 for(
int dS2=0; dS2<dDtCounter; dS2++)
584 if( (eDtStub[2][eS2]==dDtStub[2][dS2]) && (eDtStub[3][eS2]==dDtStub[3][dS2]) && (eDtStub[4][eS2]==dDtStub[4][dS2]) )
587 if( ((dDtStub[7][eS2]==-55) || (dDtStub[7][dS2]>(
abs(eDtStub[0][eS2]-dDtStub[0][dS2]))) ) && (eDtStub[6][eS2]!=1) && (dDtStub[6][dS2]!=1) )
590 dDtStub[5][dS2] = eS2;
592 eDtStub[5][eS2] = dS2;
594 dDtStub[7][dS2] =
abs(eDtStub[0][eS2]-dDtStub[0][dS2]);
601 bool dtSMulti =
false;
604 for(
int eS3a=0; eS3a<eDtCounter-1; eS3a++)
605 for(
int eS3b=eS3a+1; eS3b<eDtCounter; eS3b++)
607 if( eDtStub[5][eS3a]==eDtStub[5][eS3b] ) dtSMulti=
true;
608 if( eDtStub[5][eS3a]==-55 || eDtStub[5][eS3b]==-55 ) dtUnmap++;
612 for(
int dS3a=0; dS3a<dDtCounter-1; dS3a++)
613 for(
int dS3b=dS3a+1; dS3b<dDtCounter; dS3b++)
615 if( dDtStub[5][dS3a]==dDtStub[5][dS3b] ) dtSMulti=
true;
616 if( dDtStub[5][dS3a]==-55||dDtStub[5][dS3b]==-55 ) dtUnmap++;
623 if(dtSMulti==
false && dtUnmap==0)
625 for(
int phil=0; phil<eDtCounter; phil++)
627 if(eDtStub[6][phil]==1 || eDtStub[6][phil]==2)
629 int indexFil = eDtStub[3][phil]*2+eDtStub[4][phil]-1;
630 dtStubPhi->
Fill(eDtStub[0][phil], dDtStub[0][ eDtStub[5][phil] ]);
632 if( eDtStub[0][phil] != dDtStub[0][ eDtStub[5][phil] ])
MonitorElement * etaComp_1d
MonitorElement * qualComp
std::vector< T > get() const
std::unique_ptr< CSCTFDTReceiver > my_dtrc
edm::EDGetTokenT< L1CSCTrackCollection > emulTrackProducer
edm::EDGetTokenT< L1CSCTrackCollection > dataTrackProducer
MonitorElement * qualComp_1d
MonitorElement * dtStubPhi
MonitorElement * phiComp_1d
edm::EDGetTokenT< L1MuDTChambPhContainer > emulStubProducer
MonitorElement * ptComp_1d
Abs< T >::type abs(const T &t)
void push_many(const std::vector< T > &data)
MonitorElement * pt1Comp_1d
edm::EDGetTokenT< CSCTriggerContainer< csctf::TrackStub > > dataStubProducer
MonitorElement * occComp_1d
MonitorElement * dtStubPhi_1d
T const * product() const
MonitorElement * pt5Comp_1d
MonitorElement * badDtStubSector
MonitorElement * pt4Comp_1d
MonitorElement * pt2Comp_1d
bool isUninitialized() const
MonitorElement * pt3Comp_1d