emtf::phase2::CSCTPConverter Class Reference

#include <CSCTPConverter.h>

Inheritance diagram for emtf::phase2::CSCTPConverter:

Public Member Functions
void	convert (const TriggerPrimitive &, const TPInfo &, EMTFHit &) const final
	CSCTPConverter (const EMTFContext &, const int &, const int &)
	~CSCTPConverter () override=default
Public Member Functions inherited from emtf::phase2::TPConverter
	TPConverter ()=default
virtual	~TPConverter ()=default

Private Attributes
const EMTFContext &	context_
int	endcap_
int	sector_

Detailed Description

Definition at line 13 of file CSCTPConverter.h.

Constructor & Destructor Documentation

CSCTPConverter()

CSCTPConverter::CSCTPConverter ( const EMTFContext &  context, const int &  endcap, const int &  sector  )

explicit

Definition at line 18 of file CSCTPConverter.cc.

    : context_(context), endcap_(endcap), sector_(sector) {}

~CSCTPConverter()

emtf::phase2::CSCTPConverter::~CSCTPConverter ( )

overridedefault

Member Function Documentation

convert()

void CSCTPConverter::convert ( const TriggerPrimitive &  tp, const TPInfo &  tp_info, EMTFHit &  hit  ) const

finalvirtual

Implements emtf::phase2::TPConverter.

Definition at line 21 of file CSCTPConverter.cc.

References emtf::phase2::TPInfo::bx, emtf::phase2::tp::calcPhiInt(), emtf::phase2::tp::calcThetaInt(), emtf::phase2::TPInfo::chamber, ALPAKA_ACCELERATOR_NAMESPACE::brokenline::constexpr(), context_, emtf::phase2::TPInfo::csc_facing, emtf::phase2::TPInfo::csc_first_wire, emtf::phase2::TPInfo::csc_id, emtf::phase2::TPInfo::csc_second_wire, emtf_assert, emtf::phase2::TPInfo::endcap_pm, emtf::phase2::TPInfo::flag_substitute, emtf::phase2::EMTFContext::geometry_translator_, L1TMuon::GeometryTranslator::getGlobalPoint(), emtf::phase2::data::ZoneLut::getZones(), emtf::phase2::TPInfo::hit_id, emtf::phase2::EMTFContext::host_lut_, emtf::phase2::TPInfo::ilink, L1TMuon::kCSC, emtf::phase2::kNeighbor, emtf::phase2::csc::kRear, emtf::phase2::TPInfo::layer, emtf::phase2::data::HostLut::lookup(), emtf::phase2::data::SiteLut::lookup(), PV3DBase< T, PVType, FrameType >::perp(), PV3DBase< T, PVType, FrameType >::phi(), emtf::phase2::tp::radToDeg(), emtf::phase2::TPInfo::ring, sector_, emtf::phase2::TPInfo::segment_id, emtf::phase2::TPInfo::selection, emtf::phase2::EMTFContext::site_lut_, emtf::phase2::TPInfo::station, emtf::phase2::TPInfo::subsector, PV3DBase< T, PVType, FrameType >::theta(), cmsswSequenceInfo::tp, Geom::Phi< T1, Range >::value(), PV3DBase< T, PVType, FrameType >::z(), and emtf::phase2::EMTFContext::zone_lut_.

                                                                                                  {
  // Unpack Id
  const auto& tp_hit_id = tp_info.hit_id;
  const auto& tp_segment_id = tp_info.segment_id;

  // Unpack trigger primitive
  const auto& tp_det_id = tp.detId<CSCDetId>();
  const auto& tp_data = tp.getCSCData();

  // Unpack detector info
  const auto tp_subsystem = L1TMuon::kCSC;

  const int tp_raw_id = tp_det_id.rawId();

  const int tp_endcap_pm = tp_info.endcap_pm;
  const int tp_subsector = tp_info.subsector;
  const int tp_chamber = tp_info.chamber;
  const int tp_station = tp_info.station;
  const int tp_ring = tp_info.ring;
  const int tp_layer = tp_info.layer;

  const int tp_csc_id = tp_info.csc_id;
  const auto tp_csc_facing = tp_info.csc_facing;

  // Unpack data
  const int tp_strip = tp_data.strip;
  const int tp_strip_quart_bit = tp_data.strip_quart_bit;
  const int tp_strip_eighth_bit = tp_data.strip_eighth_bit;
  const int tp_strip_quart = tp_data.strip_quart;
  const int tp_strip_eighth = tp_data.strip_eighth;

  const int tp_wire1 = tp_info.csc_first_wire;
  const int tp_wire2 = tp_info.csc_second_wire;

  int tp_bend;
  const int tp_slope = tp_data.slope;

  const int tp_bx = tp_info.bx;
  const int tp_subbx = 0;   // no fine resolution timing
  const float tp_time = 0;  // no fine resolution timing. Note: Check with Efe, Jia Fu gets this from digi directly.

  const auto tp_selection = tp_info.selection;

  const int tp_pattern = tp_data.pattern;
  const int tp_quality = 6;

  // Apply CSC Run 2 pattern -> bend conversion
  // Override tp_bend
  constexpr int tp_bend_lut_size = 11;
  constexpr int tp_bend_lut[tp_bend_lut_size] = {-5, 5, -4, 4, -3, 3, -2, 2, -1, 1, 0};
  emtf_assert(tp_pattern < tp_bend_lut_size);
  tp_bend = tp_bend_lut[tp_pattern];
  tp_bend *= tp_endcap_pm;  // sign flip depending on endcap

  // ID scheme used in FW
  const int tp_ilink = tp_info.ilink;

  // Get Global Coordinates
  const GlobalPoint& gp_w1 = this->context_.geometry_translator_.getGlobalPoint(tp);
  const float glob_phi_w1 = tp::radToDeg(gp_w1.phi().value());
  const float glob_theta_w1 = tp::radToDeg(gp_w1.theta().value());
  const double glob_rho_w1 = gp_w1.perp();
  const double glob_z_w1 = gp_w1.z();

  // Calculate EMTF Values
  const int emtf_phi_w1 = tp::calcPhiInt(sector_, glob_phi_w1);
  const int emtf_bend_w1 = std::clamp(tp_bend * 4, -16, 15);  // 5-bit, signed
  const int emtf_theta_w1 = tp::calcThetaInt(tp_endcap_pm, glob_theta_w1);
  const int emtf_qual_w1 = std::clamp(tp_quality, 0, 15);  // 4-bit, unsigned
  const int emtf_site_w1 = context_.site_lut_.lookup({tp_subsystem, tp_station, tp_ring});
  const int emtf_host_w1 = context_.host_lut_.lookup({tp_subsystem, tp_station, tp_ring});
  const int emtf_zones_w1 = context_.zone_lut_.getZones(emtf_host_w1, emtf_theta_w1);

  // Calculated Ambiguous Info
  int emtf_theta_w2 = 0;
  int emtf_qual_w2 = tp_pattern;

  if (tp_wire2 > -1) {
    auto tp_w2 = tp;

    tp_w2.accessCSCData().keywire = tp_wire2;

    const GlobalPoint& gp_w2 = this->context_.geometry_translator_.getGlobalPoint(tp_w2);
    const double glob_theta_w2 = tp::radToDeg(gp_w2.theta().value());

    emtf_theta_w2 = tp::calcThetaInt(tp_endcap_pm, glob_theta_w2);
  }

  emtf_assert((0 <= emtf_phi_w1) and (emtf_phi_w1 < 5040));
  emtf_assert((1 <= emtf_theta_w1) and (emtf_theta_w1 < 128));
  emtf_assert((0 <= emtf_theta_w2) and (emtf_theta_w2 < 128));

  // Get flags
  const bool tp_flag_neighbor = (tp_selection == TPSelection::kNeighbor);
  const bool tp_flag_substitute = tp_info.flag_substitute;
  const bool tp_flag_valid = tp_data.valid;

  // Set properties
  hit.setId(tp_hit_id);

  hit.setRawDetId(tp_raw_id);
  hit.setSubsystem(L1TMuon::kCSC);
  hit.setEndcap(tp_endcap_pm);
  hit.setSector(sector_);
  hit.setSubsector(tp_subsector);
  hit.setStation(tp_station);
  hit.setRing(tp_ring);
  hit.setLayer(tp_layer);
  hit.setChamber(tp_chamber);

  hit.setCscId(tp_csc_id);
  hit.setCscFR(tp_csc_facing == csc::Facing::kRear);

  hit.setStrip(tp_strip);
  hit.setStripLo(tp_strip);
  hit.setStripHi(tp_strip);
  hit.setStripQuart(tp_strip_quart);
  hit.setStripEighth(tp_strip_eighth);
  hit.setStripQuartBit(tp_strip_quart_bit);
  hit.setStripEighthBit(tp_strip_eighth_bit);

  hit.setWire1(tp_wire1);
  hit.setWire2(tp_wire2);

  hit.setBend(tp_bend);
  hit.setSlope(tp_slope);

  hit.setBx(tp_bx);
  hit.setSubbx(tp_subbx);

  hit.setQuality(tp_quality);
  hit.setPattern(tp_pattern);

  hit.setGlobPhi(glob_phi_w1);
  hit.setGlobTheta(glob_theta_w1);
  hit.setGlobPerp(glob_rho_w1);
  hit.setGlobZ(glob_z_w1);
  hit.setGlobTime(tp_time);

  hit.setEmtfChamber(tp_ilink);
  hit.setEmtfSegment(tp_segment_id);
  hit.setEmtfPhi(emtf_phi_w1);
  hit.setEmtfBend(emtf_bend_w1);
  hit.setEmtfTheta1(emtf_theta_w1);
  hit.setEmtfTheta2(emtf_theta_w2);
  hit.setEmtfQual1(emtf_qual_w1);
  hit.setEmtfQual2(emtf_qual_w2);
  hit.setEmtfSite(emtf_site_w1);
  hit.setEmtfHost(emtf_host_w1);
  hit.setEmtfZones(emtf_zones_w1);

  hit.setFlagNeighbor(tp_flag_neighbor);
  hit.setFlagSubstitute(tp_flag_substitute);
  hit.setFlagValid(tp_flag_valid);
}

Member Data Documentation

context_

const EMTFContext& emtf::phase2::CSCTPConverter::context_

private

Definition at line 22 of file CSCTPConverter.h.

Referenced by convert().

endcap_

int emtf::phase2::CSCTPConverter::endcap_

private

Definition at line 24 of file CSCTPConverter.h.

sector_

int emtf::phase2::CSCTPConverter::sector_

private

Definition at line 24 of file CSCTPConverter.h.

Referenced by convert().