MicroGMTConverter.cc

Go to the documentation of this file.

#include "L1Trigger/L1TMuonEndCap/interface/MicroGMTConverter.h"

MicroGMTConverter::MicroGMTConverter() {}

MicroGMTConverter::~MicroGMTConverter() {}

void MicroGMTConverter::convert(const int global_event_BX,
                                const EMTFTrack& in_track,
                                l1t::RegionalMuonCand& out_cand) const {
  l1t::tftype tftype = (in_track.Endcap() == 1) ? l1t::tftype::emtf_pos : l1t::tftype::emtf_neg;
  int sector = in_track.Sector() - 1;

  out_cand.setHwPt(in_track.GMT_pt());
  out_cand.setHwPtUnconstrained(in_track.GMT_pt_dxy());
  out_cand.setHwDXY(in_track.GMT_dxy());
  out_cand.setHwPhi(in_track.GMT_phi());
  out_cand.setHwEta(in_track.GMT_eta());
  out_cand.setHwSign(in_track.GMT_charge());
  out_cand.setHwSignValid(in_track.GMT_charge_valid());
  out_cand.setHwQual(in_track.GMT_quality());
  out_cand.setHwHF(false);  // EMTF: halo -> 1
  out_cand.setTFIdentifiers(sector, tftype);

  const EMTFPtLUT& ptlut_data = in_track.PtLUT();

  // Form track sub addresses
  int me1_ch_id =
      (ptlut_data.bt_vi[0] == 0 && ptlut_data.bt_vi[1] != 0) ? ptlut_data.bt_ci[1] + 16 : ptlut_data.bt_ci[0];
  int me2_ch_id = ptlut_data.bt_ci[2];
  int me3_ch_id = ptlut_data.bt_ci[3];
  int me4_ch_id = ptlut_data.bt_ci[4];

  int me1_seg_id = (ptlut_data.bt_vi[0] == 0 && ptlut_data.bt_vi[1] != 0) ? ptlut_data.bt_si[1] : ptlut_data.bt_si[0];
  int me2_seg_id = ptlut_data.bt_si[2];
  int me3_seg_id = ptlut_data.bt_si[3];
  int me4_seg_id = ptlut_data.bt_si[4];

  auto get_gmt_chamber_me1 = [](int ch) {
    int gmt_ch = 0;
    if (ch == 10)
      gmt_ch = 1;
    else if (ch == 11)
      gmt_ch = 2;
    else if (ch == 12)
      gmt_ch = 3;
    else if (ch == 3 + 16)
      gmt_ch = 4;
    else if (ch == 6 + 16)
      gmt_ch = 5;
    else if (ch == 9 + 16)
      gmt_ch = 6;
    return gmt_ch;
  };

  auto get_gmt_chamber = [](int ch) {
    int gmt_ch = 0;
    if (ch == 10)
      gmt_ch = 1;
    else if (ch == 11)
      gmt_ch = 2;
    else if (ch == 3)
      gmt_ch = 3;
    else if (ch == 9)
      gmt_ch = 4;
    return gmt_ch;
  };

  // Truncate kBX to 11 bits and offset by 25 from global event BX in EMTF firmware
  int EMTF_kBX = (abs(global_event_BX) % 2048) - 25 + in_track.BX();
  if (EMTF_kBX < 0)
    EMTF_kBX += 2048;

  out_cand.setTrackSubAddress(l1t::RegionalMuonCand::kME1Seg, me1_seg_id);
  out_cand.setTrackSubAddress(l1t::RegionalMuonCand::kME1Ch, get_gmt_chamber_me1(me1_ch_id));
  out_cand.setTrackSubAddress(l1t::RegionalMuonCand::kME2Seg, me2_seg_id);
  out_cand.setTrackSubAddress(l1t::RegionalMuonCand::kME2Ch, get_gmt_chamber(me2_ch_id));
  out_cand.setTrackSubAddress(l1t::RegionalMuonCand::kME3Seg, me3_seg_id);
  out_cand.setTrackSubAddress(l1t::RegionalMuonCand::kME3Ch, get_gmt_chamber(me3_ch_id));
  out_cand.setTrackSubAddress(l1t::RegionalMuonCand::kME4Seg, me4_seg_id);
  out_cand.setTrackSubAddress(l1t::RegionalMuonCand::kME4Ch, get_gmt_chamber(me4_ch_id));
  out_cand.setTrackSubAddress(l1t::RegionalMuonCand::kTrkNum, in_track.Track_num());
  out_cand.setTrackSubAddress(l1t::RegionalMuonCand::kBX, EMTF_kBX);
}

void MicroGMTConverter::convert_all(const edm::Event& iEvent,
                                    const EMTFTrackCollection& in_tracks,
                                    l1t::RegionalMuonCandBxCollection& out_cands) const {
  int gmtMinBX = -2;
  int gmtMaxBX = +2;

  out_cands.clear();
  out_cands.setBXRange(gmtMinBX, gmtMaxBX);

  for (const auto& in_track : in_tracks) {
    int bx = in_track.BX();

    if (gmtMinBX <= bx && bx <= gmtMaxBX) {
      l1t::RegionalMuonCand out_cand;

      convert(iEvent.bunchCrossing(), in_track, out_cand);
      out_cands.push_back(bx, out_cand);
    }
  }

  // Sort by processor to match uGMT unpacked order
  emtf::sort_uGMT_muons(out_cands);
}

namespace emtf {

  void sort_uGMT_muons(l1t::RegionalMuonCandBxCollection& cands) {
    int minBX = cands.getFirstBX();
    int maxBX = cands.getLastBX();
    int emtfMinProc = 0;   // ME+ sector 1
    int emtfMaxProc = 11;  // ME- sector 6

    // New collection, sorted by processor to match uGMT unpacked order
    auto sortedCands = std::make_unique<l1t::RegionalMuonCandBxCollection>();
    sortedCands->clear();
    sortedCands->setBXRange(minBX, maxBX);
    for (int iBX = minBX; iBX <= maxBX; ++iBX) {
      for (int proc = emtfMinProc; proc <= emtfMaxProc; proc++) {
        for (l1t::RegionalMuonCandBxCollection::const_iterator cand = cands.begin(iBX); cand != cands.end(iBX);
             ++cand) {
          int cand_proc = cand->processor();
          if (cand->trackFinderType() == l1t::tftype::emtf_neg)
            cand_proc += 6;
          if (cand_proc != proc)
            continue;
          sortedCands->push_back(iBX, *cand);
        }
      }
    }

    // Return sorted collection
    std::swap(cands, *sortedCands);
    sortedCands.reset();
  }

}  // End namespace emtf