#include <TkGluedMeasurementDet.h>

Inheritance diagram for TkGluedMeasurementDet:

Classes
class	HitCollectorForFastMeasurements

class	HitCollectorForRecHits

class	HitCollectorForSimpleHits

Public Member Functions
bool	hasBadComponents (const TrajectoryStateOnSurface &tsos, const MeasurementTrackerEvent &data) const override
	return TRUE if at least one of the mono and stereo components has badChannels More...

void	init (const MeasurementDet monoDet, const MeasurementDet stereoDet, const TrackerTopology *tTopo)

bool	isActive (const MeasurementTrackerEvent &data) const override
	return TRUE if both mono and stereo components are active More...

bool	measurements (const TrajectoryStateOnSurface &stateOnThisDet, const MeasurementEstimator &est, const MeasurementTrackerEvent &data, TempMeasurements &result) const override

const TkStripMeasurementDet *	monoDet () const

template<typename Collector >
void	projectOnGluedDet (Collector &collector, const RecHitContainer &hits, const GlobalVector &gdir) const

template<typename Collector >
void	projectOnGluedDet (Collector &collector, std::vector< SiStripRecHit2D > const &hits, const GlobalVector &gdir) const

unsigned int	rawId () const

RecHitContainer	recHits (const TrajectoryStateOnSurface &, const MeasurementTrackerEvent &data) const override

bool	recHits (SimpleHitContainer &result, const TrajectoryStateOnSurface &stateOnThisDet, const MeasurementEstimator &, const MeasurementTrackerEvent &data) const override

const GluedGeomDet &	specificGeomDet () const

const TkStripMeasurementDet *	stereoDet () const

	TkGluedMeasurementDet (const GluedGeomDet gdet, const SiStripRecHitMatcher matcher, const StripClusterParameterEstimator *cpe)

Public Member Functions inherited from MeasurementDet
const GeomDet &	fastGeomDet () const

std::vector< TrajectoryMeasurement >	fastMeasurements (const TrajectoryStateOnSurface &stateOnThisDet, const TrajectoryStateOnSurface &, const Propagator &, const MeasurementEstimator &est, const MeasurementTrackerEvent &data) const

virtual const GeomDet &	geomDet () const

	MeasurementDet (const GeomDet *gdet)

const Surface::PositionType &	position () const

virtual bool	recHits (const TrajectoryStateOnSurface &stateOnThisDet, const MeasurementEstimator &, const MeasurementTrackerEvent &data, RecHitContainer &result, std::vector< float > &) const

const Surface &	surface () const

virtual	~MeasurementDet ()=default

Private Member Functions
void	checkHitProjection (const TrackingRecHit &hit, const TrajectoryStateOnSurface &ts, const GeomDet &det) const

void	checkProjection (const TrajectoryStateOnSurface &ts, const RecHitContainer &monoHits, const RecHitContainer &stereoHits) const

template<typename Collector >
void	collectRecHits (const TrajectoryStateOnSurface &, const MeasurementTrackerEvent &data, Collector &coll) const

template<typename Collector >
void	doubleMatch (const TrajectoryStateOnSurface &ts, const MeasurementTrackerEvent &data, Collector &collector) const

RecHitContainer	projectOnGluedDet (const std::vector< SiStripRecHit2D > &hits, const TrajectoryStateOnSurface &ts) const

template<typename HitCollector >
void	projectOnGluedDet (HitCollector &collector, const std::vector< SiStripRecHit2D > &hits, const GlobalVector &gdir) const

RecHitContainer	projectOnGluedDet (const RecHitContainer &hits, const TrajectoryStateOnSurface &ts) const

template<typename HitCollector >
void	projectOnGluedDet (HitCollector &collector, const RecHitContainer &hits, const GlobalVector &gdir) const

bool	testStrips (const TrajectoryStateOnSurface &tsos, const BoundPlane &gluedPlane, const TkStripMeasurementDet &mdet) const
	Test the strips on one of the two dets with projection. More...

Private Attributes
const StripClusterParameterEstimator *	theCPE

const SiStripRecHitMatcher *	theMatcher

const TkStripMeasurementDet *	theMonoDet

const TkStripMeasurementDet *	theStereoDet

const TrackerTopology *	theTopology

Additional Inherited Members
Public Types inherited from MeasurementDet
typedef TrackingRecHit::ConstRecHitContainer	RecHitContainer

using	SimpleHitContainer = std::vector< BaseTrackerRecHit * >

typedef tracking::TempMeasurements	TempMeasurements

Protected Attributes inherited from MeasurementDet
TrackingRecHit::ConstRecHitPointer	theInactiveHit

TrackingRecHit::ConstRecHitPointer	theMissingHit

Detailed Description

Definition at line 17 of file TkGluedMeasurementDet.h.

Constructor & Destructor Documentation

◆ TkGluedMeasurementDet()

TkGluedMeasurementDet::TkGluedMeasurementDet	(	const GluedGeomDet *	gdet,
		const SiStripRecHitMatcher *	matcher,
		const StripClusterParameterEstimator *	cpe
	)

Definition at line 62 of file TkGluedMeasurementDet.cc.

     : MeasurementDet(gdet),
       theMatcher(matcher),
       theCPE(cpe),
       theMonoDet(nullptr),
       theStereoDet(nullptr),
       theTopology(nullptr) {}

Member Function Documentation

◆ checkHitProjection()

void TkGluedMeasurementDet::checkHitProjection	(	const TrackingRecHit &	hit,
		const TrajectoryStateOnSurface &	ts,
		const GeomDet &	det
	)		const

private

Definition at line 347 of file TkGluedMeasurementDet.cc.

References gather_cfg::cout, EcalSCDynamicDPhiParametersESProducer_cfi::cutoff, HLT_2022v15_cff::eps, TrajectoryStateOnSurface::globalParameters(), TrajectoryStateOnSurface::localError(), TrajectoryStateOnSurface::localPosition(), mag(), SiStripPI::max, GlobalTrajectoryParameters::momentum(), LocalTrajectoryError::positionError(), RecHitPropagator::propagate(), GeomDet::surface(), theCPE, LocalError::xx(), LocalError::xy(), and LocalError::yy().

Referenced by checkProjection().

                                                                          {
   auto&& vl = projectedPos(hit, det, ts.globalParameters().momentum(), theCPE);
   ProjectedSiStripRecHit2D projectedHit(vl.first, vl.second, det, static_cast<SiStripRecHit2D const&>(hit));
 
   RecHitPropagator prop;
   TrajectoryStateOnSurface propState = prop.propagate(hit, det.surface(), ts);
 
   if ((projectedHit.localPosition() - propState.localPosition()).mag() > 0.0001f) {
     std::cout << "PROBLEM: projected and propagated hit positions differ by "
               << (projectedHit.localPosition() - propState.localPosition()).mag() << std::endl;
   }
 
   LocalError le1 = projectedHit.localPositionError();
   LocalError le2 = propState.localError().positionError();
   double eps = 1.e-5;
   double cutoff = 1.e-4;  // if element below cutoff, use absolute instead of relative accuracy
   double maxdiff = std::max(
       std::max(fabs(le1.xx() - le2.xx()) / (cutoff + le1.xx()), fabs(le1.xy() - le2.xy()) / (cutoff + fabs(le1.xy()))),
       fabs(le1.yy() - le2.yy()) / (cutoff + le1.xx()));
   if (maxdiff > eps) {
     std::cout << "PROBLEM: projected and propagated hit errors differ by " << maxdiff << std::endl;
   }
 }

◆ checkProjection()

void TkGluedMeasurementDet::checkProjection	(	const TrajectoryStateOnSurface &	ts,
		const RecHitContainer &	monoHits,
		const RecHitContainer &	stereoHits
	)		const

private

Definition at line 336 of file TkGluedMeasurementDet.cc.

References checkHitProjection(), MeasurementDet::fastGeomDet(), and mps_fire::i.

                                                                                      {
   for (RecHitContainer::const_iterator i = monoHits.begin(); i != monoHits.end(); ++i) {
     checkHitProjection(**i, ts, fastGeomDet());
   }
   for (RecHitContainer::const_iterator i = stereoHits.begin(); i != stereoHits.end(); ++i) {
     checkHitProjection(**i, ts, fastGeomDet());
   }
 }

◆ collectRecHits()

template<typename Collector >

void TkGluedMeasurementDet::collectRecHits	(	const TrajectoryStateOnSurface &	ts,
		const MeasurementTrackerEvent &	data,
		Collector &	coll
	)		const

private

Definition at line 170 of file TkGluedMeasurementDet.cc.

References data, TrajectoryStateOnSurface::globalParameters(), TrackingRecHit::hit(), TrajectoryStateOnSurface::isValid(), TrajectoryStateOnSurface::localDirection(), SiStripRecHitMatcher::match(), GlobalTrajectoryParameters::momentum(), MeasurementDet::position(), projectOnGluedDet(), TkStripMeasurementDet::recHits(), TkStripMeasurementDet::simpleRecHits(), specificGeomDet(), MeasurementDet::surface(), theMatcher, theMonoDet, theStereoDet, GloballyPositioned< T >::toLocal(), and HcalDetIdTransform::transform().

Referenced by measurements(), and recHits().

                                                                        {
   //------ WARNING: here ts is used as it is on the mono/stereo surface.
   //-----           A further propagation is necessary.
   //-----           To limit the problem, the SimpleCPE should be used
   RecHitContainer monoHits = theMonoDet->recHits(ts, data);
   GlobalVector glbDir = (ts.isValid() ? ts.globalParameters().momentum() : position() - GlobalPoint(0, 0, 0));
 
   //edm::LogWarning("TkGluedMeasurementDet::recHits") << "Query-for-detid-" << theGeomDet->geographicalId().rawId();
 
   //checkProjection(ts, monoHits, stereoHits);
 
   if (monoHits.empty()) {
     // make stereo TTRHs and project them
     projectOnGluedDet(collector, theStereoDet->recHits(ts, data), glbDir);
   } else {
     // collect simple stereo hits
     std::vector<SiStripRecHit2D> simpleSteroHitsByValue;
     theStereoDet->simpleRecHits(ts, data, simpleSteroHitsByValue);
 
     if (simpleSteroHitsByValue.empty()) {
       projectOnGluedDet(collector, monoHits, glbDir);
     } else {
       LocalVector tkDir = (ts.isValid() ? ts.localDirection() : surface().toLocal(position() - GlobalPoint(0, 0, 0)));
       SiStripRecHitMatcher::SimpleHitCollection vsStereoHits;
       vsStereoHits.resize(simpleSteroHitsByValue.size());
       std::transform(
           simpleSteroHitsByValue.begin(), simpleSteroHitsByValue.end(), vsStereoHits.begin(), take_address());
 
       // convert mono hits to type expected by matcher
       for (RecHitContainer::const_iterator monoHit = monoHits.begin(); monoHit != monoHits.end(); ++monoHit) {
         const TrackingRecHit* tkhit = (**monoHit).hit();
         const SiStripRecHit2D* verySpecificMonoHit = reinterpret_cast<const SiStripRecHit2D*>(tkhit);
         theMatcher->match(verySpecificMonoHit,
                           vsStereoHits.begin(),
                           vsStereoHits.end(),
                           collector.collector(),
                           &specificGeomDet(),
                           tkDir);
 
         if (collector.hasNewMatchedHits()) {
           collector.clearNewMatchedHitsFlag();
         } else {
           collector.addProjected(**monoHit, glbDir);
         }
       }  // loop on mono hit
     }
     //GIO// std::cerr << "TkGluedMeasurementDet hits " << monoHits.size() << "/" << stereoHits.size() << " => " << result.size() << std::endl;
   }
 }

◆ doubleMatch()

template<typename Collector >

void TkGluedMeasurementDet::doubleMatch	(	const TrajectoryStateOnSurface &	ts,
		const MeasurementTrackerEvent &	data,
		Collector &	collector
	)		const

private

◆ hasBadComponents()

bool TkGluedMeasurementDet::hasBadComponents	(	const TrajectoryStateOnSurface &	tsos,
		const MeasurementTrackerEvent &	data
	)		const

inlineoverridevirtual

return TRUE if at least one of the mono and stereo components has badChannels

Implements MeasurementDet.

Definition at line 50 of file TkGluedMeasurementDet.h.

References data, and MeasurementDet::hasBadComponents().

                                                                                                                   {
     return (monoDet()->hasBadComponents(tsos, data) || stereoDet()->hasBadComponents(tsos, data));
   }

◆ init()

void TkGluedMeasurementDet::init	(	const MeasurementDet *	monoDet,
		const MeasurementDet *	stereoDet,
		const TrackerTopology *	tTopo
	)

Definition at line 72 of file TkGluedMeasurementDet.cc.

References monoDet(), stereoDet(), theMonoDet, theStereoDet, and theTopology.

Referenced by MeasurementTrackerImpl::initGluedDet().

                                                                {
   theMonoDet = dynamic_cast<const TkStripMeasurementDet*>(monoDet);
   theStereoDet = dynamic_cast<const TkStripMeasurementDet*>(stereoDet);
   theTopology = tTopo;
 
   if ((theMonoDet == nullptr) || (theStereoDet == nullptr)) {
     throw MeasurementDetException(
         "TkGluedMeasurementDet ERROR: Trying to glue a det which is not a TkStripMeasurementDet");
   }
 }

◆ isActive()

bool TkGluedMeasurementDet::isActive ( const MeasurementTrackerEvent & data ) const

inlineoverridevirtual

return TRUE if both mono and stereo components are active

Implements MeasurementDet.

Definition at line 45 of file TkGluedMeasurementDet.h.

References data, and MeasurementDet::isActive().

                                                                     {
     return monoDet()->isActive(data) && stereoDet()->isActive(data);
   }

◆ measurements()

bool TkGluedMeasurementDet::measurements	(	const TrajectoryStateOnSurface &	stateOnThisDet,
		const MeasurementEstimator &	est,
		const MeasurementTrackerEvent &	data,
		TempMeasurements &	result
	)		const

overridevirtual

Implements MeasurementDet.

Definition at line 107 of file TkGluedMeasurementDet.cc.

References collectRecHits(), data, F(), MeasurementDet::fastGeomDet(), GeomDet::geographicalId(), MeasurementDet::geomDet(), TrajectoryStateOnSurface::globalMomentum(), TkStripMeasurementDet::hasAllGoodChannels(), TrajectoryStateOnSurface::hasError(), TkStripMeasurementDet::isActive(), cmsLHEtoEOSManager::l, MeasurementEstimator::minPt2ForHitRecoveryInGluedDet(), PV3DBase< T, PVType, FrameType >::perp2(), mps_fire::result, DetId::subdetId(), GeomDet::surface(), testStrips(), theCPE, MeasurementDet::theInactiveHit, theMatcher, MeasurementDet::theMissingHit, theMonoDet, theStereoDet, theTopology, TrackerTopology::tobLayer(), and UNLIKELY.

                                                                          {
   if UNLIKELY ((!theMonoDet->isActive(data)) && (!theStereoDet->isActive(data))) {
     //     LogDebug("TkStripMeasurementDet") << " DetID " << geomDet().geographicalId().rawId() << " (glued) fully inactive";
     result.add(theInactiveHit, 0.F);
     return true;
   }
 
   auto oldSize = result.size();
 
   HitCollectorForFastMeasurements collector(&fastGeomDet(), theMatcher, theCPE, stateOnThisDet, est, result);
   collectRecHits(stateOnThisDet, data, collector);
 
   if (result.size() > oldSize)
     return true;
 
   auto id = geomDet().geographicalId().subdetId() - 3;
   auto l = theTopology->tobLayer(geomDet().geographicalId());
   bool killHIP = (1 == l) && (2 == id);  //TOB1
   killHIP &= stateOnThisDet.globalMomentum().perp2() > est.minPt2ForHitRecoveryInGluedDet();
   if (killHIP) {
     result.add(theInactiveHit, 0.F);
     return true;
   }
 
   //LogDebug("TkStripMeasurementDet") << "No hit found on TkGlued. Testing strips...  ";
   const BoundPlane& gluedPlane = geomDet().surface();
   if (  // sorry for the big IF, but I want to exploit short-circuiting of logic
       stateOnThisDet.hasError() &&
       (/* do this only if the state has uncertainties, otherwise it will throw 
                      (states without uncertainties are passed to this code from seeding */
        (theMonoDet->isActive(data) &&
         (theMonoDet->hasAllGoodChannels() || testStrips(stateOnThisDet, gluedPlane, *theMonoDet))) /*Mono OK*/
        &&                                                                                          // was ||
        (theStereoDet->isActive(data) &&
         (theStereoDet->hasAllGoodChannels() || testStrips(stateOnThisDet, gluedPlane, *theStereoDet))) /*Stereo OK*/
        ) /* State has errors */
   ) {
     result.add(theMissingHit, 0.F);
     return false;
   }
   result.add(theInactiveHit, 0.F);
   return true;
 }

◆ monoDet()

const TkStripMeasurementDet* TkGluedMeasurementDet::monoDet ( ) const

inline

Definition at line 39 of file TkGluedMeasurementDet.h.

Referenced by init().

39 { return theMonoDet; }

TkGluedMeasurementDet::theMonoDet

const TkStripMeasurementDet * theMonoDet

Definition: TkGluedMeasurementDet.h:57

◆ projectOnGluedDet() [1/6]

TkGluedMeasurementDet::RecHitContainer TkGluedMeasurementDet::projectOnGluedDet	(	const std::vector< SiStripRecHit2D > &	hits,
		const TrajectoryStateOnSurface &	ts
	)		const

private

Definition at line 316 of file TkGluedMeasurementDet.cc.

References runTheMatrix::const, MeasurementDet::fastGeomDet(), TrajectoryStateOnSurface::globalParameters(), hfClusterShapes_cfi::hits, GlobalTrajectoryParameters::momentum(), eostools::move(), mps_fire::result, and theCPE.

Referenced by collectRecHits().

                                                                                       {
   RecHitContainer result;
   for (auto const& hit : hits) {
     auto&& vl = projectedPos(hit, fastGeomDet(), ts.globalParameters().momentum(), theCPE);
     auto&& phit = std::make_shared<ProjectedSiStripRecHit2D>(
         vl.first, vl.second, fastGeomDet(), static_cast<SiStripRecHit2D const&>(hit));
     result.push_back(std::move(phit));
   }
   return result;
 }

◆ projectOnGluedDet() [2/6]

template<typename HitCollector >

void TkGluedMeasurementDet::projectOnGluedDet	(	HitCollector &	collector,
		const std::vector< SiStripRecHit2D > &	hits,
		const GlobalVector &	gdir
	)		const

private

◆ projectOnGluedDet() [3/6]

TkGluedMeasurementDet::RecHitContainer TkGluedMeasurementDet::projectOnGluedDet	(	const RecHitContainer &	hits,
		const TrajectoryStateOnSurface &	ts
	)		const

private

Definition at line 295 of file TkGluedMeasurementDet.cc.

References runTheMatrix::const, MeasurementDet::fastGeomDet(), TrajectoryStateOnSurface::globalParameters(), hfClusterShapes_cfi::hits, GlobalTrajectoryParameters::momentum(), eostools::move(), mps_fire::result, and theCPE.

                                                                            {
   RecHitContainer result;
   for (auto const& hit : hits) {
     auto&& vl = projectedPos(*hit, fastGeomDet(), ts.globalParameters().momentum(), theCPE);
     auto&& phit = std::make_shared<ProjectedSiStripRecHit2D>(
         vl.first, vl.second, fastGeomDet(), static_cast<SiStripRecHit2D const&>(*hit));
     result.push_back(std::move(phit));
   }
   return result;
 }

◆ projectOnGluedDet() [4/6]

template<typename HitCollector >

void TkGluedMeasurementDet::projectOnGluedDet	(	HitCollector &	collector,
		const RecHitContainer &	hits,
		const GlobalVector &	gdir
	)		const

private

◆ projectOnGluedDet() [5/6]

template<typename Collector >

void TkGluedMeasurementDet::projectOnGluedDet	(	Collector &	collector,
		const RecHitContainer &	hits,
		const GlobalVector &	gdir
	)		const

Definition at line 308 of file TkGluedMeasurementDet.cc.

References hfClusterShapes_cfi::hits.

                                                                               {
   for (RecHitContainer::const_iterator ihit = hits.begin(); ihit != hits.end(); ihit++) {
     collector.addProjected(**ihit, gdir);
   }
 }

◆ projectOnGluedDet() [6/6]

template<typename Collector >

void TkGluedMeasurementDet::projectOnGluedDet	(	Collector &	collector,
		std::vector< SiStripRecHit2D > const &	hits,
		const GlobalVector &	gdir
	)		const

Definition at line 329 of file TkGluedMeasurementDet.cc.

References hfClusterShapes_cfi::hits.

                                                                               {
   for (auto const& hit : hits)
     collector.addProjected(hit, gdir);
 }

◆ rawId()

unsigned int TkGluedMeasurementDet::rawId ( ) const

inline

Definition at line 42 of file TkGluedMeasurementDet.h.

References MeasurementDet::fastGeomDet(), and GeomDet::geographicalId().

42 { return fastGeomDet().geographicalId(); }

GeomDet::geographicalId

DetId geographicalId() const

The label of this GeomDet.

Definition: GeomDet.h:64

MeasurementDet::fastGeomDet

const GeomDet & fastGeomDet() const

Definition: MeasurementDet.h:80

◆ recHits() [1/2]

TkGluedMeasurementDet::RecHitContainer TkGluedMeasurementDet::recHits	(	const TrajectoryStateOnSurface &	ts,
		const MeasurementTrackerEvent &	data
	)		const

overridevirtual

Implements MeasurementDet.

Definition at line 85 of file TkGluedMeasurementDet.cc.

References collectRecHits(), data, MeasurementDet::fastGeomDet(), mps_fire::result, theCPE, and theMatcher.

                                                                                                                  {
   RecHitContainer result;
   HitCollectorForRecHits collector(&fastGeomDet(), theMatcher, theCPE, result);
   collectRecHits(ts, data, collector);
   return result;
 }

◆ recHits() [2/2]

bool TkGluedMeasurementDet::recHits	(	SimpleHitContainer &	result,
		const TrajectoryStateOnSurface &	stateOnThisDet,
		const MeasurementEstimator &	est,
		const MeasurementTrackerEvent &	data
	)		const

overridevirtual

Reimplemented from MeasurementDet.

Definition at line 94 of file TkGluedMeasurementDet.cc.

References collectRecHits(), data, MeasurementDet::fastGeomDet(), TkStripMeasurementDet::isActive(), mps_fire::result, theCPE, theMatcher, theMonoDet, theStereoDet, and UNLIKELY.

                                                                                {
   if UNLIKELY ((!theMonoDet->isActive(data)) && (!theStereoDet->isActive(data)))
     return false;
   auto oldSize = result.size();
   HitCollectorForSimpleHits collector(&fastGeomDet(), theMatcher, theCPE, stateOnThisDet, est, result);
   collectRecHits(stateOnThisDet, data, collector);
 
   return result.size() > oldSize;
 }

◆ specificGeomDet()

const GluedGeomDet& TkGluedMeasurementDet::specificGeomDet ( ) const

inline

Definition at line 32 of file TkGluedMeasurementDet.h.

References runTheMatrix::const, and MeasurementDet::fastGeomDet().

Referenced by collectRecHits(), and MeasurementTrackerImpl::initGluedDet().

32 { return static_cast<GluedGeomDet const&>(fastGeomDet()); }

runTheMatrix.const

const

Definition: runTheMatrix.py:339

MeasurementDet::fastGeomDet

const GeomDet & fastGeomDet() const

Definition: MeasurementDet.h:80

GluedGeomDet

Definition: GluedGeomDet.h:7

◆ stereoDet()

const TkStripMeasurementDet* TkGluedMeasurementDet::stereoDet ( ) const

inline

Definition at line 40 of file TkGluedMeasurementDet.h.

Referenced by init().

40 { return theStereoDet; }

TkGluedMeasurementDet::theStereoDet

const TkStripMeasurementDet * theStereoDet

Definition: TkGluedMeasurementDet.h:58

◆ testStrips()

bool TkGluedMeasurementDet::testStrips	(	const TrajectoryStateOnSurface &	tsos,
		const BoundPlane &	gluedPlane,
		const TkStripMeasurementDet &	mdet
	)		const

private

Test the strips on one of the two dets with projection.

Definition at line 373 of file TkGluedMeasurementDet.cc.

References dumpMFGeometry_cfg::delta, submitPVResolutionJobs::err, MeasurementDet::fastGeomDet(), TrajectoryStateOnSurface::globalParameters(), TrajectoryStateOnSurface::globalPosition(), TrajectoryStateOnSurface::localError(), Topology::measurementError(), Topology::measurementPosition(), GlobalTrajectoryParameters::momentum(), LocalTrajectoryError::positionError(), TkStripMeasurementDet::specificGeomDet(), StripGeomDetUnit::specificTopology(), mathSSE::sqrt(), GeomDet::surface(), TkStripMeasurementDet::testStrips(), MeasurementError::uu(), PV2DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

Referenced by measurements().

                                                                                 {
   // from TrackingRecHitProjector
   const GeomDet& det = mdet.fastGeomDet();
   const BoundPlane& stripPlane = det.surface();
 
   //LocalPoint glp = tsos.localPosition();
   LocalError err = tsos.localError().positionError();
   /*LogDebug("TkStripMeasurementDet") << 
     "Testing local pos glued: " << glp << 
     " local err glued: " << tsos.localError().positionError() << 
       " in? " << gluedPlane.bounds().inside(glp) <<
       " in(3s)? " << gluedPlane.bounds().inside(glp, err, 3.0f);*/
 
   GlobalVector gdir = tsos.globalParameters().momentum();
 
   LocalPoint slp = stripPlane.toLocal(tsos.globalPosition());
   LocalVector sld = stripPlane.toLocal(gdir);
 
   double delta = stripPlane.localZ(tsos.globalPosition());
   LocalPoint pos = slp - sld * delta / sld.z();
 
   // now the error
   LocalVector hitXAxis = stripPlane.toLocal(gluedPlane.toGlobal(LocalVector(1, 0, 0)));
   if (stripPlane.normalVector().dot(gluedPlane.normalVector()) < 0) {
     // the two planes are inverted, and the correlation element must change sign
     err = LocalError(err.xx(), -err.xy(), err.yy());
   }
   LocalError rotatedError = err.rotate(hitXAxis.x(), hitXAxis.y());
 
   /* // This is probably meaningless 
       LogDebug("TkStripMeasurementDet") << 
       "Testing local pos on strip (SLP): " << slp << 
       " in? :" << stripPlane.bounds().inside(slp) <<
       " in(3s)? :" << stripPlane.bounds().inside(slp, rotatedError, 3.0f);
       // but it helps to test bugs in the formula for POS */
   /*LogDebug("TkStripMeasurementDet") << 
      "Testing local pos strip: " << pos << 
      " in? " << stripPlane.bounds().inside(pos) <<
      " in(3s)? " << stripPlane.bounds().inside(pos, rotatedError, 3.0f);*/
 
   // now we need to convert to MeasurementFrame
   const StripTopology& topo = mdet.specificGeomDet().specificTopology();
   float utraj = topo.measurementPosition(pos).x();
   float uerr = std::sqrt(topo.measurementError(pos, rotatedError).uu());
   return mdet.testStrips(utraj, uerr);
 }

Member Data Documentation

◆ theCPE

const StripClusterParameterEstimator* TkGluedMeasurementDet::theCPE

private

Definition at line 56 of file TkGluedMeasurementDet.h.

Referenced by checkHitProjection(), measurements(), projectOnGluedDet(), and recHits().

◆ theMatcher

const SiStripRecHitMatcher* TkGluedMeasurementDet::theMatcher

private

Definition at line 55 of file TkGluedMeasurementDet.h.

Referenced by collectRecHits(), measurements(), and recHits().

◆ theMonoDet

const TkStripMeasurementDet* TkGluedMeasurementDet::theMonoDet

private

Definition at line 57 of file TkGluedMeasurementDet.h.

Referenced by collectRecHits(), init(), measurements(), and recHits().

◆ theStereoDet

const TkStripMeasurementDet* TkGluedMeasurementDet::theStereoDet

private

Definition at line 58 of file TkGluedMeasurementDet.h.

Referenced by collectRecHits(), init(), measurements(), and recHits().

◆ theTopology

const TrackerTopology* TkGluedMeasurementDet::theTopology

private

Definition at line 59 of file TkGluedMeasurementDet.h.

Referenced by init(), and measurements().

Classes

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ TkGluedMeasurementDet()

Member Function Documentation

◆ checkHitProjection()

◆ checkProjection()

◆ collectRecHits()

◆ doubleMatch()

◆ hasBadComponents()

◆ init()

◆ isActive()

◆ measurements()

◆ monoDet()

◆ projectOnGluedDet() [1/6]

◆ projectOnGluedDet() [2/6]

◆ projectOnGluedDet() [3/6]

◆ projectOnGluedDet() [4/6]

◆ projectOnGluedDet() [5/6]

◆ projectOnGluedDet() [6/6]

◆ rawId()

◆ recHits() [1/2]

◆ recHits() [2/2]

◆ specificGeomDet()

◆ stereoDet()

◆ testStrips()

Member Data Documentation

◆ theCPE

◆ theMatcher

◆ theMonoDet

◆ theStereoDet

◆ theTopology