trackerDTC::Stub Class Reference

#include <Stub.h>

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

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

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

Detailed Description

Definition at line 12 of file Stub.h.

Constructor & Destructor Documentation

Stub()

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

Definition at line 13 of file Stub.cc.

      : setup_(&setup), 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->sin();
    // stub z in cm
    z_ = digi(sm->z() + y * sm->cos(), setup.baseZ());
 
    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.basePhi());
    // 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.baseR());
 
    // radial (cylindrical) component of sensor separation
    const double dr = sm->sep() / (sm->cos() - sm->sin() * z_ / d_);
    // converts bend into qOverPt in 1/cm
    const double qOverPtOverBend = sm->pitchRow() / dr / d_;
    // qOverPt in 1/cm
    const double qOverPt = bend_ * setup.baseBend() * qOverPtOverBend;
    // qOverPt uncertainty in 1/cm
    const double dQoverPt = setup.bendCut() * qOverPtOverBend;
    const double minPt = hybrid_ ? setup.hybridMinPt() : setup.minPt();
    const double maxQoverPt = setup.invPtToDphi() / minPt - setup.dtcBaseQoverPt() / 2.;
    double qOverPtMin = digi(qOverPt - dQoverPt, setup.dtcBaseQoverPt());
    double qOverPtMax = digi(qOverPt + dQoverPt, setup.dtcBaseQoverPt());
    if (qOverPtMin >= maxQoverPt || qOverPtMax < -maxQoverPt)
      // did not pass pt cut
      valid_ = false;
    else {
      qOverPtMin = max(qOverPtMin, -maxQoverPt);
      qOverPtMax = min(qOverPtMax, maxQoverPt);
    }
    // range of stub qOverPt in 1/cm
    qOverPt_ = {qOverPtMin, qOverPtMax};
 
    // 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_ * qOverPt_.first;
    phiT_.second = phi_ + r_ * qOverPt_.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;
  }

References funct::abs(), bend_, c, c_, col_, trackerDTC::SensorModule::cos(), cot_, d_, digi(), trackerDTC::SensorModule::Disk2S, HLT_FULL_cff::distance, flavorHistoryFilter_cfi::dr, trackerDTC::SensorModule::encodedR(), spr::find(), hybrid_, createfilelist::int, visualization-live-secondInstance_cfg::m, m_, SiStripPI::max, min(), beam_dqm_sourceclient-live_cfg::minPt, trackerDTC::SensorModule::numColumns(), trackerDTC::SensorModule::numRows(), trackerDTC::SensorModule::offsetR(), trackerDTC::SensorModule::offsetZ(), trackerDTC::SensorModule::phi(), phi_, phiT_, trackerDTC::SensorModule::pitchCol(), trackerDTC::SensorModule::pitchRow(), funct::pow(), trackerDTC::SensorModule::psModule(), qOverPt_, trackerDTC::SensorModule::r(), r_, regions_, row_, rowLUT_, rowSub_, trackerDTC::SensorModule::sep(), singleTopDQM_cfi::setup, trackerDTC::SensorModule::side(), trackerDTC::SensorModule::signBend(), trackerDTC::SensorModule::signCol(), trackerDTC::SensorModule::signRow(), trackerDTC::SensorModule::sin(), mathSSE::sqrt(), edm::swap(), ttStubRef(), trackerDTC::SensorModule::type(), valid_, trackerDTC::SensorModule::windowSize(), PV2DBase< T, PVType, FrameType >::x(), testProducerWithPsetDescEmpty_cfi::x1, PV2DBase< T, PVType, FrameType >::y(), trackerDTC::SensorModule::z(), and z_.

~Stub()

trackerDTC::Stub::~Stub ( )

inline

Definition at line 15 of file Stub.h.

{}

Member Function Documentation

bend()

int trackerDTC::Stub::bend ( ) const

inline

Definition at line 22 of file Stub.h.

{ return bend_; }

References bend_.

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

digi()

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

private

Definition at line 146 of file Stub.cc.

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

References boostedTaus_cff::precision.

Referenced by Stub().

formatHybrid()

TTDTC::BV trackerDTC::Stub::formatHybrid ( int  region ) const

private

Definition at line 149 of file Stub.cc.

                                             {
    const SensorModule::Type type = sm_->type();
    // stub phi w.r.t. processing region centre in rad
    const double phi = phi_ - (region - .5) * setup_->baseRegion() +
                       0.5 * setup_->hybridBasePhi(type) * (1 << setup_->hybridWidthPhi(type));
 
    // 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(sm_->encodedLayerId(), setup_->hybridWidthLayer());
    const TTBV hwGap(0, setup_->hybridNumUnusedBits(type));
    const TTBV hwValid(1, 1);
    // assemble final bitset
    return TTDTC::BV(hwGap.str() + hwR.str() + hwZ.str() + hwPhi.str() + hwAlpha.str() + hwBend.str() + hwLayer.str() +
                     hwValid.str());
  }

References trackerDTC::Setup::baseRegion(), bend_, trackerDTC::SensorModule::encodedLayerId(), trackerDTC::Setup::hybridBaseAlpha(), trackerDTC::Setup::hybridBasePhi(), trackerDTC::Setup::hybridBaseR(), trackerDTC::Setup::hybridBaseZ(), trackerDTC::Setup::hybridNumUnusedBits(), trackerDTC::Setup::hybridWidthAlpha(), trackerDTC::Setup::hybridWidthBend(), trackerDTC::Setup::hybridWidthLayer(), trackerDTC::Setup::hybridWidthPhi(), trackerDTC::Setup::hybridWidthR(), trackerDTC::Setup::hybridWidthZ(), phi_, r_, HLT_FULL_cff::region, row_, setup_, sm_, trackerDTC::SensorModule::type(), and z_.

Referenced by frame().

formatTMTT()

TTDTC::BV trackerDTC::Stub::formatTMTT ( int  region ) const

private

Definition at line 169 of file Stub.cc.

                                           {
    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> setcorEta({0, setup_->numSectorsEta() - 1});
    for (int bin = 0; bin < setup_->numSectorsEta(); bin++)
      if (asinh(cot_.first) < setup_->boundarieEta(bin + 1)) {
        setcorEta.first = bin;
        break;
      }
    for (int bin = setcorEta.first; bin < setup_->numSectorsEta(); bin++)
      if (asinh(cot_.second) < setup_->boundarieEta(bin + 1)) {
        setcorEta.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_->dtcNumUnusedBits());
    const TTBV hwLayer(layer, setup_->widthLayer());
    const TTBV hwSectorEtaMin(setcorEta.first, setup_->widthSectorEta());
    const TTBV hwSectorEtaMax(setcorEta.second, setup_->widthSectorEta());
    const TTBV hwR(r_, setup_->baseR(), setup_->widthR(), true);
    const TTBV hwPhi(phi, setup_->basePhi(), setup_->widthPhiDTC(), true);
    const TTBV hwZ(z_, setup_->baseZ(), setup_->widthZ(), true);
    const TTBV hwQoverPtMin(qOverPt_.first, setup_->htBaseQoverPt(), setup_->htWidthQoverPt(), true);
    const TTBV hwQoverPtMax(qOverPt_.second, setup_->htBaseQoverPt(), setup_->htWidthQoverPt(), true);
    TTBV hwSectorPhis(0, setup_->numSectorsPhi());
    for (int sectorPhi = 0; sectorPhi < setup_->numSectorsPhi(); sectorPhi++)
      hwSectorPhis[sectorPhi] = sectorsPhi[region * setup_->numSectorsPhi() + sectorPhi];
    // assemble final bitsetTTDTC::BV(hwGap.str() + hwValid.str() + hwR.str() + hwPhi.str() + hwZ.str() + hwQoverPtMin.str() +
    return TTDTC::BV(hwGap.str() + hwValid.str() + hwR.str() + hwPhi.str() + hwZ.str() + hwQoverPtMin.str() +
                     hwQoverPtMax.str() + hwSectorEtaMin.str() + hwSectorEtaMax.str() + hwSectorPhis.str() +
                     hwLayer.str());
  }

References trackerDTC::Setup::basePhi(), trackerDTC::Setup::baseR(), trackerDTC::Setup::baseRegion(), trackerDTC::Setup::baseSector(), trackerDTC::Setup::baseZ(), newFWLiteAna::bin, trackerDTC::Setup::boundarieEta(), cot_, trackerDTC::Setup::dtcNumUnusedBits(), trackerDTC::Setup::htBaseQoverPt(), trackerDTC::Setup::htWidthQoverPt(), phase1PixelTopology::layer, trackerDTC::SensorModule::layerId(), trackerDTC::Setup::numOverlappingRegions(), trackerDTC::Setup::numSectorsEta(), trackerDTC::Setup::numSectorsPhi(), phi_, phiT_, qOverPt_, r_, HLT_FULL_cff::region, TTBV::set(), setup_, sm_, TTBV::str(), trackerDTC::Setup::widthLayer(), trackerDTC::Setup::widthPhiDTC(), trackerDTC::Setup::widthR(), trackerDTC::Setup::widthSectorEta(), trackerDTC::Setup::widthZ(), and z_.

Referenced by frame().

frame()

TTDTC::BV trackerDTC::Stub::frame

(

int

region

)

const

Definition at line 140 of file Stub.cc.

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

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

inRegion()

bool trackerDTC::Stub::inRegion

(

int

region

)

const

Definition at line 143 of file Stub.cc.

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

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

ttStubRef()

TTStubRef trackerDTC::Stub::ttStubRef ( ) const

inline

Definition at line 18 of file Stub.h.

{ return ttStubRef_; }

References ttStubRef_.

Referenced by Stub().

valid()

bool trackerDTC::Stub::valid ( ) const

inline

Definition at line 20 of file Stub.h.

{ return valid_; }

References valid_.

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

Member Data Documentation

bend_

int trackerDTC::Stub::bend_

private

Definition at line 51 of file Stub.h.

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

c_

double trackerDTC::Stub::c_

private

Definition at line 65 of file Stub.h.

Referenced by Stub().

col_

int trackerDTC::Stub::col_

private

Definition at line 47 of file Stub.h.

Referenced by Stub().

cot_

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

private

Definition at line 71 of file Stub.h.

Referenced by formatTMTT(), and Stub().

d_

double trackerDTC::Stub::d_

private

Definition at line 67 of file Stub.h.

Referenced by Stub().

hybrid_

bool trackerDTC::Stub::hybrid_

private

Definition at line 43 of file Stub.h.

Referenced by frame(), and Stub().

m_

double trackerDTC::Stub::m_

private

Definition at line 63 of file Stub.h.

Referenced by Stub().

phi_

double trackerDTC::Stub::phi_

private

Definition at line 59 of file Stub.h.

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

phiT_

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

private

Definition at line 73 of file Stub.h.

Referenced by formatTMTT(), and Stub().

qOverPt_

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

private

Definition at line 69 of file Stub.h.

Referenced by formatTMTT(), and Stub().

r_

double trackerDTC::Stub::r_

private

Definition at line 57 of file Stub.h.

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

regions_

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

private

Definition at line 75 of file Stub.h.

Referenced by inRegion(), and Stub().

row_

int trackerDTC::Stub::row_

private

Definition at line 49 of file Stub.h.

Referenced by formatHybrid(), and Stub().

rowLUT_

int trackerDTC::Stub::rowLUT_

private

Definition at line 53 of file Stub.h.

Referenced by Stub().

rowSub_

int trackerDTC::Stub::rowSub_

private

Definition at line 55 of file Stub.h.

Referenced by Stub().

setup_

const Setup* trackerDTC::Stub::setup_

private

Definition at line 37 of file Stub.h.

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

sm_

SensorModule* trackerDTC::Stub::sm_

private

Definition at line 39 of file Stub.h.

Referenced by formatHybrid(), and formatTMTT().

ttStubRef_

TTStubRef trackerDTC::Stub::ttStubRef_

private

Definition at line 41 of file Stub.h.

Referenced by ttStubRef().

valid_

bool trackerDTC::Stub::valid_

private

Definition at line 45 of file Stub.h.

Referenced by Stub(), and valid().

z_

double trackerDTC::Stub::z_

private

Definition at line 61 of file Stub.h.

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