L1TMuon::GeometryTranslator Class Reference

#include <GeometryTranslator.h>

Public Member Functions
double	calculateBendAngle (const TriggerPrimitive &) const
double	calculateGlobalEta (const TriggerPrimitive &) const
double	calculateGlobalPhi (const TriggerPrimitive &) const
void	checkAndUpdateGeometry (const edm::EventSetup &)
	GeometryTranslator (edm::ConsumesCollector)
const CSCGeometry &	getCSCGeometry () const
const DTGeometry &	getDTGeometry () const
const GEMGeometry &	getGEMGeometry () const
GlobalPoint	getGlobalPoint (const TriggerPrimitive &) const
const MagneticField &	getMagneticField () const
const ME0Geometry &	getME0Geometry () const
const RPCGeometry &	getRPCGeometry () const
	~GeometryTranslator ()

Private Member Functions
double	calcCSCSpecificBend (const TriggerPrimitive &) const
double	calcCSCSpecificEta (const TriggerPrimitive &) const
double	calcCSCSpecificPhi (const TriggerPrimitive &) const
double	calcDTSpecificBend (const TriggerPrimitive &) const
double	calcDTSpecificEta (const TriggerPrimitive &) const
double	calcDTSpecificPhi (const TriggerPrimitive &) const
GlobalPoint	calcDTSpecificPoint (const TriggerPrimitive &) const
double	calcGEMSpecificBend (const TriggerPrimitive &) const
double	calcGEMSpecificEta (const TriggerPrimitive &) const
double	calcGEMSpecificPhi (const TriggerPrimitive &) const
double	calcME0SpecificBend (const TriggerPrimitive &) const
double	calcME0SpecificEta (const TriggerPrimitive &) const
double	calcME0SpecificPhi (const TriggerPrimitive &) const
double	calcRPCSpecificBend (const TriggerPrimitive &) const
double	calcRPCSpecificEta (const TriggerPrimitive &) const
double	calcRPCSpecificPhi (const TriggerPrimitive &) const
GlobalPoint	getCSCSpecificPoint (const TriggerPrimitive &) const
GlobalPoint	getGEMSpecificPoint (const TriggerPrimitive &) const
GlobalPoint	getME0SpecificPoint (const TriggerPrimitive &) const
GlobalPoint	getRPCSpecificPoint (const TriggerPrimitive &) const
bool	isCSCCounterClockwise (const std::unique_ptr< const CSCLayer > &) const

Private Attributes
edm::ESHandle< CSCGeometry >	_geocsc
edm::ESHandle< DTGeometry >	_geodt
edm::ESHandle< GEMGeometry >	_geogem
unsigned long long	_geom_cache_id
edm::ESHandle< ME0Geometry >	_geome0
edm::ESHandle< RPCGeometry >	_georpc
edm::ESHandle< MagneticField >	_magfield
unsigned long long	_magfield_cache_id
edm::ESGetToken< CSCGeometry, MuonGeometryRecord >	geocscToken_
edm::ESGetToken< DTGeometry, MuonGeometryRecord >	geodtToken_
edm::ESGetToken< GEMGeometry, MuonGeometryRecord >	geogemToken_
edm::ESGetToken< ME0Geometry, MuonGeometryRecord >	geome0Token_
edm::ESGetToken< RPCGeometry, MuonGeometryRecord >	georpcToken_
edm::ESGetToken< MagneticField, IdealMagneticFieldRecord >	magfieldToken_

Detailed Description

Definition at line 44 of file GeometryTranslator.h.

Constructor & Destructor Documentation

GeometryTranslator()

GeometryTranslator::GeometryTranslator

(

edm::ConsumesCollector

iC

)

Definition at line 25 of file GeometryTranslator.cc.

    : _geom_cache_id(0ULL),
      geodtToken_(iC.esConsumes()),
      geocscToken_(iC.esConsumes()),
      georpcToken_(iC.esConsumes()),
      geogemToken_(iC.esConsumes()),
      geome0Token_(iC.esConsumes()),
      _magfield_cache_id(0ULL),
      magfieldToken_(iC.esConsumes()) {}

~GeometryTranslator()

GeometryTranslator::~GeometryTranslator

(

)

Definition at line 35 of file GeometryTranslator.cc.

{}

Member Function Documentation

calcCSCSpecificBend()

double GeometryTranslator::calcCSCSpecificBend ( const TriggerPrimitive &  tp ) const

private

Definition at line 331 of file GeometryTranslator.cc.

References cmsswSequenceInfo::tp.

Referenced by calculateBendAngle().

{ return tp.getCSCData().bend; }

calcCSCSpecificEta()

double GeometryTranslator::calcCSCSpecificEta ( const TriggerPrimitive &  tp ) const

private

Definition at line 323 of file GeometryTranslator.cc.

References PV3DBase< T, PVType, FrameType >::eta(), getCSCSpecificPoint(), and cmsswSequenceInfo::tp.

Referenced by calculateGlobalEta().

                                                                              {
  return getCSCSpecificPoint(tp).eta();
}

calcCSCSpecificPhi()

double GeometryTranslator::calcCSCSpecificPhi ( const TriggerPrimitive &  tp ) const

private

Definition at line 327 of file GeometryTranslator.cc.

References getCSCSpecificPoint(), PV3DBase< T, PVType, FrameType >::phi(), and cmsswSequenceInfo::tp.

Referenced by calculateGlobalPhi().

                                                                              {
  return getCSCSpecificPoint(tp).phi();
}

calcDTSpecificBend()

double GeometryTranslator::calcDTSpecificBend ( const TriggerPrimitive &  tp ) const

private

Definition at line 379 of file GeometryTranslator.cc.

References cmsswSequenceInfo::tp.

Referenced by calculateBendAngle().

                                                                              {
  int bend = tp.getDTData().bendingAngle;
  double bendf = bend / 512.0;
  return bendf;
}

calcDTSpecificEta()

double GeometryTranslator::calcDTSpecificEta ( const TriggerPrimitive &  tp ) const

private

Definition at line 374 of file GeometryTranslator.cc.

References calcDTSpecificPoint(), PV3DBase< T, PVType, FrameType >::eta(), and cmsswSequenceInfo::tp.

Referenced by calculateGlobalEta().

{ return calcDTSpecificPoint(tp).eta(); }

calcDTSpecificPhi()

double GeometryTranslator::calcDTSpecificPhi ( const TriggerPrimitive &  tp ) const

private

Definition at line 376 of file GeometryTranslator.cc.

References calcDTSpecificPoint(), PV3DBase< T, PVType, FrameType >::phi(), and cmsswSequenceInfo::tp.

Referenced by calculateGlobalPhi().

{ return calcDTSpecificPoint(tp).phi(); }

calcDTSpecificPoint()

GlobalPoint GeometryTranslator::calcDTSpecificPoint ( const TriggerPrimitive &  tp ) const

private

Definition at line 342 of file GeometryTranslator.cc.

References _geodt, DTGeometry::chamber(), M_PI, PV3DBase< T, PVType, FrameType >::mag(), DTChamberId::station(), PV3DBase< T, PVType, FrameType >::theta(), and cmsswSequenceInfo::tp.

Referenced by calcDTSpecificEta(), calcDTSpecificPhi(), and getGlobalPoint().

                                                                                    {
  const DTChamberId baseid(tp.detId<DTChamberId>());
  // do not use this pointer for anything other than creating a trig geom
  std::unique_ptr<DTChamber> chamb(const_cast<DTChamber*>(_geodt->chamber(baseid)));
  std::unique_ptr<DTTrigGeom> trig_geom(new DTTrigGeom(chamb.get(), false));
  chamb.release();  // release it here so no one gets funny ideas
  // super layer one is the theta superlayer in a DT chamber
  // station 4 does not have a theta super layer
  // the BTI index from the theta trigger is an OR of some BTI outputs
  // so, we choose the BTI that's in the middle of the group
  // as the BTI that we get theta from
  // TODO:::::>>> need to make sure this ordering doesn't flip under wheel sign
  const int NBTI_theta = ((baseid.station() != 4) ? trig_geom->nCell(2) : trig_geom->nCell(3));
  const int bti_group = tp.getDTData().theta_bti_group;
  const unsigned bti_actual = bti_group * NBTI_theta / 7 + NBTI_theta / 14 + 1;
  DTBtiId thetaBTI;
  if (baseid.station() != 4 && bti_group != -1) {
    thetaBTI = DTBtiId(baseid, 2, bti_actual);
  } else {
    // since this is phi oriented it'll give us theta in the middle
    // of the chamber
    thetaBTI = DTBtiId(baseid, 3, 1);
  }
  const GlobalPoint& theta_gp = trig_geom->CMSPosition(thetaBTI);

  // local phi in sector -> global phi
  double phi = static_cast<double>(tp.getDTData().radialAngle) / 4096.0;  // 12 bits for 1 radian
  phi += tp.getDTData().sector * M_PI / 6.0;                              // add sector offset, sector is [0,11]

  return GlobalPoint(GlobalPoint::Polar(theta_gp.theta(), phi, theta_gp.mag()));
}

calcGEMSpecificBend()

double GeometryTranslator::calcGEMSpecificBend ( const TriggerPrimitive &  tp ) const

private

Definition at line 233 of file GeometryTranslator.cc.

Referenced by calculateBendAngle().

{ return 0.0; }

calcGEMSpecificEta()

double GeometryTranslator::calcGEMSpecificEta ( const TriggerPrimitive &  tp ) const

private

Definition at line 225 of file GeometryTranslator.cc.

References PV3DBase< T, PVType, FrameType >::eta(), getGEMSpecificPoint(), and cmsswSequenceInfo::tp.

Referenced by calculateGlobalEta().

                                                                              {
  return getGEMSpecificPoint(tp).eta();
}

calcGEMSpecificPhi()

double GeometryTranslator::calcGEMSpecificPhi ( const TriggerPrimitive &  tp ) const

private

Definition at line 229 of file GeometryTranslator.cc.

References getGEMSpecificPoint(), PV3DBase< T, PVType, FrameType >::phi(), and cmsswSequenceInfo::tp.

Referenced by calculateGlobalPhi().

                                                                              {
  return getGEMSpecificPoint(tp).phi();
}

calcME0SpecificBend()

double GeometryTranslator::calcME0SpecificBend ( const TriggerPrimitive &  tp ) const

private

Definition at line 207 of file GeometryTranslator.cc.

References cmsswSequenceInfo::tp.

Referenced by calculateBendAngle().

                                                                               {
  return tp.getME0Data().deltaphi * (tp.getME0Data().bend == 0 ? 1 : -1);
}

calcME0SpecificEta()

double GeometryTranslator::calcME0SpecificEta ( const TriggerPrimitive &  tp ) const

private

Definition at line 199 of file GeometryTranslator.cc.

References PV3DBase< T, PVType, FrameType >::eta(), getME0SpecificPoint(), and cmsswSequenceInfo::tp.

Referenced by calculateGlobalEta().

                                                                              {
  return getME0SpecificPoint(tp).eta();
}

calcME0SpecificPhi()

double GeometryTranslator::calcME0SpecificPhi ( const TriggerPrimitive &  tp ) const

private

Definition at line 203 of file GeometryTranslator.cc.

References getME0SpecificPoint(), PV3DBase< T, PVType, FrameType >::phi(), and cmsswSequenceInfo::tp.

Referenced by calculateGlobalPhi().

                                                                              {
  return getME0SpecificPoint(tp).phi();
}

calcRPCSpecificBend()

double GeometryTranslator::calcRPCSpecificBend ( const TriggerPrimitive &  tp ) const

private

Definition at line 259 of file GeometryTranslator.cc.

Referenced by calculateBendAngle().

{ return 0.0; }

calcRPCSpecificEta()

double GeometryTranslator::calcRPCSpecificEta ( const TriggerPrimitive &  tp ) const

private

Definition at line 249 of file GeometryTranslator.cc.

References PV3DBase< T, PVType, FrameType >::eta(), getRPCSpecificPoint(), and cmsswSequenceInfo::tp.

Referenced by calculateGlobalEta().

                                                                              {
  return getRPCSpecificPoint(tp).eta();
}

calcRPCSpecificPhi()

double GeometryTranslator::calcRPCSpecificPhi ( const TriggerPrimitive &  tp ) const

private

Definition at line 253 of file GeometryTranslator.cc.

References getRPCSpecificPoint(), PV3DBase< T, PVType, FrameType >::phi(), and cmsswSequenceInfo::tp.

Referenced by calculateGlobalPhi().

                                                                              {
  return getRPCSpecificPoint(tp).phi();
}

calculateBendAngle()

double GeometryTranslator::calculateBendAngle

(

const TriggerPrimitive &

tp

)

const

Definition at line 83 of file GeometryTranslator.cc.

References calcCSCSpecificBend(), calcDTSpecificBend(), calcGEMSpecificBend(), calcME0SpecificBend(), calcRPCSpecificBend(), L1TMuon::kCSC, L1TMuon::kDT, L1TMuon::kGEM, L1TMuon::kME0, L1TMuon::kRPC, dqmiodatasetharvest::nan, and cmsswSequenceInfo::tp.

                                                                              {
  switch (tp.subsystem()) {
    case L1TMuon::kDT:
      return calcDTSpecificBend(tp);
      break;
    case L1TMuon::kCSC:
      return calcCSCSpecificBend(tp);
      break;
    case L1TMuon::kRPC:
      return calcRPCSpecificBend(tp);
      break;
    case L1TMuon::kGEM:
      return calcGEMSpecificBend(tp);
      break;
    case L1TMuon::kME0:
      return calcME0SpecificBend(tp);
      break;
    default:
      return std::nan("Invalid TP type!");
      break;
  }
}

calculateGlobalEta()

double GeometryTranslator::calculateGlobalEta

(

const TriggerPrimitive &

tp

)

const

Definition at line 37 of file GeometryTranslator.cc.

References calcCSCSpecificEta(), calcDTSpecificEta(), calcGEMSpecificEta(), calcME0SpecificEta(), calcRPCSpecificEta(), L1TMuon::kCSC, L1TMuon::kDT, L1TMuon::kGEM, L1TMuon::kME0, L1TMuon::kRPC, dqmiodatasetharvest::nan, and cmsswSequenceInfo::tp.

                                                                              {
  switch (tp.subsystem()) {
    case L1TMuon::kDT:
      return calcDTSpecificEta(tp);
      break;
    case L1TMuon::kCSC:
      return calcCSCSpecificEta(tp);
      break;
    case L1TMuon::kRPC:
      return calcRPCSpecificEta(tp);
      break;
    case L1TMuon::kGEM:
      return calcGEMSpecificEta(tp);
      break;
    case L1TMuon::kME0:
      return calcME0SpecificEta(tp);
      break;
    default:
      return std::nan("Invalid TP type!");
      break;
  }
}

calculateGlobalPhi()

double GeometryTranslator::calculateGlobalPhi

(

const TriggerPrimitive &

tp

)

const

Definition at line 60 of file GeometryTranslator.cc.

References calcCSCSpecificPhi(), calcDTSpecificPhi(), calcGEMSpecificPhi(), calcME0SpecificPhi(), calcRPCSpecificPhi(), L1TMuon::kCSC, L1TMuon::kDT, L1TMuon::kGEM, L1TMuon::kME0, L1TMuon::kRPC, dqmiodatasetharvest::nan, and cmsswSequenceInfo::tp.

                                                                              {
  switch (tp.subsystem()) {
    case L1TMuon::kDT:
      return calcDTSpecificPhi(tp);
      break;
    case L1TMuon::kCSC:
      return calcCSCSpecificPhi(tp);
      break;
    case L1TMuon::kRPC:
      return calcRPCSpecificPhi(tp);
      break;
    case L1TMuon::kGEM:
      return calcGEMSpecificPhi(tp);
      break;
    case L1TMuon::kME0:
      return calcME0SpecificPhi(tp);
      break;
    default:
      return std::nan("Invalid TP type!");
      break;
  }
}

checkAndUpdateGeometry()

void GeometryTranslator::checkAndUpdateGeometry

(

const edm::EventSetup &

es

)

Definition at line 131 of file GeometryTranslator.cc.

References _geocsc, _geodt, _geogem, _geom_cache_id, _geome0, _georpc, _magfield, _magfield_cache_id, geocscToken_, geodtToken_, geogemToken_, relativeConstraints::geom, geome0Token_, georpcToken_, edm::EventSetup::get(), volumeBasedMagneticField_160812_cfi::magfield, and magfieldToken_.

Referenced by Phase2L1TGMTStubProducer::produce(), EMTFSetup::reload(), and emtf::phase2::EMTFContext::update().

                                                                       {
  const MuonGeometryRecord& geom = es.get<MuonGeometryRecord>();
  unsigned long long geomid = geom.cacheIdentifier();
  if (_geom_cache_id != geomid) {
    _geodt = geom.getHandle(geodtToken_);
    _geocsc = geom.getHandle(geocscToken_);
    _georpc = geom.getHandle(georpcToken_);
    _geogem = geom.getHandle(geogemToken_);
    _geome0 = geom.getHandle(geome0Token_);
    _geom_cache_id = geomid;
  }

  const IdealMagneticFieldRecord& magfield = es.get<IdealMagneticFieldRecord>();
  unsigned long long magfieldid = magfield.cacheIdentifier();
  if (_magfield_cache_id != magfieldid) {
    _magfield = magfield.getHandle(magfieldToken_);
    _magfield_cache_id = magfieldid;
  }
}

getCSCGeometry()

const CSCGeometry& L1TMuon::GeometryTranslator::getCSCGeometry ( ) const

inline

Definition at line 58 of file GeometryTranslator.h.

References _geocsc.

{ return *_geocsc; }

getCSCSpecificPoint()

GlobalPoint GeometryTranslator::getCSCSpecificPoint ( const TriggerPrimitive &  tp ) const

private

Definition at line 268 of file GeometryTranslator.cc.

References _geocsc, cms::cuda::assert(), CSCGeometry::chamber(), CSCPatternBank::getLegacyPosition(), l1ctLayer2EG_cff::id, isCSCCounterClockwise(), CSCConstants::KEY_ALCT_LAYER, PV3DBase< T, PVType, FrameType >::mag(), hltrates_dqm_sourceclient-live_cfg::offset, topSingleLeptonDQM_PU_cfi::pattern, PV3DBase< T, PVType, FrameType >::phi(), nano_mu_digi_cff::strip, PV3DBase< T, PVType, FrameType >::theta(), cmsswSequenceInfo::tp, and Geom::Phi< T1, Range >::value().

Referenced by calcCSCSpecificEta(), calcCSCSpecificPhi(), and getGlobalPoint().

                                                                                    {
  const CSCDetId id(tp.detId<CSCDetId>());
  // we should change this to weak_ptrs at some point
  // requires introducing std::shared_ptrs to geometry
  std::unique_ptr<const CSCChamber> chamb(_geocsc->chamber(id));
  assert(chamb != nullptr);  // failed to get CSC chamber
  std::unique_ptr<const CSCLayerGeometry> layer_geom(chamb->layer(CSCConstants::KEY_ALCT_LAYER)->geometry());
  std::unique_ptr<const CSCLayer> layer(chamb->layer(CSCConstants::KEY_ALCT_LAYER));

  const uint16_t halfstrip = tp.getCSCData().strip;
  const uint16_t pattern = tp.getCSCData().pattern;
  const uint16_t keyWG = tp.getCSCData().keywire;
  //const unsigned maxStrips = layer_geom->numberOfStrips();

  // so we can extend this later
  // assume TMB2007 half-strips only as baseline
  double offset = 0.0;
  switch (1) {
    case 1:
      offset = CSCPatternBank::getLegacyPosition(pattern);
  }
  const unsigned halfstrip_offs = static_cast<unsigned>(0.5 + halfstrip + offset);
  const unsigned strip = halfstrip_offs / 2 + 1;  // geom starts from 1

  // the rough location of the hit at the ALCT key layer
  // we will refine this using the half strip information
  const LocalPoint& coarse_lp = layer_geom->stripWireGroupIntersection(strip, keyWG);
  const GlobalPoint& coarse_gp = layer->surface().toGlobal(coarse_lp);

  // the strip width/4.0 gives the offset of the half-strip
  // center with respect to the strip center
  const double hs_offset = layer_geom->stripPhiPitch() / 4.0;

  // determine handedness of the chamber
  const bool ccw = isCSCCounterClockwise(layer);
  // we need to subtract the offset of even half strips and add the odd ones
  const double phi_offset = ((halfstrip_offs % 2 ? 1 : -1) * (ccw ? -hs_offset : hs_offset));

  // the global eta calculation uses the middle of the strip
  // so no need to increment it
  const GlobalPoint final_gp(
      GlobalPoint::Polar(coarse_gp.theta(), (coarse_gp.phi().value() + phi_offset), coarse_gp.mag()));

  // We need to add in some notion of the 'error' on trigger primitives
  // like the width of the wire group by the width of the strip
  // or something similar

  // release ownership of the pointers
  chamb.release();
  layer_geom.release();
  layer.release();

  return final_gp;
}

getDTGeometry()

const DTGeometry& L1TMuon::GeometryTranslator::getDTGeometry ( ) const

inline

Definition at line 57 of file GeometryTranslator.h.

References _geodt.

{ return *_geodt; }

getGEMGeometry()

const GEMGeometry& L1TMuon::GeometryTranslator::getGEMGeometry ( ) const

inline

Definition at line 60 of file GeometryTranslator.h.

References _geogem.

{ return *_geogem; }

getGEMSpecificPoint()

GlobalPoint GeometryTranslator::getGEMSpecificPoint ( const TriggerPrimitive &  tp ) const

private

Definition at line 213 of file GeometryTranslator.cc.

References _geogem, cms::cuda::assert(), GEMGeometry::etaPartition(), l1ctLayer2EG_cff::id, nano_mu_digi_cff::roll, and cmsswSequenceInfo::tp.

Referenced by calcGEMSpecificEta(), calcGEMSpecificPhi(), and getGlobalPoint().

                                                                                    {
  const GEMDetId id(tp.detId<GEMDetId>());
  const GEMEtaPartition* roll = _geogem->etaPartition(id);
  assert(roll != nullptr);  // failed to get GEM roll
  //const uint16_t pad = tp.getGEMData().pad;
  // Use half-pad precision, + 0.5 at the end to get the center of the pad (pad starts from 0)
  const float pad = (0.5 * static_cast<float>(tp.getGEMData().pad_low + tp.getGEMData().pad_hi)) + 0.5f;
  const LocalPoint& lp = roll->centreOfPad(pad);
  const GlobalPoint& gp = roll->surface().toGlobal(lp);
  return gp;
}

getGlobalPoint()

GlobalPoint GeometryTranslator::getGlobalPoint

(

const TriggerPrimitive &

tp

)

const

Definition at line 106 of file GeometryTranslator.cc.

References calcDTSpecificPoint(), getCSCSpecificPoint(), getGEMSpecificPoint(), getME0SpecificPoint(), getRPCSpecificPoint(), L1TMuon::kCSC, L1TMuon::kDT, L1TMuon::kGEM, L1TMuon::kME0, L1TMuon::kRPC, dqmiodatasetharvest::nan, runTheMatrix::ret, and cmsswSequenceInfo::tp.

Referenced by L1TPhase2GMTEndcapStubProcessor::buildCSCOnlyStub(), L1TPhase2GMTEndcapStubProcessor::buildRPCOnlyStub(), emtf::phase2::RPCTPConverter::convert(), emtf::phase2::GE0TPConverter::convert(), emtf::phase2::ME0TPConverter::convert(), emtf::phase2::GEMTPConverter::convert(), emtf::phase2::CSCTPConverter::convert(), PrimitiveConversion::convert_csc(), PrimitiveConversion::convert_dt(), PrimitiveConversion::convert_gem(), PrimitiveConversion::convert_me0(), and PrimitiveConversion::convert_rpc().

                                                                               {
  switch (tp.subsystem()) {
    case L1TMuon::kDT:
      return calcDTSpecificPoint(tp);
      break;
    case L1TMuon::kCSC:
      return getCSCSpecificPoint(tp);
      break;
    case L1TMuon::kRPC:
      return getRPCSpecificPoint(tp);
      break;
    case L1TMuon::kGEM:
      return getGEMSpecificPoint(tp);
      break;
    case L1TMuon::kME0:
      return getME0SpecificPoint(tp);
      break;
    default:
      GlobalPoint ret(
          GlobalPoint::Polar(std::nan("Invalid TP type!"), std::nan("Invalid TP type!"), std::nan("Invalid TP type!")));
      return ret;
      break;
  }
}

getMagneticField()

const MagneticField& L1TMuon::GeometryTranslator::getMagneticField ( ) const

inline

Definition at line 63 of file GeometryTranslator.h.

References _magfield.

{ return *_magfield; }

getME0Geometry()

const ME0Geometry& L1TMuon::GeometryTranslator::getME0Geometry ( ) const

inline

Definition at line 61 of file GeometryTranslator.h.

References _geome0.

{ return *_geome0; }

getME0SpecificPoint()

GlobalPoint GeometryTranslator::getME0SpecificPoint ( const TriggerPrimitive &  tp ) const

private

Definition at line 153 of file GeometryTranslator.cc.

References _geogem, _geome0, ME0Geometry::chamber(), relativeConstraints::chamber, ME0Layer::etaPartition(), GEMChamber::etaPartition(), Exception, MuonSubdetId::GEM, l1ctLayer2EG_cff::id, l1ctLayer1_patternWriters_cff::partition, nano_mu_digi_cff::roll, DetId::subdetId(), GEMGeometry::superChamber(), and cmsswSequenceInfo::tp.

Referenced by calcME0SpecificEta(), calcME0SpecificPhi(), and getGlobalPoint().

                                                                                    {
  if (tp.detId<DetId>().subdetId() == MuonSubdetId::GEM) {  // GE0
    const GEMDetId id(tp.detId<GEMDetId>());
    const GEMSuperChamber* chamber = _geogem->superChamber(id);
    const GEMChamber* keylayer = chamber->chamber(3);  // GEM key layer is layer 3
    int partition = tp.getME0Data().partition;         // 'partition' is in half-roll unit
    int iroll = (partition >> 1) + 1;
    const GEMEtaPartition* roll = keylayer->etaPartition(iroll);

    // failed to get ME0 roll
    if (roll == nullptr) {
      throw cms::Exception("Invalid GE0 Roll") << "Failed to get GEM roll for GE0" << std::endl;
    }

    // See L1Trigger/ME0Trigger/src/ME0TriggerPseudoBuilder.cc
    int phiposition = tp.getME0Data().phiposition;  // 'phiposition' is in half-strip unit
    int istrip = (phiposition >> 1);
    int phiposition2 = (phiposition & 0x1);
    float centreOfStrip = istrip + 0.25 + phiposition2 * 0.5;
    const LocalPoint& lp = roll->centreOfStrip(centreOfStrip);
    const GlobalPoint& gp = roll->surface().toGlobal(lp);
    return gp;
  } else {  // ME0
    const ME0DetId id(tp.detId<ME0DetId>());
    const ME0Chamber* chamber = _geome0->chamber(id);
    const ME0Layer* keylayer = chamber->layer(3);  // ME0 key layer is layer 3
    int partition = tp.getME0Data().partition;     // 'partition' is in half-roll unit
    int iroll = (partition >> 1) + 1;
    const ME0EtaPartition* roll = keylayer->etaPartition(iroll);

    // failed to get ME0 roll
    if (roll == nullptr) {
      throw cms::Exception("Invalid ME0 Roll") << "Failed to get ME0 roll" << std::endl;
    }

    // See L1Trigger/ME0Trigger/src/ME0TriggerPseudoBuilder.cc
    int phiposition = tp.getME0Data().phiposition;  // 'phiposition' is in half-strip unit
    int istrip = (phiposition >> 1);
    int phiposition2 = (phiposition & 0x1);
    float centreOfStrip = istrip + 0.25 + phiposition2 * 0.5;
    const LocalPoint& lp = roll->centreOfStrip(centreOfStrip);
    const GlobalPoint& gp = roll->surface().toGlobal(lp);
    return gp;
  }
}

getRPCGeometry()

const RPCGeometry& L1TMuon::GeometryTranslator::getRPCGeometry ( ) const

inline

Definition at line 59 of file GeometryTranslator.h.

References _georpc.

Referenced by emtf::phase2::RPCTPConverter::convert(), PrimitiveConversion::convert_rpc(), and EMTFSubsystemCollector::extractPrimitives().

{ return *_georpc; }

getRPCSpecificPoint()

GlobalPoint GeometryTranslator::getRPCSpecificPoint ( const TriggerPrimitive &  tp ) const

private

Definition at line 237 of file GeometryTranslator.cc.

References _georpc, cms::cuda::assert(), l1ctLayer2EG_cff::id, nano_mu_digi_cff::roll, RPCGeometry::roll(), nano_mu_digi_cff::strip, and cmsswSequenceInfo::tp.

Referenced by calcRPCSpecificEta(), calcRPCSpecificPhi(), and getGlobalPoint().

                                                                                    {
  const RPCDetId id(tp.detId<RPCDetId>());
  const RPCRoll* roll = _georpc->roll(id);
  assert(roll != nullptr);  // failed to get RPC roll
  //const int strip = static_cast<int>(tp.getRPCData().strip);
  // Use half-strip precision, - 0.5 at the end to get the center of the strip (strip starts from 1)
  const float strip = (0.5 * static_cast<float>(tp.getRPCData().strip_low + tp.getRPCData().strip_hi)) - 0.5f;
  const LocalPoint& lp = roll->centreOfStrip(strip);
  const GlobalPoint& gp = roll->surface().toGlobal(lp);
  return gp;
}

isCSCCounterClockwise()

bool GeometryTranslator::isCSCCounterClockwise ( const std::unique_ptr< const CSCLayer > &  layer ) const

private

Definition at line 333 of file GeometryTranslator.cc.

References funct::abs(), M_PI, and me0TriggerPseudoDigis_cff::nStrips.

Referenced by getCSCSpecificPoint().

                                                                                               {
  const int nStrips = layer->geometry()->numberOfStrips();
  const double phi1 = layer->centerOfStrip(1).phi();
  const double phiN = layer->centerOfStrip(nStrips).phi();
  return ((std::abs(phi1 - phiN) < M_PI && phi1 >= phiN) || (std::abs(phi1 - phiN) >= M_PI && phi1 < phiN));
}

Member Data Documentation

_geocsc

edm::ESHandle<CSCGeometry> L1TMuon::GeometryTranslator::_geocsc

private

Definition at line 68 of file GeometryTranslator.h.

Referenced by checkAndUpdateGeometry(), getCSCGeometry(), and getCSCSpecificPoint().

_geodt

edm::ESHandle<DTGeometry> L1TMuon::GeometryTranslator::_geodt

private

Definition at line 67 of file GeometryTranslator.h.

Referenced by calcDTSpecificPoint(), checkAndUpdateGeometry(), and getDTGeometry().

_geogem

edm::ESHandle<GEMGeometry> L1TMuon::GeometryTranslator::_geogem

private

Definition at line 70 of file GeometryTranslator.h.

Referenced by checkAndUpdateGeometry(), getGEMGeometry(), getGEMSpecificPoint(), and getME0SpecificPoint().

_geom_cache_id

unsigned long long L1TMuon::GeometryTranslator::_geom_cache_id

private

Definition at line 66 of file GeometryTranslator.h.

Referenced by checkAndUpdateGeometry().

_geome0

edm::ESHandle<ME0Geometry> L1TMuon::GeometryTranslator::_geome0

private

Definition at line 71 of file GeometryTranslator.h.

Referenced by checkAndUpdateGeometry(), getME0Geometry(), and getME0SpecificPoint().

_georpc

edm::ESHandle<RPCGeometry> L1TMuon::GeometryTranslator::_georpc

private

Definition at line 69 of file GeometryTranslator.h.

Referenced by checkAndUpdateGeometry(), getRPCGeometry(), and getRPCSpecificPoint().

_magfield

edm::ESHandle<MagneticField> L1TMuon::GeometryTranslator::_magfield

private

Definition at line 80 of file GeometryTranslator.h.

Referenced by checkAndUpdateGeometry(), and getMagneticField().

_magfield_cache_id

unsigned long long L1TMuon::GeometryTranslator::_magfield_cache_id

private

Definition at line 79 of file GeometryTranslator.h.

Referenced by checkAndUpdateGeometry().

geocscToken_

edm::ESGetToken<CSCGeometry, MuonGeometryRecord> L1TMuon::GeometryTranslator::geocscToken_

private

Definition at line 74 of file GeometryTranslator.h.

Referenced by checkAndUpdateGeometry().

geodtToken_

edm::ESGetToken<DTGeometry, MuonGeometryRecord> L1TMuon::GeometryTranslator::geodtToken_

private

Definition at line 73 of file GeometryTranslator.h.

Referenced by checkAndUpdateGeometry().

geogemToken_

edm::ESGetToken<GEMGeometry, MuonGeometryRecord> L1TMuon::GeometryTranslator::geogemToken_

private

Definition at line 76 of file GeometryTranslator.h.

Referenced by checkAndUpdateGeometry().

geome0Token_

edm::ESGetToken<ME0Geometry, MuonGeometryRecord> L1TMuon::GeometryTranslator::geome0Token_

private

Definition at line 77 of file GeometryTranslator.h.

Referenced by checkAndUpdateGeometry().

georpcToken_

edm::ESGetToken<RPCGeometry, MuonGeometryRecord> L1TMuon::GeometryTranslator::georpcToken_

private

Definition at line 75 of file GeometryTranslator.h.

Referenced by checkAndUpdateGeometry().

magfieldToken_

edm::ESGetToken<MagneticField, IdealMagneticFieldRecord> L1TMuon::GeometryTranslator::magfieldToken_

private

Definition at line 81 of file GeometryTranslator.h.

Referenced by checkAndUpdateGeometry().