L1MuGMTLFSortRankEtaQLUT.cc

Go to the documentation of this file.

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

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

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

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

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

void L1MuGMTLFSortRankEtaQLUT::InitParameters() {
  m_version = L1MuGMTConfig::getVersionSortRankEtaQLUT();
}

//------------------------------------------------------------------------------------
// Sort Rank LUT, Eta-Quality part 
// 
// This LUT determines the dependency of the sort rank on Eta and Quality. 
// Inputs:  Eta(6 bit) and Quality(3 bit)
// Outputs: Very-low-quality bits (2 bit)
//          Rank contribution 2-bit
//
// The very-low quality bits are used to suppress unmatched muons in certain 
// trigger-algorithms. By default vlq(1)(MSB) means: not to be used in di-muon trigger
// and vlq(0)(LSB) means not to be used in single-muon trigger. 
//
// The Rank-contribution from Eta-Quality is currently not used. Its influence on the
// sort-rank can be determined in the Sort-Rank Combine LUT.
// 
// Explanation of the meaning of m_version:
//
// m_version split to 4bit fields:
// bits 0-3: qCSC=2 code
//           - 0 - do nothing
//           - 1 - all to qGMT=3
//           - 2 - |eta|<1.5 || |eta|>1.8 to qGMT=3
//           - 3 - 1.2<|eta|<1.5 || |eta|>1.8 to qGMT=3
//           - >3 - do nothing
// bits 4-7: qCSC=1 code
//           - 0 - |eta|>1.2 to qGMT=3
//           - 1 - |eta|>1.3 to qGMT=3
//           - >1 - nothing done
// bits 8-11: qRPC<3 code
//           - 0 - qRPC=0 1.04<|eta|<1.24 || |eta|>1.48 to qGMT=2
//                 qRPC=1 0.83<|eta|<1.04 || 1.14<|eta|<1.24 || |eta|>1.36 to qGMT=2
//                 qRPC=2 0.83<|eta|<0.93 to qGMT=2
//           - 1 - do nothing
//           - 2 - remove qRPC=0 in endcap from single muon
//           - 3 - remove qRPC=0 from single muon
//           - 4 - remove qRPC=0 in endcap from all
//           - 5 - remove qRPC=0 from all
//           - >0 - do nothing
//
// m_version=1: TDR
// m_version=2: 2010 data
// m_version>2: 2011
//
//------------------------------------------------------------------------------------

unsigned L1MuGMTLFSortRankEtaQLUT::TheLookupFunction (int idx, unsigned eta, unsigned q) const {
  // idx is DT, BRPC, CSC, FRPC
  // INPUTS:  eta(6) q(3)
  // OUTPUTS: vlq(2) rank_etaq(2) 

  const L1MuTriggerScales* theTriggerScales = L1MuGMTConfig::getTriggerScales();

  int isRPC = idx % 2;
  //  int isFWD = idx / 2;

  float etaValue = fabs ( theTriggerScales->getRegionalEtaScale(idx)->getCenter( eta ) );

  //
  // very-low-quality flags
  //
  // 0 .. no VLQ set
  // 1 .. output muon will be quality 2 (disable in single and di-muon trigger)
  // 2 .. output muon will be quality 3 (disable in single-muon trigger only)
  // 3 .. output muon will be quality 4 (disable in di-muon trigger only)

  unsigned vlq = 0;
  
  int vCSC2 = (m_version) & 0xf;
  int vCSC1 = (m_version>>4) & 0xf;
  int vRPC =  (m_version>>8) & 0xf;
  
  // RPC selection
  if (isRPC) {
    if(vRPC == 0) {
      if ( ( q == 0 && ( ( etaValue > 1.04 && etaValue < 1.24 ) || // Q0, high rate, high noise
          ( etaValue > 1.48 ) ) ) ||                // Q0, high rate, high noise
          ( q == 1 && ( ( etaValue > 0.83 && etaValue < 1.04 ) || // Q1, high rate
              ( etaValue > 1.14 && etaValue < 1.24 ) || // Q1, high noise
              ( etaValue > 1.36 ) ) ) ||                // Q1, high rate
              ( q == 2 && ( etaValue > 0.83 && etaValue < 0.93 ) ) )  // Q2, high rate
        vlq = 1;
    }
// reserve vRPC == 1 for doing nothing for bwd compatibility
    if(vRPC == 2) {
      if(q == 0 && etaValue > 1.04) vlq = 2;
    }
    if(vRPC == 3) {
      if(q == 0) vlq = 2;
    }
    if(vRPC == 4) {
      if(q == 0 && etaValue > 1.04) vlq = 1;
    }
    if(vRPC == 5) {
      if(q == 0) vlq = 1;
    }
  }                         
      
  // CSC selection
  if ( idx == 2 ) { // CSC
    if (q == 2) {
      if(vCSC2 == 1) vlq = 2;
      if(vCSC2 == 2) {
        if(etaValue < 1.5 || etaValue > 1.8) vlq = 2; // disable in single-muon trigger only
      }
      if(vCSC2 == 3) {
        if( (etaValue > 1.2 && etaValue < 1.5) || etaValue > 1.8) vlq = 2; // disable in single-muon trigger only
      }
    }
    
    if (q == 1) {
      if(vCSC1 == 0) {
        if(etaValue > 1.2) vlq = 2;   // disable in single-muon trigger only
      }
      if(vCSC1 == 1) {
        if(etaValue > 1.3) vlq = 2;   // disable in single-muon trigger only
      }
    }
  }
    
  if ( L1MuGMTConfig::getDoOvlRpcAnd() ) {
    if ( idx == 0 ) { // DT
      if ( etaValue > 0.91 ) vlq = 1;
    }
    if ( idx == 2 ) { // CSC
      if ( etaValue <  1.06 ) vlq = 1;  
    }
  }

  //
  // Rank contribution from eta and quality
  //

  // by default return maximum 
  // LUT can later be used to reduce the sort rank of certain regions
  unsigned rank_etaq = 3;

  return (vlq << 2)  | rank_etaq;
}