#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.

89 { 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.

       : 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;
     }
   }

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_.

Member Function Documentation

◆ alpha()

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

inline

Definition at line 119 of file Stub.h.

119 { return alpha_; }

References alpha_.

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

◆ approxB()

double tmtt::Stub::approxB ( )

private

Definition at line 304 of file Stub.cc.

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

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

Referenced by calcDphiOverBend().

◆ assocTP()

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

inline

Definition at line 168 of file Stub.h.

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

References assocTP_.

◆ assocTPofCluster()

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

inline

Definition at line 176 of file Stub.h.

                                                     {
       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()).

References assocTPofCluster_.

◆ assocTPs()

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

inline

Definition at line 164 of file Stub.h.

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

References assocTPs_.

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

◆ barrel()

bool tmtt::Stub::barrel ( ) const

inline

Definition at line 201 of file Stub.h.

201 { return barrel_; }

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().

◆ bend()

float tmtt::Stub::bend ( ) const

inline

Definition at line 126 of file Stub.h.

126 { return bend_; }

References bend_.

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

◆ bendCut()

float tmtt::Stub::bendCut ( ) const

inline

Definition at line 128 of file Stub.h.

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

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

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

◆ bendCutInFrontend()

float tmtt::Stub::bendCutInFrontend ( ) const

inline

Definition at line 124 of file Stub.h.

124 { return settings_->bendCut(); }

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

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

◆ bendInFrontend()

float tmtt::Stub::bendInFrontend ( ) const

inline

Definition at line 123 of file Stub.h.

123 { return bendInFrontend_; }

References bendInFrontend_.

Referenced by setFrontend().

◆ calcDphiOverBend()

void tmtt::Stub::calcDphiOverBend ( )

private

Definition at line 314 of file Stub.cc.

                               {
     // 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_;
     }
   }

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().

◆ calcQoverPtrange()

void tmtt::Stub::calcQoverPtrange ( )

Definition at line 163 of file Stub.cc.

                               {
     // 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;
   }

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

Referenced by digitize(), and Stub().

◆ degradeResolution()

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

private

Definition at line 248 of file Stub.cc.

                                                                                        {
     // 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);
   }

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().

◆ digitalStub()

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

inline

Definition at line 96 of file Stub.h.

96 { return digitalStub_.get(); }

References digitalStub_.

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

◆ digitize()

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

Definition at line 203 of file Stub.cc.

                                                                   {
     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;
       }
     }
   }

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

◆ dphiOverBend()

float tmtt::Stub::dphiOverBend ( ) const

inline

Definition at line 146 of file Stub.h.

146 { return dphiOverBend_; }

References dphiOverBend_.

◆ eta()

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

inline

Definition at line 111 of file Stub.h.

111 { return asinh(z_ / r_); }

References r_, and z_.

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

◆ 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.

                                                                                  {
     //--- 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.
         }
       }
     }
   }

References assocTP_, assocTPofCluster_, assocTPs_, genuine(), genuineCluster(), settings_, tmtt::Settings::stubMatchStrict(), and ttStubRef_.

◆ frontendPass()

bool tmtt::Stub::frontendPass ( ) const

inline

Definition at line 157 of file Stub.h.

157 { return frontendPass_; }

References frontendPass_.

◆ genuine()

bool tmtt::Stub::genuine ( ) const

inline

Definition at line 167 of file Stub.h.

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

References assocTPs_.

Referenced by fillTruth().

◆ genuineCluster()

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

inline

Definition at line 173 of file Stub.h.

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

References assocTPofCluster_.

Referenced by fillTruth().

◆ index()

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

inline

Definition at line 101 of file Stub.h.

101 { return index_in_vStubs_; }

References index_in_vStubs_.

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

◆ iphi()

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

inline

Definition at line 115 of file Stub.h.

115 { return iphi_; }

References iphi_.

◆ layerId()

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

inline

Definition at line 198 of file Stub.h.

198 { return layerId_; }

References layerId_.

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

◆ layerIdReduced()

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

inline

Definition at line 200 of file Stub.h.

200 { return layerIdReduced_; }

References layerIdReduced_.

◆ localU_cluster()

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

inline

Definition at line 153 of file Stub.h.

153 { return localU_cluster_; }

References localU_cluster_.

◆ localV_cluster()

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

inline

Definition at line 154 of file Stub.h.

154 { return localV_cluster_; }

References localV_cluster_.

◆ max_qOverPt_bin()

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

inline

Definition at line 136 of file Stub.h.

136 { return max_qOverPt_bin_; }

References max_qOverPt_bin_.

◆ min_qOverPt_bin()

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

inline

Definition at line 135 of file Stub.h.

135 { return min_qOverPt_bin_; }

References min_qOverPt_bin_.

◆ nStrips()

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

inline

Definition at line 209 of file Stub.h.

209 { return nStrips_; }

References nStrips_.

Referenced by Stub().

◆ numMergedBend()

float tmtt::Stub::numMergedBend ( ) const

inline

Definition at line 130 of file Stub.h.

130 { return numMergedBend_; }

References numMergedBend_.

◆ operator==()

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

inline

Definition at line 72 of file Stub.h.

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

References index().

◆ phi()

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

inline

Definition at line 107 of file Stub.h.

107 { return phi_; }

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().

◆ phiDiff()

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

inline

Definition at line 138 of file Stub.h.

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

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

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

◆ psModule()

bool tmtt::Stub::psModule ( ) const

inline

Definition at line 196 of file Stub.h.

196 { return psModule_; }

References psModule_.

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

◆ qOverPt()

float tmtt::Stub::qOverPt ( ) const

inline

Definition at line 132 of file Stub.h.

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

References bend(), and qOverPtOverBend().

◆ qOverPtcut()

float tmtt::Stub::qOverPtcut ( ) const

inline

Definition at line 133 of file Stub.h.

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

References bendCut(), and qOverPtOverBend().

◆ qOverPtOverBend()

float tmtt::Stub::qOverPtOverBend ( ) const

inline

Definition at line 148 of file Stub.h.

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

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

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

◆ r()

float tmtt::Stub::r ( ) const

inline

Definition at line 108 of file Stub.h.

108 { return r_; }

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().

◆ setDigitizeWarningsOn()

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

inline

Definition at line 93 of file Stub.h.

93 { digitizeWarningsOn_ = newVal; }

References digitizeWarningsOn_.

◆ setFrontend()

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

private

Definition at line 267 of file Stub.cc.

                                                      {
     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;
     }
   }

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().

◆ 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.

191 { return invRoot12 * stripLength_; }

References invRoot12, and stripLength_.

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

◆ sigmaPerp()

float tmtt::Stub::sigmaPerp ( ) const

inline

Definition at line 189 of file Stub.h.

189 { return invRoot12 * stripPitch_; }

References invRoot12, and stripPitch_.

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

◆ sigmaR()

float tmtt::Stub::sigmaR ( ) const

inline

Definition at line 186 of file Stub.h.

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

References barrel(), and sigmaPar().

◆ sigmaZ()

float tmtt::Stub::sigmaZ ( ) const

inline

Definition at line 187 of file Stub.h.

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

References barrel(), and sigmaPar().

◆ stripLength()

float tmtt::Stub::stripLength ( ) const

inline

Definition at line 207 of file Stub.h.

207 { return stripLength_; }

References stripLength_.

◆ stripPitch()

float tmtt::Stub::stripPitch ( ) const

inline

Definition at line 205 of file Stub.h.

205 { return stripPitch_; }

References stripPitch_.

◆ stubFailedDegradeWindow()

bool tmtt::Stub::stubFailedDegradeWindow ( ) const

inline

Definition at line 159 of file Stub.h.

159 { return stubFailedDegradeWindow_; }

References stubFailedDegradeWindow_.

◆ theta()

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

inline

Definition at line 110 of file Stub.h.

110 { return atan2(r_, z_); }

References r_, and z_.

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

◆ tiltAngle()

float tmtt::Stub::tiltAngle ( ) const

inline

Definition at line 184 of file Stub.h.

184 { return tiltAngle_; }

References tiltAngle_.

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

◆ tiltedBarrel()

bool tmtt::Stub::tiltedBarrel ( ) const

inline

Definition at line 203 of file Stub.h.

203 { return tiltedBarrel_; }

References tiltedBarrel_.

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

◆ trackerModule()

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

inline

Definition at line 78 of file Stub.h.

78 { return trackerModule_; }

References trackerModule_.

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

◆ trkPhiAtR()

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

inline

Definition at line 140 of file Stub.h.

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

References bend_, dphiOverBend_, phi_, and r_.

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

◆ trkPhiAtRcut()

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

inline

Definition at line 142 of file Stub.h.

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

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

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

◆ ttStubRef()

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

inline

Definition at line 75 of file Stub.h.

75 { return ttStubRef_; }

References ttStubRef_.

◆ z()

float tmtt::Stub::z ( ) const

inline

Definition at line 109 of file Stub.h.

109 { return z_; }

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().

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().

Public Types

Public Member Functions

Private Member Functions

Private Attributes

Detailed Description

Member Enumeration Documentation

◆ DigiStage

Constructor & Destructor Documentation

◆ Stub() [1/2]

◆ Stub() [2/2]

Member Function Documentation

◆ alpha()

◆ approxB()

◆ assocTP()

◆ assocTPofCluster()

◆ assocTPs()

◆ barrel()

◆ bend()

◆ bendCut()

◆ bendCutInFrontend()

◆ bendInFrontend()

◆ calcDphiOverBend()

◆ calcQoverPtrange()

◆ degradeResolution()

◆ digitalStub()

◆ digitize()

◆ dphiOverBend()

◆ eta()

◆ fillTruth()

◆ frontendPass()

◆ genuine()

◆ genuineCluster()

◆ index()

◆ iphi()

◆ layerId()

◆ layerIdReduced()

◆ localU_cluster()

◆ localV_cluster()

◆ max_qOverPt_bin()

◆ min_qOverPt_bin()

◆ nStrips()

◆ numMergedBend()

◆ operator==()

◆ phi()

◆ phiDiff()

◆ psModule()

◆ qOverPt()

◆ qOverPtcut()

◆ qOverPtOverBend()

◆ r()

◆ setDigitizeWarningsOn()

◆ setFrontend()

◆ setTrackerModule()

◆ sigmaPar()

◆ sigmaPerp()

◆ sigmaR()

◆ sigmaZ()

◆ stripLength()

◆ stripPitch()

◆ stubFailedDegradeWindow()

◆ theta()

◆ tiltAngle()

◆ tiltedBarrel()

◆ trackerModule()

◆ trkPhiAtR()

◆ trkPhiAtRcut()

◆ ttStubRef()

◆ z()

Member Data Documentation

◆ alpha_

◆ assocTP_

◆ assocTPofCluster_

◆ assocTPs_

◆ barrel_

◆ bend_

◆ bendInFrontend_

◆ degradeBend_

◆ digitalStub_

◆ digitizeWarningsOn_

◆ dphiOverBend_