CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
List of all members | Public Member Functions | Private Member Functions | Private Attributes
l1t::MicroGMTCancelOutUnit Class Reference

#include <MicroGMTCancelOutUnit.h>

Public Member Functions

void initialise (L1TMuonGlobalParamsHelper *)
 Initialisation from ES record. More...
 
 MicroGMTCancelOutUnit ()
 
void setCancelOutBits (GMTInternalWedges &, tftype, cancelmode)
 Cancel out between sectors/wedges in one track finder. More...
 
void setCancelOutBitsOverlapBarrel (GMTInternalWedges &, GMTInternalWedges &, cancelmode)
 Cancel-out between overlap and barrel track finders. More...
 
void setCancelOutBitsOverlapEndcap (GMTInternalWedges &, GMTInternalWedges &, cancelmode)
 Cancel-out between overlap and endcap track finders. More...
 
virtual ~MicroGMTCancelOutUnit ()
 

Private Member Functions

void getCoordinateCancelBits (std::vector< std::shared_ptr< GMTInternalMuon >> &, std::vector< std::shared_ptr< GMTInternalMuon >> &)
 Compares all muons from coll1 with all muons from coll2 and sets the cancel-bits based on eta/phi coordinates. More...
 
void getTrackAddrCancelBits (std::vector< std::shared_ptr< GMTInternalMuon >> &, std::vector< std::shared_ptr< GMTInternalMuon >> &)
 Compares all muons from coll1 with all muons from coll2 and sets the cancel-bits based on track addresses. More...
 

Private Attributes

std::shared_ptr
< MicroGMTMatchQualLUT > 		m_boNegMatchQualLUT
 
std::shared_ptr
< MicroGMTMatchQualLUT > 		m_boPosMatchQualLUT
 
std::shared_ptr
< MicroGMTMatchQualLUT > 		m_brlSingleMatchQualLUT
 
std::shared_ptr
< MicroGMTMatchQualLUT > 		m_foNegMatchQualLUT
 
std::shared_ptr
< MicroGMTMatchQualLUT > 		m_foPosMatchQualLUT
 
std::shared_ptr
< MicroGMTMatchQualLUT > 		m_fwdNegSingleMatchQualLUT
 
std::shared_ptr
< MicroGMTMatchQualLUT > 		m_fwdPosSingleMatchQualLUT
 
std::map< int, std::shared_ptr
< MicroGMTMatchQualLUT > > 		m_lutDict
 
std::shared_ptr
< MicroGMTMatchQualLUT > 		m_ovlNegSingleMatchQualLUT
 
std::shared_ptr
< MicroGMTMatchQualLUT > 		m_ovlPosSingleMatchQualLUT
 

Detailed Description

Definition at line 15 of file MicroGMTCancelOutUnit.h.

Constructor & Destructor Documentation

l1t::MicroGMTCancelOutUnit::MicroGMTCancelOutUnit ( )

Definition at line 6 of file MicroGMTCancelOutUnit.cc.

7 {
8 }
l1t::MicroGMTCancelOutUnit::~MicroGMTCancelOutUnit ( )
virtual

Definition at line 10 of file MicroGMTCancelOutUnit.cc.

11 {
12 
13 }

Member Function Documentation

void l1t::MicroGMTCancelOutUnit::getCoordinateCancelBits ( std::vector< std::shared_ptr< GMTInternalMuon >> &  coll1,
std::vector< std::shared_ptr< GMTInternalMuon >> &  coll2 
)
private

Compares all muons from coll1 with all muons from coll2 and sets the cancel-bits based on eta/phi coordinates.

Definition at line 138 of file MicroGMTCancelOutUnit.cc.

References funct::abs(), l1t::bmtf, dPhi(), l1t::emtf_neg, l1t::emtf_pos, l1t::MicroGMTMatchQualLUT::getDeltaEtaWidth(), l1t::MicroGMTMatchQualLUT::getDeltaPhiWidth(), l1t::MicroGMTMatchQualLUT::lookup(), m_lutDict, and match().

Referenced by setCancelOutBits(), setCancelOutBitsOverlapBarrel(), and setCancelOutBitsOverlapEndcap().

139 {
140  if (coll1.size() == 0 || coll2.size() == 0) {
141  return;
142  }
143  tftype coll1TfType = (*coll1.begin())->trackFinderType();
144  tftype coll2TfType = (*coll2.begin())->trackFinderType();
145  if (coll2TfType != tftype::bmtf && coll1TfType % 2 != coll2TfType % 2) {
146  edm::LogError("Detector side mismatch") << "Overlap-Endcap cancel out between positive and negative detector side attempted. Check eta assignment. OMTF candidate: TF type: " << coll1TfType << ", hwEta: " << (*coll1.begin())->hwEta() << ". EMTF candidate: TF type: " << coll2TfType << ", hwEta: " << (*coll2.begin())->hwEta() << ". TF type even: pos. side; odd: neg. side." << std::endl;
147  return;
148  }
149 
150  MicroGMTMatchQualLUT* matchLUT = m_lutDict.at(coll1TfType+coll2TfType*10).get();
151 
152  for (auto mu_w1 = coll1.begin(); mu_w1 != coll1.end(); ++mu_w1) {
153  int etaFine1 = (*mu_w1)->hwHF();
154  // for EMTF muons set eta fine bit to true since hwHF is the halo bit
155  if (coll1TfType == tftype::emtf_pos || coll1TfType == tftype::emtf_neg) {
156  etaFine1 = 1;
157  }
158  for (auto mu_w2 = coll2.begin(); mu_w2 != coll2.end(); ++mu_w2) {
159  int etaFine2 = (*mu_w2)->hwHF();
160  // for EMTF muons set eta fine bit to true since hwHF is the halo bit
161  if (coll2TfType == tftype::emtf_pos || coll2TfType == tftype::emtf_neg) {
162  etaFine2 = 1;
163  }
164  // both muons must have the eta fine bit set in order to use the eta fine part of the LUT
165  int etaFine = (int)(etaFine1 > 0 && etaFine2 > 0);
166 
167  // The LUT for cancellation takes reduced width phi and eta, we need the LSBs
168  int dPhiMask = (1 << matchLUT->getDeltaPhiWidth()) - 1;
169  int dEtaMask = (1 << matchLUT->getDeltaEtaWidth()) - 1;
170 
171  int dPhi = (*mu_w1)->hwGlobalPhi() - (*mu_w2)->hwGlobalPhi();
172  dPhi = std::abs(dPhi);
173  if (dPhi > 338) dPhi -= 576; // shifts dPhi to [-pi, pi) in integer scale
174  dPhi = std::abs(dPhi);
175  int dEta = std::abs((*mu_w1)->hwEta() - (*mu_w2)->hwEta());
176  // check first if the delta is within the LSBs that the LUT takes, otherwise the distance
177  // is greater than what we want to cancel -> e.g. 31(int) is max => 31*0.01 = 0.31 (rad)
178  // LUT takes 5 LSB for dEta and 3 LSB for dPhi
179  if (dEta <= dEtaMask && dPhi <= dPhiMask) {
180  int match = matchLUT->lookup(etaFine, dEta & dEtaMask, dPhi & dPhiMask);
181  if (match == 1) {
182  if((*mu_w1)->hwQual() > (*mu_w2)->hwQual()) {
183  (*mu_w2)->setHwCancelBit(1);
184  } else {
185  (*mu_w1)->setHwCancelBit(1);
186  }
187  }
188  }
189  }
190  }
191 }
dPhi
double dPhi(double phi1, double phi2)
Definition: JetUtil.h:30
funct::abs
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
match
std::pair< typename Association::data_type::first_type, double > match(Reference key, Association association, bool bestMatchByMaxValue)
Generic matching function.
Definition: Utils.h:10
l1t::MicroGMTCancelOutUnit::m_lutDict
std::map< int, std::shared_ptr< MicroGMTMatchQualLUT > > m_lutDict
Definition: MicroGMTCancelOutUnit.h:43
void l1t::MicroGMTCancelOutUnit::getTrackAddrCancelBits ( std::vector< std::shared_ptr< GMTInternalMuon >> &  coll1,
std::vector< std::shared_ptr< GMTInternalMuon >> &  coll2 
)
private

Compares all muons from coll1 with all muons from coll2 and sets the cancel-bits based on track addresses.

Definition at line 194 of file MicroGMTCancelOutUnit.cc.

References l1t::bmtf, and i.

Referenced by setCancelOutBits(), setCancelOutBitsOverlapBarrel(), and setCancelOutBitsOverlapEndcap().

195 {
196  if (coll1.size() == 0 || coll2.size() == 0) {
197  return;
198  }
199  // Address based cancel out is implemented for BMTF only
200  if ((*coll1.begin())->trackFinderType() == tftype::bmtf && (*coll2.begin())->trackFinderType() == tftype::bmtf) {
201  for (auto mu_w1 = coll1.begin(); mu_w1 != coll1.end(); ++mu_w1) {
202  std::map<int, int> trkAddr_w1 = (*mu_w1)->origin().trackAddress();
203  int wheelNum_w1 = trkAddr_w1[l1t::RegionalMuonCand::bmtfAddress::kWheelNum];
204  int wheelSide_w1 = trkAddr_w1[l1t::RegionalMuonCand::bmtfAddress::kWheelSide];
205  std::vector<int> stations_w1;
206  stations_w1.push_back(trkAddr_w1[l1t::RegionalMuonCand::bmtfAddress::kStat1]);
207  stations_w1.push_back(trkAddr_w1[l1t::RegionalMuonCand::bmtfAddress::kStat2]);
208  stations_w1.push_back(trkAddr_w1[l1t::RegionalMuonCand::bmtfAddress::kStat3]);
209  stations_w1.push_back(trkAddr_w1[l1t::RegionalMuonCand::bmtfAddress::kStat4]);
210  //std::cout << "Track address 1: wheelSide (1 == negative side): " << wheelSide_w1 << ", wheelNum: " << wheelNum_w1 << ", stations1234: 0x" << hex << stations_w1[0] << stations_w1[1] << stations_w1[2] << stations_w1[3] << dec << std::endl;
211 
212  for (auto mu_w2 = coll2.begin(); mu_w2 != coll2.end(); ++mu_w2) {
213  std::map<int, int> trkAddr_w2 = (*mu_w2)->origin().trackAddress();
214  int wheelNum_w2 = trkAddr_w2[l1t::RegionalMuonCand::bmtfAddress::kWheelNum];
215  int wheelSide_w2 = trkAddr_w2[l1t::RegionalMuonCand::bmtfAddress::kWheelSide];
216  std::vector<int> stations_w2;
217  stations_w2.push_back(trkAddr_w2[l1t::RegionalMuonCand::bmtfAddress::kStat1]);
218  stations_w2.push_back(trkAddr_w2[l1t::RegionalMuonCand::bmtfAddress::kStat2]);
219  stations_w2.push_back(trkAddr_w2[l1t::RegionalMuonCand::bmtfAddress::kStat3]);
220  stations_w2.push_back(trkAddr_w2[l1t::RegionalMuonCand::bmtfAddress::kStat4]);
221  //std::cout << "Track address 2: wheelSide (1 == negative side): " << wheelSide_w2 << ", wheelNum: " << wheelNum_w2 << ", stations1234: 0x" << hex << stations_w2[0] << stations_w2[1] << stations_w2[2] << stations_w2[3] << dec << std::endl;
222 
223  int nMatchedStations = 0;
224  // search for duplicates in stations 2-4
225  for (int i = 1; i < 4; ++i) {
226  if (wheelSide_w1 == wheelSide_w2) { // both tracks are on the same detector side
227  if (wheelNum_w1 == wheelNum_w2) { // both tracks have the same reference wheel
228  if ((stations_w1[i] == 0x0 && stations_w2[i] == 0x2) ||
229  (stations_w1[i] == 0x1 && stations_w2[i] == 0x3) ||
230  (stations_w1[i] == 0x4 && stations_w2[i] == 0x0) ||
231  (stations_w1[i] == 0x5 && stations_w2[i] == 0x1) ||
232  (stations_w1[i] == 0x8 && stations_w2[i] == 0xA) ||
233  (stations_w1[i] == 0x9 && stations_w2[i] == 0xB) ||
234  (stations_w1[i] == 0xC && stations_w2[i] == 0x8) ||
235  (stations_w1[i] == 0xD && stations_w2[i] == 0x9))
236  {
237  ++nMatchedStations;
238  }
239  } else if (wheelNum_w1 == wheelNum_w2 - 1) { // track 2 is one wheel higher than track 1
240  if ((stations_w1[i] == 0x0 && stations_w2[i] == 0xA) ||
241  (stations_w1[i] == 0x1 && stations_w2[i] == 0xB) ||
242  (stations_w1[i] == 0x4 && stations_w2[i] == 0x8) ||
243  (stations_w1[i] == 0x5 && stations_w2[i] == 0x9))
244  {
245  ++nMatchedStations;
246  }
247  } else if (wheelNum_w1 == wheelNum_w2 + 1) { // track 2 is one wheel lower than track 1
248  if ((stations_w1[i] == 0x8 && stations_w2[i] == 0x2) ||
249  (stations_w1[i] == 0x9 && stations_w2[i] == 0x3) ||
250  (stations_w1[i] == 0xC && stations_w2[i] == 0x0) ||
251  (stations_w1[i] == 0xD && stations_w2[i] == 0x1))
252  {
253  ++nMatchedStations;
254  }
255  }
256  } else {
257  if (wheelNum_w1 == 0 && wheelNum_w2 == 0) { // both tracks are on either side of the central wheel (+0 and -0)
258  if ((stations_w1[i] == 0x8 && stations_w2[i] == 0xA) ||
259  (stations_w1[i] == 0x9 && stations_w2[i] == 0xB) ||
260  (stations_w1[i] == 0xC && stations_w2[i] == 0x8) ||
261  (stations_w1[i] == 0xD && stations_w2[i] == 0x9))
262  {
263  ++nMatchedStations;
264  }
265  }
266  }
267  }
268  //std::cout << "Shared hits found: " << nMatchedStations << std::endl;
269  if (nMatchedStations > 0) {
270  if ((*mu_w1)->origin().hwQual() >= (*mu_w2)->origin().hwQual()) {
271  (*mu_w2)->setHwCancelBit(1);
272  } else {
273  (*mu_w1)->setHwCancelBit(1);
274  }
275  }
276  }
277  }
278  } else {
279  edm::LogError("Cancel out not implemented") << "Address based cancel out is currently only implemented for the barrel track finder.";
280  }
281 }
i
int i
Definition: DBlmapReader.cc:9
void l1t::MicroGMTCancelOutUnit::initialise ( L1TMuonGlobalParamsHelper *  microGMTParamsHelper)

Initialisation from ES record.

Definition at line 16 of file MicroGMTCancelOutUnit.cc.

References l1t::bmtf, L1TMuonGlobalParamsHelper::bONegMatchQualLUT(), L1TMuonGlobalParamsHelper::bOPosMatchQualLUT(), l1t::MicroGMTMatchQualLUTFactory::create(), l1t::emtf_emtf_neg, l1t::emtf_emtf_pos, l1t::emtf_neg, l1t::emtf_pos, L1TMuonGlobalParamsHelper::fONegMatchQualLUT(), L1TMuonGlobalParamsHelper::fOPosMatchQualLUT(), L1TMuonGlobalParamsHelper::fwdNegSingleMatchQualLUT(), L1TMuonGlobalParamsHelper::fwdPosSingleMatchQualLUT(), stage2BMTFBufferRaw_cfi::fwVersion, L1TMuonGlobalParamsHelper::fwVersion(), m_boNegMatchQualLUT, m_boPosMatchQualLUT, m_foNegMatchQualLUT, m_foPosMatchQualLUT, m_fwdNegSingleMatchQualLUT, m_fwdPosSingleMatchQualLUT, m_lutDict, m_ovlNegSingleMatchQualLUT, m_ovlPosSingleMatchQualLUT, l1t::omtf_bmtf_neg, l1t::omtf_bmtf_pos, l1t::omtf_emtf_neg, l1t::omtf_emtf_pos, l1t::omtf_neg, l1t::omtf_omtf_neg, l1t::omtf_omtf_pos, l1t::omtf_pos, L1TMuonGlobalParamsHelper::ovlNegSingleMatchQualLUT(), and L1TMuonGlobalParamsHelper::ovlPosSingleMatchQualLUT().

Referenced by L1TMuonProducer::beginRun().

16  {
17  int fwVersion = microGMTParamsHelper->fwVersion();
18  m_boPosMatchQualLUT = l1t::MicroGMTMatchQualLUTFactory::create(microGMTParamsHelper->bOPosMatchQualLUT(), cancel_t::omtf_bmtf_pos, fwVersion);
19  m_boNegMatchQualLUT = l1t::MicroGMTMatchQualLUTFactory::create(microGMTParamsHelper->bONegMatchQualLUT(), cancel_t::omtf_bmtf_neg, fwVersion);
20  m_foPosMatchQualLUT = l1t::MicroGMTMatchQualLUTFactory::create(microGMTParamsHelper->fOPosMatchQualLUT(), cancel_t::omtf_emtf_pos, fwVersion);
21  m_foNegMatchQualLUT = l1t::MicroGMTMatchQualLUTFactory::create(microGMTParamsHelper->fONegMatchQualLUT(), cancel_t::omtf_emtf_neg, fwVersion);
22  m_ovlPosSingleMatchQualLUT = l1t::MicroGMTMatchQualLUTFactory::create(microGMTParamsHelper->ovlPosSingleMatchQualLUT(), cancel_t::omtf_omtf_pos, fwVersion);
23  m_ovlNegSingleMatchQualLUT = l1t::MicroGMTMatchQualLUTFactory::create(microGMTParamsHelper->ovlNegSingleMatchQualLUT(), cancel_t::omtf_omtf_neg, fwVersion);
24  m_fwdPosSingleMatchQualLUT = l1t::MicroGMTMatchQualLUTFactory::create(microGMTParamsHelper->fwdPosSingleMatchQualLUT(), cancel_t::emtf_emtf_pos, fwVersion);
25  m_fwdNegSingleMatchQualLUT = l1t::MicroGMTMatchQualLUTFactory::create(microGMTParamsHelper->fwdNegSingleMatchQualLUT(), cancel_t::emtf_emtf_neg, fwVersion);
26 
27  m_lutDict[tftype::omtf_neg+tftype::bmtf*10] = m_boNegMatchQualLUT;
28  m_lutDict[tftype::omtf_pos+tftype::bmtf*10] = m_boPosMatchQualLUT;
29  m_lutDict[tftype::omtf_pos+tftype::omtf_pos*10] = m_ovlPosSingleMatchQualLUT;
30  m_lutDict[tftype::omtf_neg+tftype::omtf_neg*10] = m_ovlNegSingleMatchQualLUT;
31  m_lutDict[tftype::emtf_pos+tftype::emtf_pos*10] = m_fwdPosSingleMatchQualLUT;
32  m_lutDict[tftype::emtf_neg+tftype::emtf_neg*10] = m_fwdNegSingleMatchQualLUT;
33  m_lutDict[tftype::omtf_pos+tftype::emtf_pos*10] = m_foPosMatchQualLUT;
34  m_lutDict[tftype::omtf_neg+tftype::emtf_neg*10] = m_foNegMatchQualLUT;
35 }
l1t::MicroGMTCancelOutUnit::m_ovlNegSingleMatchQualLUT
std::shared_ptr< MicroGMTMatchQualLUT > m_ovlNegSingleMatchQualLUT
Definition: MicroGMTCancelOutUnit.h:40
l1t::MicroGMTCancelOutUnit::m_foPosMatchQualLUT
std::shared_ptr< MicroGMTMatchQualLUT > m_foPosMatchQualLUT
Definition: MicroGMTCancelOutUnit.h:36
l1t::MicroGMTCancelOutUnit::m_fwdNegSingleMatchQualLUT
std::shared_ptr< MicroGMTMatchQualLUT > m_fwdNegSingleMatchQualLUT
Definition: MicroGMTCancelOutUnit.h:42
l1t::MicroGMTCancelOutUnit::m_fwdPosSingleMatchQualLUT
std::shared_ptr< MicroGMTMatchQualLUT > m_fwdPosSingleMatchQualLUT
Definition: MicroGMTCancelOutUnit.h:41
l1t::MicroGMTCancelOutUnit::m_boPosMatchQualLUT
std::shared_ptr< MicroGMTMatchQualLUT > m_boPosMatchQualLUT
Definition: MicroGMTCancelOutUnit.h:34
stage2BMTFBufferRaw_cfi.fwVersion
tuple fwVersion
Definition: stage2BMTFBufferRaw_cfi.py:62
l1t::MicroGMTCancelOutUnit::m_boNegMatchQualLUT
std::shared_ptr< MicroGMTMatchQualLUT > m_boNegMatchQualLUT
Definition: MicroGMTCancelOutUnit.h:35
l1t::MicroGMTCancelOutUnit::m_ovlPosSingleMatchQualLUT
std::shared_ptr< MicroGMTMatchQualLUT > m_ovlPosSingleMatchQualLUT
Definition: MicroGMTCancelOutUnit.h:39
l1t::MicroGMTCancelOutUnit::m_foNegMatchQualLUT
std::shared_ptr< MicroGMTMatchQualLUT > m_foNegMatchQualLUT
Definition: MicroGMTCancelOutUnit.h:37
l1t::MicroGMTMatchQualLUTFactory::create
static ReturnType create(const std::string &filename, const double maxDR, const double fEta, const double fEtaCoarse, const double fPhi, cancel_t cancelType, const int fwVersion)
Definition: MicroGMTMatchQualLUTFactory.cc:18
l1t::MicroGMTCancelOutUnit::m_lutDict
std::map< int, std::shared_ptr< MicroGMTMatchQualLUT > > m_lutDict
Definition: MicroGMTCancelOutUnit.h:43
void l1t::MicroGMTCancelOutUnit::setCancelOutBits ( GMTInternalWedges &  wedges,
tftype  trackFinder,
cancelmode  mode 
)

Cancel out between sectors/wedges in one track finder.

Definition at line 38 of file MicroGMTCancelOutUnit.cc.

References l1t::bmtf, l1t::coordinate, getCoordinateCancelBits(), getTrackAddrCancelBits(), and RPCpg::mu.

Referenced by L1TMuonProducer::produce().

39 {
40  std::vector<std::shared_ptr<GMTInternalMuon>> coll1;
41  coll1.reserve(3);
42  std::vector<std::shared_ptr<GMTInternalMuon>> coll2;
43  coll2.reserve(3);
44  int maxWedges = 6;
45  if (trackFinder == bmtf) {
46  maxWedges = 12;
47  }
48  for (int currentWedge = 0; currentWedge < maxWedges; ++currentWedge) {
49  for (auto mu : wedges.at(currentWedge)) {
50  coll1.push_back(mu);
51  }
52  // handle wrap around: max "wedge" has to be compared to first "wedge"
53  int neighbourWedge = (currentWedge + 1) % maxWedges;
54  for (auto mu : wedges.at(neighbourWedge)) {
55  coll2.push_back(mu);
56  }
57  if (mode == cancelmode::coordinate) {
58  getCoordinateCancelBits(coll2, coll1); // in case of a tie coll1 muon wins
59  } else {
60  getTrackAddrCancelBits(coll1, coll2);
61  }
62 
63  coll1.clear();
64  coll2.clear();
65  }
66 }
l1t::MicroGMTCancelOutUnit::getTrackAddrCancelBits
void getTrackAddrCancelBits(std::vector< std::shared_ptr< GMTInternalMuon >> &, std::vector< std::shared_ptr< GMTInternalMuon >> &)
Compares all muons from coll1 with all muons from coll2 and sets the cancel-bits based on track addre...
Definition: MicroGMTCancelOutUnit.cc:194
l1t::MicroGMTCancelOutUnit::getCoordinateCancelBits
void getCoordinateCancelBits(std::vector< std::shared_ptr< GMTInternalMuon >> &, std::vector< std::shared_ptr< GMTInternalMuon >> &)
Compares all muons from coll1 with all muons from coll2 and sets the cancel-bits based on eta/phi coo...
Definition: MicroGMTCancelOutUnit.cc:138
RPCpg::mu
const int mu
Definition: Constants.h:22
void l1t::MicroGMTCancelOutUnit::setCancelOutBitsOverlapBarrel ( GMTInternalWedges &  omtfSectors,
GMTInternalWedges &  bmtfWedges,
cancelmode  mode 
)

Cancel-out between overlap and barrel track finders.

Definition at line 69 of file MicroGMTCancelOutUnit.cc.

References l1t::coordinate, getCoordinateCancelBits(), getTrackAddrCancelBits(), and i.

Referenced by L1TMuonProducer::produce().

70 {
71  // overlap sector collection
72  std::vector<std::shared_ptr<GMTInternalMuon>> coll1;
73  coll1.reserve(3);
74  // barrel wedge collection with 4 wedges
75  std::vector<std::shared_ptr<GMTInternalMuon>> coll2;
76  coll2.reserve(12);
77 
78  for (int currentSector = 0; currentSector < 6; ++currentSector) {
79  for (auto omtfMuon : omtfSectors.at(currentSector)) {
80  coll1.push_back(omtfMuon);
81  }
82  // BMTF | 1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 | 10 | 11 | 0 |
83  // OMTF | 0 | 1 | 2 | 3 | 4 | 5 |
84  // cancel OMTF sector x with corresponding BMTF wedge + the two on either side;
85  // e.g. OMTF 0 with BMTF 0, 1, 2, 3, OMTF 2 with BMTF 4, 5, 6, 7 etc.
86  for (int i = 0; i < 4; ++i) {
87  int currentWedge = (currentSector * 2 + i) % 12;
88  for (auto bmtfMuon : bmtfWedges.at(currentWedge)) {
89  coll2.push_back(bmtfMuon);
90  }
91  }
92  if (mode == cancelmode::coordinate) {
93  getCoordinateCancelBits(coll1, coll2);
94  } else {
95  getTrackAddrCancelBits(coll1, coll2);
96  }
97  coll1.clear();
98  coll2.clear();
99  }
100 }
i
int i
Definition: DBlmapReader.cc:9
l1t::MicroGMTCancelOutUnit::getTrackAddrCancelBits
void getTrackAddrCancelBits(std::vector< std::shared_ptr< GMTInternalMuon >> &, std::vector< std::shared_ptr< GMTInternalMuon >> &)
Compares all muons from coll1 with all muons from coll2 and sets the cancel-bits based on track addre...
Definition: MicroGMTCancelOutUnit.cc:194
l1t::MicroGMTCancelOutUnit::getCoordinateCancelBits
void getCoordinateCancelBits(std::vector< std::shared_ptr< GMTInternalMuon >> &, std::vector< std::shared_ptr< GMTInternalMuon >> &)
Compares all muons from coll1 with all muons from coll2 and sets the cancel-bits based on eta/phi coo...
Definition: MicroGMTCancelOutUnit.cc:138
void l1t::MicroGMTCancelOutUnit::setCancelOutBitsOverlapEndcap ( GMTInternalWedges &  omtfSectors,
GMTInternalWedges &  emtfSectors,
cancelmode  mode 
)

Cancel-out between overlap and endcap track finders.

Definition at line 103 of file MicroGMTCancelOutUnit.cc.

References l1t::coordinate, getCoordinateCancelBits(), getTrackAddrCancelBits(), and i.

Referenced by L1TMuonProducer::produce().

104 {
105  // overlap sector collection
106  std::vector<std::shared_ptr<GMTInternalMuon>> coll1;
107  coll1.reserve(3);
108  // endcap sector collection with 3 sectors
109  std::vector<std::shared_ptr<GMTInternalMuon>> coll2;
110  coll2.reserve(9);
111 
112  for (int curOmtfSector = 0; curOmtfSector < 6; ++curOmtfSector) {
113  for (auto omtfMuon : omtfSectors.at(curOmtfSector)) {
114  coll1.push_back(omtfMuon);
115  }
116  // OMTF | 0 | 1 | 2 | 3 | 4 | 5 |
117  // EMTF | 0 | 1 | 2 | 3 | 4 | 5 |
118  // cancel OMTF sector x with corresponding EMTF sector + the ones on either side;
119  // e.g. OMTF 1 with EMTF 0, 1, 2; OMTF 0 with EMTF 5, 0, 1 etc.
120  for (int i = 0; i < 3; ++i) {
121  // handling the wrap around: adding 5 because 0 has to be compared to 5
122  int curEmtfSector = ((curOmtfSector + 5) + i) % 6;
123  for (auto emtfMuon : emtfSectors.at(curEmtfSector)) {
124  coll2.push_back(emtfMuon);
125  }
126  }
127  if (mode == cancelmode::coordinate) {
128  getCoordinateCancelBits(coll1, coll2);
129  } else {
130  getTrackAddrCancelBits(coll1, coll2);
131  }
132  coll1.clear();
133  coll2.clear();
134  }
135 }
i
int i
Definition: DBlmapReader.cc:9
l1t::MicroGMTCancelOutUnit::getTrackAddrCancelBits
void getTrackAddrCancelBits(std::vector< std::shared_ptr< GMTInternalMuon >> &, std::vector< std::shared_ptr< GMTInternalMuon >> &)
Compares all muons from coll1 with all muons from coll2 and sets the cancel-bits based on track addre...
Definition: MicroGMTCancelOutUnit.cc:194
l1t::MicroGMTCancelOutUnit::getCoordinateCancelBits
void getCoordinateCancelBits(std::vector< std::shared_ptr< GMTInternalMuon >> &, std::vector< std::shared_ptr< GMTInternalMuon >> &)
Compares all muons from coll1 with all muons from coll2 and sets the cancel-bits based on eta/phi coo...
Definition: MicroGMTCancelOutUnit.cc:138

Member Data Documentation

std::shared_ptr<MicroGMTMatchQualLUT> l1t::MicroGMTCancelOutUnit::m_boNegMatchQualLUT
private

Definition at line 35 of file MicroGMTCancelOutUnit.h.

Referenced by initialise().

std::shared_ptr<MicroGMTMatchQualLUT> l1t::MicroGMTCancelOutUnit::m_boPosMatchQualLUT
private

Definition at line 34 of file MicroGMTCancelOutUnit.h.

Referenced by initialise().

std::shared_ptr<MicroGMTMatchQualLUT> l1t::MicroGMTCancelOutUnit::m_brlSingleMatchQualLUT
private

Definition at line 38 of file MicroGMTCancelOutUnit.h.

std::shared_ptr<MicroGMTMatchQualLUT> l1t::MicroGMTCancelOutUnit::m_foNegMatchQualLUT
private

Definition at line 37 of file MicroGMTCancelOutUnit.h.

Referenced by initialise().

std::shared_ptr<MicroGMTMatchQualLUT> l1t::MicroGMTCancelOutUnit::m_foPosMatchQualLUT
private

Definition at line 36 of file MicroGMTCancelOutUnit.h.

Referenced by initialise().

std::shared_ptr<MicroGMTMatchQualLUT> l1t::MicroGMTCancelOutUnit::m_fwdNegSingleMatchQualLUT
private

Definition at line 42 of file MicroGMTCancelOutUnit.h.

Referenced by initialise().

std::shared_ptr<MicroGMTMatchQualLUT> l1t::MicroGMTCancelOutUnit::m_fwdPosSingleMatchQualLUT
private

Definition at line 41 of file MicroGMTCancelOutUnit.h.

Referenced by initialise().

std::map<int, std::shared_ptr<MicroGMTMatchQualLUT> > l1t::MicroGMTCancelOutUnit::m_lutDict
private

Definition at line 43 of file MicroGMTCancelOutUnit.h.

Referenced by getCoordinateCancelBits(), and initialise().

std::shared_ptr<MicroGMTMatchQualLUT> l1t::MicroGMTCancelOutUnit::m_ovlNegSingleMatchQualLUT
private

Definition at line 40 of file MicroGMTCancelOutUnit.h.

Referenced by initialise().

std::shared_ptr<MicroGMTMatchQualLUT> l1t::MicroGMTCancelOutUnit::m_ovlPosSingleMatchQualLUT
private

Definition at line 39 of file MicroGMTCancelOutUnit.h.

Referenced by initialise().

Generated for CMSSW Reference Manual by doxygen 1.8.5