L1MuDTSectorReceiver Class Reference

#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 54 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 64 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 304 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 318 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 218 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(), 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 96 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(), 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 91 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 75 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);
  }
}

Member Data Documentation

edm::EDGetTokenT<CSCTriggerContainer<csctf::TrackStub> > L1MuDTSectorReceiver::m_CSCTrSToken

private

Definition at line 84 of file L1MuDTSectorReceiver.h.

Referenced by receiveCSCData().

edm::EDGetTokenT<L1MuDTChambPhContainer> L1MuDTSectorReceiver::m_DTDigiToken

private

Definition at line 83 of file L1MuDTSectorReceiver.h.

Referenced by receiveDTBXData().

L1MuDTSectorProcessor& L1MuDTSectorReceiver::m_sp

private

Definition at line 82 of file L1MuDTSectorReceiver.h.

Referenced by address2sector(), address2wheel(), receiveCSCData(), receiveDTBXData(), and run().

edm::ESHandle<L1MuDTTFMasks> L1MuDTSectorReceiver::msks

private

Definition at line 87 of file L1MuDTSectorReceiver.h.

Referenced by receiveCSCData(), receiveDTBXData(), and run().

edm::ESHandle<L1MuDTTFParameters> L1MuDTSectorReceiver::pars

private

Definition at line 86 of file L1MuDTSectorReceiver.h.

Referenced by receiveCSCData(), receiveDTBXData(), and run().