trackerDTC::Stub Class Reference

Class to represent an outer tracker Stub. More...

#include <Stub.h>

Public Member Functions
int	bend () const
tt::Frame	frame (int region) const
bool	inRegion (int region) const
	Stub (const edm::ParameterSet &, const tt::Setup *, const LayerEncoding *, tt::SensorModule *, const TTStubRef &)
TTStubRef	ttStubRef () const
bool	valid () const
	~Stub ()

Private Member Functions
double	digi (double value, double precision) const
tt::Frame	formatHybrid (int region) const
tt::Frame	formatTMTT (int region) const

Private Attributes
int	bend_
double	c_
int	col_
std::pair< double, double >	cot_
double	d_
bool	hybrid_
std::pair< double, double >	inv2R_
const LayerEncoding *	layerEncoding_
double	m_
double	phi_
std::pair< double, double >	phiT_
double	r_
std::vector< int >	regions_
int	row_
int	rowLUT_
int	rowSub_
const tt::Setup *	setup_
tt::SensorModule *	sm_
TTStubRef	ttStubRef_
bool	valid_
double	z_

Detailed Description

Class to represent an outer tracker Stub.

Author

Thomas Schuh

Date

2020, Jan

Definition at line 18 of file Stub.h.

Constructor & Destructor Documentation

Stub()

trackerDTC::Stub::Stub	(	const edm::ParameterSet &	iConfig,
		const tt::Setup *	setup,
		const LayerEncoding *	layerEncoding,
		tt::SensorModule *	sm,
		const TTStubRef &	ttStubRef
	)

Definition at line 14 of file Stub.cc.

References funct::abs(), bend_, DummyCfis::c, c_, col_, cot_, d_, digi(), HLT_2024v14_cff::distance, spr::find(), hybrid_, createfilelist::int, inv2R_, visualization-live-secondInstance_cfg::m, m_, SiStripPI::max, SiStripPI::min, PV_cfg::minPt, phi_, phiT_, funct::pow(), r_, regions_, row_, rowLUT_, rowSub_, singleTopDQM_cfi::setup, ALPAKA_ACCELERATOR_NAMESPACE::ecal::reconstruction::internal::barrel::sm(), mathSSE::sqrt(), edm::swap(), ttStubRef(), valid_, PV2DBase< T, PVType, FrameType >::x(), x, PV2DBase< T, PVType, FrameType >::y(), and z_.

      : setup_(setup),
        layerEncoding_(layerEncoding),
        sm_(sm),
        ttStubRef_(ttStubRef),
        hybrid_(iConfig.getParameter<bool>("UseHybrid")),
        valid_(true) {
    regions_.reserve(setup->numOverlappingRegions());
    // get stub local coordinates
    const MeasurementPoint& mp = ttStubRef->clusterRef(0)->findAverageLocalCoordinatesCentered();

    // convert to uniformed local coordinates

    // column number in pitch units
    col_ = (int)floor(pow(-1, sm->signCol()) * (mp.y() - sm->numColumns() / 2) / setup->baseCol());
    // row number in half pitch units
    row_ = (int)floor(pow(-1, sm->signRow()) * (mp.x() - sm->numRows() / 2) / setup->baseRow());
    // bend number in quarter pitch units
    bend_ = (int)floor(pow(-1, sm->signBend()) * (ttStubRef->bendBE()) / setup->baseBend());
    // reduced row number for look up
    rowLUT_ = (int)floor((double)row_ / pow(2., setup->widthRow() - setup->dtcWidthRowLUT()));
    // sub row number inside reduced row number
    rowSub_ = row_ - (rowLUT_ + .5) * pow(2, setup->widthRow() - setup->dtcWidthRowLUT());

    // convert local to global coordinates

    const double y = (col_ + .5) * setup->baseCol() * sm->pitchCol();
    // radius of a column of strips/pixel in cm
    d_ = sm->r() + y * sm->sinTilt();
    // stub z in cm
    z_ = digi(sm->z() + y * sm->cosTilt(), setup->tmttBaseZ());

    const double x0 = rowLUT_ * setup->baseRow() * setup->dtcNumMergedRows() * sm->pitchRow();
    const double x1 = (rowLUT_ + 1) * setup->baseRow() * setup->dtcNumMergedRows() * sm->pitchRow();
    const double x = (rowLUT_ + .5) * setup->baseRow() * setup->dtcNumMergedRows() * sm->pitchRow();
    // stub r in cm
    r_ = sqrt(d_ * d_ + x * x);

    const double phi0 = sm->phi() + atan2(x0, d_);
    const double phi1 = sm->phi() + atan2(x1, d_);
    const double c = (phi0 + phi1) / 2.;
    const double m = (phi1 - phi0) / setup->dtcNumMergedRows();

    // intercept of linearized stub phi in rad
    c_ = digi(c, setup->tmttBasePhi());
    // slope of linearized stub phi in rad / strip
    m_ = digi(m, setup->dtcBaseM());

    if (hybrid_) {
      if (abs(z_ / r_) > setup->hybridMaxCot())
        // did not pass eta cut
        valid_ = false;
    } else {
      // extrapolated z at radius T assuming z0=0
      const double zT = setup->chosenRofZ() * z_ / r_;
      // extrapolated z0 window at radius T
      const double dZT = setup->beamWindowZ() * abs(1. - setup->chosenRofZ() / r_);
      double zTMin = zT - dZT;
      double zTMax = zT + dZT;
      if (zTMin >= setup->maxZT() || zTMax < -setup->maxZT())
        // did not pass "eta" cut
        valid_ = false;
      else {
        zTMin = max(zTMin, -setup->maxZT());
        zTMax = min(zTMax, setup->maxZT());
      }
      // range of stub cot(theta)
      cot_ = {zTMin / setup->chosenRofZ(), zTMax / setup->chosenRofZ()};
    }

    // stub r w.r.t. chosenRofPhi in cm
    static const double chosenRofPhi = hybrid_ ? setup->hybridChosenRofPhi() : setup->chosenRofPhi();
    r_ = digi(r_ - chosenRofPhi, setup->tmttBaseR());

    // radial (cylindrical) component of sensor separation
    const double dr = sm->sep() / (sm->cosTilt() - sm->sinTilt() * z_ / d_);
    // converts bend into inv2R in 1/cm
    const double inv2ROverBend = sm->pitchRow() / dr / d_;
    // inv2R in 1/cm
    const double inv2R = -bend_ * setup->baseBend() * inv2ROverBend;
    // inv2R uncertainty in 1/cm
    const double dInv2R = setup->bendCut() * inv2ROverBend;
    const double minPt = hybrid_ ? setup->hybridMinPtStub() : setup->minPt();
    const double maxInv2R = setup->invPtToDphi() / minPt - setup->dtcBaseInv2R() / 2.;
    double inv2RMin = digi(inv2R - dInv2R, setup->dtcBaseInv2R());
    double inv2RMax = digi(inv2R + dInv2R, setup->dtcBaseInv2R());
    if (inv2RMin > maxInv2R || inv2RMax < -maxInv2R) {
      // did not pass pt cut
      valid_ = false;
    } else {
      inv2RMin = max(inv2RMin, -maxInv2R);
      inv2RMax = min(inv2RMax, maxInv2R);
    }
    // range of stub inv2R in 1/cm
    inv2R_ = {inv2RMin, inv2RMax};

    // stub phi w.r.t. detector region centre in rad
    phi_ = c_ + rowSub_ * m_;

    // range of stub extrapolated phi to radius chosenRofPhi in rad
    phiT_.first = phi_ - r_ * inv2R_.first;
    phiT_.second = phi_ - r_ * inv2R_.second;
    if (phiT_.first > phiT_.second)
      swap(phiT_.first, phiT_.second);

    if (phiT_.first < 0.)
      regions_.push_back(0);
    if (phiT_.second >= 0.)
      regions_.push_back(1);

    // apply data format specific manipulations
    if (!hybrid_)
      return;

    // stub r w.r.t. an offset in cm
    r_ -= sm->offsetR() - chosenRofPhi;
    // stub z w.r.t. an offset in cm
    z_ -= sm->offsetZ();
    if (sm->type() == SensorModule::Disk2S) {
      // encoded r
      r_ = sm->encodedR() + (sm->side() ? -col_ : (col_ + sm->numColumns() / 2));
      r_ = (r_ + 0.5) * setup->hybridBaseR(sm->type());
    }

    // encode bend
    const vector<double>& encodingBend = setup->encodingBend(sm->windowSize(), sm->psModule());
    const auto pos = find(encodingBend.begin(), encodingBend.end(), abs(ttStubRef->bendBE()));
    const int uBend = distance(encodingBend.begin(), pos);
    bend_ = pow(-1, signbit(bend_)) * uBend;
  }

~Stub()

trackerDTC::Stub::~Stub ( )

inline

Definition at line 21 of file Stub.h.

{}

Member Function Documentation

bend()

int trackerDTC::Stub::bend ( ) const

inline

Definition at line 27 of file Stub.h.

References bend_.

Referenced by trackerDTC::DTC::DTC().

{ return bend_; }

digi()

double trackerDTC::Stub::digi ( double  value, double  precision  ) const

private

Definition at line 156 of file Stub.cc.

References hcalRecHitTable_cff::precision.

Referenced by Stub().

{ return (floor(value / precision) + .5) * precision; }

formatHybrid()

Frame trackerDTC::Stub::formatHybrid ( int  region ) const

private

Definition at line 159 of file Stub.cc.

References tt::Setup::baseRegion(), bend_, trackerDTC::LayerEncoding::decode(), l1trig_cff::hwPhi, tt::Setup::hybridBaseAlpha(), tt::Setup::hybridBasePhi(), tt::Setup::hybridBaseR(), tt::Setup::hybridBaseZ(), tt::Setup::hybridNumUnusedBits(), tt::Setup::hybridRangePhi(), tt::Setup::hybridWidthAlpha(), tt::Setup::hybridWidthBend(), tt::Setup::hybridWidthLayerId(), tt::Setup::hybridWidthPhi(), tt::Setup::hybridWidthR(), tt::Setup::hybridWidthZ(), layerEncoding_, phi_, r_, nano_mu_digi_cff::region, row_, setup_, sm_, tt::SensorModule::type(), and z_.

Referenced by frame().

                                           {
    const SensorModule::Type type = sm_->type();
    // layer encoding
    const int decodedLayerId = layerEncoding_->decode(sm_);
    // stub phi w.r.t. processing region border in rad
    double phi = phi_ - (region - .5) * setup_->baseRegion() + setup_->hybridRangePhi() / 2.;
    if (phi >= setup_->hybridRangePhi())
      phi = setup_->hybridRangePhi() - setup_->hybridBasePhi(type) / 2.;
    // convert stub variables into bit vectors
    const TTBV hwR(r_, setup_->hybridBaseR(type), setup_->hybridWidthR(type), true);
    const TTBV hwPhi(phi, setup_->hybridBasePhi(type), setup_->hybridWidthPhi(type), true);
    const TTBV hwZ(z_, setup_->hybridBaseZ(type), setup_->hybridWidthZ(type), true);
    const TTBV hwAlpha(row_, setup_->hybridBaseAlpha(type), setup_->hybridWidthAlpha(type), true);
    const TTBV hwBend(bend_, setup_->hybridWidthBend(type), true);
    const TTBV hwLayer(decodedLayerId, setup_->hybridWidthLayerId());
    const TTBV hwGap(0, setup_->hybridNumUnusedBits(type));
    const TTBV hwValid(1, 1);
    // assemble final bitset
    return Frame(hwGap.str() + hwR.str() + hwZ.str() + hwPhi.str() + hwAlpha.str() + hwBend.str() + hwLayer.str() +
                 hwValid.str());
  }

formatTMTT()

Frame trackerDTC::Stub::formatTMTT ( int  region ) const

private

Definition at line 181 of file Stub.cc.

References tt::Setup::baseRegion(), tt::Setup::baseSector(), newFWLiteAna::bin, tt::Setup::boundarieEta(), cot_, l1trig_cff::hwPhi, inv2R_, tt::SensorModule::layerId(), tt::Setup::numOverlappingRegions(), tt::Setup::numSectorsEta(), tt::Setup::numSectorsPhi(), phi_, phiT_, r_, nano_mu_digi_cff::region, setup_, sm_, TTBV::str(), tt::Setup::tmttBaseInv2R(), tt::Setup::tmttBasePhi(), tt::Setup::tmttBaseR(), tt::Setup::tmttBaseZ(), tt::Setup::tmttNumUnusedBits(), tt::Setup::tmttWidthInv2R(), tt::Setup::tmttWidthLayer(), tt::Setup::tmttWidthPhi(), tt::Setup::tmttWidthR(), tt::Setup::tmttWidthSectorEta(), tt::Setup::tmttWidthZ(), and z_.

Referenced by frame().

                                         {
    int layerM = sm_->layerId();
    // convert unique layer id [1-6,11-15] into reduced layer id [0-6]
    // a fiducial track may not cross more then 7 detector layers, for stubs from a given track the reduced layer id is actually unique
    int layer(-1);
    if (layerM == 1)
      layer = 0;
    else if (layerM == 2)
      layer = 1;
    else if (layerM == 6 || layerM == 11)
      layer = 2;
    else if (layerM == 5 || layerM == 12)
      layer = 3;
    else if (layerM == 4 || layerM == 13)
      layer = 4;
    else if (layerM == 14)
      layer = 5;
    else if (layerM == 3 || layerM == 15)
      layer = 6;
    // assign stub to phi sectors within a processing region, to be generalized
    TTBV sectorsPhi(0, setup_->numOverlappingRegions() * setup_->numSectorsPhi());
    if (phiT_.first < 0.) {
      if (phiT_.first < -setup_->baseSector())
        sectorsPhi.set(0);
      else
        sectorsPhi.set(1);
      if (phiT_.second < 0. && phiT_.second >= -setup_->baseSector())
        sectorsPhi.set(1);
    }
    if (phiT_.second >= 0.) {
      if (phiT_.second < setup_->baseSector())
        sectorsPhi.set(2);
      else
        sectorsPhi.set(3);
      if (phiT_.first >= 0. && phiT_.first < setup_->baseSector())
        sectorsPhi.set(2);
    }
    // assign stub to eta sectors within a processing region
    pair<int, int> sectorEta({0, setup_->numSectorsEta() - 1});
    for (int bin = 0; bin < setup_->numSectorsEta(); bin++)
      if (asinh(cot_.first) < setup_->boundarieEta(bin + 1)) {
        sectorEta.first = bin;
        break;
      }
    for (int bin = sectorEta.first; bin < setup_->numSectorsEta(); bin++)
      if (asinh(cot_.second) < setup_->boundarieEta(bin + 1)) {
        sectorEta.second = bin;
        break;
      }
    // stub phi w.r.t. processing region centre in rad
    const double phi = phi_ - (region - .5) * setup_->baseRegion();
    // convert stub variables into bit vectors
    const TTBV hwValid(1, 1);
    const TTBV hwGap(0, setup_->tmttNumUnusedBits());
    const TTBV hwLayer(layer, setup_->tmttWidthLayer());
    const TTBV hwSectorEtaMin(sectorEta.first, setup_->tmttWidthSectorEta());
    const TTBV hwSectorEtaMax(sectorEta.second, setup_->tmttWidthSectorEta());
    const TTBV hwR(r_, setup_->tmttBaseR(), setup_->tmttWidthR(), true);
    const TTBV hwPhi(phi, setup_->tmttBasePhi(), setup_->tmttWidthPhi(), true);
    const TTBV hwZ(z_, setup_->tmttBaseZ(), setup_->tmttWidthZ(), true);
    const TTBV hwInv2RMin(inv2R_.first, setup_->tmttBaseInv2R(), setup_->tmttWidthInv2R(), true);
    const TTBV hwInv2RMax(inv2R_.second, setup_->tmttBaseInv2R(), setup_->tmttWidthInv2R(), true);
    TTBV hwSectorPhis(0, setup_->numSectorsPhi());
    for (int sectorPhi = 0; sectorPhi < setup_->numSectorsPhi(); sectorPhi++)
      hwSectorPhis[sectorPhi] = sectorsPhi[region * setup_->numSectorsPhi() + sectorPhi];
    // assemble final bitset
    return Frame(hwGap.str() + hwValid.str() + hwR.str() + hwPhi.str() + hwZ.str() + hwLayer.str() +
                 hwSectorPhis.str() + hwSectorEtaMin.str() + hwSectorEtaMax.str() + hwInv2RMin.str() +
                 hwInv2RMax.str());
  }

frame()

Frame trackerDTC::Stub::frame

(

int

region

)

const

Definition at line 150 of file Stub.cc.

References formatHybrid(), formatTMTT(), hybrid_, and nano_mu_digi_cff::region.

{ return hybrid_ ? formatHybrid(region) : formatTMTT(region); }

inRegion()

bool trackerDTC::Stub::inRegion

(

int

region

)

const

Definition at line 153 of file Stub.cc.

References spr::find(), nano_mu_digi_cff::region, and regions_.

{ return find(regions_.begin(), regions_.end(), region) != regions_.end(); }

ttStubRef()

TTStubRef trackerDTC::Stub::ttStubRef ( ) const

inline

Definition at line 23 of file Stub.h.

References ttStubRef_.

Referenced by Stub().

{ return ttStubRef_; }

valid()

bool trackerDTC::Stub::valid ( ) const

inline

Definition at line 25 of file Stub.h.

References valid_.

Referenced by trackerDTC::DTC::DTC().

{ return valid_; }

Member Data Documentation

bend_

int trackerDTC::Stub::bend_

private

Definition at line 57 of file Stub.h.

Referenced by bend(), formatHybrid(), and Stub().

c_

double trackerDTC::Stub::c_

private

Definition at line 71 of file Stub.h.

Referenced by Stub().

col_

int trackerDTC::Stub::col_

private

Definition at line 53 of file Stub.h.

Referenced by Stub().

cot_

std::pair<double, double> trackerDTC::Stub::cot_

private

Definition at line 77 of file Stub.h.

Referenced by formatTMTT(), and Stub().

d_

double trackerDTC::Stub::d_

private

Definition at line 73 of file Stub.h.

Referenced by Stub().

hybrid_

bool trackerDTC::Stub::hybrid_

private

Definition at line 49 of file Stub.h.

Referenced by frame(), and Stub().

inv2R_

std::pair<double, double> trackerDTC::Stub::inv2R_

private

Definition at line 75 of file Stub.h.

Referenced by formatTMTT(), and Stub().

layerEncoding_

const LayerEncoding* trackerDTC::Stub::layerEncoding_

private

Definition at line 43 of file Stub.h.

Referenced by formatHybrid().

m_

double trackerDTC::Stub::m_

private

Definition at line 69 of file Stub.h.

Referenced by Stub().

phi_

double trackerDTC::Stub::phi_

private

Definition at line 65 of file Stub.h.

Referenced by formatHybrid(), formatTMTT(), and Stub().

phiT_

std::pair<double, double> trackerDTC::Stub::phiT_

private

Definition at line 79 of file Stub.h.

Referenced by formatTMTT(), and Stub().

r_

double trackerDTC::Stub::r_

private

Definition at line 63 of file Stub.h.

Referenced by formatHybrid(), formatTMTT(), and Stub().

regions_

std::vector<int> trackerDTC::Stub::regions_

private

Definition at line 81 of file Stub.h.

Referenced by inRegion(), and Stub().

row_

int trackerDTC::Stub::row_

private

Definition at line 55 of file Stub.h.

Referenced by formatHybrid(), and Stub().

rowLUT_

int trackerDTC::Stub::rowLUT_

private

Definition at line 59 of file Stub.h.

Referenced by Stub().

rowSub_

int trackerDTC::Stub::rowSub_

private

Definition at line 61 of file Stub.h.

Referenced by Stub().

setup_

const tt::Setup* trackerDTC::Stub::setup_

private

Definition at line 41 of file Stub.h.

Referenced by formatHybrid(), formatTMTT(), and upgradeWorkflowComponents.UpgradeWorkflow::setup().

sm_

tt::SensorModule* trackerDTC::Stub::sm_

private

Definition at line 45 of file Stub.h.

Referenced by formatHybrid(), and formatTMTT().

ttStubRef_

TTStubRef trackerDTC::Stub::ttStubRef_

private

Definition at line 47 of file Stub.h.

Referenced by ttStubRef().

valid_

bool trackerDTC::Stub::valid_

private

Definition at line 51 of file Stub.h.

Referenced by Stub(), and valid().

z_

double trackerDTC::Stub::z_

private

Definition at line 67 of file Stub.h.

Referenced by formatHybrid(), formatTMTT(), and Stub().