TTPrimitiveConversion.cc

Go to the documentation of this file.

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

#include "L1Trigger/L1TMuonEndCap/interface/SectorProcessorLUT.h"
#include "L1Trigger/L1TMuonEndCap/interface/TrackTools.h"


void TTPrimitiveConversion::configure(
    const TTGeometryTranslator* tp_ttgeom,
    const SectorProcessorLUT* lut,
    int verbose, int endcap, int sector, int bx
) {
  assert(tp_ttgeom != nullptr);
  assert(lut != nullptr);

  tp_ttgeom_ = tp_ttgeom;
  lut_       = lut;  // not used

  verbose_ = verbose;
  endcap_  = endcap; // 1 for ME+, 2 for ME-
  sector_  = sector;
  bx_      = bx;
}

void TTPrimitiveConversion::process(
    const std::map<int, TTTriggerPrimitiveCollection>& selected_ttprim_map,
    EMTFHitCollection& conv_hits
) const {

  for (const auto& map_tp_it : selected_ttprim_map) {
    for (const auto& tp_it : map_tp_it.second) {
      EMTFHit conv_hit;
      convert_tt(tp_it, conv_hit);
      conv_hits.push_back(conv_hit);
    }
  }
}

void TTPrimitiveConversion::process_no_prim_sel(
    const TTTriggerPrimitiveCollection& ttmuon_primitives,
    EMTFHitCollection& conv_hits
) const {

  for (const auto& tp_it : ttmuon_primitives) {
    if (endcap_ == 1 && sector_ == 1 && bx_ == tp_it.getTTData().bx) {  //FIXME: stupidly put everything into sector +1, to be fixed.
      EMTFHit conv_hit;
      convert_tt(tp_it, conv_hit);
      conv_hits.push_back(conv_hit);
    }
  }
}


// _____________________________________________________________________________
// TT functions
void TTPrimitiveConversion::convert_tt(
    const TTTriggerPrimitive& ttmuon_primitive,
    EMTFHit& conv_hit
) const {
  //const DetId&   tp_detId = ttmuon_primitive.detId();
  const TTData&  tp_data  = ttmuon_primitive.getTTData();

  int tp_region    = tp_ttgeom_->region(ttmuon_primitive);     // 0 for Barrel, +/-1 for +/- Endcap
  int tp_endcap    = (tp_region == -1) ? 2 : tp_region;
  int tp_station   = tp_ttgeom_->layer(ttmuon_primitive);
  int tp_ring      = tp_ttgeom_->ring(ttmuon_primitive);
  int tp_chamber   = tp_ttgeom_->module(ttmuon_primitive);
  int tp_sector    = 1;  //FIXME
  int tp_subsector = 0;  //FIXME

  const bool is_neighbor = false;  //FIXME

  // Set properties
  //conv_hit.SetTTDetId      ( tp_detId );

  conv_hit.set_endcap      ( (tp_endcap == 2) ? -1 : tp_endcap );
  conv_hit.set_station     ( tp_station );
  conv_hit.set_ring        ( tp_ring );
  //conv_hit.set_roll        ( tp_roll );
  conv_hit.set_chamber     ( tp_chamber );
  conv_hit.set_sector      ( tp_sector );
  conv_hit.set_subsector   ( tp_subsector );
  //conv_hit.set_csc_ID      ( tp_csc_ID );
  //conv_hit.set_csc_nID     ( csc_nID );
  //conv_hit.set_track_num   ( tp_data.trknmb );
  //conv_hit.set_sync_err    ( tp_data.syncErr );

  conv_hit.set_bx          ( tp_data.bx );
  conv_hit.set_subsystem   ( TTTriggerPrimitive::kTT );
  conv_hit.set_is_CSC      ( false );
  conv_hit.set_is_RPC      ( false );
  conv_hit.set_is_GEM      ( false );

  //conv_hit.set_pc_sector   ( pc_sector );
  //conv_hit.set_pc_station  ( pc_station );
  //conv_hit.set_pc_chamber  ( pc_chamber );
  //conv_hit.set_pc_segment  ( pc_segment );

  conv_hit.set_valid       ( true );
  conv_hit.set_strip       ( static_cast<int>(tp_data.row_f) );
  //conv_hit.set_strip_low   ( tp_data.strip_low );
  //conv_hit.set_strip_hi    ( tp_data.strip_hi );
  conv_hit.set_wire        ( static_cast<int>(tp_data.col_f) );
  //conv_hit.set_quality     ( tp_data.quality );
  //conv_hit.set_pattern     ( tp_data.pattern );
  conv_hit.set_bend        ( tp_data.bend );
  //conv_hit.set_time        ( tp_data.time );

  conv_hit.set_neighbor    ( is_neighbor );
  conv_hit.set_sector_idx  ( (endcap_ == 1) ? sector_ - 1 : sector_ + 5 );

  // Add coordinates from fullsim
  {
    const GlobalPoint& gp = tp_ttgeom_->getGlobalPoint(ttmuon_primitive);
    double glob_phi   = emtf::rad_to_deg(gp.phi().value());
    double glob_theta = emtf::rad_to_deg(gp.theta());
    double glob_eta   = gp.eta();
    double glob_rho   = gp.perp();
    double glob_z     = gp.z();

    conv_hit.set_phi_sim   ( glob_phi );
    conv_hit.set_theta_sim ( glob_theta );
    conv_hit.set_eta_sim   ( glob_eta );
    conv_hit.set_rho_sim   ( glob_rho );
    conv_hit.set_z_sim     ( glob_z );
  }
}