#include <Stub.h>

Inheritance diagram for l1tVertexFinder::Stub:

Public Member Functions
const TrackingParticlePtr	assocTP () const

const std::set< TrackingParticlePtr > &	assocTPs () const

void	fillTruth (edm::Handle< TTStubAssMap > mcTruthTTStubHandle, edm::Handle< TTClusterAssMap > mcTruthTTClusterHandle)

unsigned int	layerId () const

float	phi () const

bool	psModule () const

float	r () const

	Stub ()

	Stub (const TTStubRef &ttStubRef, const AnalysisSettings &settings, const TrackerGeometry trackerGeometry, const TrackerTopology trackerTopology)

float	z () const

	~Stub ()

Public Member Functions inherited from edm::Ref< C, T, F >
T const *	get () const
	Returns C++ pointer to the item. More...

bool	hasProductCache () const
	Returns true if container referenced by the Ref has been cached. More...

ProductID	id () const
	Accessor for product ID. More...

key_type	index () const

bool	isAvailable () const

bool	isNonnull () const
	Checks for non-null. More...

bool	isNull () const
	Checks for null. More...

bool	isTransient () const
	Checks if this ref is transient (i.e. not persistable). More...

key_type	key () const
	Accessor for product key. More...

bool	operator! () const
	Checks for null. More...

T const &	operator* () const
	Dereference operator. More...

T const *	operator-> () const
	Member dereference operator. More...

EDProductGetter const *	productGetter () const
	Accessor for product getter. More...

	Ref ()
	Default constructor needed for reading from persistent store. Not for direct use. More...

	Ref (Handle< C > const &handle, key_type itemKey, bool setNow=true)
	General purpose constructor from handle. More...

	Ref (OrphanHandle< C > const &handle, key_type itemKey, bool setNow=true)
	General purpose constructor from orphan handle. More...

	Ref (C const *product, key_type itemKey, bool setNow=true)
	Constructors for ref to object that is not in an event. More...

	Ref (TestHandle< C > const &handle, key_type itemKey, bool setNow=true)
	Constructor from test handle. More...

	Ref (ProductID const &productID, key_type itemKey, EDProductGetter const *prodGetter)

	Ref (ProductID const &iProductID, T const item, key_type itemKey, C const )
	Constructor for use in the various X::fillView(...) functions. More...

	Ref (ProductID const &iProductID, T const *item, key_type itemKey)

	Ref (ProductID const &iProductID, T const *item, key_type itemKey, bool transient)

	Ref (ProductID const &iId)

	Ref (RefProd< C > const &refProd, key_type itemKey)
	Constructor from RefProd<C> and key. More...

	Ref (RefCore const &iRefCore, key_type const &iKey)

RefCore const &	refCore () const

	~Ref ()
	Destructor. More...

Private Member Functions
void	setModuleInfo (const TrackerGeometry trackerGeometry, const TrackerTopology trackerTopology, const DetId &detId)

Private Attributes
TrackingParticlePtr	assocTP_

std::set< TrackingParticlePtr >	assocTPs_

bool	barrel_

unsigned int	endcapRing_

unsigned int	idDet_

unsigned int	layerId_

float	moduleMaxPhi_

float	moduleMaxR_

float	moduleMaxZ_

float	moduleMinPhi_

float	moduleMinR_

float	moduleMinZ_

unsigned int	nStrips_

bool	outerModuleAtSmallerR_

float	phi_

bool	psModule_

float	r_

float	sensorWidth_

const AnalysisSettings *	settings_

float	sigmaPar_

float	sigmaPerp_

float	stripLength_

float	stripPitch_

float	z_

Additional Inherited Members
Public Types inherited from edm::Ref< C, T, F >
typedef boost::binary_traits< F >::second_argument_type	argument_type

typedef T const	element_type

typedef F	finder_type

typedef std::remove_cv< typename std::remove_reference< argument_type >::type >::type	key_type

typedef C	product_type
	for export More...

typedef T	value_type

Static Public Member Functions inherited from edm::Ref< C, T, F >
static short	Class_Version ()

static key_type	invalidKey ()

Detailed Description

Definition at line 37 of file Stub.h.

Constructor & Destructor Documentation

◆ Stub() [1/2]

l1tVertexFinder::Stub::Stub ( )

Definition at line 7 of file Stub.cc.

       : settings_(nullptr),
         phi_(0.),
         r_(0.),
         z_(0.),
         idDet_(0),
         moduleMinR_(0.),
         moduleMaxR_(0.),
         moduleMinPhi_(0.),
         moduleMaxPhi_(0.),
         moduleMinZ_(0.),
         moduleMaxZ_(0.),
         psModule_(false),
         layerId_(0),
         endcapRing_(0),
         barrel_(false),
         sigmaPerp_(0.),
         sigmaPar_(0.),
         stripPitch_(0.),
         stripLength_(0.),
         nStrips_(0),
         sensorWidth_(0.),
         outerModuleAtSmallerR_(false) {}

◆ Stub() [2/2]

l1tVertexFinder::Stub::Stub	(	const TTStubRef &	ttStubRef,
		const AnalysisSettings &	settings,
		const TrackerGeometry *	trackerGeometry,
		const TrackerTopology *	trackerTopology
	)

Definition at line 31 of file Stub.cc.

References funct::abs(), Exception, TrackerGeometry::idToDet(), phi_, r_, setModuleInfo(), settings_, l1tVertexFinder::AlgoSettings::trackerHalfLength(), l1tVertexFinder::AlgoSettings::trackerOuterRadius(), and z_.

       : TTStubRef(ttStubRef), settings_(&settings) {
     auto geoDetId = trackerGeometry->idToDet(ttStubRef->clusterRef(0)->getDetId())->geographicalId();
     auto theGeomDet = trackerGeometry->idToDet(geoDetId);
     auto measurementPoint = ttStubRef->clusterRef(0)->findAverageLocalCoordinatesCentered();
     auto pos = theGeomDet->surface().toGlobal(theGeomDet->topology().localPosition(measurementPoint));
 
     phi_ = pos.phi();
     r_ = pos.perp();
     z_ = pos.z();
 
     if (r_ < settings_->trackerInnerRadius() || r_ > settings_->trackerOuterRadius() ||
         std::abs(z_) > settings_->trackerHalfLength()) {
       throw cms::Exception(
           "Stub: Stub found outside assumed tracker volume. Please update tracker dimensions specified in Settings.h!")
           << " r=" << r_ << " z=" << z_ << " " << ttStubRef->getDetId().subdetId() << std::endl;
     }
 
     // Set info about the module this stub is in
     this->setModuleInfo(trackerGeometry, trackerTopology, geoDetId);
   }

◆ ~Stub()

l1tVertexFinder::Stub::~Stub ( )

inline

Definition at line 45 of file Stub.h.

45 {}

Member Function Documentation

◆ assocTP()

const TrackingParticlePtr l1tVertexFinder::Stub::assocTP ( ) const

inline

Definition at line 74 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.

◆ assocTPs()

const std::set<TrackingParticlePtr>& l1tVertexFinder::Stub::assocTPs ( ) const

inline

Definition at line 71 of file Stub.h.

References assocTPs_.

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

◆ fillTruth()

void l1tVertexFinder::Stub::fillTruth	(	edm::Handle< TTStubAssMap >	mcTruthTTStubHandle,
		edm::Handle< TTClusterAssMap >	mcTruthTTClusterHandle
	)

Definition at line 57 of file Stub.cc.

References assocTP_, assocTPs_, TTStubAssociationMap< T >::findTrackingParticlePtr(), TTClusterAssociationMap< T >::findTrackingParticlePtrs(), TTStubAssociationMap< T >::isGenuine(), edm::Ptr< T >::isNull(), settings_, and l1tVertexFinder::AnalysisSettings::stubMatchStrict().

Referenced by TPStubValueMapProducer::produce().

                                                                           {
     const TTStubRef& ttStubRef(*this);  // Cast to base class
 
     //--- Fill assocTP_ info. If both clusters in this stub were produced by the same single tracking particle, find out which one it was.
 
     // Require same TP contributed to both clusters.
     if (mcTruthTTStubHandle->isGenuine(ttStubRef)) {
       assocTP_ = mcTruthTTStubHandle->findTrackingParticlePtr(ttStubRef);
     }
 
     // Fill assocTPs_ info.
     if (settings_->stubMatchStrict()) {
       // We consider only stubs in which this TP contributed to both clusters.
       if (!assocTP_.isNull())
         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.
         std::vector<edm::Ptr<TrackingParticle>> vecTpPtr =
             mcTruthTTClusterHandle->findTrackingParticlePtrs(ttClusterRef);
 
         for (const edm::Ptr<TrackingParticle>& tpPtr : vecTpPtr) {
           assocTPs_.insert(tpPtr);
         }
       }
     }
   }

◆ layerId()

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

inline

Definition at line 66 of file Stub.h.

References layerId_.

66 { return layerId_; }

l1tVertexFinder::Stub::layerId_

unsigned int layerId_

Definition: Stub.h:100

◆ phi()

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

inline

Definition at line 57 of file Stub.h.

References phi_.

Referenced by Particle.Particle::__str__(), datamodel.Object::DeltaR(), datamodel.Object::p4(), and ntupleDataFormat.Track::phiPull().

57 { return phi_; }

l1tVertexFinder::Stub::phi_

float phi_

Definition: Stub.h:87

◆ psModule()

bool l1tVertexFinder::Stub::psModule ( ) const

inline

Definition at line 64 of file Stub.h.

References psModule_.

64 { return psModule_; }

l1tVertexFinder::Stub::psModule_

bool psModule_

Definition: Stub.h:99

◆ r()

float l1tVertexFinder::Stub::r ( ) const

inline

Definition at line 58 of file Stub.h.

References r_.

58 { return r_; }

l1tVertexFinder::Stub::r_

float r_

Definition: Stub.h:88

◆ setModuleInfo()

void l1tVertexFinder::Stub::setModuleInfo	(	const TrackerGeometry *	trackerGeometry,
		const TrackerTopology *	trackerTopology,
		const DetId &	detId
	)

private

Definition at line 90 of file Stub.cc.

References barrel_, Surface::bounds(), endcapRing_, TrackerGeometry::getDetectorType(), idDet_, TrackerGeometry::idToDetUnit(), TrackerTopology::layer(), layerId_, SiStripPI::max, SiStripPI::min, moduleMaxPhi_, moduleMaxR_, moduleMaxZ_, moduleMinPhi_, moduleMinR_, moduleMinZ_, PixelTopology::nrows(), nStrips_, outerModuleAtSmallerR_, TrackerTopology::partnerDetId(), PV3DBase< T, PVType, FrameType >::perp(), TrackerGeometry::Ph2PSP, PV3DBase< T, PVType, FrameType >::phi(), PixelTopology::pitch(), GeomDet::position(), psModule_, HLT_2023v12_cff::R0, sensorWidth_, TrackerTopology::side(), sigmaPar_, sigmaPerp_, mathSSE::sqrt(), stripLength_, stripPitch_, DetId::subdetId(), GeomDet::surface(), StripSubdetector::TIB, TrackerTopology::tidRing(), TrackerTopology::tidWheel(), StripSubdetector::TOB, unit(), PV3DBase< T, PVType, FrameType >::z(), and Z0.

Referenced by Stub().

                                                {
     idDet_ = detId();
 
     // Get min & max (r,phi,z) coordinates of the centre of the two sensors containing this stub.
     const GeomDetUnit* det0 = trackerGeometry->idToDetUnit(detId);
     const GeomDetUnit* det1 = trackerGeometry->idToDetUnit(trackerTopology->partnerDetId(detId));
 
     float R0 = det0->position().perp();
     float R1 = det1->position().perp();
     float PHI0 = det0->position().phi();
     float PHI1 = det1->position().phi();
     float Z0 = det0->position().z();
     float Z1 = det1->position().z();
     moduleMinR_ = std::min(R0, R1);
     moduleMaxR_ = std::max(R0, R1);
     moduleMinPhi_ = std::min(PHI0, PHI1);
     moduleMaxPhi_ = std::max(PHI0, PHI1);
     moduleMinZ_ = std::min(Z0, Z1);
     moduleMaxZ_ = std::max(Z0, Z1);
 
     // Note if module is PS or 2S, and whether in barrel or endcap.
     psModule_ =
         trackerGeometry->getDetectorType(detId) ==
         TrackerGeometry::ModuleType::
             Ph2PSP;  // From https://github.com/cms-sw/cmssw/blob/CMSSW_8_1_X/Geometry/TrackerGeometryBuilder/README.md
     barrel_ = detId.subdetId() == StripSubdetector::TOB || detId.subdetId() == StripSubdetector::TIB;
 
     // Encode layer ID.
     if (barrel_) {
       layerId_ = trackerTopology->layer(detId);  // barrel layer 1-6 encoded as 1-6
     } else {
       // layerId_ = 10*detId.iSide() + detId.iDisk(); // endcap layer 1-5 encoded as 11-15 (endcap A) or 21-25 (endcapB)
       // EJC This seems to give the same encoding
       layerId_ = 10 * trackerTopology->side(detId) + trackerTopology->tidWheel(detId);
     }
 
     // Note module ring in endcap
     // endcapRing_ = barrel_  ?  0  :  detId.iRing();
     endcapRing_ = barrel_ ? 0 : trackerTopology->tidRing(detId);
 
     // Get sensor strip or pixel pitch using innermost sensor of pair.
 
     const PixelGeomDetUnit* unit = reinterpret_cast<const PixelGeomDetUnit*>(det0);
     const PixelTopology& topo = unit->specificTopology();
     const Bounds& bounds = det0->surface().bounds();
 
     std::pair<float, float> pitch = topo.pitch();
     stripPitch_ = pitch.first;      // Strip pitch (or pixel pitch along shortest axis)
     stripLength_ = pitch.second;    //  Strip length (or pixel pitch along longest axis)
     nStrips_ = topo.nrows();        // No. of strips in sensor
     sensorWidth_ = bounds.width();  // Width of sensitive region of sensor (= stripPitch * nStrips).
 
     outerModuleAtSmallerR_ = false;
     if (barrel_ && det0->position().perp() > det1->position().perp()) {
       outerModuleAtSmallerR_ = true;
     }
 
     sigmaPerp_ = stripPitch_ / sqrt(12.);  // resolution perpendicular to strip (or to longest pixel axis)
     sigmaPar_ = stripLength_ / sqrt(12.);  // resolution parallel to strip (or to longest pixel axis)
   }

◆ z()

float l1tVertexFinder::Stub::z ( ) const

inline

Definition at line 59 of file Stub.h.

References z_.

Referenced by geometryXMLparser.Alignable::pos(), and ntupleDataFormat._HitObject::r3D().

59 { return z_; }

l1tVertexFinder::Stub::z_

float z_

Definition: Stub.h:89

Member Data Documentation

◆ assocTP_

TrackingParticlePtr l1tVertexFinder::Stub::assocTP_

private

Definition at line 111 of file Stub.h.

Referenced by assocTP(), and fillTruth().

◆ assocTPs_

std::set<TrackingParticlePtr> l1tVertexFinder::Stub::assocTPs_

private

Definition at line 112 of file Stub.h.

Referenced by assocTPs(), and fillTruth().

◆ barrel_

bool l1tVertexFinder::Stub::barrel_

private

Definition at line 102 of file Stub.h.

Referenced by setModuleInfo().

◆ endcapRing_

unsigned int l1tVertexFinder::Stub::endcapRing_

private

Definition at line 101 of file Stub.h.

Referenced by setModuleInfo().

◆ idDet_

unsigned int l1tVertexFinder::Stub::idDet_

private

Definition at line 92 of file Stub.h.

Referenced by setModuleInfo().

◆ layerId_

unsigned int l1tVertexFinder::Stub::layerId_

private

Definition at line 100 of file Stub.h.

Referenced by layerId(), and setModuleInfo().

◆ moduleMaxPhi_

float l1tVertexFinder::Stub::moduleMaxPhi_

private

Definition at line 96 of file Stub.h.

Referenced by setModuleInfo().

◆ moduleMaxR_

float l1tVertexFinder::Stub::moduleMaxR_

private

Definition at line 94 of file Stub.h.

Referenced by setModuleInfo().

◆ moduleMaxZ_

float l1tVertexFinder::Stub::moduleMaxZ_

private

Definition at line 98 of file Stub.h.

Referenced by setModuleInfo().

◆ moduleMinPhi_

float l1tVertexFinder::Stub::moduleMinPhi_

private

Definition at line 95 of file Stub.h.

Referenced by setModuleInfo().

◆ moduleMinR_

float l1tVertexFinder::Stub::moduleMinR_

private

Definition at line 93 of file Stub.h.

Referenced by setModuleInfo().

◆ moduleMinZ_

float l1tVertexFinder::Stub::moduleMinZ_

private

Definition at line 97 of file Stub.h.

Referenced by setModuleInfo().

◆ nStrips_

unsigned int l1tVertexFinder::Stub::nStrips_

private

Definition at line 107 of file Stub.h.

Referenced by setModuleInfo().

◆ outerModuleAtSmallerR_

bool l1tVertexFinder::Stub::outerModuleAtSmallerR_

private

Definition at line 109 of file Stub.h.

Referenced by setModuleInfo().

◆ phi_

float l1tVertexFinder::Stub::phi_

private

Definition at line 87 of file Stub.h.

Referenced by phi(), and Stub().

◆ psModule_

bool l1tVertexFinder::Stub::psModule_

private

Definition at line 99 of file Stub.h.

Referenced by psModule(), and setModuleInfo().

◆ r_

float l1tVertexFinder::Stub::r_

private

Definition at line 88 of file Stub.h.

Referenced by r(), and Stub().

◆ sensorWidth_

float l1tVertexFinder::Stub::sensorWidth_

private

Definition at line 108 of file Stub.h.

Referenced by setModuleInfo().

◆ settings_

const AnalysisSettings* l1tVertexFinder::Stub::settings_

private

Definition at line 85 of file Stub.h.

Referenced by fillTruth(), and Stub().

◆ sigmaPar_

float l1tVertexFinder::Stub::sigmaPar_

private

Definition at line 104 of file Stub.h.

Referenced by setModuleInfo().

◆ sigmaPerp_

float l1tVertexFinder::Stub::sigmaPerp_

private

Definition at line 103 of file Stub.h.

Referenced by setModuleInfo().

◆ stripLength_

float l1tVertexFinder::Stub::stripLength_

private

Definition at line 106 of file Stub.h.

Referenced by setModuleInfo().

◆ stripPitch_

float l1tVertexFinder::Stub::stripPitch_

private

Definition at line 105 of file Stub.h.

Referenced by setModuleInfo().

◆ z_

float l1tVertexFinder::Stub::z_

private

Definition at line 89 of file Stub.h.

Referenced by Stub(), and z().

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ Stub() [1/2]

◆ Stub() [2/2]

◆ ~Stub()

Member Function Documentation

◆ assocTP()

◆ assocTPs()

◆ fillTruth()

◆ layerId()

◆ phi()

◆ psModule()

◆ r()

◆ setModuleInfo()

◆ z()

Member Data Documentation

◆ assocTP_

◆ assocTPs_

◆ barrel_

◆ endcapRing_

◆ idDet_

◆ layerId_

◆ moduleMaxPhi_

◆ moduleMaxR_

◆ moduleMaxZ_

◆ moduleMinPhi_

◆ moduleMinR_

◆ moduleMinZ_

◆ nStrips_

◆ outerModuleAtSmallerR_

◆ phi_

◆ psModule_

◆ r_

◆ sensorWidth_

◆ settings_

◆ sigmaPar_

◆ sigmaPerp_

◆ stripLength_

◆ stripPitch_

◆ z_