L1MuGMTLFPhiProLUT.cc

Go to the documentation of this file.

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

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

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

//#include <iostream>

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

#include "L1Trigger/GlobalMuonTrigger/src/L1MuGMTPhiLUT.h"

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

void L1MuGMTLFPhiProLUT::InitParameters() {
}

//------------------------
// The Lookup Function  --
//------------------------

//
// Phi projection LUT to project the phi coordinates to the vertex
// 
// The output is 9 bits signed and should be in the range of -32 <= dphi < 32
// (+/- 80 deg)
// 
// The input eta is converted from six to four bits as in the MIP/ISO AU chips
//
// The same parameterization as in the MIP/ISO AU chips can be used (proj. to vertex for ISO).

unsigned L1MuGMTLFPhiProLUT::TheLookupFunction (int idx, unsigned eta, unsigned pt, unsigned charge) const {
  // idx is DT, BRPC, CSC, FRPC
  // INPUTS:  eta(4) pt(5) charge(1)
  // OUTPUTS: dphi(9) 

//  const L1MuTriggerScales* theTriggerScales = L1MuGMTConfig::getTriggerScales();
  const L1MuTriggerPtScale* theTriggerPtScale = L1MuGMTConfig::getTriggerPtScale();

  //  static bool doProjection = SimpleConfigurable<bool> (false, "L1GlobalMuonTrigger:PropagatePhi" );
  static bool doProjection = L1MuGMTConfig::getPropagatePhi();

  if (!doProjection) return 0;

  int isRPC = idx % 2;
  int isFWD = idx / 2;
      
  int isys = isFWD + 2 * isRPC; // DT, CSC, BRPC, FRPC
  int ch_idx = (charge == 0) ? 1 : 0; // positive charge is 0 (but idx 1)

  // currently only support 3-bit eta (3 lower bits); ignore 4th bit
  if (eta>7) eta -= 8;

  float dphi =  L1MuGMTPhiLUT::dphi (isys, 1, ch_idx, (int) eta, 
     theTriggerPtScale->getPtScale()->getLowEdge(pt) );  // use old LUT, here
  // theTriggerScales->getPtScale()->getLowEdge(pt) );  // use old LUT, here
  
  int dphi_int = (int) ( (-dphi + 1.25 / 180. * M_PI + 2* M_PI ) / ( 2.5 / 180. * M_PI ) ) - 144;
    
  L1MuSignedPacking<9> PhiPacking;
  return PhiPacking.packedFromIdx(dphi_int);
}