OmtfAngleConverter Class Reference

#include <OmtfAngleConverter.h>

Inheritance diagram for OmtfAngleConverter:

Public Member Functions
virtual int	getGlobalEta (const DTChamberId dTChamberId, const L1MuDTChambThContainer *dtThDigis, int bxNum) const
virtual int	getGlobalEta (unsigned int rawid, const CSCCorrelatedLCTDigi &aDigi, float &r) const
	Convert local eta coordinate to global digital microGMT scale. More...
virtual EtaValue	getGlobalEta (unsigned int rawid, const unsigned int &aDigi) const
virtual std::vector< EtaValue >	getGlobalEta (const L1MuDTChambThContainer *dtThDigis, int bxFrom, int bxTo) const
virtual EtaValue	getGlobalEta (const CSCDetId &detId, const CSCCorrelatedLCTDigi &aDigi) const
	Convert local eta coordinate to global digital microGMT scale. More...
virtual void	getGlobalEta (const L1MuDTChambThDigi &thetaDigi, std::vector< EtaValue > &etaSegments) const
virtual int	getGlobalEtaRpc (unsigned int rawid, const unsigned int &aDigi, float &r) const
	Convert local eta coordinate to global digital microGMT scale. More...
	OmtfAngleConverter ()
	~OmtfAngleConverter () override
Public Member Functions inherited from AngleConverterBase
	AngleConverterBase ()
virtual void	checkAndUpdateGeometry (const edm::EventSetup &, const ProcConfigurationBase *config, const MuonGeometryTokens &muonGeometryTokens)
	Update the Geometry with current Event Setup. More...
float	cscChamberEtaSize (const CSCDetId &id) const
virtual void	getGlobalEta (const L1MuDTChambThDigi &thetaDigi, std::vector< EtaValue > &etaSegments) const
virtual std::vector< EtaValue >	getGlobalEta (const L1MuDTChambThContainer *dtThDigis, int bxFrom, int bxTo) const
virtual EtaValue	getGlobalEta (const CSCDetId &detId, const CSCCorrelatedLCTDigi &aDigi) const
	Convert local eta coordinate to global digital microGMT scale. More...
virtual EtaValue	getGlobalEta (unsigned int rawid, const unsigned int &aDigi) const
virtual EtaValue	getGlobalEtaCsc (const CSCDetId &detId) const
	returns the eta position of the CSC chamber More...
virtual EtaValue	getGlobalEtaDt (const DTChamberId &detId) const
virtual int	getProcessorPhi (int phiZero, l1t::tftype part, int dtScNum, int dtPhi) const
virtual int	getProcessorPhi (int phiZero, l1t::tftype part, const CSCDetId &csc, const CSCCorrelatedLCTDigi &digi, unsigned int iInput) const
virtual int	getProcessorPhi (unsigned int iProcessor, l1t::tftype part, const RPCDetId &rollId, const unsigned int &digi) const
virtual int	getProcessorPhi (int phiZero, l1t::tftype part, const RPCDetId &rollId, const unsigned int &digi1, const unsigned int &digi2) const
virtual	~AngleConverterBase ()

Additional Inherited Members
Protected Member Functions inherited from AngleConverterBase
virtual const int	findBTIgroup (const L1MuDTChambPhDigi &aDigi, const L1MuDTChambThContainer *dtThDigis)
	Find BTI group. More...
virtual bool	isCSCCounterClockwise (const CSCLayer *layer) const
	Check orientation of strips in given CSC chamber. More...
Protected Attributes inherited from AngleConverterBase
edm::ESHandle< CSCGeometry >	_geocsc
edm::ESHandle< DTGeometry >	_geodt
unsigned long long	_geom_cache_id = 0
edm::ESHandle< RPCGeometry >	_georpc
const ProcConfigurationBase *	config = nullptr
edm::ESWatcher< MuonGeometryRecord >	muonGeometryRecordWatcher
unsigned int	nPhiBins = 0
	Number of phi bins along 2Pi. More...

Detailed Description

Definition at line 13 of file OmtfAngleConverter.h.

Constructor & Destructor Documentation

OmtfAngleConverter()

OmtfAngleConverter::OmtfAngleConverter ( )

inline

Definition at line 15 of file OmtfAngleConverter.h.

: AngleConverterBase() {}

~OmtfAngleConverter()

OmtfAngleConverter::~OmtfAngleConverter ( )

override

Definition at line 101 of file OmtfAngleConverter.cc.

{}

Member Function Documentation

getGlobalEta() [1/6]

int OmtfAngleConverter::getGlobalEta ( const DTChamberId  dTChamberId, const L1MuDTChambThContainer *  dtThDigis, int  bxNum  ) const

virtual

Convert local eta coordinate to global digital microGMT scale. theta is returned only if in the dtThDigis is only one hit, otherwise eta = 95 or middle of the chamber

Definition at line 105 of file OmtfAngleConverter.cc.

References AngleConverterBase::_geodt, funct::abs(), ProcConfigurationBase::bits, DTGeometry::chamber(), L1MuDTChambThContainer::chThetaSegm(), PV3DBase< T, PVType, FrameType >::eta(), OMTFConfiguration::eta2Bits(), mps_fire::i, l1tPhase2CaloJetEmulator_cfi::iEta, L1MuDTChambThDigi::position(), DTChamberId::sector(), FWPFMaths::sgn(), DTChamberId::station(), ProcConfigurationBase::valueP1Scale, and DTChamberId::wheel().

Referenced by CscDigiToStubsConverterOmtf::addCSCstubs(), DtDigiToStubsConverterOmtf::addDTphiDigi(), and DtPhase2DigiToStubsConverterOmtf::addDTphiDigi().

                                                      {
  //const DTChamberId dTChamberId(aDigi.whNum(),aDigi.stNum(),aDigi.scNum()+1);
  DTTrigGeom trig_geom(_geodt->chamber(dTChamberId), false);

  // 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 = ((dTChamberId.station() != 4) ? trig_geom.nCell(2) : trig_geom.nCell(3));

  const L1MuDTChambThDigi *theta_segm =
      dtThDigis->chThetaSegm(dTChamberId.wheel(), dTChamberId.station(), dTChamberId.sector() - 1, bxNum);

  int bti_group = -1;
  if (theta_segm) {
    for (unsigned int i = 0; i < 7; ++i)
      if (theta_segm->position(i) && bti_group < 0)
        bti_group = i;
      else if (theta_segm->position(i) && bti_group > -1)
        bti_group = 511;
  }

  //TODO why bti_group == 511 (meaning there is more then one bit fired) is converted into 95, but not into middle of the chamber?
  int iEta = 0;
  if (bti_group == 511)
    iEta = 95;
  else if (bti_group == -1 && dTChamberId.station() == 1)
    iEta = 92;
  else if (bti_group == -1 && dTChamberId.station() == 2)
    iEta = 79;
  else if (bti_group == -1 && dTChamberId.station() == 3)
    iEta = 75;
  else if (dTChamberId.station() != 4 && bti_group >= 0) {
    unsigned bti_actual = bti_group * NBTI_theta / 7 + NBTI_theta / 14 + 1;
    DTBtiId thetaBTI = DTBtiId(dTChamberId, 2, bti_actual);
    GlobalPoint theta_gp = trig_geom.CMSPosition(thetaBTI);
    iEta = etaVal2Code(fabs(theta_gp.eta()));
  }
  int signEta = sgn(dTChamberId.wheel());
  iEta *= signEta;

  if (config->getStubEtaEncoding() == ProcConfigurationBase::StubEtaEncoding::bits)
    return OMTFConfiguration::eta2Bits(abs(iEta));
  else if (config->getStubEtaEncoding() == ProcConfigurationBase::StubEtaEncoding::valueP1Scale)
    return abs(iEta);

  return 0;
}

getGlobalEta() [2/6]

int OmtfAngleConverter::getGlobalEta ( unsigned int  rawid, const CSCCorrelatedLCTDigi &  aDigi, float &  r  ) const

virtual

Convert local eta coordinate to global digital microGMT scale.

Code taken from GeometryTranslator. Will be replaced by direct CSC phi local to global scale transformation as used in FPGA implementation

Definition at line 160 of file OmtfAngleConverter.cc.

References AngleConverterBase::_geocsc, funct::abs(), ProcConfigurationBase::bits, CSCGeometry::chamber(), PV3DBase< T, PVType, FrameType >::eta(), OMTFConfiguration::eta2Bits(), CSCLayer::geometry(), CSCCorrelatedLCTDigi::getKeyWG(), CSCCorrelatedLCTDigi::getStrip(), EcalPhiSymFlatTableProducers_cfi::id, AngleConverterBase::isCSCCounterClockwise(), CSCConstants::KEY_ALCT_LAYER, nano_mu_digi_cff::layer, CSCChamber::layer(), PV3DBase< T, PVType, FrameType >::mag(), HLT_IsoTrack_cff::offset, PV3DBase< T, PVType, FrameType >::perp(), PV3DBase< T, PVType, FrameType >::phi(), alignCSCRings::r, nano_mu_digi_cff::strip, PV3DBase< T, PVType, FrameType >::theta(), Geom::Phi< T1, Range >::value(), and ProcConfigurationBase::valueP1Scale.

                                                                                                          {

  // a lot of this is transcription and consolidation of the CSC
  // global phi calculation code
  // this works directly with the geometry
  // rather than using the old phi luts
  const CSCDetId id(rawid);
  // we should change this to weak_ptrs at some point
  // requires introducing std::shared_ptrs to geometry
  auto chamb = _geocsc->chamber(id);
  auto layer_geom = chamb->layer(CSCConstants::KEY_ALCT_LAYER)->geometry();
  auto layer = chamb->layer(CSCConstants::KEY_ALCT_LAYER);

  const uint16_t halfstrip = aDigi.getStrip();
  //const uint16_t pattern = aDigi.getPattern();
  const uint16_t keyWG = aDigi.getKeyWG();
  //const unsigned maxStrips = layer_geom->numberOfStrips();

  // so we can extend this later
  // assume TMB2007 half-strips only as baseline
  double offset = 0.0;
  //K.B. CSCPatternLUT is removed since CMSSW_11_2_1, but it looks that this offset is effectively not needed here
  //offset = CSCPatternLUT::get2007Position(pattern);

  const unsigned halfstrip_offs = 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()));

  //TODO for phase 2, add firmware like, fixed point conversion from keyWG to eta and r
  const LocalPoint lpWg = layer_geom->localCenterOfWireGroup(keyWG);
  const GlobalPoint gpWg = layer->surface().toGlobal(lpWg);
  r = round(gpWg.perp());

  if (config->getStubEtaEncoding() == ProcConfigurationBase::StubEtaEncoding::bits)
    return OMTFConfiguration::eta2Bits(abs(etaKeyWG2Code(id, keyWG)));
  else if (config->getStubEtaEncoding() == ProcConfigurationBase::StubEtaEncoding::valueP1Scale) {
    return config->etaToHwEta(abs(gpWg.eta()));
  } else {
    return 0;
  }
}

getGlobalEta() [3/6]

EtaValue AngleConverterBase::getGlobalEta

Convert local eta coordinate to global digital microGMT scale. EtaValue::etaSigma is half of the strip

Definition at line 386 of file AngleConverterBase.cc.

                                                                                             {
  const RPCDetId id(rawid);

  auto roll = _georpc->roll(id);
  const LocalPoint lp = roll->centreOfStrip((int)strip);
  const GlobalPoint gp = roll->toGlobal(lp);

  int neighbRoll = 1;  //neighbor roll in eta
  //roll->chamber()->nrolls() does not work
  if (id.region() == 0) {  //barel
    if (id.station() == 2 && ((std::abs(id.ring()) == 2 && id.layer() == 2) ||
                              (std::abs(id.ring()) != 2 && id.layer() == 1))) {  //three-roll chamber
      if (id.roll() == 2)
        neighbRoll = 1;
      else {
        neighbRoll = 2;
      }
    } else  //two-roll chamber
      neighbRoll = (id.roll() == 1 ? 3 : 1);
  } else {  //endcap
    neighbRoll = id.roll() + (id.roll() == 1 ? +1 : -1);
  }

  const RPCDetId idNeigh =
      RPCDetId(id.region(), id.ring(), id.station(), id.sector(), id.layer(), id.subsector(), neighbRoll);
  //LogTrace("l1tOmtfEventPrint")<<__FUNCTION__<<":"<<__LINE__<<" rpc "<<id<<std::endl;
  //LogTrace("l1tOmtfEventPrint")<<__FUNCTION__<<":"<<__LINE__<<" rpc "<<idNeigh<<std::endl;

  auto rollNeigh = _georpc->roll(idNeigh);
  const LocalPoint lpNeigh = rollNeigh->centreOfStrip((int)strip);
  const GlobalPoint gpNeigh = rollNeigh->toGlobal(lpNeigh);

  EtaValue etaValue = {config->etaToHwEta(gp.eta()),
                       config->etaToHwEta(std::abs(gp.eta() - gpNeigh.eta())) /
                           2,  //half of the size of the strip in eta - not precise, but OK
                       0};

  //LogTrace("l1tOmtfEventPrint")<<__FUNCTION__<<":"<<__LINE__<<" rpc "<<id<<" eta "<<gp.eta()<<" etaHw "<<etaValue.eta<<" etaSigma "<<etaValue.etaSigma<<std::endl;
  return etaValue;
}

getGlobalEta() [4/6]

std::vector< EtaValue > AngleConverterBase::getGlobalEta

Definition at line 287 of file AngleConverterBase.cc.

                                                                       {
  //LogTrace("l1tOmtfEventPrint")<<__FUNCTION__<<":"<<__LINE__<<" dtThDigis size "<<dtThDigis->getContainer()->size()<<std::endl;

  std::vector<EtaValue> etaSegments;

  for (auto& thetaDigi : (*(dtThDigis->getContainer()))) {
    if (thetaDigi.bxNum() >= bxFrom && thetaDigi.bxNum() <= bxTo) {
      getGlobalEta(thetaDigi, etaSegments);
    }
  }
  return etaSegments;
}

getGlobalEta() [5/6]

void AngleConverterBase::getGlobalEta

Definition at line 256 of file AngleConverterBase.cc.

                                                                                                                {
  const DTChamberId baseid(thetaDigi.whNum(), thetaDigi.stNum(), thetaDigi.scNum() + 1);
  DTTrigGeom trig_geom(_geodt->chamber(baseid), false);

  // super layer 2 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 = trig_geom.nCell(2);
  for (unsigned int btiGroup = 0; btiGroup < 7; ++btiGroup) {
    if (thetaDigi.position(btiGroup)) {
      unsigned btiActual = btiGroup * NBTI_theta / 7 + NBTI_theta / 14 + 1;
      DTBtiId thetaBTI = DTBtiId(baseid, 2, btiActual);
      GlobalPoint theta_gp = trig_geom.CMSPosition(thetaBTI);

      EtaValue etaValue = {
          config->etaToHwEta(theta_gp.eta()),
          0,
          thetaDigi.quality(btiGroup),
          thetaDigi.bxNum(),
          0  //TODO what about sub-bx timing???
      };
      etaSegments.emplace_back(etaValue);

      //LogTrace("l1tOmtfEventPrint")<<__FUNCTION__<<":"<<__LINE__<<" bx "<<thetaDigi.bxNum()<<" baseid "<<baseid<<" btiGroup "<<btiGroup<<" quality "<<thetaDigi.quality(btiGroup)<<" theta_gp.eta() "<<theta_gp.eta()<<" eta "<<etaValue.eta<<" etaSigma "<<etaValue.etaSigma<<std::endl;
    }
  }
}

getGlobalEta() [6/6]

EtaValue AngleConverterBase::getGlobalEta

Convert local eta coordinate to global digital microGMT scale.

Code taken from GeometryTranslator. Will be replaced by direct CSC phi local to global scale transformation as used in FPGA implementation

Definition at line 333 of file AngleConverterBase.cc.

                                                                                                        {

  // alot of this is transcription and consolidation of the CSC
  // global phi calculation code
  // this works directly with the geometry
  // rather than using the old phi luts

  auto chamb = _geocsc->chamber(detId);
  auto layer_geom = chamb->layer(CSCConstants::KEY_ALCT_LAYER)->geometry();
  auto layer = chamb->layer(CSCConstants::KEY_ALCT_LAYER);

  const uint16_t keyWG = aDigi.getKeyWG();

  const LocalPoint lpWg = layer_geom->localCenterOfWireGroup(keyWG);
  const GlobalPoint gpWg = layer->surface().toGlobal(lpWg);

  EtaValue etaSegment = {
      config->etaToHwEta(gpWg.eta()),
      0,  //config->etaToHwEta(cscChamberEtaSize(id) ),
      0,
      aDigi.getBX(),
      0  //tming???
  };

  //LogTrace("l1tOmtfEventPrint")<<__FUNCTION__<<":"<<__LINE__<<" csc "<<detId<<" eta "<<gpWg.eta()<<" etaHw "<<etaSegment.eta<<" etaSigma "<<etaSegment.etaSigma<<std::endl;
  return etaSegment;
}

getGlobalEtaRpc()

int OmtfAngleConverter::getGlobalEtaRpc ( unsigned int  rawid, const unsigned int &  aDigi, float &  r  ) const

virtual

Convert local eta coordinate to global digital microGMT scale.

Definition at line 224 of file OmtfAngleConverter.cc.

References AngleConverterBase::_georpc, funct::abs(), ProcConfigurationBase::bits, OMTFConfiguration::eta2Bits(), runTauDisplay::gp, EcalPhiSymFlatTableProducers_cfi::id, alignCSCRings::r, nano_mu_digi_cff::region, nano_mu_digi_cff::roll, RPCGeometry::roll(), nano_mu_digi_cff::strip, and ProcConfigurationBase::valueP1Scale.

Referenced by RpcDigiToStubsConverterOmtf::addRPCstub().

                                                                                                     {
  const RPCDetId id(rawid);
  auto roll = _georpc->roll(id);
  const LocalPoint lp = roll->centreOfStrip((int)strip);
  const GlobalPoint gp = roll->toGlobal(lp);

  if (id.region() != 0) {  //outside barrel
    r = gp.perp();
  }

  if (config->getStubEtaEncoding() == ProcConfigurationBase::StubEtaEncoding::bits)
    return OMTFConfiguration::eta2Bits(abs(etaVal2Code(gp.eta())));
  else if (config->getStubEtaEncoding() == ProcConfigurationBase::StubEtaEncoding::valueP1Scale)
    return abs(config->etaToHwEta((gp.eta())));

  return 0;
}