#include <L1MuDTSectorReceiver.h>

Public Member Functions
	L1MuDTSectorReceiver (L1MuDTSectorProcessor &, edm::ConsumesCollector &&iC)
	constructor More...

void	reset ()
	clear Sector Receiver More...

void	run (int bx, const edm::Event &e, const edm::EventSetup &c)
	receive track segment data from the DTBX and CSC chamber triggers More...

virtual	~L1MuDTSectorReceiver ()
	destructor More...

Private Member Functions
int	address2sector (int adr) const
	find the right sector for a given address More...

int	address2wheel (int adr) const
	find the right wheel for a given address More...

void	receiveCSCData (int bx, const edm::Event &e, const edm::EventSetup &c)
	receive track segment data from CSC chamber trigger More...

void	receiveDTBXData (int bx, const edm::Event &e, const edm::EventSetup &c)
	receive track segment data from DTBX chamber trigger More...

Private Attributes
edm::EDGetTokenT< CSCTriggerContainer< csctf::TrackStub > >	m_CSCTrSToken

edm::EDGetTokenT< L1MuDTChambPhContainer >	m_DTDigiToken

L1MuDTSectorProcessor &	m_sp

edm::ESHandle< L1MuDTTFMasks >	msks

edm::ESHandle< L1MuDTTFParameters >	pars

Detailed Description

Sector Receiver:

The Sector Receiver receives track segment data from the DTBX and CSC chamber triggers and stores it into the Data Buffer

N. Neumeister CERN EP J. Troconiz UAM Madrid

Definition at line 53 of file L1MuDTSectorReceiver.h.

Constructor & Destructor Documentation

L1MuDTSectorReceiver::L1MuDTSectorReceiver	(	L1MuDTSectorProcessor &	sp,
		edm::ConsumesCollector &&	iC
	)

constructor

Definition at line 56 of file L1MuDTSectorReceiver.cc.

                                                                                                 : 
         m_sp(sp),
         m_DTDigiToken(iC.consumes<L1MuDTChambPhContainer>(m_sp.tf().config()->getDTDigiInputTag())),
         m_CSCTrSToken(iC.mayConsume<CSCTriggerContainer<csctf::TrackStub> >(m_sp.tf().config()->getCSCTrSInputTag())) {
 }

L1MuDTSectorReceiver::~L1MuDTSectorReceiver ( )

virtual

destructor

Definition at line 66 of file L1MuDTSectorReceiver.cc.

                                             { 
 
 //  reset();
   
 }

Member Function Documentation

int L1MuDTSectorReceiver::address2sector ( int adr ) const

private

find the right sector for a given address

Definition at line 294 of file L1MuDTSectorReceiver.cc.

References L1MuDTSectorProcessor::id(), m_sp, and L1MuDTSecProcId::sector().

Referenced by receiveDTBXData().

                                                       {
 
   int sector = m_sp.id().sector();
 
   if ( adr >= 4 && adr <= 7  ) sector = (sector+13)%12; // +1
   if ( adr >= 8 && adr <= 11 ) sector = (sector+11)%12; // -1
 
   return sector;
 
 }

int L1MuDTSectorReceiver::address2wheel ( int adr ) const

private

find the right wheel for a given address

Definition at line 309 of file L1MuDTSectorReceiver.cc.

References L1MuDTSectorProcessor::id(), L1MuDTSecProcId::locwheel(), m_sp, L1MuDTSecProcId::wheel(), and makeMuonMisalignmentScenario::wheel.

Referenced by receiveDTBXData().

                                                      {
 
   int wheel = m_sp.id().locwheel();
 
   // for 2, 3, 6, 7, 10, 11
   if ( (adr/2)%2 == 1 ) wheel = m_sp.id().wheel();
 
   return wheel;
 
 }

void L1MuDTSectorReceiver::receiveCSCData	(	int	bx,
		const edm::Event &	e,
		const edm::EventSetup &	c
	)

private

receive track segment data from CSC chamber trigger

Definition at line 223 of file L1MuDTSectorReceiver.cc.

References L1MuDTDataBuffer::addTSphi(), L1MuDTTrackFinder::config(), L1MuDTSectorProcessor::data(), L1MuDTTFMasks::get_inrec_chdis_csc(), L1MuDTTFParameters::get_soc_csc_etacanc(), edm::Event::getByToken(), L1MuDTTFConfig::getCSCTrSInputTag(), L1MuDTTFConfig::getEtaCanc(), L1MuDTSectorProcessor::id(), diffTwoXMLs::label, m_CSCTrSToken, m_sp, msks, pars, phi, L1MuDTSecProcId::sector(), relativeConstraints::station, GlobalPosition_Frontier_DevDB_cff::tag, L1MuDTSectorProcessor::tf(), L1MuDTSecProcId::wheel(), and makeMuonMisalignmentScenario::wheel.

Referenced by run().

                                                                                            {
   
   if ( (m_sp.tf().config()->getCSCTrSInputTag()).label() == "none" ) return;
 
   if ( bx < -6 || bx > 6 ) return;
 
   edm::Handle<CSCTriggerContainer<csctf::TrackStub> > csctrig;
   e.getByToken(m_CSCTrSToken,csctrig);
 
   const int bxCSC = 6;
   
   vector<csctf::TrackStub> csc_list;
   vector<csctf::TrackStub>::const_iterator csc_iter;  
   
   int station = 1; // only ME13
   int wheel = m_sp.id().wheel();
   int side = ( wheel == 3 ) ? 1 : 2;
   int sector = m_sp.id().sector();
   int csc_sector = ( sector == 0 ) ? 6 : (sector+1)/2;
   int subsector = ( sector%2 == 0 ) ? 2 : 1;
 
   csc_list = csctrig->get(side,station,csc_sector,subsector,bxCSC+bx);
   int ncsc = 0;
   for ( csc_iter = csc_list.begin(); csc_iter != csc_list.end(); csc_iter++ ) {
     bool etaFlag = ( csc_iter->etaPacked() > 17 ); 
     int qualCSC = csc_iter->getQuality();
            
     // convert CSC quality code to DTBX quality code
     unsigned int qual = 7;
     if ( qualCSC ==  2 ) qual = 0;
     if ( qualCSC ==  6 ) qual = 1;
     if ( qualCSC ==  7 ) qual = 2;
     if ( qualCSC ==  8 ) qual = 2;
     if ( qualCSC ==  9 ) qual = 3;
     if ( qualCSC == 10 ) qual = 3;
     if ( qualCSC == 11 ) qual = 4;
     if ( qualCSC == 12 ) qual = 5;
     if ( qualCSC == 13 ) qual = 5;
     if ( qualCSC == 14 ) qual = 6;
     if ( qualCSC == 15 ) qual = 6;
     if ( qual == 7 ) continue;
 
     // convert CSC phi to DTBX phi
     int phi = csc_iter->phiPacked();
     if ( phi > 2047 ) phi -= 4096; 
     if ( phi < -2048 || phi > 2047 ) continue; 
 
     if ( msks->get_inrec_chdis_csc(m_sp.id().wheel(), m_sp.id().sector()) ) continue;
     if ( qual < pars->get_soc_qual_csc(m_sp.id().wheel(), m_sp.id().sector()) ) continue;
     if ( pars->get_soc_csc_etacanc(m_sp.id().wheel(), m_sp.id().sector()) && etaFlag ) continue;
     if ( m_sp.tf().config()->getEtaCanc() && etaFlag ) continue;
 
     if ( ncsc < 2 ) {
       int address = 16 + ncsc;
       bool tag = (ncsc == 1 ) ? true : false;
       L1MuDTTrackSegPhi tmpts(wheel,sector,station+2,phi,0,
                               static_cast<L1MuDTTrackSegPhi::TSQuality>(qual),
                               tag,bx,etaFlag);
       m_sp.data()->addTSphi(address,tmpts);
       ncsc++;
     }
     //    else cout << "too many CSC track segments!" << endl;
   }  
 
 }

void L1MuDTSectorReceiver::receiveDTBXData	(	int	bx,
		const edm::Event &	e,
		const edm::EventSetup &	c
	)

private

receive track segment data from DTBX chamber trigger

Definition at line 107 of file L1MuDTSectorReceiver.cc.

References funct::abs(), address2sector(), address2wheel(), L1MuDTDataBuffer::addTSphi(), L1MuDTChambPhContainer::chPhiSegm1(), L1MuDTChambPhContainer::chPhiSegm2(), L1MuDTChambPhDigi::code(), L1MuDTTrackFinder::config(), L1MuDTSectorProcessor::data(), L1MuDTTFMasks::get_inrec_chdis_st1(), L1MuDTTFMasks::get_inrec_chdis_st2(), L1MuDTTFMasks::get_inrec_chdis_st3(), L1MuDTTFMasks::get_inrec_chdis_st4(), L1MuDTTFParameters::get_soc_nbx_del(), edm::Event::getByToken(), L1MuDTTFConfig::getNbitsExtPhi(), L1MuDTTFConfig::getTSOutOfTimeFilter(), L1MuDTTFConfig::getTSOutOfTimeWindow(), L1MuDTSectorProcessor::id(), m_DTDigiToken, m_sp, msks, L1MuDTSectorProcessor::ovl(), pars, phi, L1MuDTChambPhDigi::phi(), L1MuDTChambPhDigi::phiB(), L1MuDTSecProcId::sector(), relativeConstraints::station, GlobalPosition_Frontier_DevDB_cff::tag, L1MuDTSectorProcessor::tf(), RecoTauCleanerPlugins::tolerance, L1MuDTSecProcId::wheel(), and makeMuonMisalignmentScenario::wheel.

Referenced by run().

                                                                                             {
 
   edm::Handle<L1MuDTChambPhContainer> dttrig;
   e.getByToken(m_DTDigiToken,dttrig);
 
   L1MuDTChambPhDigi const* ts=nullptr;
 
   // const int bx_offset = dttrig->correctBX();
   int bx_offset=0;
   bx = bx + bx_offset;
 
   // get DTBX phi track segments  
   int address = 0;
   for ( int station = 1; station <= 4; station++ ) {
     int max_address = (station == 1) ? 2 : 12;
     for (int reladr =0; reladr < max_address; reladr++) {
       address++;
       if ( m_sp.ovl() && (reladr/2)%2 != 0 ) continue;
       int wheel  = address2wheel(reladr);
       int sector = address2sector(reladr);     
       if ( reladr%2 == 0 ) ts = dttrig->chPhiSegm1(wheel,station,sector,bx);
       if ( reladr%2 == 1 ) ts = dttrig->chPhiSegm2(wheel,station,sector,bx);
       if ( ts ) {
         int phi  = ts->phi();
         int phib = ts->phiB();
         int qual = ts->code();
         bool tag = (reladr%2 == 1) ? true : false;
 
         int lwheel = m_sp.id().wheel();
         lwheel = abs(lwheel)/lwheel*(abs(wheel)+1);
 
         if ( station == 1 ) {
           if ( msks->get_inrec_chdis_st1(lwheel, sector) ) continue;
           if ( qual < pars->get_inrec_qual_st1(lwheel, sector) ) continue;
         } 
         else if ( station == 2 ) {
           if ( msks->get_inrec_chdis_st2(lwheel, sector) ) continue;
           if ( qual < pars->get_inrec_qual_st2(lwheel, sector) ) continue;
           } 
         else if ( station == 3 ) {
           if ( msks->get_inrec_chdis_st3(lwheel, sector) ) continue;
           if ( qual < pars->get_inrec_qual_st3(lwheel, sector) ) continue;
         } 
         else if ( station == 4 ) {
           if ( msks->get_inrec_chdis_st4(lwheel, sector) ) continue;
           if ( qual < pars->get_inrec_qual_st4(lwheel, sector) ) continue;
         } 
 
         if ( reladr/2 == 1 && qual < pars->get_soc_stdis_n(m_sp.id().wheel(), m_sp.id().sector())  ) continue;
         if ( reladr/2 == 2 && qual < pars->get_soc_stdis_wl(m_sp.id().wheel(), m_sp.id().sector()) ) continue;
         if ( reladr/2 == 3 && qual < pars->get_soc_stdis_zl(m_sp.id().wheel(), m_sp.id().sector()) ) continue;
         if ( reladr/2 == 4 && qual < pars->get_soc_stdis_wr(m_sp.id().wheel(), m_sp.id().sector()) ) continue;
         if ( reladr/2 == 5 && qual < pars->get_soc_stdis_zr(m_sp.id().wheel(), m_sp.id().sector()) ) continue;
 
         //
         // out-of-time TS filter (compare TS at +-1 bx)
         // 
         bool skipTS = false;
 
         bool nbx_del = pars->get_soc_nbx_del(m_sp.id().wheel(), m_sp.id().sector());
         if ( m_sp.tf().config()->getTSOutOfTimeFilter() || nbx_del ) {
  
           int sh_phi = 12 - m_sp.tf().config()->getNbitsExtPhi();
           int tolerance = m_sp.tf().config()->getTSOutOfTimeWindow();
 
           L1MuDTChambPhDigi const * tsPreviousBX_1 = dttrig->chPhiSegm1(wheel,station,sector,bx-1);
           if ( tsPreviousBX_1 ) {
             int phiBX  = tsPreviousBX_1->phi();
             int qualBX = tsPreviousBX_1->code();
             if ( abs( (phi >> sh_phi) - (phiBX >> sh_phi) ) <= tolerance &&
                  qualBX > qual ) skipTS = true;
           }
           
           L1MuDTChambPhDigi const * tsPreviousBX_2 = dttrig->chPhiSegm2(wheel,station,sector,bx-1);
           if ( tsPreviousBX_2 ) {
             int phiBX  = tsPreviousBX_2->phi();
             int qualBX = tsPreviousBX_2->code();
             if ( abs( (phi >> sh_phi) - (phiBX >> sh_phi) ) <= tolerance &&
                  qualBX > qual ) skipTS = true;
           }
      
           L1MuDTChambPhDigi const * tsNextBX_1 = dttrig->chPhiSegm1(wheel,station,sector,bx+1);
           if ( tsNextBX_1 ) {
             int phiBX  = tsNextBX_1->phi();
             int qualBX = tsNextBX_1->code();
             if ( abs( (phi >> sh_phi) - (phiBX >> sh_phi) ) <= tolerance &&
                  qualBX > qual ) skipTS = true;
           }
 
           L1MuDTChambPhDigi const * tsNextBX_2 = dttrig->chPhiSegm2(wheel,station,sector,bx+1);
           if ( tsNextBX_2 ) {
             int phiBX  = tsNextBX_2->phi();
             int qualBX = tsNextBX_2->code();
             if ( abs( (phi >> sh_phi) - (phiBX >> sh_phi) ) <= tolerance &&
                  qualBX > qual ) skipTS = true;
           }
 
         }
 
         if ( !skipTS ) { 
           L1MuDTTrackSegPhi tmpts(wheel,sector,station,phi,phib,
                                   static_cast<L1MuDTTrackSegPhi::TSQuality>(qual),
                                   tag,bx-bx_offset);
           m_sp.data()->addTSphi(address-1,tmpts);
         }
       }
 
     }
   }
 
 }

void L1MuDTSectorReceiver::reset ( void )

clear Sector Receiver

Definition at line 99 of file L1MuDTSectorReceiver.cc.

Referenced by L1MuDTSectorProcessor::reset().

                                  {
 
 }

void L1MuDTSectorReceiver::run	(	int	bx,
		const edm::Event &	e,
		const edm::EventSetup &	c
	)

receive track segment data from the DTBX and CSC chamber triggers

Definition at line 80 of file L1MuDTSectorReceiver.cc.

References L1MuDTTrackFinder::config(), edm::EventSetup::get(), m_sp, msks, L1MuDTTFConfig::overlap(), L1MuDTSectorProcessor::ovl(), pars, receiveCSCData(), receiveDTBXData(), and L1MuDTSectorProcessor::tf().

Referenced by L1MuDTSectorProcessor::run().

                                                                                 {
 
   c.get< L1MuDTTFParametersRcd >().get( pars );
   c.get< L1MuDTTFMasksRcd >().get( msks );
 
   // get track segments from DTBX chamber trigger
   receiveDTBXData(bx, e, c);
   
   // get track segments from CSC chamber trigger
   if ( m_sp.tf().config()->overlap() && m_sp.ovl() ) { 
     receiveCSCData(bx, e, c);
   }
 
 }