tmtt::Stub Class Reference

#include <Stub.h>

Public Types
enum	DigiStage {   DigiStage::NONE, DigiStage::GP, DigiStage::HT, DigiStage::SF,   DigiStage::TF }

Public Member Functions
float	alpha () const
const TP *	assocTP () const
std::array< const TP *, 2 >	assocTPofCluster () const
const std::set< const TP * > &	assocTPs () const
bool	barrel () const
float	bend () const
float	bendCut () const
float	bendCutInFrontend () const
float	bendInFrontend () const
void	calcQoverPtrange ()
const DigitalStub *	digitalStub () const
void	digitize (unsigned int iPhiSec, DigiStage digiStep)
float	dphiOverBend () const
float	eta () const
void	fillTruth (const std::map< edm::Ptr< TrackingParticle >, const TP *> &translateTP, const edm::Handle< TTStubAssMap > &mcTruthTTStubHandle, const edm::Handle< TTClusterAssMap > &mcTruthTTClusterHandle)
bool	frontendPass () const
bool	genuine () const
std::array< bool, 2 >	genuineCluster () const
unsigned int	index () const
unsigned int	iphi () const
unsigned int	layerId () const
unsigned int	layerIdReduced () const
std::array< float, 2 >	localU_cluster () const
std::array< float, 2 >	localV_cluster () const
unsigned int	max_qOverPt_bin () const
unsigned int	min_qOverPt_bin () const
unsigned int	nStrips () const
float	numMergedBend () const
bool	operator== (const Stub &stubOther)
float	phi () const
float	phiDiff (float rad, float Pt) const
bool	psModule () const
float	qOverPt () const
float	qOverPtcut () const
float	qOverPtOverBend () const
float	r () const
void	setDigitizeWarningsOn (bool newVal)
float	sigmaPar () const
float	sigmaPerp () const
float	sigmaR () const
float	sigmaZ () const
float	stripLength () const
float	stripPitch () const
	Stub (const Settings *settings, unsigned int idStub, double phi, double r, double z, double bend, unsigned int iphi, double alpha, unsigned int layerId, unsigned int iPhiSec, bool psModule, bool barrel, bool tiltedBarrel, float stripPitch, float stripLength, unsigned int nStrips)
	Stub (const TTStubRef &ttStubRef, unsigned int index_in_vStubs, const Settings *settings, const TrackerTopology *trackerTopology, const TrackerModule *trackerModule, const DegradeBend *degradeBend, const StubKiller *stubKiller)
bool	stubFailedDegradeWindow () const
float	theta () const
float	tiltAngle () const
bool	tiltedBarrel () const
const TrackerModule *	trackerModule () const
float	trkPhiAtR (float rad) const
float	trkPhiAtRcut (float rad) const
const TTStubRef &	ttStubRef () const
float	z () const

Private Member Functions
double	approxB ()
void	calcDphiOverBend ()
void	degradeResolution (float bend, float &degradedBend, unsigned int &num) const
void	setFrontend (const StubKiller *stubKiller)
void	setTrackerModule (const TrackerGeometry *trackerGeometry, const TrackerTopology *trackerTopology, const DetId &detId)

Private Attributes
float	alpha_
const TP *	assocTP_
std::array< const TP *, 2 >	assocTPofCluster_
std::set< const TP * >	assocTPs_
bool	barrel_
float	bend_
float	bendInFrontend_
const DegradeBend *	degradeBend_
std::unique_ptr< DigitalStub >	digitalStub_
bool	digitizeWarningsOn_
float	dphiOverBend_
bool	frontendPass_
unsigned int	index_in_vStubs_
const float	invRoot12 = sqrt(1. / 12.)
unsigned int	iphi_
DigiStage	lastDigiStep_
unsigned int	layerId_
unsigned int	layerIdReduced_
std::array< float, 2 >	localU_cluster_
std::array< float, 2 >	localV_cluster_
unsigned int	max_qOverPt_bin_
unsigned int	min_qOverPt_bin_
unsigned int	nStrips_
unsigned int	numMergedBend_
float	phi_
bool	psModule_
float	r_
const float	rejectedStubBend_ = 99999.
const Settings *	settings_
float	stripLength_
float	stripPitch_
bool	stubFailedDegradeWindow_
float	tiltAngle_
bool	tiltedBarrel_
const TrackerModule *	trackerModule_
TTStubRef	ttStubRef_
float	z_

Detailed Description

Definition at line 43 of file Stub.h.

Member Enumeration Documentation

DigiStage

enum tmtt::Stub::DigiStage

strong

Enumerator
NONE
GP
HT
SF
TF

Definition at line 89 of file Stub.h.

{ NONE, GP, HT, SF, TF };

Constructor & Destructor Documentation

Stub() [1/2]

tmtt::Stub::Stub	(	const Settings *	settings,
		unsigned int	idStub,
		double	phi,
		double	r,
		double	z,
		double	bend,
		unsigned int	iphi,
		double	alpha,
		unsigned int	layerId,
		unsigned int	iPhiSec,
		bool	psModule,
		bool	barrel,
		bool	tiltedBarrel,
		float	stripPitch,
		float	stripLength,
		unsigned int	nStrips
	)

Definition at line 21 of file Stub.cc.

      : index_in_vStubs_(idStub),  // A unique ID to label the stub.
        phi_(phi),
        r_(r),
        z_(z),
        bend_(bend),
        iphi_(iphi),
        alpha_(alpha),
        digitalStub_(std::make_unique<DigitalStub>(settings, r, phi, z, iPhiSec)),
        layerId_(layerId),
        layerIdReduced_(TrackerModule::calcLayerIdReduced(layerId)),
        stripPitch_(stripPitch),
        stripLength_(stripLength),
        nStrips_(nStrips),
        psModule_(psModule),
        barrel_(barrel),
        tiltedBarrel_(tiltedBarrel) {}

Stub() [2/2]

tmtt::Stub::Stub	(	const TTStubRef &	ttStubRef,
		unsigned int	index_in_vStubs,
		const Settings *	settings,
		const TrackerTopology *	trackerTopology,
		const TrackerModule *	trackerModule,
		const DegradeBend *	degradeBend,
		const StubKiller *	stubKiller
	)

Definition at line 56 of file Stub.cc.

References alpha_, assocTPofCluster_, barrel(), bend_, bendInFrontend_, calcDphiOverBend(), calcQoverPtrange(), TTStub< T >::clusterRef(), tmtt::Settings::degradeBendRes(), degradeResolution(), dqmMemoryStats::float, edm::Ref< C, T, F >::get(), iphi_, Topology::localPosition(), localU_cluster_, localV_cluster_, nStrips(), numMergedBend_, tmtt::TrackerModule::outerModuleAtSmallerR(), phi_, psModule(), r_, tmtt::TrackerModule::sensorWidth(), setFrontend(), tmtt::TrackerModule::specDet(), tmtt::TrackerModule::specTopol(), GeomDet::surface(), Surface::toGlobal(), trackerModule_, ttStubRef_, and z_.

      : ttStubRef_(ttStubRef),
        settings_(settings),
        index_in_vStubs_(index_in_vStubs),
        assocTP_(nullptr),  // Initialize in case job is using no MC truth info.
        lastDigiStep_(Stub::DigiStage::NONE),
        digitizeWarningsOn_(true),
        trackerModule_(trackerModule),  // Info about tracker module containing stub
        degradeBend_(degradeBend),      // Used to degrade stub bend information.
        // Module related variables (need to be stored for Hybrid)
        layerId_(trackerModule->layerId()),
        layerIdReduced_(trackerModule->layerIdReduced()),
        tiltAngle_(trackerModule->tiltAngle()),
        stripPitch_(trackerModule->stripPitch()),
        stripLength_(trackerModule->stripLength()),
        nStrips_(trackerModule->nStrips()),
        psModule_(trackerModule->psModule()),
        barrel_(trackerModule->barrel()),
        tiltedBarrel_(trackerModule->tiltedBarrel()) {
    // Get coordinates of stub.
    const TTStub<Ref_Phase2TrackerDigi_>* ttStubP = ttStubRef_.get();

    const PixelGeomDetUnit* specDet = trackerModule_->specDet();
    const PixelTopology* specTopol = trackerModule_->specTopol();
    MeasurementPoint measurementPoint = ttStubRef_->clusterRef(0)->findAverageLocalCoordinatesCentered();
    LocalPoint clustlp = specTopol->localPosition(measurementPoint);
    GlobalPoint pos = specDet->surface().toGlobal(clustlp);

    phi_ = pos.phi();
    r_ = pos.perp();
    z_ = pos.z();

    // Get the coordinates of the two clusters that make up this stub, measured in units of strip pitch, and measured
    // in the local frame of the sensor. They have a granularity  of 0.5*pitch.
    for (unsigned int iClus = 0; iClus <= 1; iClus++) {  // Loop over two clusters in stub.
      localU_cluster_[iClus] = ttStubP->clusterRef(iClus)->findAverageLocalCoordinatesCentered().x();
      localV_cluster_[iClus] = ttStubP->clusterRef(iClus)->findAverageLocalCoordinatesCentered().y();
    }

    // Get location of stub in module in units of strip number (or pixel number along finest granularity axis).
    // Range from 0 to (nStrips - 1) inclusive.
    // N.B. Since iphi is integer, this degrades the granularity by a factor 2. This seems silly, but track fit wants it.
    iphi_ = localU_cluster_[0];  // granularity 1*strip (unclear why we want to degrade it ...)

    // Determine alpha correction for non-radial strips in endcap 2S modules.
    // (If true hit at larger r than stub r by deltaR, then stub phi needs correcting by +alpha*deltaR).
    alpha_ = 0.;
    if ((not barrel()) && (not psModule())) {
      float fracPosInModule = (float(iphi_) - 0.5 * float(nStrips())) / float(nStrips());
      float phiRelToModule = trackerModule_->sensorWidth() * fracPosInModule / r_;
      if (z_ < 0)
        phiRelToModule *= -1;
      if (trackerModule_->outerModuleAtSmallerR())
        phiRelToModule *= -1;  // Module flipped.
      // If true hit at larger r than stub r by deltaR, then stub phi needs correcting by +alpha*deltaR.
      alpha_ = -phiRelToModule / r_;
    }

    // Calculate variables giving ratio of track intercept angle to stub bend.
    this->calcDphiOverBend();

    // Get stub bend that is available in front-end electronics, where bend is displacement between
    // two hits in stubs in units of strip pitch.
    bendInFrontend_ = ttStubRef_->bendFE();
    if ((not barrel()) && pos.z() > 0)
      bendInFrontend_ *= -1;
    if (barrel())
      bendInFrontend_ *= -1;

    // Get stub bend that is available in off-detector electronics, allowing for degredation of
    // bend resolution due to bit encoding by FE chip if required.
    numMergedBend_ = 1;  // Number of bend values merged into single degraded one.
    if (settings->degradeBendRes() == 2) {
      float degradedBend;  // degraded bend
      // This returns values of degradedBend & numMergedBend_
      this->degradeResolution(bendInFrontend_, degradedBend, numMergedBend_);
      bend_ = degradedBend;
    } else if (settings->degradeBendRes() == 1) {
      bend_ = ttStubRef_->bendBE();  // Degraded bend from official CMS recipe.
      if ((not barrel()) && pos.z() > 0)
        bend_ *= -1;
      if (barrel())
        bend_ *= -1;
    } else {
      bend_ = bendInFrontend_;
    }

    // Fill frontendPass_ flag, indicating if frontend readout electronics will output this stub.
    this->setFrontend(stubKiller);

    // Calculate bin range along q/Pt axis of r-phi Hough transform array consistent with bend of this stub.
    this->calcQoverPtrange();

    // Initialize truth info to false in case job is using no MC truth info.
    for (unsigned int iClus = 0; iClus <= 1; iClus++) {
      assocTPofCluster_[iClus] = nullptr;
    }
  }

Member Function Documentation

alpha()

float tmtt::Stub::alpha ( void  ) const

inline

Definition at line 119 of file Stub.h.

References alpha_.

Referenced by tmtt::KFParamsComb::matrixV(), and tmtt::KFbase::residual().

{ return alpha_; }

approxB()

double tmtt::Stub::approxB ( )

private

Definition at line 304 of file Stub.cc.

References funct::abs(), tmtt::Settings::bApprox_gradient(), tmtt::Settings::bApprox_intercept(), barrel(), r_, settings_, tiltedBarrel(), and z_.

Referenced by calcDphiOverBend().

                       {
    if (tiltedBarrel()) {
      return settings_->bApprox_gradient() * std::abs(z_) / r_ + settings_->bApprox_intercept();
    } else {
      return barrel() ? 1 : std::abs(z_) / r_;
    }
  }

assocTP()

const TP* tmtt::Stub::assocTP ( ) const

inline

Definition at line 168 of file Stub.h.

References assocTP_.

                              {
      return assocTP_;
    }  // If only one TP contributed to both clusters, this tells you which TP it is. Returns nullptr if none.

assocTPofCluster()

std::array<const TP*, 2> tmtt::Stub::assocTPofCluster ( ) const

inline

Definition at line 176 of file Stub.h.

References assocTPofCluster_.

                                                    {
      return assocTPofCluster_;
    }  // Which TP made each cluster. Warning: If cluster was not produced by a single TP, then returns nullptr! (P.S. If both clusters match same TP, then this will equal assocTP()).

assocTPs()

const std::set<const TP*>& tmtt::Stub::assocTPs ( ) const

inline

Definition at line 164 of file Stub.h.

References assocTPs_.

Referenced by tmtt::KalmanState::good(), and tmtt::KFbase::printStub().

                                              {
      return assocTPs_;
    }  // Return TPs associated to this stub. (Whether only TPs contributing to both clusters are returned is determined by "StubMatchStrict" config param.)

barrel()

bool tmtt::Stub::barrel ( ) const

inline

Definition at line 201 of file Stub.h.

References barrel_.

Referenced by approxB(), tmtt::KalmanState::KalmanState(), tmtt::KFbase::kalmanUpdate(), tmtt::KFParamsComb::matrixV(), tmtt::KFbase::residual(), sigmaR(), sigmaZ(), Stub(), tmtt::TP::trkRAtStub(), and tmtt::TP::trkZAtStub().

{ return barrel_; }

bend()

float tmtt::Stub::bend ( ) const

inline

Definition at line 126 of file Stub.h.

References bend_.

Referenced by calcQoverPtrange(), degradeResolution(), qOverPt(), and setFrontend().

{ return bend_; }

bendCut()

float tmtt::Stub::bendCut ( ) const

inline

Definition at line 128 of file Stub.h.

References tmtt::Settings::bendCut(), tmtt::Settings::bendCutExtra(), numMergedBend_, and settings_.

Referenced by calcQoverPtrange(), qOverPtcut(), setFrontend(), and trkPhiAtRcut().

{ return (settings_->bendCut() + (numMergedBend_ - 1) * settings_->bendCutExtra()); }

bendCutInFrontend()

float tmtt::Stub::bendCutInFrontend ( ) const

inline

Definition at line 124 of file Stub.h.

References tmtt::Settings::bendCut(), and settings_.

Referenced by degradeResolution(), tmtt::StubWindowSuggest::process(), and setFrontend().

{ return settings_->bendCut(); }

bendInFrontend()

float tmtt::Stub::bendInFrontend ( ) const

inline

Definition at line 123 of file Stub.h.

References bendInFrontend_.

Referenced by setFrontend().

{ return bendInFrontend_; }

calcDphiOverBend()

void tmtt::Stub::calcDphiOverBend ( )

private

Definition at line 314 of file Stub.cc.

References funct::abs(), approxB(), funct::cos(), dphiOverBend_, tmtt::TrackerModule::pitchOverSep(), settings_, funct::sin(), theta(), tmtt::TrackerModule::tiltAngle(), trackerModule_, and tmtt::Settings::useApproxB().

Referenced by digitize(), and Stub().

                              {
    // Uses stub (r,z) instead of module (r,z). Logically correct but has negligable effect on results.
    if (settings_->useApproxB()) {
      float dphiOverBendCorrection_approx_ = approxB();
      dphiOverBend_ = trackerModule_->pitchOverSep() * dphiOverBendCorrection_approx_;
    } else {
      float dphiOverBendCorrection_ = std::abs(cos(this->theta() - trackerModule_->tiltAngle()) / sin(this->theta()));
      dphiOverBend_ = trackerModule_->pitchOverSep() * dphiOverBendCorrection_;
    }
  }

calcQoverPtrange()

void tmtt::Stub::calcQoverPtrange

(

)

Definition at line 163 of file Stub.cc.

References bend(), bendCut(), tmtt::Settings::houghMinPt(), tmtt::Settings::houghNbinsPt(), createfilelist::int, SiStripPI::max, max_qOverPt_bin_, SiStripPI::min, min_qOverPt_bin_, L1MuTriggerPtScale_cfi::nbinsPt, qOverPtOverBend(), settings_, tmtt::Settings::shape(), and createJobs::tmp.

Referenced by digitize(), and Stub().

                              {
    // First determine bin range along q/Pt axis of HT array
    // (Use "int" as nasty things happen if multiply "int" and "unsigned int").
    const int nbinsPt = (int)settings_->houghNbinsPt();
    const int min_array_bin = 0;
    const int max_array_bin = nbinsPt - 1;
    // Now calculate range of q/Pt bins allowed by bend filter.
    float qOverPtMin = this->qOverPtOverBend() * (this->bend() - this->bendCut());
    float qOverPtMax = this->qOverPtOverBend() * (this->bend() + this->bendCut());
    int houghNbinsPt = settings_->houghNbinsPt();
    const float houghMaxInvPt = 1. / settings_->houghMinPt();
    float qOverPtBinSize = (2. * houghMaxInvPt) / houghNbinsPt;
    if (settings_->shape() == 2 || settings_->shape() == 1 || settings_->shape() == 3)  // Non-square HT cells.
      qOverPtBinSize = 2. * houghMaxInvPt / (houghNbinsPt - 1);
    // Convert to bin number along q/Pt axis of HT array.
    // N.B. For square HT cells, setting "tmp = -0.5" causeas cell to be accepted if q/Pt at its centre is consistent
    // with the stub bend. Instead using "tmp = 0.0" accepts cells if q/Pt at any point in cell is consistent with bend.
    // So if you use change from -0.5 to 0.0, you have to tighten the bend cut (by ~0.05) to get similar performance.
    // Decision to set tmp = 0.0 taken in softare & GP firmware on 9th August 2016.

    float tmp = (settings_->shape() == 2 || settings_->shape() == 1 || settings_->shape() == 3) ? 1. : 0.;
    int min_bin = std::floor(-tmp + (qOverPtMin + houghMaxInvPt) / qOverPtBinSize);
    int max_bin = std::floor(tmp + (qOverPtMax + houghMaxInvPt) / qOverPtBinSize);

    // Limit it to range of HT array.
    min_bin = max(min_bin, min_array_bin);
    max_bin = min(max_bin, max_array_bin);
    // If min_bin > max_bin at this stage, it means that the Pt estimated from the bend is below the cutoff for track-finding.
    // Keep min_bin > max_bin, so such stubs can be rejected, but set both variables to values inside the HT bin range.
    if (min_bin > max_bin) {
      min_bin = max_array_bin;
      max_bin = min_array_bin;
    }
    min_qOverPt_bin_ = (unsigned int)min_bin;
    max_qOverPt_bin_ = (unsigned int)max_bin;
  }

degradeResolution()

void tmtt::Stub::degradeResolution ( float  bend, float &  degradedBend, unsigned int &  num  ) const

private

Definition at line 248 of file Stub.cc.

References funct::abs(), bend(), bendCutInFrontend(), tmtt::DegradeBend::degrade(), degradeBend_, tmtt::TrackerModule::detId(), tmtt::Settings::houghMinPt(), tmtt::Settings::killLowPtStubs(), EgammaValidation_cff::num, psModule(), qOverPtOverBend(), rejectedStubBend_, settings_, and trackerModule_.

Referenced by Stub().

                                                                                       {
    // If TMTT code is tightening official CMS FE stub window cuts, then calculate TMTT stub windows.
    float windowFE;
    if (settings_->killLowPtStubs()) {
      // Window size corresponding to Pt cut used for tracking.
      float invPtMax = 1. / (settings_->houghMinPt());
      windowFE = invPtMax / std::abs(this->qOverPtOverBend());
      // Increase half-indow size to allow for resolution in bend.
      windowFE += this->bendCutInFrontend();
    } else {
      windowFE = rejectedStubBend_;  // TMTT is not tightening windows.
    }

    degradeBend_->degrade(bend, psModule(), trackerModule_->detId(), windowFE, degradedBend, num);
  }

digitalStub()

const DigitalStub* tmtt::Stub::digitalStub ( ) const

inline

Definition at line 96 of file Stub.h.

References digitalStub_.

Referenced by tmtt::Sector::insideEtaSubSecs().

{ return digitalStub_.get(); }

digitize()

void tmtt::Stub::digitize	(	unsigned int	iPhiSec,
		Stub::DigiStage	digiStep
	)

Definition at line 203 of file Stub.cc.

References bend_, calcDphiOverBend(), calcQoverPtrange(), digitalStub_, digitizeWarningsOn_, tmtt::Settings::enableDigitize(), GP, HT, lastDigiStep_, max_qOverPt_bin_, min_qOverPt_bin_, phi_, r_, settings_, and z_.

                                                                  {
    if (settings_->enableDigitize()) {
      bool updated = true;
      if (not digitalStub_) {
        // Digitize stub if not yet done.
        digitalStub_ =
            std::make_unique<DigitalStub>(settings_, phi_, r_, z_, min_qOverPt_bin_, max_qOverPt_bin_, bend_, iPhiSec);
      } else {
        // If digitization already done, redo phi digi if phi sector has changed.
        updated = digitalStub_->changePhiSec(iPhiSec);
      }

      // Save CPU by only updating if something has changed.
      if (updated || digiStep != lastDigiStep_) {
        lastDigiStep_ = digiStep;

        // Replace stub coords with those degraded by digitization process.
        if (digiStep == DigiStage::GP) {
          phi_ = digitalStub_->phi_GP();
        } else {
          phi_ = digitalStub_->phi_HT_TF();
        }
        if (digiStep == DigiStage::GP || digiStep == DigiStage::HT) {
          r_ = digitalStub_->r_GP_HT();
        } else {
          r_ = digitalStub_->r_SF_TF();
        }
        z_ = digitalStub_->z();
        bend_ = digitalStub_->bend();

        // Update data members that depend on updated coords.
        // (Logically part of digitisation, so disable warnings)
        digitizeWarningsOn_ = false;
        if (digiStep == DigiStage::GP)
          this->calcDphiOverBend();
        if (digiStep == DigiStage::HT)
          this->calcQoverPtrange();
        digitizeWarningsOn_ = true;
      }
    }
  }

dphiOverBend()

float tmtt::Stub::dphiOverBend ( ) const

inline

Definition at line 146 of file Stub.h.

References dphiOverBend_.

{ return dphiOverBend_; }

eta()

float tmtt::Stub::eta ( void  ) const

inline

Definition at line 111 of file Stub.h.

References r_, and z_.

Referenced by Particle.Particle::__str__(), Jet.Jet::jetID(), Jet.Jet::puJetId(), and setFrontend().

{ return asinh(z_ / r_); }

fillTruth()

void tmtt::Stub::fillTruth	(	const std::map< edm::Ptr< TrackingParticle >, const TP *> &	translateTP,
		const edm::Handle< TTStubAssMap > &	mcTruthTTStubHandle,
		const edm::Handle< TTClusterAssMap > &	mcTruthTTClusterHandle
	)

Definition at line 328 of file Stub.cc.

References assocTP_, assocTPofCluster_, assocTPs_, TTClusterAssociationMap< T >::findTrackingParticlePtr(), TTStubAssociationMap< T >::findTrackingParticlePtr(), TTClusterAssociationMap< T >::findTrackingParticlePtrs(), genuine(), genuineCluster(), TTClusterAssociationMap< T >::isGenuine(), TTStubAssociationMap< T >::isGenuine(), settings_, tmtt::Settings::stubMatchStrict(), and ttStubRef_.

                                                                                 {
    //--- Fill assocTP_ info. If both clusters in this stub were produced by the same single tracking particle, find out which one it was.

    bool genuine = mcTruthTTStubHandle->isGenuine(ttStubRef_);  // Same TP contributed to both clusters?
    assocTP_ = nullptr;

    // Require same TP contributed to both clusters.
    if (genuine) {
      edm::Ptr<TrackingParticle> tpPtr = mcTruthTTStubHandle->findTrackingParticlePtr(ttStubRef_);
      if (translateTP.find(tpPtr) != translateTP.end()) {
        assocTP_ = translateTP.at(tpPtr);
        // N.B. Since not all tracking particles are stored in InputData::vTPs_, sometimes no match will be found.
      }
    }

    // Fill assocTPs_ info.

    if (settings_->stubMatchStrict()) {
      // We consider only stubs in which this TP contributed to both clusters.
      if (assocTP_ != nullptr)
        assocTPs_.insert(assocTP_);

    } else {
      // We consider stubs in which this TP contributed to either cluster.

      for (unsigned int iClus = 0; iClus <= 1; iClus++) {  // Loop over both clusters that make up stub.
        const TTClusterRef& ttClusterRef = ttStubRef_->clusterRef(iClus);

        // Now identify all TP's contributing to either cluster in stub.
        vector<edm::Ptr<TrackingParticle> > vecTpPtr = mcTruthTTClusterHandle->findTrackingParticlePtrs(ttClusterRef);

        for (const edm::Ptr<TrackingParticle>& tpPtr : vecTpPtr) {
          if (translateTP.find(tpPtr) != translateTP.end()) {
            assocTPs_.insert(translateTP.at(tpPtr));
            // N.B. Since not all tracking particles are stored in InputData::vTPs_, sometimes no match will be found.
          }
        }
      }
    }

    //--- Also note which tracking particles produced the two clusters that make up the stub

    for (unsigned int iClus = 0; iClus <= 1; iClus++) {  // Loop over both clusters that make up stub.
      const TTClusterRef& ttClusterRef = ttStubRef_->clusterRef(iClus);

      bool genuineCluster = mcTruthTTClusterHandle->isGenuine(ttClusterRef);  // Only 1 TP made cluster?
      assocTPofCluster_[iClus] = nullptr;

      // Only consider clusters produced by just one TP.
      if (genuineCluster) {
        edm::Ptr<TrackingParticle> tpPtr = mcTruthTTClusterHandle->findTrackingParticlePtr(ttClusterRef);

        if (translateTP.find(tpPtr) != translateTP.end()) {
          assocTPofCluster_[iClus] = translateTP.at(tpPtr);
          // N.B. Since not all tracking particles are stored in InputData::vTPs_, sometimes no match will be found.
        }
      }
    }
  }

frontendPass()

bool tmtt::Stub::frontendPass ( ) const

inline

Definition at line 157 of file Stub.h.

References frontendPass_.

{ return frontendPass_; }

genuine()

bool tmtt::Stub::genuine ( ) const

inline

Definition at line 167 of file Stub.h.

References assocTPs_.

Referenced by fillTruth().

{ return (not assocTPs_.empty()); }  // Did stub match at least one TP?

genuineCluster()

std::array<bool, 2> tmtt::Stub::genuineCluster ( ) const

inline

Definition at line 173 of file Stub.h.

References assocTPofCluster_.

Referenced by fillTruth().

                                             {
      return std::array<bool, 2>{{(assocTPofCluster_[0] != nullptr), (assocTPofCluster_[1] != nullptr)}};
    }  // Was cluster produced by a single TP?

index()

unsigned int tmtt::Stub::index ( ) const

inline

Definition at line 101 of file Stub.h.

References index_in_vStubs_.

Referenced by operator==(), and tmtt::KFbase::printStub().

{ return index_in_vStubs_; }

iphi()

unsigned int tmtt::Stub::iphi ( ) const

inline

Definition at line 115 of file Stub.h.

References iphi_.

{ return iphi_; }

layerId()

unsigned int tmtt::Stub::layerId ( ) const

inline

Definition at line 198 of file Stub.h.

References layerId_.

Referenced by tmtt::KFbase::printStub().

{ return layerId_; }

layerIdReduced()

unsigned int tmtt::Stub::layerIdReduced ( ) const

inline

Definition at line 200 of file Stub.h.

References layerIdReduced_.

{ return layerIdReduced_; }

localU_cluster()

std::array<float, 2> tmtt::Stub::localU_cluster ( ) const

inline

Definition at line 153 of file Stub.h.

References localU_cluster_.

{ return localU_cluster_; }

localV_cluster()

std::array<float, 2> tmtt::Stub::localV_cluster ( ) const

inline

Definition at line 154 of file Stub.h.

References localV_cluster_.

{ return localV_cluster_; }

max_qOverPt_bin()

unsigned int tmtt::Stub::max_qOverPt_bin ( ) const

inline

Definition at line 136 of file Stub.h.

References max_qOverPt_bin_.

{ return max_qOverPt_bin_; }

min_qOverPt_bin()

unsigned int tmtt::Stub::min_qOverPt_bin ( ) const

inline

Definition at line 135 of file Stub.h.

References min_qOverPt_bin_.

{ return min_qOverPt_bin_; }

nStrips()

unsigned int tmtt::Stub::nStrips ( ) const

inline

Definition at line 209 of file Stub.h.

References nStrips_.

Referenced by Stub().

{ return nStrips_; }

numMergedBend()

float tmtt::Stub::numMergedBend ( ) const

inline

Definition at line 130 of file Stub.h.

References numMergedBend_.

{ return numMergedBend_; }

operator==()

bool tmtt::Stub::operator== ( const Stub &  stubOther )

inline

Definition at line 72 of file Stub.h.

References index().

{ return (this->index() == stubOther.index()); }

phi()

float tmtt::Stub::phi ( void  ) const

inline

Definition at line 107 of file Stub.h.

References phi_.

Referenced by Particle.Particle::__str__(), tmtt::Sector::insidePhi(), tmtt::HTrphi::iPhiRange(), ntupleDataFormat.Track::phiPull(), tmtt::KFbase::printStub(), tmtt::HTrphi::store(), and tmtt::KFParamsComb::vectorM().

{ return phi_; }

phiDiff()

float tmtt::Stub::phiDiff ( float  rad, float  Pt  ) const

inline

Definition at line 138 of file Stub.h.

References funct::abs(), tmtt::Settings::invPtToDphi(), r_, and settings_.

Referenced by tmtt::Sector::insidePhi().

{ return std::abs(r_ - rad) * (settings_->invPtToDphi()) / Pt; }

psModule()

bool tmtt::Stub::psModule ( ) const

inline

Definition at line 196 of file Stub.h.

References psModule_.

Referenced by degradeResolution(), tmtt::KFParamsComb::matrixV(), tmtt::KFbase::residual(), and Stub().

{ return psModule_; }

qOverPt()

float tmtt::Stub::qOverPt ( ) const

inline

Definition at line 132 of file Stub.h.

References bend(), and qOverPtOverBend().

{ return (this->qOverPtOverBend() * this->bend()); }

qOverPtcut()

float tmtt::Stub::qOverPtcut ( ) const

inline

Definition at line 133 of file Stub.h.

References bendCut(), and qOverPtOverBend().

{ return (this->qOverPtOverBend() * this->bendCut()); }

qOverPtOverBend()

float tmtt::Stub::qOverPtOverBend ( ) const

inline

Definition at line 148 of file Stub.h.

References dphiOverBend_, tmtt::Settings::invPtToDphi(), r_, and settings_.

Referenced by calcQoverPtrange(), degradeResolution(), tmtt::StubWindowSuggest::process(), qOverPt(), qOverPtcut(), and setFrontend().

{ return dphiOverBend_ / (r_ * settings_->invPtToDphi()); }

r()

float tmtt::Stub::r ( ) const

inline

Definition at line 108 of file Stub.h.

References r_.

Referenced by tmtt::Sector::insideEtaRange(), tmtt::HTrphi::iPhiRange(), tmtt::KalmanState::KalmanState(), tmtt::KFParamsComb::matrixH(), tmtt::KFParamsComb::matrixV(), tmtt::KFbase::printStub(), tmtt::KFbase::residual(), tmtt::HTrphi::store(), tmtt::TP::trkRAtStub(), and tmtt::TP::trkZAtStub().

{ return r_; }

setDigitizeWarningsOn()

void tmtt::Stub::setDigitizeWarningsOn ( bool  newVal )

inline

Definition at line 93 of file Stub.h.

References digitizeWarningsOn_.

{ digitizeWarningsOn_ = newVal; }

setFrontend()

void tmtt::Stub::setFrontend ( const StubKiller *  stubKiller )

private

Definition at line 267 of file Stub.cc.

References funct::abs(), bend(), bendCut(), bendCutInFrontend(), bendInFrontend(), eta(), Exception, frontendPass_, edm::Ref< C, T, F >::get(), tmtt::Settings::houghMinPt(), tmtt::Settings::killLowPtStubs(), tmtt::Settings::killScenario(), tmtt::StubKiller::killStub(), tmtt::Settings::maxStubEta(), tmtt::StubKiller::none, qOverPtOverBend(), rejectedStubBend_, settings_, stubFailedDegradeWindow_, and ttStubRef_.

Referenced by Stub().

                                                     {
    frontendPass_ = true;              // Did stub pass cuts applied in front-end chip
    stubFailedDegradeWindow_ = false;  // Did it only fail cuts corresponding to windows encoded in DegradeBend.h?
    // Don't use stubs at large eta, since it is impossible to form L1 tracks from them, so they only contribute to combinatorics.
    if (std::abs(this->eta()) > settings_->maxStubEta())
      frontendPass_ = false;
    // Don't use stubs whose Pt is significantly below the Pt cut used in the L1 tracking, allowing for uncertainty in q/Pt due to stub bend resolution.
    const float qOverPtCut = 1. / settings_->houghMinPt();
    if (settings_->killLowPtStubs()) {
      // Apply this cut in the front-end electronics.
      if (std::abs(this->bendInFrontend()) - this->bendCutInFrontend() > qOverPtCut / this->qOverPtOverBend())
        frontendPass_ = false;
    }

    if (frontendPass_ && this->bend() == rejectedStubBend_) {
      throw cms::Exception(
          "BadConfig: FE stub bend window sizes provided in cfg ES source are tighter than those to make the stubs. "
          "Please fix them");
    }

    if (settings_->killLowPtStubs()) {
      // Reapply the same cut using the degraded bend information available in the off-detector electronics.
      // The reason is  that the bend degredation can move the Pt below the Pt cut, making the stub useless to the off-detector electronics.
      if (std::abs(this->bend()) - this->bendCut() > qOverPtCut / this->qOverPtOverBend())
        frontendPass_ = false;
    }

    // Emulate stubs in dead tracker regions..
    StubKiller::KillOptions killScenario = static_cast<StubKiller::KillOptions>(settings_->killScenario());
    if (killScenario != StubKiller::KillOptions::none) {
      bool kill = stubKiller->killStub(ttStubRef_.get());
      if (kill)
        frontendPass_ = false;
    }
  }

setTrackerModule()

void tmtt::Stub::setTrackerModule ( const TrackerGeometry *  trackerGeometry, const TrackerTopology *  trackerTopology, const DetId &  detId  )

private

sigmaPar()

float tmtt::Stub::sigmaPar ( ) const

inline

Definition at line 191 of file Stub.h.

References invRoot12, and stripLength_.

Referenced by tmtt::KFParamsComb::matrixV(), tmtt::KFbase::printStub(), sigmaR(), and sigmaZ().

{ return invRoot12 * stripLength_; }

sigmaPerp()

float tmtt::Stub::sigmaPerp ( ) const

inline

Definition at line 189 of file Stub.h.

References invRoot12, and stripPitch_.

Referenced by tmtt::KFParamsComb::matrixV(), and tmtt::KFbase::printStub().

{ return invRoot12 * stripPitch_; }

sigmaR()

float tmtt::Stub::sigmaR ( ) const

inline

Definition at line 186 of file Stub.h.

References barrel(), and sigmaPar().

{ return (barrel() ? 0. : sigmaPar()); }

sigmaZ()

float tmtt::Stub::sigmaZ ( ) const

inline

Definition at line 187 of file Stub.h.

References barrel(), and sigmaPar().

{ return (barrel() ? sigmaPar() : 0.); }

stripLength()

float tmtt::Stub::stripLength ( ) const

inline

Definition at line 207 of file Stub.h.

References stripLength_.

{ return stripLength_; }

stripPitch()

float tmtt::Stub::stripPitch ( ) const

inline

Definition at line 205 of file Stub.h.

References stripPitch_.

{ return stripPitch_; }

stubFailedDegradeWindow()

bool tmtt::Stub::stubFailedDegradeWindow ( ) const

inline

Definition at line 159 of file Stub.h.

References stubFailedDegradeWindow_.

{ return stubFailedDegradeWindow_; }

theta()

float tmtt::Stub::theta ( void  ) const

inline

Definition at line 110 of file Stub.h.

References r_, and z_.

Referenced by calcDphiOverBend(), and Tau.Tau::zImpact().

{ return atan2(r_, z_); }

tiltAngle()

float tmtt::Stub::tiltAngle ( ) const

inline

Definition at line 184 of file Stub.h.

References tiltAngle_.

Referenced by tmtt::KFParamsComb::matrixV().

{ return tiltAngle_; }

tiltedBarrel()

bool tmtt::Stub::tiltedBarrel ( ) const

inline

Definition at line 203 of file Stub.h.

References tiltedBarrel_.

Referenced by approxB(), and tmtt::KFParamsComb::matrixV().

{ return tiltedBarrel_; }

trackerModule()

const TrackerModule* tmtt::Stub::trackerModule ( ) const

inline

Definition at line 78 of file Stub.h.

References trackerModule_.

Referenced by tmtt::StubWindowSuggest::updateStoredWindow().

{ return trackerModule_; }

trkPhiAtR()

float tmtt::Stub::trkPhiAtR ( float  rad ) const

inline

Definition at line 140 of file Stub.h.

References bend_, dphiOverBend_, phi_, and r_.

Referenced by tmtt::Sector::insidePhi().

{ return phi_ + (bend_ * dphiOverBend_) * (1. - rad / r_); }

trkPhiAtRcut()

float tmtt::Stub::trkPhiAtRcut ( float  rad ) const

inline

Definition at line 142 of file Stub.h.

References funct::abs(), bendCut(), dphiOverBend_, and r_.

Referenced by tmtt::Sector::insidePhi().

{ return (bendCut() * dphiOverBend_) * std::abs(1. - rad / r_); }

ttStubRef()

const TTStubRef& tmtt::Stub::ttStubRef ( ) const

inline

Definition at line 75 of file Stub.h.

References ttStubRef_.

{ return ttStubRef_; }

z()

float tmtt::Stub::z ( ) const

inline

Definition at line 109 of file Stub.h.

References z_.

Referenced by tmtt::Sector::insideEtaRange(), tmtt::KalmanState::KalmanState(), tmtt::KFParamsComb::matrixV(), geometryXMLparser.Alignable::pos(), tmtt::KFbase::printStub(), ntupleDataFormat._HitObject::r3D(), tmtt::KFbase::residual(), tmtt::TP::trkRAtStub(), tmtt::TP::trkZAtStub(), and tmtt::KFParamsComb::vectorM().

{ return z_; }

Member Data Documentation

alpha_

float tmtt::Stub::alpha_

private

Definition at line 252 of file Stub.h.

Referenced by alpha(), and Stub().

assocTP_

const TP* tmtt::Stub::assocTP_

private

Definition at line 264 of file Stub.h.

Referenced by assocTP(), and fillTruth().

assocTPofCluster_

std::array<const TP*, 2> tmtt::Stub::assocTPofCluster_

private

Definition at line 267 of file Stub.h.

Referenced by assocTPofCluster(), fillTruth(), genuineCluster(), and Stub().

assocTPs_

std::set<const TP*> tmtt::Stub::assocTPs_

private

Definition at line 265 of file Stub.h.

Referenced by assocTPs(), fillTruth(), and genuine().

barrel_

bool tmtt::Stub::barrel_

private

Definition at line 288 of file Stub.h.

Referenced by barrel().

bend_

float tmtt::Stub::bend_

private

Definition at line 242 of file Stub.h.

Referenced by bend(), digitize(), Stub(), and trkPhiAtR().

bendInFrontend_

float tmtt::Stub::bendInFrontend_

private

Definition at line 259 of file Stub.h.

Referenced by bendInFrontend(), and Stub().

degradeBend_

const DegradeBend* tmtt::Stub::degradeBend_

private

Definition at line 277 of file Stub.h.

Referenced by degradeResolution().

digitalStub_

std::unique_ptr<DigitalStub> tmtt::Stub::digitalStub_

private

Definition at line 269 of file Stub.h.

Referenced by digitalStub(), and digitize().

digitizeWarningsOn_

bool tmtt::Stub::digitizeWarningsOn_

private

Definition at line 271 of file Stub.h.

Referenced by digitize(), and setDigitizeWarningsOn().

dphiOverBend_

float tmtt::Stub::dphiOverBend_

private

Definition at line 243 of file Stub.h.

Referenced by calcDphiOverBend(), dphiOverBend(), qOverPtOverBend(), trkPhiAtR(), and trkPhiAtRcut().

frontendPass_

bool tmtt::Stub::frontendPass_

private

Definition at line 255 of file Stub.h.

Referenced by frontendPass(), and setFrontend().

index_in_vStubs_

unsigned int tmtt::Stub::index_in_vStubs_

private

Definition at line 235 of file Stub.h.

Referenced by index().

invRoot12

const float tmtt::Stub::invRoot12 = sqrt(1. / 12.)

private

Definition at line 293 of file Stub.h.

Referenced by sigmaPar(), and sigmaPerp().

iphi_

unsigned int tmtt::Stub::iphi_

private

Definition at line 251 of file Stub.h.

Referenced by iphi(), and Stub().

lastDigiStep_

DigiStage tmtt::Stub::lastDigiStep_

private

Definition at line 270 of file Stub.h.

Referenced by digitize().

layerId_

unsigned int tmtt::Stub::layerId_

private

Definition at line 281 of file Stub.h.

Referenced by layerId().

layerIdReduced_

unsigned int tmtt::Stub::layerIdReduced_

private

Definition at line 282 of file Stub.h.

Referenced by layerIdReduced().

localU_cluster_

std::array<float, 2> tmtt::Stub::localU_cluster_

private

Definition at line 248 of file Stub.h.

Referenced by localU_cluster(), and Stub().

localV_cluster_

std::array<float, 2> tmtt::Stub::localV_cluster_

private

Definition at line 249 of file Stub.h.

Referenced by localV_cluster(), and Stub().

max_qOverPt_bin_

unsigned int tmtt::Stub::max_qOverPt_bin_

private

Definition at line 245 of file Stub.h.

Referenced by calcQoverPtrange(), digitize(), and max_qOverPt_bin().

min_qOverPt_bin_

unsigned int tmtt::Stub::min_qOverPt_bin_

private

Definition at line 244 of file Stub.h.

Referenced by calcQoverPtrange(), digitize(), and min_qOverPt_bin().

nStrips_

unsigned int tmtt::Stub::nStrips_

private

Definition at line 286 of file Stub.h.

Referenced by nStrips().

numMergedBend_

unsigned int tmtt::Stub::numMergedBend_

private

Definition at line 261 of file Stub.h.

Referenced by bendCut(), numMergedBend(), and Stub().

phi_

float tmtt::Stub::phi_

private

Definition at line 239 of file Stub.h.

Referenced by digitize(), phi(), Stub(), and trkPhiAtR().

psModule_

bool tmtt::Stub::psModule_

private

Definition at line 287 of file Stub.h.

Referenced by psModule().

r_

float tmtt::Stub::r_

private

Definition at line 240 of file Stub.h.

Referenced by approxB(), digitize(), eta(), phiDiff(), qOverPtOverBend(), r(), Stub(), theta(), trkPhiAtR(), and trkPhiAtRcut().

rejectedStubBend_

const float tmtt::Stub::rejectedStubBend_ = 99999.

private

Definition at line 291 of file Stub.h.

Referenced by degradeResolution(), and setFrontend().

settings_

const Settings* tmtt::Stub::settings_

private

Definition at line 233 of file Stub.h.

Referenced by approxB(), bendCut(), bendCutInFrontend(), calcDphiOverBend(), calcQoverPtrange(), degradeResolution(), digitize(), fillTruth(), phiDiff(), qOverPtOverBend(), and setFrontend().

stripLength_

float tmtt::Stub::stripLength_

private

Definition at line 285 of file Stub.h.

Referenced by sigmaPar(), and stripLength().

stripPitch_

float tmtt::Stub::stripPitch_

private

Definition at line 284 of file Stub.h.

Referenced by sigmaPerp(), and stripPitch().

stubFailedDegradeWindow_

bool tmtt::Stub::stubFailedDegradeWindow_

private

Definition at line 257 of file Stub.h.

Referenced by setFrontend(), and stubFailedDegradeWindow().

tiltAngle_

float tmtt::Stub::tiltAngle_

private

Definition at line 283 of file Stub.h.

Referenced by tiltAngle().

tiltedBarrel_

bool tmtt::Stub::tiltedBarrel_

private

Definition at line 289 of file Stub.h.

Referenced by tiltedBarrel().

trackerModule_

const TrackerModule* tmtt::Stub::trackerModule_

private

Definition at line 274 of file Stub.h.

Referenced by calcDphiOverBend(), degradeResolution(), Stub(), and trackerModule().

ttStubRef_

TTStubRef tmtt::Stub::ttStubRef_

private

Definition at line 231 of file Stub.h.

Referenced by fillTruth(), setFrontend(), Stub(), and ttStubRef().

z_

float tmtt::Stub::z_

private

Definition at line 241 of file Stub.h.

Referenced by approxB(), digitize(), eta(), Stub(), theta(), and z().