emtf::phase2::RPCTPConverter Class Reference

#include <RPCTPConverter.h>

Inheritance diagram for emtf::phase2::RPCTPConverter:

Public Member Functions
void	convert (const TriggerPrimitive &, const TPInfo &, EMTFHit &) const final
	RPCTPConverter (const EMTFContext &, const int &, const int &)
	~RPCTPConverter () 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 11 of file RPCTPConverter.h.

Constructor & Destructor Documentation

RPCTPConverter()

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

explicit

Definition at line 20 of file RPCTPConverter.cc.

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

~RPCTPConverter()

emtf::phase2::RPCTPConverter::~RPCTPConverter ( )

overridedefault

Member Function Documentation

convert()

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

finalvirtual

Implements emtf::phase2::TPConverter.

Definition at line 23 of file RPCTPConverter.cc.

References funct::abs(), emtf::phase2::TPInfo::bx, emtf::phase2::tp::calcPhiInt(), emtf::phase2::tp::calcThetaInt(), emtf::phase2::TPInfo::chamber, context_, emtf::phase2::TPInfo::csc_facing, emtf::phase2::TPInfo::csc_id, MillePedeFileConverter_cfg::e, emtf_assert, emtf::phase2::TPInfo::endcap_pm, emtf::phase2::TPInfo::flag_substitute, emtf::phase2::EMTFContext::geometry_translator_, L1TMuon::GeometryTranslator::getGlobalPoint(), L1TMuon::GeometryTranslator::getRPCGeometry(), emtf::phase2::data::ZoneLut::getZones(), emtf::phase2::TPInfo::hit_id, emtf::phase2::EMTFContext::host_lut_, emtf::phase2::TPInfo::ilink, emtf::phase2::rpc::kiRPC, emtf::phase2::kNeighbor, emtf::phase2::csc::kRear, L1TMuon::kRPC, emtf::phase2::TPInfo::layer, emtf::phase2::data::HostLut::lookup(), emtf::phase2::data::SiteLut::lookup(), or, PV3DBase< T, PVType, FrameType >::perp(), PV3DBase< T, PVType, FrameType >::phi(), emtf::phase2::tp::radToDeg(), emtf::phase2::TPInfo::ring, emtf::phase2::TPInfo::roll, nano_mu_digi_cff::roll, RPCGeometry::roll(), emtf::phase2::TPInfo::rpc_type, 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<RPCDetId>();
  const auto& tp_data = tp.getRPCData();

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

  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_station = tp_info.station;
  const int tp_ring = tp_info.ring;
  const int tp_roll = tp_info.roll;
  const int tp_layer = tp_info.layer;
  const int tp_chamber = tp_info.chamber;

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

  // Unpack data
  const int tp_strip = (tp_data.strip_low + tp_data.strip_hi) / 2;  // in full-strip unit
  const int tp_strip_lo = tp_data.strip_low;
  const int tp_strip_hi = tp_data.strip_hi;
  const int tp_clus_width = (tp_strip_hi - tp_strip_lo + 1);

  int tp_bend = 0;  // not applicable

  const int tp_bx = tp_info.bx;
  const float tp_time = tp_data.time;
  float tp_subbx_f32 = tp_time - (std::round(tp_time / 25.) * 25.);  // reduce range to [-12.5,12.5] ns
  int tp_subbx = static_cast<int>(std::round(tp_subbx_f32 * 16. / 25.));
  tp_subbx = std::clamp(tp_subbx, -8, 7);  // 4-bit, signed
  int tp_bx_check = static_cast<int>(std::round(tp_time / 25.));

  // Not sure why sometimes digi.time() returns 0?
  emtf_assert(((not(std::abs(tp_time) < 1e-6)) and (tp_bx == tp_bx_check)) or (std::abs(tp_time) < 1e-6));

  const auto tp_selection = tp_info.selection;

  // Use cluster width as quality.
  int tp_quality;

  if (tp_rpc_type == rpc::Type::kiRPC) {
    tp_quality = tp_clus_width;
  } else {
    tp_quality = tp_clus_width * 3 / 2;  // RPC strip pitch is 1.5 times the iRPC strip pitch.
  }

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

  // Get Global Coordinates
  float glob_phi;
  float glob_theta;
  double glob_rho;
  double glob_z;

  if (tp_rpc_type == rpc::Type::kiRPC) {
    // Handle iRPC Coordinates
    const RPCRoll* roll =
        dynamic_cast<const RPCRoll*>(this->context_.geometry_translator_.getRPCGeometry().roll(tp_det_id));
    const GlobalPoint& irpc_gp = roll->surface().toGlobal(LocalPoint(tp_data.x, tp_data.y, 0));

    glob_phi = tp::radToDeg(irpc_gp.phi().value());
    glob_theta = tp::radToDeg(irpc_gp.theta().value());
    glob_rho = irpc_gp.perp();
    glob_z = irpc_gp.z();
  } else {
    // Handle RPC Coordinates
    const GlobalPoint& gp = this->context_.geometry_translator_.getGlobalPoint(tp);
    glob_phi = tp::radToDeg(gp.phi().value());
    glob_theta = tp::radToDeg(gp.theta().value());
    glob_rho = gp.perp();
    glob_z = gp.z();
  }

  // Calculate EMTF Values
  const int emtf_phi = tp::calcPhiInt(sector_, glob_phi);
  const int emtf_bend = 0;
  const int emtf_theta = tp::calcThetaInt(tp_endcap_pm, glob_theta);
  const int emtf_qual = 0;
  const int emtf_site = context_.site_lut_.lookup({tp_subsystem, tp_station, tp_ring});
  const int emtf_host = context_.host_lut_.lookup({tp_subsystem, tp_station, tp_ring});
  const int emtf_zones = context_.zone_lut_.getZones(emtf_host, emtf_theta);

  emtf_assert((0 <= emtf_phi) and (emtf_phi < 5040));
  emtf_assert((1 <= emtf_theta) and (emtf_theta < 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 all the variables
  hit.setId(tp_hit_id);

  hit.setRawDetId(tp_raw_id);
  hit.setSubsystem(tp_subsystem);
  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_lo);
  hit.setStripHi(tp_strip_hi);

  hit.setWire1(tp_roll);
  hit.setWire2(0);

  hit.setBend(tp_bend);

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

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

  hit.setGlobPhi(glob_phi);
  hit.setGlobTheta(glob_theta);
  hit.setGlobPerp(glob_rho);
  hit.setGlobZ(glob_z);
  hit.setGlobTime(tp_time);

  hit.setEmtfChamber(tp_ilink);
  hit.setEmtfSegment(tp_segment_id);
  hit.setEmtfPhi(emtf_phi);
  hit.setEmtfBend(emtf_bend);
  hit.setEmtfTheta1(emtf_theta);
  hit.setEmtfTheta2(0);
  hit.setEmtfQual1(emtf_qual);
  hit.setEmtfQual2(0);
  hit.setEmtfSite(emtf_site);
  hit.setEmtfHost(emtf_host);
  hit.setEmtfZones(emtf_zones);

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

Member Data Documentation

context_

const EMTFContext& emtf::phase2::RPCTPConverter::context_

private

Definition at line 20 of file RPCTPConverter.h.

Referenced by convert().

endcap_

int emtf::phase2::RPCTPConverter::endcap_

private

Definition at line 22 of file RPCTPConverter.h.

sector_

int emtf::phase2::RPCTPConverter::sector_

private

Definition at line 22 of file RPCTPConverter.h.

Referenced by convert().