#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
	return TRUE if at least one of the mono and stereo components has badChannels More...

void	init (const MeasurementDet monoDet, const MeasurementDet stereoDet)

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

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

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

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

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

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

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

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 19 of file TkGluedMeasurementDet.h.

Constructor & Destructor Documentation

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

Definition at line 69 of file TkGluedMeasurementDet.cc.

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

Member Function Documentation

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

private

Definition at line 333 of file TkGluedMeasurementDet.cc.

References gather_cfg::cout, prof2calltree::cutoff, TrajectoryStateOnSurface::globalParameters(), TrajectoryStateOnSurface::localError(), TrajectoryStateOnSurface::localPosition(), mag(), bookConverter::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;
   }
   
 }

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

private

Definition at line 321 of file TkGluedMeasurementDet.cc.

References checkHitProjection(), MeasurementDet::fastGeomDet(), and 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());
   }
 }

template<typename Collector >

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

private

Definition at line 171 of file TkGluedMeasurementDet.cc.

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

template<typename Collector >

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

private

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

inlinevirtual

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 MeasurementDet::hasBadComponents().

50 {

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

TkGluedMeasurementDet::monoDet

const TkStripMeasurementDet * monoDet() const

Definition: TkGluedMeasurementDet.h:40

TkGluedMeasurementDet::hasBadComponents

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

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

Definition: TkGluedMeasurementDet.h:50

TkGluedMeasurementDet::stereoDet

const TkStripMeasurementDet * stereoDet() const

Definition: TkGluedMeasurementDet.h:41

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

Definition at line 77 of file TkGluedMeasurementDet.cc.

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

Referenced by MeasurementTrackerImpl::initGluedDet().

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

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

inlinevirtual

return TRUE if both mono and stereo components are active

Implements MeasurementDet.

Definition at line 47 of file TkGluedMeasurementDet.h.

References MeasurementDet::isActive().

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

TkGluedMeasurementDet::monoDet

const TkStripMeasurementDet * monoDet() const

Definition: TkGluedMeasurementDet.h:40

TkStripMeasurementDet::isActive

bool isActive(const MeasurementTrackerEvent &data) const

Is this module active in reconstruction? It must be both 'setActiveThisEvent' and 'setActive'...

Definition: TkStripMeasurementDet.h:118

TkGluedMeasurementDet::stereoDet

const TkStripMeasurementDet * stereoDet() const

Definition: TkGluedMeasurementDet.h:41

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

virtual

Implements MeasurementDet.

Definition at line 114 of file TkGluedMeasurementDet.cc.

References tracking::TempMeasurements::add(), collectRecHits(), F(), MeasurementDet::fastGeomDet(), MeasurementDet::geomDet(), TkStripMeasurementDet::hasAllGoodChannels(), TrajectoryStateOnSurface::hasError(), TkStripMeasurementDet::isActive(), tracking::TempMeasurements::size(), GeomDet::surface(), testStrips(), theCPE, MeasurementDet::theInactiveHit, theMatcher, MeasurementDet::theMissingHit, theMonoDet, theStereoDet, 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;
    
    //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*/ || 
                      (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;
    
 }

const TkStripMeasurementDet* TkGluedMeasurementDet::monoDet ( ) const

inline

Definition at line 40 of file TkGluedMeasurementDet.h.

Referenced by init().

40 { return theMonoDet;}

TkGluedMeasurementDet::theMonoDet

const TkStripMeasurementDet * theMonoDet

Definition: TkGluedMeasurementDet.h:56

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

private

Definition at line 295 of file TkGluedMeasurementDet.cc.

References MeasurementDet::fastGeomDet(), TrajectoryStateOnSurface::globalParameters(), GlobalTrajectoryParameters::momentum(), query::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;
 }

template<typename HitCollector >

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

private

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

private

Definition at line 271 of file TkGluedMeasurementDet.cc.

References MeasurementDet::fastGeomDet(), TrajectoryStateOnSurface::globalParameters(), GlobalTrajectoryParameters::momentum(), query::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;
 }

template<typename HitCollector >

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

private

template<typename Collector >

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

Definition at line 285 of file TkGluedMeasurementDet.cc.

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

template<typename Collector >

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

Definition at line 309 of file TkGluedMeasurementDet.cc.

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

unsigned int TkGluedMeasurementDet::rawId ( ) const

inline

Definition at line 43 of file TkGluedMeasurementDet.h.

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

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

MeasurementDet::fastGeomDet

const GeomDet & fastGeomDet() const

Definition: MeasurementDet.h:77

GeomDet::geographicalId

DetId geographicalId() const

The label of this GeomDet.

Definition: GeomDet.h:77

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

virtual

Implements MeasurementDet.

Definition at line 88 of file TkGluedMeasurementDet.cc.

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

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

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

virtual

Reimplemented from MeasurementDet.

Definition at line 99 of file TkGluedMeasurementDet.cc.

References collectRecHits(), MeasurementDet::fastGeomDet(), TkStripMeasurementDet::isActive(), 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;
   
 }

const GluedGeomDet& TkGluedMeasurementDet::specificGeomDet ( ) const

inline

Definition at line 34 of file TkGluedMeasurementDet.h.

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

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

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

MeasurementDet::fastGeomDet

const GeomDet & fastGeomDet() const

Definition: MeasurementDet.h:77

compareJSON.const

string const

Definition: compareJSON.py:14

GluedGeomDet

Definition: GluedGeomDet.h:7

const TkStripMeasurementDet* TkGluedMeasurementDet::stereoDet ( ) const

inline

Definition at line 41 of file TkGluedMeasurementDet.h.

Referenced by init().

41 { return theStereoDet;}

TkGluedMeasurementDet::theStereoDet

const TkStripMeasurementDet * theStereoDet

Definition: TkGluedMeasurementDet.h:57

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 363 of file TkGluedMeasurementDet.cc.

References MeasurementDet::fastGeomDet(), TrajectoryStateOnSurface::globalParameters(), TrajectoryStateOnSurface::globalPosition(), TrajectoryStateOnSurface::localError(), Topology::measurementError(), Topology::measurementPosition(), GlobalTrajectoryParameters::momentum(), LocalTrajectoryError::positionError(), LocalError::rotate(), TkStripMeasurementDet::specificGeomDet(), StripGeomDetUnit::specificTopology(), mathSSE::sqrt(), GeomDet::surface(), TkStripMeasurementDet::testStrips(), MeasurementError::uu(), PV2DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::x(), LocalError::xx(), LocalError::xy(), PV3DBase< T, PVType, FrameType >::y(), LocalError::yy(), 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);
 } 