EMTFUnpackerTools.cc

Go to the documentation of this file.

#include "EMTFUnpackerTools.h"

namespace l1t {
  namespace stage2 {
    namespace emtf {

      void ImportME(EMTFHit& _hit, const l1t::emtf::ME _ME, const int _endcap, const int _evt_sector) {
        _hit.set_endcap(_endcap == 1 ? 1 : -1);
        _hit.set_sector_idx(_endcap == 1 ? _evt_sector - 1 : _evt_sector + 5);

        _hit.set_wire(_ME.Wire());
        _hit.set_strip(_ME.Strip());
        _hit.set_quality(_ME.Quality());
        _hit.set_pattern(_ME.CLCT_pattern());
        _hit.set_bend((_ME.LR() == 1) ? 1 : -1);
        _hit.set_valid(_ME.VP());
        _hit.set_sync_err(_ME.SE());
        _hit.set_bx(_ME.TBIN() - 3);
        _hit.set_bc0(_ME.BC0());
        _hit.set_is_CSC(true);
        _hit.set_is_RPC(false);
        _hit.set_subsystem(1);
        // _hit.set_layer();

        _hit.set_ring(L1TMuonEndCap::calc_ring(_hit.Station(), _hit.CSC_ID(), _hit.Strip()));
        _hit.set_chamber(
            L1TMuonEndCap::calc_chamber(_hit.Station(), _hit.Sector(), _hit.Subsector(), _hit.Ring(), _hit.CSC_ID()));

        _hit.SetCSCDetId(_hit.CreateCSCDetId());
        //_hit.SetCSCLCTDigi ( _hit.CreateCSCCorrelatedLCTDigi() );

        // Station, CSC_ID, Sector, Subsector, and Neighbor filled in
        // EventFilter/L1TRawToDigi/src/implementations_stage2/EMTFBlockME.cc
        // "set_layer()" is not invoked, so Layer is not yet filled - AWB 21.04.16

      }  // End ImportME

      void ImportRPC(EMTFHit& _hit, const l1t::emtf::RPC _RPC, const int _endcap, const int _evt_sector) {
        _hit.set_endcap(_endcap == 1 ? 1 : -1);
        _hit.set_sector_idx(_endcap == 1 ? _evt_sector - 1 : _evt_sector + 5);

        _hit.set_phi_fp(_RPC.Phi() * 4);      // 1/4th the precision of CSC LCTs
        _hit.set_theta_fp(_RPC.Theta() * 4);  // 1/4th the precision of CSC LCTs
        _hit.set_bx(_RPC.TBIN() - 3);
        _hit.set_valid(_RPC.VP());
        _hit.set_bc0(_RPC.BC0());
        _hit.set_is_CSC(false);
        _hit.set_is_RPC(true);
        _hit.set_subsystem(2);

        _hit.SetRPCDetId(_hit.CreateRPCDetId());
        // // Not yet implemented - AWB 15.03.17
        // _hit.SetRPCDigi  ( _hit.CreateRPCDigi() );

        // Convert integer values to degrees
        _hit.set_phi_loc(L1TMuonEndCap::calc_phi_loc_deg(_hit.Phi_fp()));
        _hit.set_phi_glob(L1TMuonEndCap::calc_phi_glob_deg(_hit.Phi_loc(), _evt_sector));
        _hit.set_theta(L1TMuonEndCap::calc_theta_deg_from_int(_hit.Theta_fp()));
        _hit.set_eta(L1TMuonEndCap::calc_eta_from_theta_deg(_hit.Theta(), _hit.Endcap()));

        // Station, Ring, Sector, Subsector, Neighbor, and PC/FS/BT_segment filled in
        // EventFilter/L1TRawToDigi/src/implementations_stage2/EMTFBlockRPC.cc - AWB 02.05.17

      }  // End ImportRPC

      void ImportSP(EMTFTrack& _track, const l1t::emtf::SP _SP, const int _endcap, const int _evt_sector) {
        _track.set_endcap((_endcap == 1) ? 1 : -1);
        _track.set_sector(_evt_sector);
        _track.set_sector_idx((_endcap == 1) ? _evt_sector - 1 : _evt_sector + 5);
        _track.set_mode(_SP.Mode());
        _track.set_mode_inv((((_SP.Mode() >> 0) & 1) << 3) | (((_SP.Mode() >> 1) & 1) << 2) |
                            (((_SP.Mode() >> 2) & 1) << 1) | (((_SP.Mode() >> 3) & 1) << 0));
        _track.set_charge((_SP.C() == 1) ? -1 : 1);  // uGMT uses opposite of physical charge (to match pdgID)
        _track.set_bx(_SP.TBIN() - 3);
        _track.set_phi_fp(_SP.Phi_full());
        _track.set_phi_loc(L1TMuonEndCap::calc_phi_loc_deg(_SP.Phi_full()));
        _track.set_phi_glob(L1TMuonEndCap::calc_phi_glob_deg(_track.Phi_loc(), _track.Sector()));
        _track.set_eta(L1TMuonEndCap::calc_eta(_SP.Eta_GMT()));
        _track.set_pt((_SP.Pt_GMT() - 1) * 0.5);

        _track.set_gmt_pt(_SP.Pt_GMT());
        _track.set_gmt_phi(_SP.Phi_GMT());
        _track.set_gmt_eta(_SP.Eta_GMT());
        _track.set_gmt_quality(_SP.Quality_GMT());
        _track.set_gmt_charge(_SP.C());
        _track.set_gmt_charge_valid(_SP.VC());

        EMTFPtLUT _lut = {};
        _lut.address = _SP.Pt_LUT_addr();
        _track.set_PtLUT(_lut);

        // First_bx, Second_bx, Track_num, Has_neighbor, All_neighbor, and Hits should be filled in
        // EventFilter/L1TRawToDigi/src/implementations_stage2/EMTFBlockSP.cc - AWB 07.03.17

      }  // End ImportSP

    }  // End namespace emtf
  }    // End namespace stage2
}  // End namespace l1t