CSCSectorReceiverMiniLUT Class Reference

#include <CSCSectorReceiverMiniLUT.h>

Static Public Member Functions
static global_eta_data	calcGlobalEtaMEMini (unsigned short endcap, unsigned short sector, unsigned short station, unsigned short subsector, unsigned theadd, const bool gangedME1a)
static global_phi_data	calcGlobalPhiMBMini (unsigned short endcap, unsigned short sector, unsigned short subsector, unsigned theadd, const bool gangedME1a)
static global_phi_data	calcGlobalPhiMEMini (unsigned short endcap, unsigned short sector, unsigned short station, unsigned short subsector, unsigned theadd, const bool gangedME1a)
static lclphidat	calcLocalPhiMini (unsigned theadd, const bool gangedME1a)

Static Private Attributes
static const unsigned short int	gbl_eta_bounds [2][6][4][2][4][9][2]
static const float	gbl_eta_params [2][6][4][2][4][9][3]
static const float	gbl_phi_mb_params [2][6][2][9][2]
static const float	gbl_phi_me_params [2][6][4][2][9][2]
static const float	lcl_phi_param0 [1<< 4]
static const float	lcl_phi_param1 = 6.4

Detailed Description

Author

Brett Jackson

Provides a new way of defining the Lookup tables used by the core. Defines the lookup tables as parameterized functions in order to save on memory usage when compared to the current standard definitions of LUTs

Definition at line 15 of file CSCSectorReceiverMiniLUT.h.

Member Function Documentation

calcGlobalEtaMEMini()

global_eta_data CSCSectorReceiverMiniLUT::calcGlobalEtaMEMini ( unsigned short  endcap, unsigned short  sector, unsigned short  station, unsigned short  subsector, unsigned  theadd, const bool  gangedME1a  )

static

Definition at line 50 of file CSCSectorReceiverMiniLUT.cc.

References data, makeMuonMisalignmentScenario::endcap, egammaIdentification::eta_max, egammaIdentification::eta_min, l1temulator_dqm_sourceclient-live_cfg::gangedME1a, gbl_eta_bounds, gbl_eta_params, dqm-mbProfile::log, CSCConstants::MAX_NUM_WIREGROUPS, nano_mu_digi_cff::sector, relativeConstraints::station, and nano_mu_digi_cff::subsector.

Referenced by CSCSectorReceiverLUT::globalEtaME().

                                                                                     {
  if (endcap < 1 || endcap > 2)
    edm::LogWarning("CSCSectorReceiverMiniLUT")
        << "+++ Value of endcap, " << endcap << ", is out of bounds, [1, 2] +++\n";
  if (sector < 1 || sector > 6)
    edm::LogWarning("CSCSectorReceiverMiniLUT")
        << "+++ Value of sector, " << sector << ", is out of bounds, [1, 6] +++\n";
  if (station < 1 || station > 4)
    edm::LogWarning("CSCSectorReceiverMiniLUT")
        << "+++ Value of station, " << station << ", is out of bounds, [1, 4] +++\n";

  gbletadat data(0);

  unsigned short int tcscid = ((theadd >> 15) & 0xf);
  unsigned short int lclPhi = ((theadd >> 6) & 0x3);
  unsigned short int WG = ((theadd >> 8) & 0x7f);
  unsigned short int bend = ((theadd)&0x3f);

  int eta_temp = 999, eta_min = 999, eta_max = 999;

  if ((tcscid > 0) && (tcscid <= 12) && (WG < CSCConstants::MAX_NUM_WIREGROUPS)) {
    unsigned short int cscid = (tcscid > 9) ? tcscid - 9 : tcscid;
    if (station == 1) {
      unsigned short int lclPhip = 0;
      if (lclPhi == 1 || lclPhi == 3)
        lclPhip = 2;
      // use only eta correction for first and last third of ME1/1 chamber since local phi scaling changed
      if (gangedME1a) {
        eta_temp =
            (gbl_eta_params[endcap - 1][sector - 1][station - 1][subsector - 1][lclPhi][cscid - 1][0] +
             gbl_eta_params[endcap - 1][sector - 1][station - 1][subsector - 1][lclPhi][cscid - 1][1] *
                 log(gbl_eta_params[endcap - 1][sector - 1][station - 1][subsector - 1][lclPhi][cscid - 1][2] + WG));
        eta_min = gbl_eta_bounds[endcap - 1][sector - 1][station - 1][subsector - 1][lclPhi][cscid - 1][0];
        eta_max = gbl_eta_bounds[endcap - 1][sector - 1][station - 1][subsector - 1][lclPhi][cscid - 1][1];
      } else {  // DA and MDG, if unganged replace "lclPhi" index with "lclPhip"
        eta_temp =
            (gbl_eta_params[endcap - 1][sector - 1][station - 1][subsector - 1][lclPhip][cscid - 1][0] +
             gbl_eta_params[endcap - 1][sector - 1][station - 1][subsector - 1][lclPhip][cscid - 1][1] *
                 log(gbl_eta_params[endcap - 1][sector - 1][station - 1][subsector - 1][lclPhip][cscid - 1][2] + WG));
        eta_min = gbl_eta_bounds[endcap - 1][sector - 1][station - 1][subsector - 1][lclPhip][cscid - 1][0];
        eta_max = gbl_eta_bounds[endcap - 1][sector - 1][station - 1][subsector - 1][lclPhip][cscid - 1][1];
      }

      // add offset to ME+11a, subtract for ME-11a (wire tilt and strip direction)1
      // only is ganged
      if (gangedME1a && (tcscid < 4) && (lclPhi == 3)) {
        if (endcap == 1)
          eta_temp += 3;
        else
          eta_temp -= 3;
      }
    } else {
      eta_temp = (gbl_eta_params[endcap - 1][sector - 1][station - 1][0][lclPhi][cscid - 1][0] +
                  gbl_eta_params[endcap - 1][sector - 1][station - 1][0][lclPhi][cscid - 1][1] *
                      log(gbl_eta_params[endcap - 1][sector - 1][station - 1][0][lclPhi][cscid - 1][2] + WG));
      eta_min = gbl_eta_bounds[endcap - 1][sector - 1][station - 1][0][lclPhi][cscid - 1][0];
      eta_max = gbl_eta_bounds[endcap - 1][sector - 1][station - 1][0][lclPhi][cscid - 1][1];
    }
  } else {
    edm::LogWarning("CSCSectorReceiverMiniLUT")
        << "+++ Value of cscid, " << tcscid << ", is out of bounds, [1, 9] -- or --"
        << " Value of wire group, " << WG << ", exceeds max allowed, " << CSCConstants::MAX_NUM_WIREGROUPS << " +++\n";
  }

  // protect from negative numbers.  If the value of eta_temp is <0, set global eta to the minimum value
  if ((eta_temp >= eta_min) && (eta_temp <= eta_max))
    data.global_eta = eta_temp;
  else if (eta_temp < eta_min)
    data.global_eta = eta_min;
  else
    data.global_eta = eta_max;

  //  data.global_bend = 0;
  // Just pass through lowest 5 bits of local bend (drop 1 MSB)
  data.global_bend = bend & 0x1F;

  return data;
}

calcGlobalPhiMBMini()

global_phi_data CSCSectorReceiverMiniLUT::calcGlobalPhiMBMini ( unsigned short  endcap, unsigned short  sector, unsigned short  subsector, unsigned  theadd, const bool  gangedME1a  )

static

Definition at line 217 of file CSCSectorReceiverMiniLUT.cc.

References data, makeMuonMisalignmentScenario::endcap, l1temulator_dqm_sourceclient-live_cfg::gangedME1a, gbl_phi_mb_params, CSCBitWidths::kGlobalPhiDataBitWidth, CSCBitWidths::kLocalPhiDataBitWidth, nano_mu_digi_cff::sector, and nano_mu_digi_cff::subsector.

                                                                                                                    {
  if (endcap < 1 || endcap > 2)
    edm::LogWarning("CSCSectorReceiverMiniLUT")
        << "+++ Value of endcap, " << endcap << ", is out of bounds, [1, 2] +++\n";
  if (sector < 1 || sector > 6)
    edm::LogWarning("CSCSectorReceiverMiniLUT")
        << "+++ Value of sector, " << sector << ", is out of bounds, [1, 6] +++\n";

  gblphidat data(0);

  unsigned short int maxPhiL = 1 << CSCBitWidths::kLocalPhiDataBitWidth;
  unsigned short int maxPhiG = 1 << CSCBitWidths::kGlobalPhiDataBitWidth;
  unsigned short int cscid = ((theadd >> 15) & 0xf);
  unsigned short int lclPhi = (theadd & 0x3ff);
  if (!gangedME1a)
    lclPhi = lclPhi / 0.625;  // DA and MDG, recover old scaling of local phi

  if ((cscid <= 3) || (cscid >= 7))
    maxPhiL =
        maxPhiL *
        (64. / 80);  // currently a hack that is in place to handle the different number of strips in ME1/1 and ME1/3

  if ((cscid > 0) && (cscid <= 9)) {
    if (lclPhi < maxPhiL)
      data.global_phi = (gbl_phi_mb_params[endcap - 1][sector - 1][subsector - 1][cscid - 1][0] +
                         gbl_phi_mb_params[endcap - 1][sector - 1][subsector - 1][cscid - 1][1] * lclPhi);
    else
      data.global_phi = (gbl_phi_mb_params[endcap - 1][sector - 1][subsector - 1][cscid - 1][0] +
                         gbl_phi_mb_params[endcap - 1][sector - 1][subsector - 1][cscid - 1][1] * (maxPhiL - 1));
  } else
    edm::LogWarning("CSCSectorReceiverMiniLUT")
        << "+++ Value of cscid, " << cscid << ", is out of bounds, [1, 9] +++\n";

  if (data.global_phi >= maxPhiG)
    edm::LogWarning("CSCSectorReceiverMiniLUT")
        << "+++ Value of global_phi, " << data.global_phi << ", exceeds max allowed, " << maxPhiG - 1 << " +++\n";

  if (data.global_phi >= maxPhiG)
    edm::LogWarning("CSCSectorReceiverMiniLUT")
        << "+++ Value of global_phi, " << data.global_phi << ", exceeds max allowed, " << maxPhiG - 1 << " +++\n";

  return data;
}

calcGlobalPhiMEMini()

global_phi_data CSCSectorReceiverMiniLUT::calcGlobalPhiMEMini ( unsigned short  endcap, unsigned short  sector, unsigned short  station, unsigned short  subsector, unsigned  theadd, const bool  gangedME1a  )

static

Definition at line 134 of file CSCSectorReceiverMiniLUT.cc.

References data, makeMuonMisalignmentScenario::endcap, l1temulator_dqm_sourceclient-live_cfg::gangedME1a, gbl_phi_me_params, CSCBitWidths::kGlobalPhiDataBitWidth, CSCBitWidths::kLocalPhiDataBitWidth, CSCConstants::MAX_NUM_STRIPS_RUN1, nano_mu_digi_cff::sector, relativeConstraints::station, nano_mu_digi_cff::strip, and nano_mu_digi_cff::subsector.

Referenced by CSCSectorReceiverLUT::globalPhiME().

                                                                                     {
  if (endcap < 1 || endcap > 2)
    edm::LogWarning("CSCSectorReceiverMiniLUT")
        << "+++ Value of endcap, " << endcap << ", is out of bounds, [1, 2] +++\n";
  if (sector < 1 || sector > 6)
    edm::LogWarning("CSCSectorReceiverMiniLUT")
        << "+++ Value of sector, " << sector << ", is out of bounds, [1, 6] +++\n";
  if (station < 1 || station > 4)
    edm::LogWarning("CSCSectorReceiverMiniLUT")
        << "+++ Value of station, " << station << ", is out of bounds, [1, 4] +++\n";

  gblphidat data(0);

  unsigned short int maxPhiL = 1 << CSCBitWidths::kLocalPhiDataBitWidth;
  unsigned short int maxPhiG = 1 << CSCBitWidths::kGlobalPhiDataBitWidth;
  unsigned short int cscid = ((theadd >> 15) & 0xf);
  unsigned short int lclPhi = (theadd & 0x3ff);
  if (!gangedME1a)
    lclPhi = lclPhi / 0.625;  // DA and MDG recover old scaling of local phi

  // 12/11/09
  // GP et DA: how to identify the strip number and isolate and shift the localPhi value
  const double binPhiL = static_cast<double>(maxPhiL) / (2 * CSCConstants::MAX_NUM_STRIPS_RUN1);

  int strip = static_cast<int>(lclPhi / binPhiL);
  if (station == 1 && (cscid <= 3) &&
      (strip >= 127 && strip < 224)) {  // 160 --> 224, change range for ME1/1a acceptance, DA and MDG
    // in this case need to redefine lclPhi in order to
    // place local phi in the middle of the 5th CFEB
    // and not on the first third of the CFEB as default

    gangedME1a
        ? lclPhi = (strip - 127 + 31) * (4 * binPhiL / 3)
        : lclPhi =
              (strip - 127) *
              (4 * binPhiL /
               3);  //DA and MDG remove offset to center of ME1/1a (no ganging), and reset ME1/1a strip number to start from 0, and scale 48 strips to match ME1/1b 64 strips
  }
  // end GP et DA

  if (station == 1 && ((cscid <= 3) || (cscid >= 7))) {
    //if ( (strip >= 127 && strip < 160) || (cscid >= 10) ) // VK: the || (cscid >= 10) for unganged ME1a
    //  maxPhiL = maxPhiL*(48./80); // GP et DA: currently a hack that is in place to handle the different number of strips in ME1/1a and ME1/3
    //else
    maxPhiL =
        maxPhiL *
        (64. / 80);  // currently a hack that is in place to handle the different number of strips in ME1/1 and ME1/3
  }

  // VK: The the unganged ME1a hack
  if (station == 1 && (cscid >= 10)) {
    lclPhi = strip * (4 * binPhiL / 3);
    cscid = cscid - 9;  // back to normal 1-9 range
  }
  // end VK

  if ((cscid > 0) && (cscid <= 9)) {
    if ((station == 1) && (lclPhi < maxPhiL))
      data.global_phi = (gbl_phi_me_params[endcap - 1][sector - 1][station - 1][subsector - 1][cscid - 1][0] +
                         gbl_phi_me_params[endcap - 1][sector - 1][station - 1][subsector - 1][cscid - 1][1] * lclPhi);
    else if ((station == 1) && (lclPhi >= maxPhiL))
      data.global_phi =
          (gbl_phi_me_params[endcap - 1][sector - 1][station - 1][subsector - 1][cscid - 1][0] +
           gbl_phi_me_params[endcap - 1][sector - 1][station - 1][subsector - 1][cscid - 1][1] * (maxPhiL - 1));
    else
      data.global_phi = (gbl_phi_me_params[endcap - 1][sector - 1][station - 1][0][cscid - 1][0] +
                         gbl_phi_me_params[endcap - 1][sector - 1][station - 1][0][cscid - 1][1] * lclPhi);
  } else
    edm::LogWarning("CSCSectorReceiverMiniLUT")
        << "+++ Value of cscid, " << cscid << ", is out of bounds, [1, 9] +++\n";

  if (data.global_phi >= maxPhiG)
    edm::LogWarning("CSCSectorReceiverMiniLUT")
        << "+++ Value of global_phi, " << data.global_phi << ", exceeds max allowed, " << maxPhiG - 1 << " +++\n";

  return data;
}

calcLocalPhiMini()

lclphidat CSCSectorReceiverMiniLUT::calcLocalPhiMini ( unsigned  theadd, const bool  gangedME1a  )

static

Definition at line 18 of file CSCSectorReceiverMiniLUT.cc.

References ALPAKA_ACCELERATOR_NAMESPACE::brokenline::constexpr(), data, l1temulator_dqm_sourceclient-live_cfg::gangedME1a, CSCBitWidths::kLocalPhiDataBitWidth, lcl_phi_param0, lcl_phi_param1, CSCConstants::MAX_NUM_STRIPS_RUN1, CSCConstants::NUM_CLCT_PATTERNS, topSingleLeptonDQM_PU_cfi::pattern, and nano_mu_digi_cff::strip.

Referenced by CSCSectorReceiverLUT::localPhi().

                                                                                           {
  // This method is ripped from CSCSectorReceverLUT.cc with minor changes

  lclphidat data;

  constexpr int maxPhiL = 1 << CSCBitWidths::kLocalPhiDataBitWidth;
  unsigned short int pattern = ((theadd >> 8) & 0xf);
  unsigned short int strip = (theadd & 0xff);

  if (strip < 2 * (CSCConstants::MAX_NUM_STRIPS_RUN1 * 7 / 5) &&
      pattern <
          CSCConstants::
              NUM_CLCT_PATTERNS) {  // MDG, DA and RW, for ME1 we have 7CFEBs and not just 5, so the num_strips can go up to 16 * 7 but only for ME1
    data.phi_local = gangedME1a ? static_cast<unsigned>((lcl_phi_param0[pattern] + strip) * lcl_phi_param1)
                                : static_cast<unsigned>((lcl_phi_param0[pattern] + strip) * 0.625 * lcl_phi_param1);
    //DA and MDG, rescale range of local phi so ME1/1b fits in 0-511
  } else
    edm::LogWarning("CSCSectorReceiverMiniLUT") << "+++ Value of strip, " << strip << ", exceeds max allowed, "
                                                << 2 * CSCConstants::MAX_NUM_STRIPS_RUN1 - 1 << " +++\n";

  if (data.phi_local >= maxPhiL)
    edm::LogWarning("CSCSectorReceiverMiniLUT")
        << "+++ Value of phi_local, " << data.phi_local << ", exceeds max allowed, "
        << CSCConstants::NUM_CLCT_PATTERNS - 1 << " +++\n";

  //  data.phi_bend_local = 0;
  // Just pass through all bits of pattern as bend angle (so 2 MSB unfilled)
  data.phi_bend_local = pattern & 0x3F;

  return data;
}

Member Data Documentation

gbl_eta_bounds

const unsigned short int CSCSectorReceiverMiniLUT::gbl_eta_bounds

staticprivate

Definition at line 40 of file CSCSectorReceiverMiniLUT.h.

Referenced by calcGlobalEtaMEMini().

gbl_eta_params

const float CSCSectorReceiverMiniLUT::gbl_eta_params

staticprivate

Definition at line 38 of file CSCSectorReceiverMiniLUT.h.

Referenced by calcGlobalEtaMEMini().

gbl_phi_mb_params

const float CSCSectorReceiverMiniLUT::gbl_phi_mb_params

staticprivate

Definition at line 44 of file CSCSectorReceiverMiniLUT.h.

Referenced by calcGlobalPhiMBMini().

gbl_phi_me_params

const float CSCSectorReceiverMiniLUT::gbl_phi_me_params

staticprivate

Definition at line 42 of file CSCSectorReceiverMiniLUT.h.

Referenced by calcGlobalPhiMEMini().

lcl_phi_param0

const float CSCSectorReceiverMiniLUT::lcl_phi_param0

staticprivate

Initial value:

= {
    0.50, 0.50, -0.10, 1.10, -0.14, 1.14, 0.27, 0.73, 0.29, 0.71, 0.50, 0, 0, 0, 0, 0}

Definition at line 34 of file CSCSectorReceiverMiniLUT.h.

Referenced by calcLocalPhiMini().

lcl_phi_param1

const float CSCSectorReceiverMiniLUT::lcl_phi_param1 = 6.4

staticprivate

Definition at line 35 of file CSCSectorReceiverMiniLUT.h.

Referenced by calcLocalPhiMini().