#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

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

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

l1tVertexFinder::Stub::~Stub ( )

inline

Definition at line 45 of file Stub.h.

45 {}

Member Function Documentation

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.

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

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

Definition at line 57 of file Stub.cc.

References assocTP_, assocTPs_, 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);
         }
       }
     }
   }

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

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

inline

Definition at line 57 of file Stub.h.

References phi_.

Referenced by Particle.Particle::__str__(), and ntupleDataFormat.Track::phiPull().

57 { return phi_; }

l1tVertexFinder::Stub::phi_

float phi_

Definition: Stub.h:87

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

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

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, 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_FULL_cff::R0, sensorWidth_, TrackerTopology::side(), sigmaPar_, sigmaPerp_, PixelGeomDetUnit::specificTopology(), mathSSE::sqrt(), stripLength_, stripPitch_, DetId::subdetId(), GeomDet::surface(), StripSubdetector::TIB, TrackerTopology::tidRing(), TrackerTopology::tidWheel(), StripSubdetector::TOB, unit(), Bounds::width(), 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)
   }