L1MuTriggerScales.h

Go to the documentation of this file.

//-------------------------------------------------
//
//   \class L1MuTriggerScales
//
//
//   $Date: 2008/11/24 10:26:31 $
//   $Revision: 1.5 $
//
//   Original Author :
//   Hannes Sakulin      HEPHY / Vienna
//
//   Migrated to CMSSW:
//   I. Mikulec
//
//--------------------------------------------------
#ifndef CondFormatsL1TObjects_L1MuTriggerScales_h
#define CondFormatsL1TObjects_L1MuTriggerScales_h

#include "CondFormats/Serialization/interface/Serializable.h"

#include <cmath>
#include <iostream>
#include <vector>

#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "CondFormats/L1TObjects/interface/L1MuScale.h"

class L1MuTriggerScales {
public:
  L1MuTriggerScales() {}

  L1MuTriggerScales(int nbitPackingDTEta,
                    bool signedPackingDTEta,
                    int nbinsDTEta,
                    float minDTEta,
                    float maxDTEta,
                    int offsetDTEta,

                    int nbitPackingCSCEta,
                    int nbinsCSCEta,
                    float minCSCEta,
                    float maxCSCEta,

                    const std::vector<double>& scaleRPCEta,
                    int nbitPackingBrlRPCEta,
                    bool signedPackingBrlRPCEta,
                    int nbinsBrlRPCEta,
                    int offsetBrlRPCEta,
                    int nbitPackingFwdRPCEta,
                    bool signedPackingFwdRPCEta,
                    int nbinsFwdRPCEta,
                    int offsetFwdRPCEta,

                    int nbitPackingGMTEta,
                    int nbinsGMTEta,
                    const std::vector<double>& scaleGMTEta,

                    int nbitPackingPhi,
                    bool signedPackingPhi,
                    int nbinsPhi,
                    float minPhi,
                    float maxPhi

                    /*            int nbitPackingPt, */
                    /*            bool signedPackingPt, */
                    /*            int nbinsPt, */
                    /*            const std::vector<double>& scalePt */

  ) {
    //
    // Regional Muon Trigger Eta Scales
    //
    /*     const float rpcetabins[34]= { */
    /*       -2.10, -1.97, -1.85, -1.73, -1.61, -1.48, */
    /*       -1.36, -1.24, -1.14, -1.04, -0.93, -0.83,  */
    /*       -0.72, -0.58, -0.44, -0.27, -0.07,      */
    /*               0.07,  0.27,  0.44,  0.58,  0.72, */
    /*        0.83,  0.93,  1.04,  1.14,  1.24,  1.36, */
    /*        1.48,  1.61,  1.73,  1.85,  1.97,  2.10}; */

    // DT
    //m_RegionalEtaScale[0] = L1MuBinnedScale( 6, true, 64, -1.2, 1.2, 32);
    m_RegionalEtaScale[0] =
        L1MuBinnedScale(nbitPackingDTEta, signedPackingDTEta, nbinsDTEta, minDTEta, maxDTEta, offsetDTEta);

    // RPC index -16 .. 16, brl RPC
    // m_RegionalEtaScale[1] = L1MuBinnedScale (6, true, 33, rpcetabins, 16);
    m_RegionalEtaScale[1] =
        L1MuBinnedScale(nbitPackingBrlRPCEta, signedPackingBrlRPCEta, nbinsBrlRPCEta, scaleRPCEta, offsetBrlRPCEta);

    // CSC
    m_RegionalEtaScale[2] = L1MuBinnedScale();
    // // m_RegionalEtaScale[2] = L1MuSymmetricBinnedScale ( 6, 32, 0.9, 2.5);
    //    m_RegionalEtaScaleCSC = L1MuSymmetricBinnedScale ( 6, 32, 0.9, 2.5);
    m_RegionalEtaScaleCSC = L1MuSymmetricBinnedScale(nbitPackingCSCEta, nbinsCSCEta, minCSCEta, maxCSCEta);

    // RPC index -16 .. 16, fwd RPC
    // m_RegionalEtaScale[3] = L1MuBinnedScale (6, true, 33, rpcetabins, 16);
    m_RegionalEtaScale[3] =
        L1MuBinnedScale(nbitPackingFwdRPCEta, signedPackingFwdRPCEta, nbinsFwdRPCEta, scaleRPCEta, offsetFwdRPCEta);

    //
    // Eta scale at GMT output
    //

    /*     const float gmt_outputetascale[32] = {  0.00, */
    /*             0.10,  0.20,  0.30,  0.40,  0.50,  0.60,  0.70,  0.80,  */
    /*             0.90,  1.00,  1.10,  1.20,  1.30,  1.40,  1.50,  1.60, */
    /*             1.70,  1.75,  1.80,  1.85,  1.90,  1.95,  2.00,  2.05, */
    /*             2.10,  2.15,  2.20,  2.25,  2.30,  2.35,  2.40 }; */

    // m_GMTEtaScale = L1MuSymmetricBinnedScale (6, 31, gmt_outputetascale);
    m_GMTEtaScale = L1MuSymmetricBinnedScale(nbitPackingGMTEta, nbinsGMTEta, scaleGMTEta);

    //
    // Phi Scale. Common to all Regioanl Muon Triggers and GMT
    //

    // m_PhiScale = L1MuBinnedScale (8, false, 144, 0., 2. * M_PI);
    m_PhiScale = L1MuBinnedScale(nbitPackingPhi, signedPackingPhi, nbinsPhi, minPhi, maxPhi);

    //
    // Pt Scale. Common to all Regioanl Muon Triggers and GMT
    //

    // pt scale in GeV
    // low edges of pt bins
    /*     const float ptscale[33] = {  */
    /*       -1.,   0.0,   1.5,   2.0,   2.5,   3.0,   3.5,   4.0, */
    /*       4.5,   5.0,   6.0,   7.0,   8.0,  10.0,  12.0,  14.0,   */
    /*       16.0,  18.0,  20.0,  25.0,  30.0,  35.0,  40.0,  45.0,  */
    /*       50.0,  60.0,  70.0,  80.0,  90.0, 100.0, 120.0, 140.0, 1.E6 }; */

    // m_PtScale = L1MuBinnedScale ( 5, false, 32, ptscale) ;
    /*     m_PtScale = L1MuBinnedScale ( nbitPackingPt, */
    /*                signedPackingPt, */
    /*                nbinsPt, */
    /*                scalePt ) ; */
  };

  virtual ~L1MuTriggerScales(){
      //     for (int i=0; i<4; i++)
      //       delete m_RegionalEtaScale[i];

      //     delete m_GMTEtaScale;
      //     delete m_PhiScale;
      //     delete m_PtScale;
  };

  const L1MuScale* getRegionalEtaScale(int isys) const {
    if (isys < 0 || isys > 3)
      edm::LogWarning("ScaleRangeViolation") << "Error in L1MuTriggerScales:: isys out of range: " << isys;
    if (isys == 2) {
      return &m_RegionalEtaScaleCSC;
    } else {
      return &(m_RegionalEtaScale[isys]);
    }
  };

  const L1MuScale* getGMTEtaScale() const { return &m_GMTEtaScale; };

  void setGMTEtaScale(const L1MuSymmetricBinnedScale& scale) { m_GMTEtaScale = scale; };

  const L1MuScale* getPhiScale() const { return &m_PhiScale; };

  void setPhiScale(const L1MuBinnedScale& scale) { m_PhiScale = scale; };

  /*   /// get the Pt scale */
  /*   const L1MuScale* getPtScale() const { return &m_PtScale;}; */

private:
  L1MuBinnedScale m_RegionalEtaScale[4];  // 2=csc will be empty
  L1MuSymmetricBinnedScale m_RegionalEtaScaleCSC;
  L1MuSymmetricBinnedScale m_GMTEtaScale;
  L1MuBinnedScale m_PhiScale;
  //  L1MuBinnedScale m_PtScale;

  COND_SERIALIZABLE;
};

#endif