L1MuGMTLFCOUDeltaEtaLUT.cc

Go to the documentation of this file.

//-------------------------------------------------
//
//   Class: L1MuGMTLFCOUDeltaEtaLUT
//
// 
//
//   Author :
//   H. Sakulin            HEPHY Vienna
//
//   Migrated to CMSSW:
//   I. Mikulec
//
//--------------------------------------------------

//-----------------------
// This Class's Header --
//-----------------------
#include "L1Trigger/GlobalMuonTrigger/src/L1MuGMTLFCOUDeltaEtaLUT.h"

//---------------
// C++ Headers --
//---------------

//#include <iostream>

//-------------------------------
// Collaborating Class Headers --
//-------------------------------
#include "L1Trigger/GlobalMuonTrigger/src/L1MuGMTConfig.h"
#include "CondFormats/L1TObjects/interface/L1MuGMTScales.h"
#include "CondFormats/L1TObjects/interface/L1MuPacking.h"

//-------------------
// InitParameters  --
//-------------------

void L1MuGMTLFCOUDeltaEtaLUT::InitParameters() {
}

//---------------------------------------------------------------------------
// Cancel-Out-Unit Delta-Eta LUT
// =============================
// 
// Calculate delta-eta between two eta values coded in a special 4-bit scale 
// used in the overlap region
//
// This scale is used to save I/O pins and in order to be able to do the delta-eta 
// calculateion in a distributed RAM LUT
//
//---------------------------------------------------------------------------

unsigned L1MuGMTLFCOUDeltaEtaLUT::TheLookupFunction (int idx, unsigned eta1, unsigned eta2) const {
  // idx is DTCSC, CSCDT, bRPCCSC, fRPCDT
  // INPUTS:  eta1(4) eta2(4)
  // OUTPUTS: delta_eta(4) 

  const L1MuGMTScales* theGMTScales = L1MuGMTConfig::getGMTScales();

  // check out of range in inputs
  L1MuSignedPacking<4> pack;
  unsigned delta_eta_OOR = pack.packedFromIdx (-8);

  if (eta1==7 || eta2 ==7) return delta_eta_OOR;

  int type1=0, type2=0;
  switch (idx) {
    case 0: type1 = 0; type2 = 2; break; // DT, CSC
    case 1: type1 = 2; type2 = 0; break; // CSC, DT
    case 2: type1 = 1; type2 = 2; break; // bRPC, CSC
    case 3: type1 = 3; type2 = 0; break; // fRPC, DT
  }

  float etaValue1 = theGMTScales->getOvlEtaScale(type1)->getCenter(eta1);
  float etaValue2 = theGMTScales->getOvlEtaScale(type2)->getCenter(eta2); 

  float delta_eta = etaValue1 - etaValue2;

  unsigned delta_eta_4bit = 0;

  // check out of range
  if (delta_eta < theGMTScales->getDeltaEtaScale(idx+2)->getScaleMin() ||
      delta_eta > theGMTScales->getDeltaEtaScale(idx+2)->getScaleMax()) {
    delta_eta_4bit = delta_eta_OOR;
  }
  else {
    delta_eta_4bit = theGMTScales->getDeltaEtaScale(idx+2)->getPacked( delta_eta );
  }
    
  return delta_eta_4bit;
}