PtAssignmentEngineAux.cc

Go to the documentation of this file.

#include "L1Trigger/L1TMuonEndCap/interface/PtAssignmentEngineAux.h"

// _____________________________________________________________________________
static const int GMT_eta_from_theta[128] = {
    239, 235, 233, 230, 227, 224, 222, 219, 217, 214, 212, 210, 207, 205, 203, 201, 199, 197, 195, 193, 191, 189,
    187, 186, 184, 182, 180, 179, 177, 176, 174, 172, 171, 169, 168, 166, 165, 164, 162, 161, 160, 158, 157, 156,
    154, 153, 152, 151, 149, 148, 147, 146, 145, 143, 142, 141, 140, 139, 138, 137, 136, 135, 134, 133, 132, 131,
    130, 129, 128, 127, 126, 125, 124, 123, 122, 121, 120, 119, 118, 117, 116, 116, 115, 114, 113, 112, 111, 110,
    110, 109, 108, 107, 106, 106, 105, 104, 103, 102, 102, 101, 100, 99,  99,  98,  97,  96,  96,  95,  94,  93,
    93,  92,  91,  91,  90,  89,  89,  88,  87,  87,  86,  85,  84,  84,  83,  83,  82,  81};

int PtAssignmentEngineAux::getGMTPt(float pt) const {
  // compressed pt = pt*2 (scale) + 1 (pt = 0 is empty candidate)
  int gmt_pt = (pt * 2) + 1;
  gmt_pt = (gmt_pt > 511) ? 511 : gmt_pt;
  return gmt_pt;
}

int PtAssignmentEngineAux::getGMTPtDxy(float pt) const {
  // compressed pt = pt*1 (scale) + 1 (pt = 0 is empty candidate)
  int gmt_pt_dxy = (pt * 1) + 1;
  gmt_pt_dxy = (gmt_pt_dxy > 255) ? 255 : gmt_pt_dxy;
  return gmt_pt_dxy;
}

int PtAssignmentEngineAux::getGMTDxy(float dxy) const {
  int gmt_dxy = 0;
  if (std::abs(dxy) < 32.) {
    gmt_dxy = 0;
  } else if (std::abs(dxy) < 64.) {
    gmt_dxy = 1;
  } else if (std::abs(dxy) < 96.) {
    gmt_dxy = 2;
  } else {
    gmt_dxy = 3;
  }
  return gmt_dxy;
}

float PtAssignmentEngineAux::getPtFromGMTPt(int gmt_pt) const {
  float pt = (gmt_pt <= 0) ? 0 : 0.5 * (gmt_pt - 1);
  return pt;
}

float PtAssignmentEngineAux::getPtFromGMTPtDxy(int gmt_pt_dxy) const {
  float pt = (gmt_pt_dxy <= 0) ? 0 : 1.0 * (gmt_pt_dxy - 1);
  return pt;
}

int PtAssignmentEngineAux::getGMTPhi(int phi) const {
  // convert phi into gmt scale according to DN15-017
  // full scale is -16 to 100, or 116 values, covers range -10 to 62.5 deg
  // my internal ph scale is 0..5000, covers from -22 to 63.333 deg
  // converted to GMT scale it is from -35 to 95
  // bt_phi * 107.01/4096, equivalent to bt_phi * 6849/0x40000
  phi *= 6849;
  phi >>= 18;  // divide by 0x40000
  phi -= 35;   // offset of -22 deg
  return phi;
}

int PtAssignmentEngineAux::getGMTPhiV2(int phi) const {
  // convert phi into gmt scale according to DN15-017
  phi *= 6991;
  phi >>= 18;  // divide by 0x40000
  phi -= 35;   // offset of -22 deg
  return phi;
}

int PtAssignmentEngineAux::getGMTEta(int theta, int endcap) const {  // [-1,+1]
  if (theta < 0)
    return 0;
  if (endcap == -1 && theta > 127)
    return -240;
  if (endcap == +1 && theta > 127)
    return 239;

  int eta = GMT_eta_from_theta[theta];
  if (endcap == -1)
    eta = -eta;
  return eta;
}

int PtAssignmentEngineAux::getGMTQuality(int mode, int theta, bool promoteMode7, int version) const {
  int quality = 0;

  if (version == 2) {  // 2018 mode <--> quality mapping
    switch (mode) {
      case 15:
        quality = 15;
        break;
      case 14:
        quality = 14;
        break;
      case 13:
        quality = 13;
        break;
      case 12:
        quality = 7;
        break;
      case 11:
        quality = 12;
        break;
      case 10:
        quality = 10;
        break;
      case 9:
        quality = 9;
        break;
      case 7:
        quality = 11;
        break;
      case 6:
        quality = 6;
        break;
      case 5:
        quality = 5;
        break;
      case 3:
        quality = 4;
        break;
      default:
        quality = 0;
        break;
    }

    if (theta > 88) {  // if (eta < 1.19)
      quality /= 4;
    }
  }  // End 2018 mode <--> quality mapping

  else {               // 2016 and 2017 mode <--> quality mapping
    if (theta > 87) {  // if (eta < 1.2)
      switch (mode) {
        case 15:
          quality = 8;
          break;
        case 14:
          quality = 4;
          break;
        case 13:
          quality = 4;
          break;
        case 12:
          quality = 4;
          break;
        case 11:
          quality = 4;
          break;
        default:
          quality = 4;
          break;
      }
    } else {
      switch (mode) {
        case 15:
          quality = 12;
          break;
        case 14:
          quality = 12;
          break;
        case 13:
          quality = 12;
          break;
        case 12:
          quality = 8;
          break;
        case 11:
          quality = 12;
          break;
        case 10:
          quality = 8;
          break;
        case 7:
          quality = 8;
          break;
        default:
          quality = 4;
          break;
      }
    }
    quality |= (mode & 3);
  }  // End 2016 and 2017 mode <--> quality mapping

  // Fix for missing CSC LCTs in ME1/1, including dead "water-leak" chambers ME1/1/34 and 35
  if (promoteMode7 && mode == 7 && theta <= 50)
    quality = 12;

  return quality;
}

std::pair<int, int> PtAssignmentEngineAux::getGMTCharge(int mode, const std::vector<int>& phidiffs) const {
  // -1 = postive physical charge to match pdgId code (i.e. -13 is positive, anti-muon). +1 = negative physical charge.
  // Also matches DN-2015/017 format for track finder --> uGMT interface format, where 0 indicates positive, 1 negative.
  int emuCharge = 0;

  // Note: sign_ph[0] == 1 in firmware actually translates to phidiffs[0] >= 0 (instead of phidiffs[0] > 0 in the old emulator)
  // The effect needs to be checked

  switch (mode) {
    case 15:                 // 1-2-3-4
      if (phidiffs[0] >= 0)  // 1-2 (should use > 0)
        emuCharge = 1;
      else if (phidiffs[0] == 0 && phidiffs[1] < 0)  // 1-3
        emuCharge = 1;
      else if (phidiffs[1] == 0 && phidiffs[2] < 0)  // 1-4
        emuCharge = 1;
      else
        emuCharge = -1;
      break;

    case 14:                // 1-2-3
      if (phidiffs[0] < 0)  // 1-2
        emuCharge = -1;
      else if (phidiffs[0] == 0 && phidiffs[1] < 0)  // 1-3
        emuCharge = -1;
      else
        emuCharge = 1;
      break;

    case 13:                 // 1-2-4
      if (phidiffs[0] >= 0)  // 1-2 (should use > 0)
        emuCharge = 1;
      else if (phidiffs[0] == 0 && phidiffs[2] < 0)  // 1-4
        emuCharge = 1;
      else
        emuCharge = -1;
      break;

    case 12:                 // 1-2
      if (phidiffs[0] >= 0)  // 1-2
        emuCharge = 1;
      else
        emuCharge = -1;
      break;

    case 11:                 // 1-3-4
      if (phidiffs[1] >= 0)  // 1-3 (should use > 0)
        emuCharge = 1;
      else if (phidiffs[1] == 0 && phidiffs[2] < 0)  // 1-4
        emuCharge = 1;
      else
        emuCharge = -1;
      break;

    case 10:                 // 1-3
      if (phidiffs[1] >= 0)  // 1-3
        emuCharge = 1;
      else
        emuCharge = -1;
      break;

    case 9:                  // 1-4
      if (phidiffs[2] >= 0)  // 1-4
        emuCharge = 1;
      else
        emuCharge = -1;
      break;

    case 7:                  // 2-3-4
      if (phidiffs[3] >= 0)  // 2-3 (should use > 0)
        emuCharge = 1;
      else if (phidiffs[3] == 0 && phidiffs[4] < 0)  // 2-4
        emuCharge = 1;
      else
        emuCharge = -1;
      break;

    case 6:                  // 2-3
      if (phidiffs[3] >= 0)  // 2-3
        emuCharge = 1;
      else
        emuCharge = -1;
      break;

    case 5:                  // 2-4
      if (phidiffs[4] >= 0)  // 2-4
        emuCharge = 1;
      else
        emuCharge = -1;
      break;

    case 3:                  // 3-4
      if (phidiffs[5] >= 0)  // 3-4
        emuCharge = 1;
      else
        emuCharge = -1;
      break;

    default:
      //emuCharge = -1;
      emuCharge = 0;
      break;
  }

  int charge = 0;
  if (emuCharge == 1)
    charge = 1;

  int charge_valid = 1;
  if (emuCharge == 0)
    charge_valid = 0;
  return std::make_pair(charge, charge_valid);
}