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
edm::ESGetToken< L1MuDTTFMasks, L1MuDTTFMasksRcd >	m_msksToken
edm::ESGetToken< L1MuDTTFParameters, L1MuDTTFParametersRcd >	m_parsToken
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 56 of file L1MuDTSectorReceiver.h.

Constructor & Destructor Documentation

L1MuDTSectorReceiver()

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())),
      m_parsToken(iC.esConsumes()),
      m_msksToken(iC.esConsumes()) {}

~L1MuDTSectorReceiver()

L1MuDTSectorReceiver::~L1MuDTSectorReceiver ( )

virtual

destructor

Definition at line 66 of file L1MuDTSectorReceiver.cc.

                                            {
  //  reset();
}

Member Function Documentation

address2sector()

int L1MuDTSectorReceiver::address2sector ( int  adr ) const

private

find the right sector for a given address

Definition at line 306 of file L1MuDTSectorReceiver.cc.

References L1MuDTSectorProcessor::id(), m_sp, hgcalTBTopologyTester_cfi::sector, 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;
}

address2wheel()

int L1MuDTSectorReceiver::address2wheel ( int  adr ) const

private

find the right wheel for a given address

Definition at line 320 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;
}

receiveCSCData()

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 220 of file L1MuDTSectorReceiver.cc.

References L1MuDTDataBuffer::addTSphi(), simKBmtfDigis_cfi::bx, L1MuDTTrackFinder::config(), L1MuDTSectorProcessor::data(), MillePedeFileConverter_cfg::e, L1MuDTTFMasks::get_inrec_chdis_csc(), L1MuDTTFParameters::get_soc_csc_etacanc(), L1MuDTTFConfig::getCSCTrSInputTag(), L1MuDTTFConfig::getEtaCanc(), L1MuDTSectorProcessor::id(), label, m_CSCTrSToken, m_sp, msks, pars, phi, hgcalTBTopologyTester_cfi::sector, L1MuDTSecProcId::sector(), relativeConstraints::station, makeGlobalPositionRcd_cfg::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;
  }
}

receiveDTBXData()

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 98 of file L1MuDTSectorReceiver.cc.

References funct::abs(), address2sector(), address2wheel(), L1MuDTDataBuffer::addTSphi(), simKBmtfDigis_cfi::bx, L1MuDTChambPhContainer::chPhiSegm1(), L1MuDTChambPhContainer::chPhiSegm2(), L1MuDTChambPhDigi::code(), L1MuDTTrackFinder::config(), L1MuDTSectorProcessor::data(), MillePedeFileConverter_cfg::e, 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(), L1MuDTTFConfig::getNbitsExtPhi(), L1MuDTTFConfig::getTSOutOfTimeFilter(), L1MuDTTFConfig::getTSOutOfTimeWindow(), L1MuDTSectorProcessor::id(), m_DTDigiToken, m_sp, msks, L1MuDTSectorProcessor::ovl(), pars, phi, L1MuDTChambPhDigi::phi(), L1MuDTChambPhDigi::phiB(), hgcalTBTopologyTester_cfi::sector, L1MuDTSecProcId::sector(), relativeConstraints::station, makeGlobalPositionRcd_cfg::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);
        }
      }
    }
  }
}

reset()

void L1MuDTSectorReceiver::reset

(

void

)

clear Sector Receiver

Definition at line 93 of file L1MuDTSectorReceiver.cc.

Referenced by L1MuDTSectorProcessor::reset().

{}

run()

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 77 of file L1MuDTSectorReceiver.cc.

References simKBmtfDigis_cfi::bx, HltBtagPostValidation_cff::c, L1MuDTTrackFinder::config(), MillePedeFileConverter_cfg::e, m_msksToken, m_parsToken, m_sp, msks, L1MuDTTFConfig::overlap(), L1MuDTSectorProcessor::ovl(), pars, receiveCSCData(), receiveDTBXData(), and L1MuDTSectorProcessor::tf().

Referenced by L1MuDTSectorProcessor::run().

                                                                                {
  pars = c.getHandle(m_parsToken);
  msks = c.getHandle(m_msksToken);

  // 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

m_CSCTrSToken

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

private

Definition at line 86 of file L1MuDTSectorReceiver.h.

Referenced by receiveCSCData().

m_DTDigiToken

edm::EDGetTokenT<L1MuDTChambPhContainer> L1MuDTSectorReceiver::m_DTDigiToken

private

Definition at line 85 of file L1MuDTSectorReceiver.h.

Referenced by receiveDTBXData().

m_msksToken

edm::ESGetToken<L1MuDTTFMasks, L1MuDTTFMasksRcd> L1MuDTSectorReceiver::m_msksToken

private

Definition at line 89 of file L1MuDTSectorReceiver.h.

Referenced by run().

m_parsToken

edm::ESGetToken<L1MuDTTFParameters, L1MuDTTFParametersRcd> L1MuDTSectorReceiver::m_parsToken

private

Definition at line 88 of file L1MuDTSectorReceiver.h.

Referenced by run().

m_sp

L1MuDTSectorProcessor& L1MuDTSectorReceiver::m_sp

private

Definition at line 84 of file L1MuDTSectorReceiver.h.

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

msks

edm::ESHandle<L1MuDTTFMasks> L1MuDTSectorReceiver::msks

private

Definition at line 91 of file L1MuDTSectorReceiver.h.

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

pars

edm::ESHandle<L1MuDTTFParameters> L1MuDTSectorReceiver::pars

private

Definition at line 90 of file L1MuDTSectorReceiver.h.

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