tmtt::Sector Class Reference

#include <Sector.h>

Public Member Functions
float	etaMax () const
float	etaMin () const
unsigned int	iEtaReg () const
bool	inside (const Stub *stub) const
bool	insideEta (const Stub *stub) const
bool	insideEtaReg (const TP &tp) const
std::vector< bool >	insideEtaSubSecs (const Stub *stub) const
bool	insidePhi (const Stub *stub) const
bool	insidePhiSec (const TP &tp) const
unsigned int	iPhiSec () const
void	numStubsInside (const TP &tp, unsigned int &nStubsInsideEtaPhi, unsigned int &nStubsInsideEta, unsigned int &nStubsInsidePhi) const
float	phiCentre () const
	Sector (const Settings *settings, unsigned int iPhiSec, unsigned int iEtaSec)
float	sectorHalfWidth () const
std::unordered_map< const Stub *, std::pair< bool, bool > >	stubsInside (const TP &tp) const
float	zAtChosenR_Max () const
float	zAtChosenR_Min () const

Private Member Functions
int64_t	forceBitWidth (const float value, const UInt_t nBits) const
bool	insideEtaRange (const Stub *stub, float zRangeMin, float zRangeMax) const
std::vector< bool >	subEtaFwCalc (const int rT, const int z) const

Private Attributes
float	assumedPhiTrkRes_
float	beamWindowZ_
bool	calcPhiTrkRes_
float	chosenRofPhi_
float	chosenRofZ_
float	etaMax_
float	etaMin_
unsigned int	iEtaReg_
unsigned int	iPhiSec_
float	minPt_
unsigned int	numSubSecsEta_
float	phiCentre_
float	sectorHalfWidth_
const Settings *	settings_
bool	useStubPhi_
bool	useStubPhiTrk_
float	zOuterMax_
std::vector< float >	zOuterMaxSub_
float	zOuterMin_
std::vector< float >	zOuterMinSub_

Detailed Description

Definition at line 14 of file Sector.h.

Constructor & Destructor Documentation

Sector()

tmtt::Sector::Sector	(	const Settings *	settings,
		unsigned int	iPhiSec,
		unsigned int	iEtaSec
	)

Definition at line 13 of file Sector.cc.

References chosenRofZ_, etaMax_, etaMin_, JetChargeProducer_cfi::exp, dqmMemoryStats::float, mps_fire::i, iPhiSec(), M_PI, tmtt::Settings::numPhiNonants(), tmtt::Settings::numPhiSectors(), numSubSecsEta_, phiCentre_, sectorHalfWidth_, funct::tan(), zOuterMaxSub_, and zOuterMinSub_.

      : settings_(settings),
        // Sector number
        iPhiSec_(iPhiSec),
        iEtaReg_(iEtaReg),

        beamWindowZ_(settings->beamWindowZ()),  // Assumed half-length of beam-spot

        //===  Characteristics of this eta region.
        // Using lines of specified rapidity drawn from centre of CMS, determine the z coords at which
        // they cross the radius chosenRofZ_.
        etaMin_(settings->etaRegions()[iEtaReg]),
        etaMax_(settings->etaRegions()[iEtaReg + 1]),
        chosenRofZ_(settings->chosenRofZ()),
        // Get range in z of tracks covered by this sector at chosen radius from beam-line
        zOuterMin_(chosenRofZ_ / tan(2. * atan(exp(-etaMin_)))),
        zOuterMax_(chosenRofZ_ / tan(2. * atan(exp(-etaMax_)))),

        //=== Characteristics of this phi region.
        chosenRofPhi_(settings->chosenRofPhi()),
        minPt_(settings->houghMinPt()),  // Min Pt covered by  HT array.
        assumedPhiTrkRes_(settings->assumedPhiTrkRes()),
        useStubPhi_(settings->useStubPhi()),
        useStubPhiTrk_(settings->useStubPhiTrk()),
        calcPhiTrkRes_(settings->calcPhiTrkRes()),
        //=== Check if subsectors in eta are being used within each sector.
        numSubSecsEta_(settings->numSubSecsEta()) {
    // Centre of phi (tracking) nonant zero must be along x-axis to be consistent with tracker cabling map.
    // Define phi sector zero  to start at lower end of phi range in nonant 0.
    float phiCentreSec0 = -M_PI / float(settings->numPhiNonants()) + M_PI / float(settings->numPhiSectors());
    // Centre of sector in phi
    phiCentre_ = 2. * M_PI * float(iPhiSec) / float(settings->numPhiSectors()) + phiCentreSec0;
    sectorHalfWidth_ = M_PI / float(settings->numPhiSectors());  // Sector half width excluding overlaps.

    // If eta subsectors have equal width in rapidity, do this.
    float subSecWidth = (etaMax_ - etaMin_) / float(numSubSecsEta_);
    for (unsigned int i = 0; i < numSubSecsEta_; i++) {
      float subSecEtaMin = etaMin_ + i * subSecWidth;
      float subSecEtaMax = subSecEtaMin + subSecWidth;
      float subSecZmin = chosenRofZ_ / tan(2. * atan(exp(-subSecEtaMin)));
      float subSecZmax = chosenRofZ_ / tan(2. * atan(exp(-subSecEtaMax)));
      zOuterMinSub_.push_back(subSecZmin);
      zOuterMaxSub_.push_back(subSecZmax);
    }
  }

Member Function Documentation

etaMax()

float tmtt::Sector::etaMax ( ) const

inline

Definition at line 31 of file Sector.h.

References etaMax_.

Referenced by tmtt::L1track3D::cheat().

{ return etaMax_; }        // Eta range covered by this sector.

etaMin()

float tmtt::Sector::etaMin ( ) const

inline

Definition at line 30 of file Sector.h.

References etaMin_.

Referenced by tmtt::L1track3D::cheat().

{ return etaMin_; }        // Eta range covered by this sector.

forceBitWidth()

int64_t tmtt::Sector::forceBitWidth ( const float  value, const UInt_t  nBits  ) const

private

Definition at line 192 of file Sector.cc.

References funct::abs(), gpuClustering::pixelStatus::mask, mps_fire::result, and Validation_hcalonly_cfi::sign.

Referenced by subEtaFwCalc().

                                                                           {
    // slightly hand-waving treatment of 2s complement
    int64_t sign = 1;
    if (value < 0)
      sign = -1;
    int64_t iValue = int64_t(std::abs(value));
    int64_t mask = (int64_t(1) << nBits) - int64_t(1);
    int64_t result = sign * (iValue & mask);
    if (std::abs(result - value) > 1)
      throw cms::Exception("LogicError")
          << "Sector::forceBitWidth is messing up by using too few bits to digitize number"
          << " nBits=" << nBits << " Input float=" << value << " Output digi = " << result;
    return result;
    // Check that result is compatible with value. Throw error if not.
  }

iEtaReg()

unsigned int tmtt::Sector::iEtaReg ( ) const

inline

Definition at line 28 of file Sector.h.

References iEtaReg_.

{ return iEtaReg_; }

inside()

bool tmtt::Sector::inside ( const Stub *  stub ) const

inline

Definition at line 20 of file Sector.h.

References insideEta(), and insidePhi().

Referenced by insideEta(), insideEtaRange(), insideEtaSubSecs(), and stubsInside().

{ return (this->insideEta(stub) && this->insidePhi(stub)); }

insideEta()

bool tmtt::Sector::insideEta

(

const Stub *

stub

)

const

Definition at line 61 of file Sector.cc.

References inside(), insideEtaRange(), zOuterMax_, and zOuterMin_.

Referenced by inside(), numStubsInside(), and stubsInside().

                                               {
    // Lower edge of this eta region defined by line from (r,z) = (0,-beamWindowZ) to (chosenRofZ_, zOuterMin_).
    // Upper edge of this eta region defined by line from (r,z) = (0, beamWindowZ) to (chosenRofZ_, zOuterMax_).

    bool inside = this->insideEtaRange(stub, zOuterMin_, zOuterMax_);
    return inside;
  }

insideEtaRange()

bool tmtt::Sector::insideEtaRange ( const Stub *  stub, float  zRangeMin, float  zRangeMax  ) const

private

Definition at line 93 of file Sector.cc.

References funct::abs(), beamWindowZ_, chosenRofZ_, inside(), tmtt::Stub::r(), tmtt::Stub::z(), photonAnalyzer_cfi::zMax, and photonAnalyzer_cfi::zMin.

Referenced by insideEta(), and insideEtaSubSecs().

                                                                                      {
    // Lower edge of this eta region defined by line from (r,z) = (0,-beamWindowZ) to (chosenRofZ_, zRangeMin).
    // Upper edge of this eta region defined by line from (r,z) = (0, beamWindowZ) to (chosenRofZ_, zRangeMax).

    float zMin, zMax;
    bool inside;

    // Calculate z coordinate of lower edge of this eta region, evaluated at radius of stub.
    zMin = (zRangeMin * stub->r() - beamWindowZ_ * std::abs(stub->r() - chosenRofZ_)) / chosenRofZ_;
    // Calculate z coordinate of upper edge of this eta region, evaluated at radius of stub.
    zMax = (zRangeMax * stub->r() + beamWindowZ_ * std::abs(stub->r() - chosenRofZ_)) / chosenRofZ_;

    inside = (stub->z() > zMin && stub->z() < zMax);
    return inside;
  }

insideEtaReg()

bool tmtt::Sector::insideEtaReg ( const TP &  tp ) const

inline

Definition at line 56 of file Sector.h.

References chosenRofZ_, cmsswSequenceInfo::tp, zOuterMax_, and zOuterMin_.

                                          {
      return (tp.trkZAtR(chosenRofZ_) > zOuterMin_ && tp.trkZAtR(chosenRofZ_) < zOuterMax_);
    }

insideEtaSubSecs()

vector< bool > tmtt::Sector::insideEtaSubSecs

(

const Stub *

stub

)

const

Definition at line 71 of file Sector.cc.

References tmtt::Stub::digitalStub(), tmtt::Settings::enableDigitize(), mps_fire::i, tmtt::DigitalStub::iDigi_Rt(), tmtt::DigitalStub::iDigi_Z(), inside(), insideEtaRange(), numSubSecsEta_, settings_, subEtaFwCalc(), zOuterMaxSub_, and zOuterMinSub_.

                                                              {
    if (settings_->enableDigitize() && numSubSecsEta_ == 2) {
      // Use (complicated) digitized firmware emulation
      return subEtaFwCalc(stub->digitalStub()->iDigi_Rt(), stub->digitalStub()->iDigi_Z());

    } else {
      // Use (simpler) floating point calculation.

      vector<bool> insideVec;

      // Loop over subsectors.
      for (unsigned int i = 0; i < numSubSecsEta_; i++) {
        bool inside = this->insideEtaRange(stub, zOuterMinSub_[i], zOuterMaxSub_[i]);
        insideVec.push_back(inside);
      }

      return insideVec;
    }
  }

insidePhi()

bool tmtt::Sector::insidePhi

(

const Stub *

stub

)

const

Definition at line 111 of file Sector.cc.

References funct::abs(), assumedPhiTrkRes_, calcPhiTrkRes_, chosenRofPhi_, delPhi(), reco::deltaPhi(), SiStripPI::min, minPt_, tmtt::Stub::phi(), phiCentre_, tmtt::Stub::phiDiff(), sectorHalfWidth_, tmtt::Stub::trkPhiAtR(), tmtt::Stub::trkPhiAtRcut(), useStubPhi_, and useStubPhiTrk_.

Referenced by inside(), numStubsInside(), and stubsInside().

                                               {
    // N.B. The logic here for preventing a stub being assigned to > 2 sectors seems overly agressive.
    // But attempts at improving it have failed ...

    bool okPhi = true;
    bool okPhiTrk = true;

    if (useStubPhi_) {
      float delPhi =
          reco::deltaPhi(stub->phi(), phiCentre_);  // Phi difference between stub & sector in range -PI to +PI.
      float tolerancePhi = stub->phiDiff(
          chosenRofPhi_, minPt_);  // How much stub phi might differ from track phi because of track curvature.
      float outsidePhi = std::abs(delPhi) - sectorHalfWidth_ -
                         tolerancePhi;  // If > 0, then stub is not compatible with being inside this sector.
      if (outsidePhi > 0)
        okPhi = false;
    }

    if (useStubPhiTrk_) {
      // Estimate either phi0 of track from stub info, or phi of the track at radius chosenRofPhi_.
      float phiTrk = stub->trkPhiAtR(chosenRofPhi_);
      // Phi difference between stub & sector in range -PI to +PI.
      float delPhiTrk = reco::deltaPhi(phiTrk, phiCentre_);
      // Set tolerance equal to nominal resolution assumed in phiTrk
      float tolerancePhiTrk = assumedPhiTrkRes_ * (2 * sectorHalfWidth_);
      if (calcPhiTrkRes_) {
        // Calculate uncertainty in phiTrk due to poor resolution in stub bend
        float phiTrkRes = stub->trkPhiAtRcut(chosenRofPhi_);
        // Reduce tolerance if this is smaller than the nominal assumed resolution.
        tolerancePhiTrk = min(tolerancePhiTrk, phiTrkRes);
      }
      // If following > 0, then stub is not compatible with being inside this sector.
      float outsidePhiTrk = std::abs(delPhiTrk) - sectorHalfWidth_ - tolerancePhiTrk;

      if (outsidePhiTrk > 0)
        okPhiTrk = false;
    }

    return (okPhi && okPhiTrk);
  }

insidePhiSec()

bool tmtt::Sector::insidePhiSec ( const TP &  tp ) const

inline

Definition at line 53 of file Sector.h.

References funct::abs(), chosenRofPhi_, phiCentre_, sectorHalfWidth_, and cmsswSequenceInfo::tp.

                                          {
      return (std::abs(tp.trkPhiAtR(chosenRofPhi_) - phiCentre_) < sectorHalfWidth_);
    }

iPhiSec()

unsigned int tmtt::Sector::iPhiSec ( ) const

inline

Definition at line 27 of file Sector.h.

References iPhiSec_.

Referenced by Sector().

{ return iPhiSec_; }  // Sector number.

numStubsInside()

void tmtt::Sector::numStubsInside	(	const TP &	tp,
		unsigned int &	nStubsInsideEtaPhi,
		unsigned int &	nStubsInsideEta,
		unsigned int &	nStubsInsidePhi
	)		const

Definition at line 171 of file Sector.cc.

References insideEta(), insidePhi(), and stubsInside().

                                                                   {
    nStubsInsideEtaPhi = 0;
    nStubsInsideEta = 0;
    nStubsInsidePhi = 0;
    for (const auto& iter : this->stubsInside(tp)) {
      bool insidePhi = iter.second.first;
      bool insideEta = iter.second.second;
      if (insidePhi && insideEta)
        nStubsInsideEtaPhi++;
      if (insideEta)
        nStubsInsideEta++;
      if (insidePhi)
        nStubsInsidePhi++;
    }
  }

phiCentre()

float tmtt::Sector::phiCentre ( ) const

inline

Definition at line 29 of file Sector.h.

References phiCentre_.

Referenced by tmtt::L1track3D::cheat().

{ return phiCentre_; }  // Return phi of centre of this sector.

sectorHalfWidth()

float tmtt::Sector::sectorHalfWidth ( ) const

inline

Definition at line 33 of file Sector.h.

References sectorHalfWidth_.

{ return sectorHalfWidth_; }  // Half width in phi of sector measured in radians.

stubsInside()

unordered_map< const Stub *, pair< bool, bool > > tmtt::Sector::stubsInside

(

const TP &

tp

)

const

Definition at line 156 of file Sector.cc.

References inside(), insideEta(), insidePhi(), and cmsswSequenceInfo::tp.

Referenced by numStubsInside().

                                                                                      {
    unordered_map<const Stub*, pair<bool, bool> > inside;
    // Loop over stubs produced by tracking particle
    const vector<const Stub*>& assStubs = tp.assocStubs();
    for (const Stub* stub : assStubs) {
      // Check if this stub is inside sector
      inside[stub] = pair<bool, bool>(this->insidePhi(stub), this->insideEta(stub));
    }
    return inside;
  }

subEtaFwCalc()

vector< bool > tmtt::Sector::subEtaFwCalc ( const int  rT, const int  z  ) const

private

Definition at line 211 of file Sector.cc.

References funct::abs(), beamWindowZ_, chosenRofPhi_, chosenRofZ_, cm_to_mm, tmtt::Settings::etaRegions(), JetChargeProducer_cfi::exp, forceBitWidth(), g, iEtaReg_, funct::pow(), tmtt::Settings::rtBits(), tmtt::Settings::rtRange(), settings_, funct::tan(), tmtt::Settings::zBits(), tmtt::Settings::zRange(), and beampixel_dqm_sourceclient-live_cfg::zRange.

Referenced by insideEtaSubSecs().

                                                                   {
    // Note number of reference bits used to digitize rT and z, used when GP authors determined some constants below.
    unsigned int rtBitsRef = 10;
    unsigned int zBitsRef = 12;

    // This replaces Kristian's hard-wired constants with configurable ones.
    unsigned int rtBits = settings_->rtBits();
    unsigned int zBits = settings_->zBits();
    float rtRange = settings_->rtRange();
    float zRange = settings_->zRange();
    constexpr float cm_to_mm = 10.;  // firwmare is in mm and CMSSW in cm.
    float zBase = cm_to_mm / (pow(2, zBits) / zRange);
    float rTBase = cm_to_mm / (pow(2, rtBits) / rtRange);

    // Number of bits used by DSP in UltraScale-Plus FPGA (where DSP does D = A*B + C)
    constexpr unsigned int nDSPa = 27;
    //constexpr unsigned int nDSPb = 18;
    constexpr unsigned int nDSPc = 48;
    constexpr unsigned int nDSPd = 48;

    // unit transformations: firmware uses mm, software uses cm
    float BeamWindow = cm_to_mm * beamWindowZ_;
    float T_rphi = cm_to_mm * chosenRofPhi_;
    float T_rz = cm_to_mm * chosenRofZ_;

    // actual algorithm as used in firmware, mostly using same variable names
    float Beam_over_T = BeamWindow / T_rz;
    // Value chosen so that number digitized below when calculating "bot" uses most of the nDSPa bits, without overflowing them. This is done assuming reference number of bits for rT and z mentioned above.
    unsigned int nShiftA = 24;
    // Guess from to keep "bot" in correct range (nDSPa) if number of digitsation bits are changed.
    nShiftA += (rtBits - rtBitsRef) - (zBits - zBitsRef);
    float Beam_over_T_base = 1. / (1 << nShiftA);
    int64_t bot = forceBitWidth(Beam_over_T * rTBase / zBase / Beam_over_T_base, nDSPa);
    int64_t bw = forceBitWidth(BeamWindow / zBase / Beam_over_T_base, nDSPc);
    float etaSecMid = (settings_->etaRegions()[iEtaReg_] + settings_->etaRegions()[iEtaReg_ + 1]) / 2.0;
    float tanlSecMid = 1.0 / tan(2.0 * atan(exp(-etaSecMid)));
    // Value chosen so that number digitized below when calculating "tanlSec_Mid" uses most of the nDSPa bits, without overflowing them. This is done assuming reference number of bits for rT and z mentioned above.
    unsigned int nShiftB = 16;
    // Guess to keep "tanlSec_Mid" in correct range (nDSPa) if number of digitsation bits are changed.
    nShiftB += (rtBits - rtBitsRef) - (zBits - zBitsRef);
    float tanlSecBase = 1. / (1 << nShiftB);
    int64_t tanlSec_Mid = forceBitWidth(int(tanlSecMid * rTBase / zBase / tanlSecBase), nDSPa);
    // Number of extra bits used to digitise r instead of rT within GP code, if both encoded as signed int.
    constexpr unsigned int nExtraBitsR = 2;
    unsigned int rBits = rtBits + nExtraBitsR;
    int64_t r = forceBitWidth(rT + T_rphi / rTBase, rBits);
    int64_t g = forceBitWidth(bot * r - bw, nDSPd);
    int64_t absg = abs(g);
    // Number of useful bits left of the nDSPd assigned to "absg" after right-shifting by nShiftA bits.
    const unsigned nBitsRemainingA = nDSPd - nShiftA;
    int64_t shift_g = forceBitWidth((absg >> nShiftA), nBitsRemainingA);
    // Number of bits is sum of those in two numbers being multiplied.
    int64_t tlsr = forceBitWidth(tanlSec_Mid * r, nDSPa + rBits);
    // Number of useful bits left of (nDSPa + rBits) assigned to "tlsr" after right-shifting by nShiftB bits.
    const unsigned nBitsRemainingB = (nDSPa + rBits) - nShiftB;
    int64_t shift_tlsr = forceBitWidth((tlsr >> nShiftB), nBitsRemainingB);

    vector<bool> insideVec;
    insideVec.push_back(z <= (shift_tlsr + shift_g));
    insideVec.push_back(z >= (shift_tlsr - shift_g));
    return insideVec;
  }

zAtChosenR_Max()

float tmtt::Sector::zAtChosenR_Max ( ) const

inline

Definition at line 37 of file Sector.h.

References zOuterMax_.

{ return zOuterMax_; }

zAtChosenR_Min()

float tmtt::Sector::zAtChosenR_Min ( ) const

inline

Definition at line 34 of file Sector.h.

References zOuterMin_.

                                 {
      return zOuterMin_;
    }  // Range in z of particle at chosen radius from beam line covered by this sector.

Member Data Documentation

assumedPhiTrkRes_

float tmtt::Sector::assumedPhiTrkRes_

private

Definition at line 88 of file Sector.h.

Referenced by insidePhi().

beamWindowZ_

float tmtt::Sector::beamWindowZ_

private

Definition at line 76 of file Sector.h.

Referenced by insideEtaRange(), and subEtaFwCalc().

calcPhiTrkRes_

bool tmtt::Sector::calcPhiTrkRes_

private

Definition at line 91 of file Sector.h.

Referenced by insidePhi().

chosenRofPhi_

float tmtt::Sector::chosenRofPhi_

private

Definition at line 86 of file Sector.h.

Referenced by insidePhi(), insidePhiSec(), and subEtaFwCalc().

chosenRofZ_

float tmtt::Sector::chosenRofZ_

private

Definition at line 79 of file Sector.h.

Referenced by insideEtaRange(), insideEtaReg(), Sector(), and subEtaFwCalc().

etaMax_

float tmtt::Sector::etaMax_

private

Definition at line 78 of file Sector.h.

Referenced by etaMax(), and Sector().

etaMin_

float tmtt::Sector::etaMin_

private

Definition at line 77 of file Sector.h.

Referenced by etaMin(), and Sector().

iEtaReg_

unsigned int tmtt::Sector::iEtaReg_

private

Definition at line 75 of file Sector.h.

Referenced by iEtaReg(), and subEtaFwCalc().

iPhiSec_

unsigned int tmtt::Sector::iPhiSec_

private

Definition at line 74 of file Sector.h.

Referenced by iPhiSec().

minPt_

float tmtt::Sector::minPt_

private

Definition at line 87 of file Sector.h.

Referenced by insidePhi().

numSubSecsEta_

unsigned int tmtt::Sector::numSubSecsEta_

private

Definition at line 94 of file Sector.h.

Referenced by insideEtaSubSecs(), and Sector().

phiCentre_

float tmtt::Sector::phiCentre_

private

Definition at line 84 of file Sector.h.

Referenced by insidePhi(), insidePhiSec(), phiCentre(), and Sector().

sectorHalfWidth_

float tmtt::Sector::sectorHalfWidth_

private

Definition at line 85 of file Sector.h.

Referenced by insidePhi(), insidePhiSec(), Sector(), and sectorHalfWidth().

settings_

const Settings* tmtt::Sector::settings_

private

Definition at line 71 of file Sector.h.

Referenced by insideEtaSubSecs(), and subEtaFwCalc().

useStubPhi_

bool tmtt::Sector::useStubPhi_

private

Definition at line 89 of file Sector.h.

Referenced by insidePhi().

useStubPhiTrk_

bool tmtt::Sector::useStubPhiTrk_

private

Definition at line 90 of file Sector.h.

Referenced by insidePhi().

zOuterMax_

float tmtt::Sector::zOuterMax_

private

Definition at line 81 of file Sector.h.

Referenced by insideEta(), insideEtaReg(), and zAtChosenR_Max().

zOuterMaxSub_

std::vector<float> tmtt::Sector::zOuterMaxSub_

private

Definition at line 96 of file Sector.h.

Referenced by insideEtaSubSecs(), and Sector().

zOuterMin_

float tmtt::Sector::zOuterMin_

private

Definition at line 80 of file Sector.h.

Referenced by insideEta(), insideEtaReg(), and zAtChosenR_Min().

zOuterMinSub_

std::vector<float> tmtt::Sector::zOuterMinSub_

private

Definition at line 95 of file Sector.h.

Referenced by insideEtaSubSecs(), and Sector().