MicroGMTConfiguration.cc

Go to the documentation of this file.

#include "L1Trigger/L1TMuon/interface/MicroGMTConfiguration.h"

unsigned
l1t::MicroGMTConfiguration::getTwosComp(const int signed_int, const int width) {
  if (signed_int >= 0) {
    return (unsigned)signed_int;
  }
  int all_one = (1 << width)-1;
  return ((-signed_int) ^ all_one) + 1;
}

int
l1t::MicroGMTConfiguration::calcGlobalPhi(int locPhi, tftype t, int proc) {
  int globPhi = 0;
  if (t == bmtf) {
      // each BMTF processor corresponds to a 30 degree wedge = 48 in int-scale
      globPhi = (proc) * 48 + locPhi;
      // first processor starts at CMS phi = -15 degrees...
      globPhi += 576-24;
      // handle wrap-around (since we add the 576-24, the value will never be negative!)
      globPhi = globPhi%576;
  } else {
      // all others correspond to 60 degree sectors = 96 in int-scale
      globPhi = (proc) * 96 + locPhi;
      // first processor starts at CMS phi = 15 degrees (24 in int)... Handle wrap-around with %. Add 576 to make sure the number is positive
      globPhi = (globPhi + 600) % 576;
  }
  return globPhi;
}

int
l1t::MicroGMTConfiguration::setOutputMuonQuality(int muQual, tftype type, int haloBit) {
  // use only the two MSBs for the muon to the uGT
  int outQual = muQual & 0xC;
  if (haloBit == 1 && (type == tftype::emtf_neg || type == tftype::emtf_pos)) {
    // set quality to 0xF if the halo bit is set
    outQual = 0xF;
  }
  return outQual;
}

int
l1t::MicroGMTConfiguration::calcMuonHwEtaExtra(const l1t::Muon& mu)
{
  return mu.hwEta() + mu.hwDEtaExtra();
}

int
l1t::MicroGMTConfiguration::calcMuonHwPhiExtra(const l1t::Muon& mu)
{
  auto hwPhiExtra = mu.hwPhi() + mu.hwDPhiExtra();
  while (hwPhiExtra < 0) {
    hwPhiExtra += 576;
  }
  while (hwPhiExtra > 575) {
    hwPhiExtra -= 576;
  }
  return hwPhiExtra;
}

double
l1t::MicroGMTConfiguration::calcMuonEtaExtra(const l1t::Muon& mu)
{
  return calcMuonHwEtaExtra(mu) * 0.010875;
}

double
l1t::MicroGMTConfiguration::calcMuonPhiExtra(const l1t::Muon& mu)
{
  //use the LorentzVector to get phi in the range -pi to +pi exactly as in the emulator
  math::PtEtaPhiMLorentzVector vec{0., 0., calcMuonHwPhiExtra(mu)*0.010908, 0.};
  return vec.phi();
}