TkGluedMeasurementDet Class Reference

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

Constructor & Destructor Documentation

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

Definition at line 71 of file TkGluedMeasurementDet.cc.

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

Member Function Documentation

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

private

Definition at line 359 of file TkGluedMeasurementDet.cc.

References gather_cfg::cout, qjetsadder_cfi::cutoff, 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;
  }
  
}

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

private

Definition at line 347 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());
  }
}

template<typename Collector >

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

private

Definition at line 187 of file TkGluedMeasurementDet.cc.

References TrajectoryStateOnSurface::globalParameters(), TrackingRecHit::hit(), TrajectoryStateOnSurface::isValid(), TrajectoryStateOnSurface::localDirection(), funct::m, match(), SiStripRecHitMatcher::match(), GlobalTrajectoryParameters::momentum(), MeasurementDet::position(), projectOnGluedDet(), TkStripMeasurementDet::recHits(), alignCSCRings::s, TkStripMeasurementDet::simpleRecHits(), specificGeomDet(), trackingPlots::stat, MeasurementDet::surface(), lumiQTWidget::t, theMatcher, theMonoDet, theStereoDet, GloballyPositioned< T >::toLocal(), and create_public_lumi_plots::transform.

Referenced by measurements(), take_address::operator()(), 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

Referenced by take_address::operator()().

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

References MeasurementDet::hasBadComponents().

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

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

Definition at line 80 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");
  }
}

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

inlineoverridevirtual

return TRUE if both mono and stereo components are active

Implements MeasurementDet.

Definition at line 48 of file TkGluedMeasurementDet.h.

References MeasurementDet::isActive().

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

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

overridevirtual

Implements MeasurementDet.

Definition at line 119 of file TkGluedMeasurementDet.cc.

References tracking::TempMeasurements::add(), collectRecHits(), F(), MeasurementDet::fastGeomDet(), GeomDet::geographicalId(), MeasurementDet::geomDet(), TrajectoryStateOnSurface::globalMomentum(), TkStripMeasurementDet::hasAllGoodChannels(), TrajectoryStateOnSurface::hasError(), TkStripMeasurementDet::isActive(), checklumidiff::l, MeasurementEstimator::minPt2ForHitRecoveryInGluedDet(), PV3DBase< T, PVType, FrameType >::perp2(), tracking::TempMeasurements::size(), 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;
   
}

const TkStripMeasurementDet* TkGluedMeasurementDet::monoDet ( ) const

inline

Definition at line 41 of file TkGluedMeasurementDet.h.

Referenced by init().

{ return theMonoDet;} 

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

private

Referenced by collectRecHits(), and projectOnGluedDet().

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

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

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

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

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

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

unsigned int TkGluedMeasurementDet::rawId ( ) const

inline

Definition at line 44 of file TkGluedMeasurementDet.h.

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

{ return  fastGeomDet().geographicalId(); }

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

overridevirtual

Implements MeasurementDet.

Definition at line 93 of file TkGluedMeasurementDet.cc.

References collectRecHits(), MeasurementDet::fastGeomDet(), mps_fire::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

overridevirtual

Reimplemented from MeasurementDet.

Definition at line 104 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 35 of file TkGluedMeasurementDet.h.

References MeasurementDet::fastGeomDet(), and MeasurementDet::measurements().

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

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

const TkStripMeasurementDet* TkGluedMeasurementDet::stereoDet ( ) const

inline

Definition at line 42 of file TkGluedMeasurementDet.h.

Referenced by init().

{ return theStereoDet;} 

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

Member Data Documentation

const StripClusterParameterEstimator* TkGluedMeasurementDet::theCPE

private

Definition at line 56 of file TkGluedMeasurementDet.h.

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

const SiStripRecHitMatcher* TkGluedMeasurementDet::theMatcher

private

Definition at line 55 of file TkGluedMeasurementDet.h.

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

const TkStripMeasurementDet* TkGluedMeasurementDet::theMonoDet

private

Definition at line 57 of file TkGluedMeasurementDet.h.

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

const TkStripMeasurementDet* TkGluedMeasurementDet::theStereoDet

private

Definition at line 58 of file TkGluedMeasurementDet.h.

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

const TrackerTopology* TkGluedMeasurementDet::theTopology

private

Definition at line 59 of file TkGluedMeasurementDet.h.

Referenced by init(), and measurements().