Constants.h

#ifndef L1TMUONPHASE2_PHASE2GMT_CONSTANTS
#define L1TMUONPHASE2_PHASE2GMT_CONSTANTS

#include "ap_int.h"

namespace Phase2L1GMT {

  // INPUT Dataformat
  // Track
  const int BITSTTCURV = 15;
  const int BITSTTCURV2 = 27;
  const int BITSTTPHI = 12;
  const int BITSTTTANL = 16;
  // Not used in the emulator now, but will in future
  const int BITSTTZ0 = 12;
  const int BITSTTD0 = 13;
  const int BITSTTCHI2 = 4;
  const int BITSTTBENDCHI2 = 3;
  const int BITSTTMASK = 7;
  const int BITSTTTRACKMVA = 3;
  const int BITSTTOTHERMVA = 6;

  // Bitwidth for common internal variables in GMT
  const int BITSPT = 13;
  const int BITSPHI = 13;
  const int BITSETA = 13;
  const int BITSZ0 = 10;
  const int BITSD0 = 12;

  //Muon ROI
  const int BITSSTUBCOORD = 10;
  const int BITSSTUBETA = 8;
  const int BITSSTUBID = 12;
  const int BITSSTUBPHIQUALITY = 4;
  const int BITSSTUBETAQUALITY = 4;
  const int BITSSTUBTIME = 8;
  const int BITSSTAMUONQUALITY = 6;
  const int BITSMUONBX = 4;

  //PreTrackMatherdMuon
  const int BITSMATCHQUALITY = 9;
  const int BITSMUONBETA = 4;

  //Track Muon Match
  const int BITSSIGMAETA = 4;
  const int BITSSIGMACOORD = 4;
  const int BITSPROPCOORD = 9;
  const int BITSPROPSIGMACOORD_A = 6;
  const int BITSPROPSIGMACOORD_B = 6;
  const int BITSPROPSIGMAETA_A = 5;
  const int BITSPROPSIGMAETA_B = 5;

  // OUTPUT Dataformat
  // Bitwidth for standalone muons to CL1 and GT
  const int BITSSAZ0 = 5;
  const int BITSSAD0 = 7;
  const int BITSSAQUAL = 4;

  // Bitwidth for dataformat to GT
  const int BITSGTPT = 16;
  const int BITSGTPHI = 13;
  const int BITSGTETA = 14;
  const int BITSGTZ0 = 10;
  const int BITSGTD0 = 10;
  const int BITSGTQUAL = 6;
  const int BITSGTQUALSC = 7;
  const int BITSGTISO = 6;
  const int BITSGTBETA = 4;

  // Bitwidth for Tau->3mu object
  const int BITSTMPT = 8;
  const int BITSTMPHI = 8;
  const int BITSTMETA = 8;
  const int BITSTMMASS2 = 8;
  const int BITSTMTYPE = 6;
  const int BITSTMIDX = 4;
  const int BITSTMQUAL = 4;

  const float maxCurv_ = 0.00855;  // 2 GeV pT Rinv is in cm
  const float maxPhi_ = 1.026;     // relative to the center of the sector
  const float maxTanl_ = 8.0;
  const float maxZ0_ = 25.6;
  const float maxD0_ = 15.36;
  // Updated barrelLimit according to Karol, https://indico.cern.ch/event/1113802/#1-phase2-gmt-performance-and-i
  const int barrelLimit0_ = 1.4 / 0.00076699039 / 8;
  const int barrelLimit1_ = 1.1 / 0.00076699039 / 8;
  const int barrelLimit2_ = 0.95 / 0.00076699039 / 8;
  const int barrelLimit3_ = 0.95 / 0.00076699039 / 8;
  const int barrelLimit4_ = 0;

  // LSB
  const float LSBpt = 0.03125;
  const float LSBphi = 2. * M_PI / pow(2, BITSPHI);
  const float LSBeta = 2. * M_PI / pow(2, BITSETA);
  const float LSBGTiso = 0.25;  // 0.25GeV, LSBphi * 8
  const float LSBGTz0 = 0.05;   // 0.5mm, in sync with GTT and Correlator
  const float LSBGTd0 = 0.03;   // from GT interface doc
  const float LSBSAz0 = 1.6;    // 0.05 * 32 cm, with range +- 25.6
  const float LSBSAd0 = 3.84;   // 0.03 * 128 cm, with range +- 245.76

  typedef ap_uint<64> wordtype;
  typedef ap_uint<1> valid_gt_t;             //valid
  typedef ap_uint<1> q_gt_t;                 //charge
  typedef ap_uint<BITSGTPT> pt_gt_t;         //pt        of tracker muon
  typedef ap_int<BITSGTPHI> phi_gt_t;        //phi       of tracker muon
  typedef ap_int<BITSGTETA> eta_gt_t;        //eta       of tracker muon
  typedef ap_int<BITSGTZ0> z0_gt_t;          //z0        of tracker muon
  typedef ap_int<BITSGTD0> d0_gt_t;          //d0        of tracker muon
  typedef ap_uint<BITSGTISO> iso_gt_t;       //isolation of tracker muon
  typedef ap_uint<BITSGTBETA> beta_gt_t;     //beta      of tracker muon
  typedef ap_uint<BITSGTQUAL> qual_gt_t;     //quality   of tracker muon
  typedef ap_int<BITSGTQUALSC> qualsc_gt_t;  //quality score of tracker muon (qual/q)

  //Standalone muon datatype
  typedef ap_uint<1> valid_sa_t;          //valid
  typedef ap_uint<BITSGTPT> pt_sa_t;      //pt      of standalone muon
  typedef ap_int<BITSGTPHI> phi_sa_t;     //phi     of standalone muon
  typedef ap_int<BITSGTETA> eta_sa_t;     //eta     of standalone muon
  typedef ap_int<BITSSAZ0> z0_sa_t;       //z0      of standalone muon
  typedef ap_int<BITSSAD0> d0_sa_t;       //d0      of standalone muon
  typedef ap_uint<1> q_sa_t;              //charge  of standalone muon
  typedef ap_uint<BITSSAQUAL> qual_sa_t;  //quality of standalone muon

  inline uint64_t twos_complement(long long int v, uint bits) {
    uint64_t mask = (1 << bits) - 1;
    if (v >= 0)
      return v & mask;
    else
      return (~(-v) + 1) & mask;
  }

  template <class T>
  inline int wordconcat(T& word, int bstart, long int input, int bitsize) {
    int bend = bstart + bitsize - 1;
    word.range(bend, bstart) = twos_complement(input, bitsize);
    return bend + 1;
  }

}  // namespace Phase2L1GMT
#endif