#include <L1MuBMEtaProcessor.h>

Public Member Functions
int	eta (int id) const
	return eta values, index [0,11] More...

bool	fineBit (int id) const
	return fine bit, index [0,11] More...

int	id () const
	return Eta Processor identifier (0-11) More...

	L1MuBMEtaProcessor (const L1MuBMTrackFinder &, int id, edm::ConsumesCollector &&iC)
	constructor More...

void	print () const
	print muon candidates found by the Eta Processor More...

virtual void	reset ()
	reset the Eta Processor More...

virtual void	run (int bx, const edm::Event &e, const edm::EventSetup &c)
	run the Eta Processor More...

const L1MuBMTrackFinder &	tf () const
	return reference to barrel MTTF More...

virtual	~L1MuBMEtaProcessor ()
	destructor More...

Private Member Functions
void	assign ()
	assign eta and etaFineBit More...

void	receiveAddresses ()
	receive addresses (from 6 Sector Processors) More...

void	receiveData (int bx, const edm::Event &e, const L1TMuonBarrelParams &params)
	receive data (eta trigger primitives) More...

void	runEtaMatchingUnit (const L1TMuonBarrelParams &params)
	run Eta Matching Unit (EMU) More...

void	runEtaTrackFinder (const L1TMuonBarrelParams &params)
	run Eta Track Finder (ETF) More...

Static Private Member Functions
static int	quality (int id, int stat)
	get quality code; id [0,26], stat [1,3] More...

Private Attributes
int	m_address [12]

edm::ESGetToken < L1TMuonBarrelParams, L1TMuonBarrelParamsRcd >	m_bmtfParamsToken

edm::EDGetTokenT < L1MuDTChambThContainer >	m_DTDigiToken

int	m_epid

int	m_eta [12]

bool	m_fine [12]

std::vector< int >	m_foundPattern

int	m_mask

int	m_pattern [12]

const L1MuBMTrackFinder &	m_tf

L1MuBMTrack *	m_TrackCand [12]

L1MuBMTrack *	m_TracKCand [12]

std::vector< const L1MuBMTrackSegEta * >	m_tseta

L1MuDTTFMasks	msks

L1MuBMTEtaPatternLut	theEtaPatternLUT

L1MuBMTQualPatternLut	theQualPatternLUT

Detailed Description

Eta Processor:

An Eta Processor consists of :

one receiver unit,
one Eta Track Finder (ETF) and
one Eta Matching Unit (EMU)

N. Neumeister CERN EP J. Troconiz UAM Madrid

Definition at line 63 of file L1MuBMEtaProcessor.h.

Constructor & Destructor Documentation

L1MuBMEtaProcessor::L1MuBMEtaProcessor	(	const L1MuBMTrackFinder &	tf,
		int	id,
		edm::ConsumesCollector &&	iC
	)

constructor

Definition at line 64 of file L1MuBMEtaProcessor.cc.

References m_tseta.

     : m_tf(tf),
       m_epid(id),
       m_foundPattern(0),
       m_tseta(15),
       m_bmtfParamsToken(iC.esConsumes()),
       m_DTDigiToken(iC.consumes<L1MuDTChambThContainer>(L1MuBMTFConfig::getBMThetaDigiInputTag())) {
   m_tseta.reserve(15);
 }

L1MuBMEtaProcessor::~L1MuBMEtaProcessor ( )

virtual

destructor

Definition at line 78 of file L1MuBMEtaProcessor.cc.

78 {}

Member Function Documentation

void L1MuBMEtaProcessor::assign ( )

private

assign eta and etaFineBit

Definition at line 418 of file L1MuBMEtaProcessor.cc.

References L1MuBMTEtaPatternLut::getPattern(), mps_fire::i, m_eta, m_fine, m_pattern, m_TrackCand, m_TracKCand, m_tseta, AlCaHLTBitMon_ParallelJobs::p, L1MuDTEtaPattern::position(), L1MuBMTrack::setEta(), L1MuBMTrack::setFineEtaBit(), L1MuBMTrack::setTSeta(), edm_modernize_messagelogger::stat, theEtaPatternLUT, and L1MuDTEtaPattern::wheel().

Referenced by run().

                                 {
   for (int i = 0; i < 12; i++) {
     if (m_TrackCand[i]) {
       if (m_eta[i] != 99) {
         m_TrackCand[i]->setEta(m_eta[i]);
         m_TracKCand[i]->setEta(m_eta[i]);
       } else {
         //        if ( i/2 != 2 ) cerr << "L1MuBMEtaProcessor: assign invalid eta" << " " << m_address[i] << endl;
       }
       if (m_fine[i]) {
         m_TrackCand[i]->setFineEtaBit();
         m_TracKCand[i]->setFineEtaBit();
         // find all contributing track segments
         const L1MuDTEtaPattern p = theEtaPatternLUT.getPattern(m_pattern[i]);
         vector<const L1MuBMTrackSegEta*> TSeta;
         const L1MuBMTrackSegEta* ts = nullptr;
         for (int stat = 0; stat < 3; stat++) {
           int wh = p.wheel(stat + 1);
           int pos = p.position(stat + 1);
           if (pos == 0)
             continue;
           ts = m_tseta[wh + 2 + stat * 5];
           TSeta.push_back(ts);
         }
         m_TrackCand[i]->setTSeta(TSeta);
         m_TracKCand[i]->setTSeta(TSeta);
       }
     }
   }
 }

int L1MuBMEtaProcessor::eta ( int id ) const

inline

return eta values, index [0,11]

Definition at line 87 of file L1MuBMEtaProcessor.h.

References id(), and m_eta.

Referenced by Particle.Particle::__str__(), Jet.Jet::jetID(), and Jet.Jet::puJetId().

87 { return m_eta[id]; }

L1MuBMEtaProcessor::m_eta

int m_eta[12]

Definition: L1MuBMEtaProcessor.h:117

L1MuBMEtaProcessor::id

int id() const

return Eta Processor identifier (0-11)

Definition: L1MuBMEtaProcessor.h:72

bool L1MuBMEtaProcessor::fineBit ( int id ) const

inline

return fine bit, index [0,11]

Definition at line 90 of file L1MuBMEtaProcessor.h.

References id(), and m_fine.

90 { return m_fine[id]; }

L1MuBMEtaProcessor::m_fine

bool m_fine[12]

Definition: L1MuBMEtaProcessor.h:118

L1MuBMEtaProcessor::id

int id() const

return Eta Processor identifier (0-11)

Definition: L1MuBMEtaProcessor.h:72

int L1MuBMEtaProcessor::id ( void ) const

inline

return Eta Processor identifier (0-11)

Definition at line 72 of file L1MuBMEtaProcessor.h.

References m_epid.

Referenced by eta(), fineBit(), and quality().

72 { return m_epid; }

L1MuBMEtaProcessor::m_epid

int m_epid

Definition: L1MuBMEtaProcessor.h:113

void L1MuBMEtaProcessor::print ( void ) const

print muon candidates found by the Eta Processor

Definition at line 132 of file L1MuBMEtaProcessor.cc.

References gather_cfg::cout, L1MuDTEtaPattern::eta(), L1MuBMTQualPatternLut::getCoarseEta(), L1MuBMTEtaPatternLut::getPattern(), mps_fire::i, L1MuDTEtaPattern::id(), dqmiolumiharvest::j, m_address, m_epid, m_eta, m_fine, m_foundPattern, m_pattern, m_tseta, AlCaHLTBitMon_ParallelJobs::p, position, L1MuDTEtaPattern::quality(), quality(), relativeConstraints::station, theEtaPatternLUT, and theQualPatternLUT.

                                      {
   bool empty1 = true;
   for (int i = 0; i < 15; i++) {
     empty1 &= (m_tseta[i] == nullptr || m_tseta[i]->empty());
   }
 
   bool empty2 = true;
   for (int i = 0; i < 12; i++) {
     empty2 &= (m_address[i] == 0);
   }
 
   if (!empty1 || !empty2) {
     cout << "Eta processor " << m_epid << " : " << endl;
 
     // print local pattern
     if (!empty1) {
       cout << "Local pattern : " << endl;
       for (int i = 0; i < 15; i++) {
         if ((i + 5) % 5 == 0)
           cout << "station " << m_tseta[i]->station() << " : ";
         bitset<7> pos(m_tseta[i]->position());
         bitset<7> qua(m_tseta[i]->quality());
         for (int j = 6; j >= 0; j--) {
           cout << pos[j] + qua[j];
         }
         cout << " ";
         if ((i + 1) % 5 == 0)
           cout << endl;
       }
       cout << "Found patterns :" << endl;
       vector<int>::const_iterator iter;
       for (iter = m_foundPattern.begin(); iter != m_foundPattern.end(); iter++) {
         const L1MuDTEtaPattern p = theEtaPatternLUT.getPattern(*iter);
         int qualitycode = p.quality();
         cout << "ID = " << setw(4) << p.id() << "  "
              << "eta = " << setw(3) << p.eta() << "  "
              << "quality = " << setw(2) << qualitycode << " (" << quality(qualitycode, 1) << " "
              << quality(qualitycode, 2) << " " << quality(qualitycode, 3) << ")";
         for (int i = 0; i < 12; i++) {
           if (m_pattern[i] == p.id())
             cout << " <--";
         }
         cout << endl;
       }
     }
 
     cout << "Received addresses : " << endl;
     for (int i = 0; i < 12; i++)
       cout << setw(3) << m_address[i] << " ";
     cout << endl;
 
     if (!empty1) {
       cout << "Matched patterns : " << endl;
       for (int i = 0; i < 12; i++) {
         if (m_fine[i]) {
           const L1MuDTEtaPattern p = theEtaPatternLUT.getPattern(m_pattern[i]);
           int fineeta = p.eta();
           int coarseeta = theQualPatternLUT.getCoarseEta(i / 2 + 1, m_address[i]);
           cout << "Index = " << setw(2) << i << ", "
                << "address = " << setw(2) << m_address[i] << " --> "
                << "pattern = " << setw(4) << m_pattern[i] << " "
                << "eta (coarse) = " << setw(3) << coarseeta << " "
                << "eta (fine) = " << setw(3) << fineeta << " "
                << "quality = " << setw(2) << p.quality() << endl;
         }
       }
     }
 
     cout << "Eta values and fine bits : " << endl;
     for (int i = 0; i < 12; i++)
       cout << setw(3) << m_eta[i] << " ";
     cout << endl;
     for (int i = 0; i < 12; i++)
       cout << setw(3) << m_fine[i] << " ";
     cout << endl;
   }
 }

int L1MuBMEtaProcessor::quality	(	int	id,
		int	stat
	)

staticprivate

get quality code; id [0,26], stat [1,3]

Definition at line 452 of file L1MuBMEtaProcessor.cc.

References id().

Referenced by print(), and runEtaTrackFinder().

                                                 {
   // quality codes as defined in CMS Note 2001/027
   // This QualityPatterns are used to have a defined Quality-Identifier for
   // all possible found tracks.
   // Therefore three integer numbers ( Station 1, 2, 3 from left to right )
   // can have numbers like 0, 1 or 2
   //    0 ... no hit is given
   //    1 ... a LTRG is given
   //    2 ... a HTRG is given
 
   const int qualcode[27][3] = {{0, 0, 0}, {1, 0, 0}, {0, 1, 0}, {0, 0, 1}, {2, 0, 0}, {0, 2, 0}, {0, 0, 2},
                                {1, 1, 0}, {1, 0, 1}, {0, 1, 1}, {2, 1, 0}, {1, 2, 0}, {2, 0, 1}, {1, 0, 2},
                                {0, 2, 1}, {0, 1, 2}, {2, 2, 0}, {2, 0, 2}, {0, 2, 2}, {1, 1, 1}, {2, 1, 1},
                                {1, 2, 1}, {1, 1, 2}, {2, 2, 1}, {2, 1, 2}, {1, 2, 2}, {2, 2, 2}};
 
   return qualcode[id][stat - 1];
 }

void L1MuBMEtaProcessor::receiveAddresses ( )

private

receive addresses (from 6 Sector Processors)

Definition at line 279 of file L1MuBMEtaProcessor.cc.

References L1MuBMTrack::address(), L1MuBMTrack::empty(), mps_fire::i, m_address, m_epid, m_tf, m_TrackCand, m_TracKCand, contentValuesFiles::number, L1MuBMTrackFinder::sp(), L1MuBMSectorProcessor::track(), L1MuBMSectorProcessor::tracK(), and L1MuBMAddressArray::trackAddressCode().

Referenced by run().

                                           {
   // get track address code of all track segments
   // 6*2 times 5 bits; valid range [1-22]
 
   int sector = m_epid;
 
   int i = 0;
   for (int wheel = -3; wheel <= 3; wheel++) {
     if (wheel == 0)
       continue;
     L1MuBMSecProcId tmpspid(wheel, sector);
     for (int number = 0; number < 2; number++) {
       const L1MuBMTrack* cand = m_tf.sp(tmpspid)->track(number);
       const L1MuBMTrack* canD = m_tf.sp(tmpspid)->tracK(number);
       if (cand) {
         m_address[i] = cand->address().trackAddressCode();
         if (!cand->empty()) {
           m_TrackCand[i] = const_cast<L1MuBMTrack*>(cand);
           m_TracKCand[i] = const_cast<L1MuBMTrack*>(canD);
         }
         i++;
       }
     }
   }
 }

void L1MuBMEtaProcessor::receiveData	(	int	bx,
		const edm::Event &	e,
		const L1TMuonBarrelParams &	params
	)

private

receive data (eta trigger primitives)

Definition at line 213 of file L1MuBMEtaProcessor.cc.

References L1TMuonBarrelParams::LUTParams::eta_lut_, L1MuDTTFMasks::get_etsoc_chdis_st1(), L1MuDTTFMasks::get_etsoc_chdis_st2(), L1MuDTTFMasks::get_etsoc_chdis_st3(), edm::Event::getByToken(), mps_fire::i, L1TMuonBarrelParams::l1mudttfmasks, L1TMuonBarrelParams::lutparams_, m_DTDigiToken, m_epid, L1MuBMTEtaPatternLut::m_lut, L1MuBMTQualPatternLut::m_lut, m_mask, m_tseta, msks, L1MuDTChambThDigi::position(), L1TMuonBarrelParams::LUTParams::qp_lut_, L1MuDTChambThDigi::quality(), edm_modernize_messagelogger::stat, theEtaPatternLUT, and theQualPatternLUT.

Referenced by run().

                                                                                                      {
   msks = bmtfParams.l1mudttfmasks;
   theEtaPatternLUT.m_lut = bmtfParams.lutparams_.eta_lut_;  //l1mudttfetaplut;  // ETF look-up table
   theQualPatternLUT.m_lut = bmtfParams.lutparams_.qp_lut_;  //l1mudttfqualplut; // EMU look-up tables
 
   edm::Handle<L1MuDTChambThContainer> dttrig;
   e.getByToken(m_DTDigiToken, dttrig);
 
   // const int bx_offset = dttrig->correctBX();
   int bx_offset = 0;
   bx = bx + bx_offset;
 
   //
   // get 5*3 eta track segments
   //
   int sector = m_epid;
   for (int stat = 1; stat <= 3; stat++) {
     for (int wheel = -2; wheel <= 2; wheel++) {
       L1MuDTChambThDigi const* tseta = dttrig->chThetaSegm(wheel, stat, sector, bx);
       bitset<7> pos;
       bitset<7> qual;
 
       int lwheel = wheel + 1;
       if (wheel < 0)
         lwheel = wheel - 1;
 
       bool masked = false;
       if (stat == 1)
         masked = msks.get_etsoc_chdis_st1(lwheel, sector);
       if (stat == 2)
         masked = msks.get_etsoc_chdis_st2(lwheel, sector);
       if (stat == 3)
         masked = msks.get_etsoc_chdis_st3(lwheel, sector);
 
       if (!masked)
         m_mask = false;
 
       if (tseta && !masked) {
         if (wheel == -2 || wheel == -1 ||
             (wheel == 0 && (sector == 0 || sector == 3 || sector == 4 || sector == 7 || sector == 8 || sector == 11))) {
           for (int i = 0; i < 7; i++) {
             if (tseta->position(i))
               pos.set(6 - i);
             if (tseta->quality(i))
               qual.set(6 - i);
           }
         } else {
           for (int i = 0; i < 7; i++) {
             if (tseta->position(i))
               pos.set(i);
             if (tseta->quality(i))
               qual.set(i);
           }
         }
       }
 
       const L1MuBMTrackSegEta* tmpts =
           new L1MuBMTrackSegEta(wheel, sector, stat, pos.to_ulong(), qual.to_ulong(), bx - bx_offset);
       m_tseta.push_back(tmpts);
     }
   }
 }

void L1MuBMEtaProcessor::reset ( void )

virtual

reset the Eta Processor

Definition at line 103 of file L1MuBMEtaProcessor.cc.

References mps_fire::i, m_address, m_eta, m_fine, m_foundPattern, m_mask, m_pattern, m_TrackCand, m_TracKCand, and m_tseta.

                                {
   vector<const L1MuBMTrackSegEta*>::iterator iter = m_tseta.begin();
   while (iter != m_tseta.end()) {
     if (*iter) {
       delete *iter;
       *iter = nullptr;
     }
     iter++;
   }
 
   m_tseta.clear();
 
   for (int i = 0; i < 12; i++) {
     m_eta[i] = 99;
     m_fine[i] = false;
     m_pattern[i] = 0;
     m_address[i] = 0;
     m_TrackCand[i] = nullptr;
     m_TracKCand[i] = nullptr;
   }
 
   m_foundPattern.clear();
 
   m_mask = true;
 }

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

virtual

run the Eta Processor

Definition at line 87 of file L1MuBMEtaProcessor.cc.

References assign(), edm::EventSetup::getData(), L1MuBMTFConfig::getEtaTF(), m_bmtfParamsToken, submitPVValidationJobs::params, receiveAddresses(), receiveData(), runEtaMatchingUnit(), and runEtaTrackFinder().

                                                                               {
   auto const& params = c.getData(m_bmtfParamsToken);
   if (L1MuBMTFConfig::getEtaTF()) {
     receiveData(bx, e, params);
     runEtaTrackFinder(params);
   }
 
   receiveAddresses();
   runEtaMatchingUnit(params);
 
   assign();
 }

void L1MuBMEtaProcessor::runEtaMatchingUnit ( const L1TMuonBarrelParams & params )

private

run Eta Matching Unit (EMU)

Definition at line 354 of file L1MuBMEtaProcessor.cc.

References L1MuDTEtaPattern::eta(), spr::find(), L1MuBMTQualPatternLut::getCoarseEta(), L1MuBMTEtaPatternLut::getPattern(), L1MuBMTQualPatternLut::getQualifiedPatterns(), mps_fire::i, L1TMuonBarrelParams::lutparams_, m_address, m_eta, m_fine, m_foundPattern, L1MuBMTQualPatternLut::m_lut, m_mask, m_pattern, AlCaHLTBitMon_ParallelJobs::p, L1TMuonBarrelParams::LUTParams::qp_lut_, theEtaPatternLUT, and theQualPatternLUT.

Referenced by run().

                                                                                  {
   theQualPatternLUT.m_lut = bmtfParams.lutparams_.qp_lut_;  //l1mudttfqualplut; // EMU look-up tables
 
   // loop over all addresses
   for (int i = 0; i < 12; i++) {
     int adr = m_address[i];
     if (adr == 0)
       continue;
     int sp = i / 2 + 1;  //sector processor [1,6]
 
     // assign coarse eta value
     if (!m_mask)
       m_eta[i] = theQualPatternLUT.getCoarseEta(sp, adr);
 
     if (m_foundPattern.empty())
       continue;
 
     // get list of qualified patterns ordered by quality
     // and compare with found patterns
     const vector<short>& qualifiedPatterns = theQualPatternLUT.getQualifiedPatterns(sp, adr);
     vector<short>::const_iterator iter;
     vector<int>::const_iterator f_iter;
     for (iter = qualifiedPatterns.begin(); iter != qualifiedPatterns.end(); iter++) {
       f_iter = find(m_foundPattern.begin(), m_foundPattern.end(), (*iter));
       // found
       if (f_iter != m_foundPattern.end()) {
         const L1MuDTEtaPattern p = theEtaPatternLUT.getPattern(*f_iter);
         // assign fine eta value
         m_fine[i] = true;
         m_eta[i] = p.eta();  // improved eta
         ;
         m_pattern[i] = (*f_iter);
         break;
       }
     }
   }
 
   // if both tracks from one sector processor deliver the same track address
   // both tracks get only a coarse eta value!
 
   // loop over sector processors
   for (int i = 0; i < 6; i++) {
     int idx1 = 2 * i;
     int idx2 = 2 * i + 1;
     int adr1 = m_address[idx1];
     int adr2 = m_address[idx2];
     if (adr1 == 0 || adr2 == 0)
       continue;
     if (adr1 == adr2 && !m_mask) {
       // both tracks get coarse (default) eta value
       m_eta[idx1] = theQualPatternLUT.getCoarseEta(i + 1, adr1);
       m_pattern[idx1] = 0;
       m_fine[idx1] = false;
       m_eta[idx2] = theQualPatternLUT.getCoarseEta(i + 1, adr2);
 
       m_pattern[idx2] = 0;
       m_fine[idx2] = false;
     }
   }
 }

void L1MuBMEtaProcessor::runEtaTrackFinder ( const L1TMuonBarrelParams & params )

private

run Eta Track Finder (ETF)

Definition at line 308 of file L1MuBMEtaProcessor.cc.

References L1MuBMTEtaPatternLut::begin(), newFWLiteAna::bin, relativeConstraints::empty, L1MuBMTEtaPatternLut::end(), L1TMuonBarrelParams::LUTParams::eta_lut_, good, mps_fire::i, L1MuDTEtaPattern::id(), L1TMuonBarrelParams::lutparams_, m_foundPattern, L1MuBMTEtaPatternLut::m_lut, m_tseta, L1MuDTEtaPattern::position(), submitPVResolutionJobs::q, L1MuDTEtaPattern::quality(), quality(), relativeConstraints::station, theEtaPatternLUT, and L1MuDTEtaPattern::wheel().

Referenced by run().

                                                                                 {
   theEtaPatternLUT.m_lut = bmtfParams.lutparams_.eta_lut_;  //l1mudttfetaplut;  // ETF look-up table
 
   // check if there are any data
   bool empty = true;
   for (int i = 0; i < 15; i++) {
     empty &= m_tseta[i]->empty();
   }
   if (empty)
     return;
 
   // Pattern comparator:
   // loop over all patterns and compare with local chamber pattern
   // result : list of valid pattern IDs ( m_foundPattern )
   L1MuBMTEtaPatternLut::ETFLut_iter it = theEtaPatternLUT.begin();
   while (it != theEtaPatternLUT.end()) {
     const L1MuDTEtaPattern pattern = (*it).second;
     int qualitycode = pattern.quality();
 
     bool good = true;
 
     for (int station = 0; station < 3; station++) {
       int q = quality(qualitycode, station + 1);
       int wheel = pattern.wheel(station + 1);
       int bin = pattern.position(station + 1);
       if (bin == 0)
         continue;
       bitset<7> pos = m_tseta[wheel + 2 + station * 5]->position();
       bitset<7> qual = m_tseta[wheel + 2 + station * 5]->quality();
       // compare position
       good &= pos.test(bin - 1);
       // compare quality
       if (q == 2)
         good &= qual.test(bin - 1);
     }
 
     if (good)
       m_foundPattern.push_back(pattern.id());
 
     it++;
   }
 }