GlobalMuonTrigger/src/L1MuGMTPhiLUT.cc

0001 //-------------------------------------------------
0002 //
0003 //   Class: L1MuGMTPhiLUT
0004 //
0005 //   Description: Look-up table for GMT Phi Projection Unit
0006 //
0007 //
0008 //
0009 //   Author :
0010 //   H. Sakulin                CERN EP
0011 //
0012 //   Migrated to CMSSW:
0013 //   I. Mikulec
0014 //
0015 //--------------------------------------------------
0016
0017 //-----------------------
0018 // This Class's Header --
0019 //-----------------------
0020
0021 #include "L1Trigger/GlobalMuonTrigger/src/L1MuGMTPhiLUT.h"
0022
0023 //---------------
0024 // C++ Headers --
0025 //---------------
0026
0027 #include <iostream>
0028 #include <vector>
0029 #include <cmath>
0030
0031 //-------------------------------
0032 // Collaborating Class Headers --
0033 //-------------------------------
0034
0035 #include "FWCore/MessageLogger/interface/MessageLogger.h"
0036
0037 // --------------------------------
0038 //       class L1MuGMTPhiLUT
0039 //---------------------------------
0040
0041 //----------------
0042 // Constructors --
0043 //----------------
0044 L1MuGMTPhiLUT::L1MuGMTPhiLUT() {}
0045
0046 //--------------
0047 // Destructor --
0048 //--------------
0049 L1MuGMTPhiLUT::~L1MuGMTPhiLUT() {}
0050
0051 //--------------
0052 // Operations --
0053 //--------------
0054
0055 int L1MuGMTPhiLUT::etabin(float eta, int isys) {
0056   int i;
0057   for (i = 0; i < (int)NETA; i++)
0058     if (eta >= etabins[isys][i] && eta < etabins[isys][i + 1])
0059       break;
0060   if (i >= (int)NETA)
0061     edm::LogWarning("LUTProblem") << "L1MuGMTPhiLUT::etabin(): could not assign eta bin ";
0062   return i;
0063 }
0064
0065 //
0066 // project muon phi to calo/vertex
0067 //
0068 float L1MuGMTPhiLUT::dphi(int isys, int isISO, int icharge, int ieta, float pt) {
0069   if ((isys == DT || isys == bRPC) && pt < 4.)
0070     pt = 4.;  // cut off
0071   if ((isys == CSC || isys == fRPC) && pt < 3.)
0072     pt = 3.;
0073
0074   float a = fitparams_phi[isISO][isys][ieta][icharge][0];
0075   float b = fitparams_phi[isISO][isys][ieta][icharge][1];
0076   float c = fitparams_phi[isISO][isys][ieta][icharge][2];
0077
0078   float dphi = a + b / pt + c / (pt * pt);
0079
0080   // dphi is difference between low edge of muon phi bin and real phi in calo/vtx
0081
0082   // correct for RPC LUTs which were calculated with center of bin phi
0083   if (isys == bRPC || isys == fRPC)
0084     dphi -= 1.25 / 180. * M_PI;
0085
0086   return dphi;
0087 }
0088
0089 //
0090 // static parameters of LUTs
0091 //
0092 float L1MuGMTPhiLUT::etabins[L1MuGMTPhiLUT::NSYS][L1MuGMTPhiLUT::NETA + 1] = {
0093     {0.00, 0.22, 0.27, 0.58, 0.77, 0.87, 0.92, 1.24, 1.24},  // DT
0094     {0.00, 1.06, 1.26, 1.46, 1.66, 1.86, 2.06, 2.26, 2.46},  // CSC
0095     {0.00, 0.06, 0.25, 0.41, 0.54, 0.70, 0.83, 0.93, 1.04},  // bRPC
0096     {1.04, 1.24, 1.36, 1.48, 1.61, 1.73, 1.85, 1.97, 2.10}   // fRPC
0097 };
0098
0099 float L1MuGMTPhiLUT::fitparams_phi[L1MuGMTPhiLUT::NRP][L1MuGMTPhiLUT::NSYS][L1MuGMTPhiLUT::NETA][2][3] = {
0100     {// projection to HCAL derived from HCAL positions retrieved in ORCA
0101      // calo phi projection parametrization for DT, projection to HCAL in (for MIP)
0102      {{{-0.001290, 0.504308, 4.549529}, {-0.023984, -0.970760, -2.263805}},
0103       {{-0.001920, 0.489318, 4.356311}, {-0.022833, -1.012401, -2.018399}},
0104       {{-0.000194, 0.552174, 4.758974}, {-0.023465, -1.057848, -2.331967}},
0105       {{-0.002269, 0.604601, 4.030573}, {-0.023426, -1.007778, -2.420637}},
0106       {{-0.000271, 0.533090, 4.488132}, {-0.025969, -0.981604, -2.149771}},
0107       {{0.004051, 0.374228, 4.799831}, {-0.031902, -0.732712, -2.896422}},
0108       {{-0.004666, 0.583151, 3.559875}, {-0.011293, -1.095987, -1.514169}},
0109       {{0.000000, 0.000000, 0.000000}, {0.000000, 0.000000, 0.000000}}},
0110      // calo phi projection parametrization for CSC, projection to HCAL in (for MIP)
0111      {{{-0.001691, 0.581263, 3.408109}, {-0.019783, -1.150412, -0.622612}},
0112       {{-0.002138, 0.554844, 2.839203}, {-0.018853, -0.986130, -1.098399}},
0113       {{0.002930, 0.300343, 2.637024}, {-0.021499, -0.705878, -1.024004}},
0114       {{0.002278, 0.217012, 1.529797}, {-0.019729, -0.581867, -0.343904}},
0115       {{0.017130, 0.166188, 0.849492}, {-0.035404, -0.321889, -0.612119}},
0116       {{0.024403, 0.133149, 0.520598}, {-0.043978, -0.156961, -0.719900}},
0117       {{0.026036, 0.161303, 0.114142}, {-0.045904, -0.096463, -0.652872}},
0118       {{0.024208, 0.146788, 0.016299}, {-0.043546, -0.023951, -0.683871}}},
0119      // calo phi projection parametrization for bRPC, projection to HCAL in (for MIP)
0120      {{{0.015837, 0.665802, 3.613197}, {-0.015709, -0.464655, -4.528647}},
0121       {{0.017926, 0.613610, 3.745701}, {-0.016081, -0.448283, -4.606902}},
0122       {{0.014551, 0.797055, 3.479120}, {-0.014036, -0.612790, -4.167947}},
0123       {{0.018614, 0.637761, 4.156929}, {-0.017071, -0.504211, -4.651026}},
0124       {{0.017844, 0.685002, 3.569384}, {-0.017014, -0.487703, -4.822967}},
0125       {{0.016250, 0.753998, 3.303425}, {-0.017553, -0.538919, -4.559156}},
0126       {{0.019015, 0.653392, 3.993684}, {-0.018937, -0.525990, -4.335586}},
0127       {{0.017044, 0.659718, 3.461996}, {-0.018881, -0.479496, -4.242701}}},
0128      // calo phi projection parametrization for fRPC, projection to HCAL in (for MIP)
0129      {{{0.015301, 0.681718, 2.707917}, {-0.015104, -0.555634, -3.052359}},
0130       {{0.018868, 0.511874, 2.826352}, {-0.015763, -0.549633, -2.301781}},
0131       {{0.022818, 0.267658, 2.445503}, {-0.018882, -0.216642, -2.313509}},
0132       {{0.021048, 0.269100, 1.373464}, {-0.018502, -0.228849, -1.458785}},
0133       {{0.019976, 0.269698, 0.896744}, {-0.021801, -0.227103, -1.032946}},
0134       {{0.039155, 0.112461, 0.965459}, {-0.043126, -0.065388, -1.062142}},
0135       {{0.040538, 0.146610, 0.512473}, {-0.045719, -0.079835, -0.681237}},
0136       {{0.040019, 0.120561, 0.511092}, {-0.044669, -0.031321, -0.745264}}}},
0137     {// calo phi projection parametrization for DT, projection to vertex (for ISO)
0138      {{{-0.019085, 2.530752, 3.112250}, {-0.023016, -2.590798, -2.725504}},
0139       {{-0.018956, 2.533785, 2.969394}, {-0.023637, -2.568775, -2.964180}},
0140       {{-0.018869, 2.533260, 2.940186}, {-0.022805, -2.598547, -2.607942}},
0141       {{-0.019732, 2.527371, 2.611957}, {-0.023624, -2.510180, -2.870058}},
0142       {{-0.019042, 2.477538, 2.628565}, {-0.022541, -2.560252, -1.989453}},
0143       {{-0.018623, 2.408255, 2.631045}, {-0.022602, -2.486282, -2.123144}},
0144       {{-0.019051, 2.424498, 2.422834}, {-0.015844, -2.610096, -2.047908}},
0145       {{0.000000, 0.000000, 0.000000}, {0.000000, 0.000000, 0.000000}}},
0146      // calo phi projection parametrization for CSC, projection to vertex (for ISO)
0147      {{{-0.017973, 2.508663, 1.548810}, {-0.020354, -2.641325, -0.938284}},
0148       {{-0.018821, 2.346581, 1.584596}, {-0.018082, -2.470272, -1.086931}},
0149       {{-0.016044, 1.919469, 1.765014}, {-0.021392, -2.039066, -1.013819}},
0150       {{-0.013416, 1.526508, 1.576988}, {-0.020571, -1.715829, -0.641736}},
0151       {{-0.011872, 1.200872, 1.248002}, {-0.019462, -1.444363, -0.276027}},
0152       {{-0.009545, 0.882696, 1.369971}, {-0.019899, -1.152519, -0.299660}},
0153       {{-0.007845, 0.657952, 1.103447}, {-0.019512, -0.927346, -0.260993}},
0154       {{-0.006395, 0.478389, 1.283519}, {-0.018998, -0.796310, -0.126124}}},
0155      // calo phi projection parametrization for bRPC, projection to vertex (for ISO)
0156      {{{0.002479, 2.626429, 2.457669}, {0.000336, -2.438207, -3.453805}},
0157       {{0.002440, 2.639809, 2.115057}, {0.000804, -2.486881, -3.044556}},
0158       {{0.002161, 2.641236, 2.264882}, {0.001000, -2.491077, -2.856622}},
0159       {{0.001901, 2.637712, 2.031976}, {-0.001031, -2.410168, -3.299684}},
0160       {{0.002487, 2.595609, 2.105496}, {0.000474, -2.422474, -3.296519}},
0161       {{0.001624, 2.614643, 2.065672}, {-0.000492, -2.468336, -2.884312}},
0162       {{0.002451, 2.566535, 2.097821}, {-0.002069, -2.451065, -2.481968}},
0163       {{0.002541, 2.517676, 2.129723}, {0.000255, -2.467283, -2.342318}}},
0164      // calo phi projection parametrization for fRPC, projection to vertex (for ISO)
0165      {{{0.000191, 2.586133, 1.009632}, {0.002492, -2.417437, -1.790964}},
0166       {{0.000549, 2.282987, 0.992959}, {0.003267, -2.192674, -1.343063}},
0167       {{0.000117, 2.081965, 0.577031}, {0.003213, -1.933657, -1.071969}},
0168       {{-0.001817, 1.899919, 0.061467}, {0.003783, -1.746777, -0.644316}},
0169       {{-0.001978, 1.613175, 0.304733}, {-0.000203, -1.567548, -0.433546}},
0170       {{-0.003000, 1.437795, 0.192639}, {-0.001441, -1.342958, -0.556960}},
0171       {{-0.003725, 1.309459, -0.159897}, {-0.002232, -1.182838, -0.440219}},
0172       {{-0.005232, 1.179625, -0.153718}, {0.000748, -1.088718, -0.151436}}}}};