TkClonerImpl Class Reference

#include <TkClonerImpl.h>

Inheritance diagram for TkClonerImpl:

Public Member Functions
virtual TrackingRecHit::ConstRecHitPointer	makeShared (SiPixelRecHit const &hit, TrajectoryStateOnSurface const &tsos) const override
virtual TrackingRecHit::ConstRecHitPointer	makeShared (SiStripRecHit2D const &hit, TrajectoryStateOnSurface const &tsos) const override
virtual TrackingRecHit::ConstRecHitPointer	makeShared (SiStripRecHit1D const &hit, TrajectoryStateOnSurface const &tsos) const override
virtual TrackingRecHit::ConstRecHitPointer	makeShared (SiStripMatchedRecHit2D const &hit, TrajectoryStateOnSurface const &tsos) const override
virtual TrackingRecHit::ConstRecHitPointer	makeShared (ProjectedSiStripRecHit2D const &hit, TrajectoryStateOnSurface const &tsos) const override
virtual std::unique_ptr < SiPixelRecHit >	operator() (SiPixelRecHit const &hit, TrajectoryStateOnSurface const &tsos) const override
virtual std::unique_ptr < SiStripRecHit2D >	operator() (SiStripRecHit2D const &hit, TrajectoryStateOnSurface const &tsos) const override
virtual std::unique_ptr < SiStripRecHit1D >	operator() (SiStripRecHit1D const &hit, TrajectoryStateOnSurface const &tsos) const override
virtual std::unique_ptr < SiStripMatchedRecHit2D >	operator() (SiStripMatchedRecHit2D const &hit, TrajectoryStateOnSurface const &tsos) const override
virtual std::unique_ptr < ProjectedSiStripRecHit2D >	operator() (ProjectedSiStripRecHit2D const &hit, TrajectoryStateOnSurface const &tsos) const override
std::unique_ptr < ProjectedSiStripRecHit2D >	project (SiStripMatchedRecHit2D const &hit, bool mono, TrajectoryStateOnSurface const &tsos) const
	TkClonerImpl ()
	TkClonerImpl (const PixelClusterParameterEstimator *ipixelCPE, const StripClusterParameterEstimator *istripCPE, const SiStripRecHitMatcher *iMatcher)
Public Member Functions inherited from TkCloner
TrackingRecHit *operator()[[cms TrackingRecHit::ConstRecHitPointer	makeShared (TrackingRecHit::ConstRecHitPointer const &hit, TrajectoryStateOnSurface const &tsos) const

Private Attributes
const PixelClusterParameterEstimator *	pixelCPE
const StripClusterParameterEstimator *	stripCPE
const SiStripRecHitMatcher *	theMatcher

Detailed Description

Definition at line 12 of file TkClonerImpl.h.

Constructor & Destructor Documentation

TkClonerImpl::TkClonerImpl ( )

inline

Definition at line 14 of file TkClonerImpl.h.

{}

TkClonerImpl::TkClonerImpl ( const PixelClusterParameterEstimator *  ipixelCPE, const StripClusterParameterEstimator *  istripCPE, const SiStripRecHitMatcher *  iMatcher  )

inline

Definition at line 15 of file TkClonerImpl.h.

            : pixelCPE(ipixelCPE), stripCPE(istripCPE), theMatcher(iMatcher){}

Member Function Documentation

TrackingRecHit::ConstRecHitPointer TkClonerImpl::makeShared ( SiPixelRecHit const &  hit, TrajectoryStateOnSurface const &  tsos  ) const

overridevirtual

Implements TkCloner.

Definition at line 47 of file TkClonerImpl.cc.

References SiPixelRecHit::cluster(), TrackingRecHit::det(), TrackerSingleRecHit::detUnit(), PixelClusterParameterEstimator::getParameters(), and pixelCPE.

Referenced by DAFTrackProducerAlgorithm::filter(), TrackingRecHitPropagator::project(), and SiTrackerMultiRecHitUpdator::update().

                                                                                                                               {
 // std::cout << "cloning " << typeid(hit).name() << std::endl;
  const SiPixelCluster& clust = *hit.cluster();  
  auto && params = pixelCPE->getParameters( clust, *hit.detUnit(), tsos);
  return std::make_shared<SiPixelRecHit>(std::get<0>(params), std::get<1>(params), std::get<2>(params), *hit.det(), hit.cluster());
}

TrackingRecHit::ConstRecHitPointer TkClonerImpl::makeShared ( SiStripRecHit2D const &  hit, TrajectoryStateOnSurface const &  tsos  ) const

overridevirtual

FIXME: this only uses the first cluster and ignores the others

Implements TkCloner.

Definition at line 54 of file TkClonerImpl.cc.

References TrackingRecHit::det(), TrackerSingleRecHit::detUnit(), StripClusterParameterEstimator::localParameters(), TrackerSingleRecHit::omniCluster(), TrackerSingleRecHit::stripCluster(), and stripCPE.

                                                                                                                                 {
  // std::cout << "cloning " << typeid(hit).name()  << std::endl;
    const SiStripCluster&  clust = hit.stripCluster();  
    StripClusterParameterEstimator::LocalValues lv = 
      stripCPE->localParameters( clust, *hit.detUnit(), tsos);
    return std::make_shared<SiStripRecHit2D>(lv.first, lv.second, *hit.det(), hit.omniCluster());
}

TrackingRecHit::ConstRecHitPointer TkClonerImpl::makeShared ( SiStripRecHit1D const &  hit, TrajectoryStateOnSurface const &  tsos  ) const

overridevirtual

FIXME: this only uses the first cluster and ignores the others

Implements TkCloner.

Definition at line 63 of file TkClonerImpl.cc.

References TrackingRecHit::det(), TrackerSingleRecHit::detUnit(), asciidump::le, StripClusterParameterEstimator::localParameters(), max(), TrackerSingleRecHit::omniCluster(), TrackerSingleRecHit::stripCluster(), and stripCPE.

                                                                                                                                 {
  // std::cout << "cloning " << typeid(hit).name()  << std::endl;
  const SiStripCluster&  clust = hit.stripCluster();  
  StripClusterParameterEstimator::LocalValues lv = 
    stripCPE->localParameters( clust, *hit.detUnit(), tsos);
  LocalError le(lv.second.xx(),0.,std::numeric_limits<float>::max()); //Correct??
  return  std::make_shared<SiStripRecHit1D>(lv.first, le, *hit.det(), hit.omniCluster());
}

TrackingRecHit::ConstRecHitPointer TkClonerImpl::makeShared ( SiStripMatchedRecHit2D const &  hit, TrajectoryStateOnSurface const &  tsos  ) const

overridevirtual

std::cout << "cloning " << typeid(hit).name() << std::endl;

Implements TkCloner.

Definition at line 132 of file TkClonerImpl.cc.

                                                                                                                                        {
  return TrackingRecHit::ConstRecHitPointer((*this)(hit,tsos));
}

TrackingRecHit::ConstRecHitPointer TkClonerImpl::makeShared ( ProjectedSiStripRecHit2D const &  hit, TrajectoryStateOnSurface const &  tsos  ) const

overridevirtual

Implements TkCloner.

Definition at line 137 of file TkClonerImpl.cc.

                                                                                                                                          {
  // std::cout << "cloning " << typeid(hit).name()  << std::endl;
   return TrackingRecHit::ConstRecHitPointer((*this)(hit,tsos));
}

std::unique_ptr< SiPixelRecHit > TkClonerImpl::operator() ( SiPixelRecHit const &  hit, TrajectoryStateOnSurface const &  tsos  ) const

overridevirtual

Implements TkCloner.

Definition at line 24 of file TkClonerImpl.cc.

References SiPixelRecHit::cluster(), TrackingRecHit::det(), TrackerSingleRecHit::detUnit(), PixelClusterParameterEstimator::getParameters(), and pixelCPE.

                                                                                                                            {
  const SiPixelCluster& clust = *hit.cluster();  
  auto && params = pixelCPE->getParameters( clust, *hit.detUnit(), tsos);
  return std::unique_ptr<SiPixelRecHit>(new SiPixelRecHit(std::get<0>(params), std::get<1>(params), std::get<2>(params), *hit.det(), hit.cluster()));
}

std::unique_ptr< SiStripRecHit2D > TkClonerImpl::operator() ( SiStripRecHit2D const &  hit, TrajectoryStateOnSurface const &  tsos  ) const

overridevirtual

FIXME: this only uses the first cluster and ignores the others

Implements TkCloner.

Definition at line 30 of file TkClonerImpl.cc.

References TrackingRecHit::det(), TrackerSingleRecHit::detUnit(), StripClusterParameterEstimator::localParameters(), TrackerSingleRecHit::omniCluster(), TrackerSingleRecHit::stripCluster(), and stripCPE.

                                                                                                                               {
    const SiStripCluster&  clust = hit.stripCluster();  
    StripClusterParameterEstimator::LocalValues lv = 
      stripCPE->localParameters( clust, *hit.detUnit(), tsos);
    return std::unique_ptr<SiStripRecHit2D>{new SiStripRecHit2D(lv.first, lv.second, *hit.det(), hit.omniCluster())};
}

std::unique_ptr< SiStripRecHit1D > TkClonerImpl::operator() ( SiStripRecHit1D const &  hit, TrajectoryStateOnSurface const &  tsos  ) const

overridevirtual

FIXME: this only uses the first cluster and ignores the others

Implements TkCloner.

Definition at line 38 of file TkClonerImpl.cc.

References TrackingRecHit::det(), TrackerSingleRecHit::detUnit(), asciidump::le, StripClusterParameterEstimator::localParameters(), max(), TrackerSingleRecHit::omniCluster(), TrackerSingleRecHit::stripCluster(), and stripCPE.

                                                                                                                               {
  const SiStripCluster&  clust = hit.stripCluster();  
  StripClusterParameterEstimator::LocalValues lv = 
    stripCPE->localParameters( clust, *hit.detUnit(), tsos);
  LocalError le(lv.second.xx(),0.,std::numeric_limits<float>::max()); //Correct??
  return std::unique_ptr<SiStripRecHit1D>{new SiStripRecHit1D(lv.first, le, *hit.det(), hit.omniCluster())};
}

std::unique_ptr< SiStripMatchedRecHit2D > TkClonerImpl::operator() ( SiStripMatchedRecHit2D const &  hit, TrajectoryStateOnSurface const &  tsos  ) const

overridevirtual

Implements TkCloner.

Definition at line 103 of file TkClonerImpl.cc.

References SiStripMatchedRecHit2D::clone(), TrackingRecHit::det(), TrajectoryStateOnSurface::isValid(), TrajectoryStateOnSurface::localDirection(), StripClusterParameterEstimator::localParameters(), SiStripRecHitMatcher::match(), SiStripMatchedRecHit2D::monoCluster(), SiStripMatchedRecHit2D::monoClusterRef(), GluedGeomDet::monoDet(), GeomDet::position(), SiStripMatchedRecHit2D::stereoCluster(), SiStripMatchedRecHit2D::stereoClusterRef(), GluedGeomDet::stereoDet(), stripCPE, GeomDet::surface(), groupFilesInBlocks::temp, theMatcher, and GloballyPositioned< T >::toLocal().

                                                                                                                                             {
    const GeomDet * det = hit.det();
    const GluedGeomDet *gdet = static_cast<const GluedGeomDet *> (det);
    LocalVector tkDir = (tsos.isValid() ? tsos.localDirection() : 
             det->surface().toLocal( det->position()-GlobalPoint(0,0,0)));

    const SiStripCluster& monoclust   = hit.monoCluster();  
    const SiStripCluster& stereoclust = hit.stereoCluster();
    
    StripClusterParameterEstimator::LocalValues lvMono = 
      stripCPE->localParameters( monoclust, *gdet->monoDet(), tsos);
    StripClusterParameterEstimator::LocalValues lvStereo = 
      stripCPE->localParameters( stereoclust, *gdet->stereoDet(), gluedToStereo(tsos, gdet));
      
    SiStripRecHit2D monoHit = SiStripRecHit2D( lvMono.first, lvMono.second, 
                           *gdet->monoDet(),
                           hit.monoClusterRef());
    SiStripRecHit2D stereoHit = SiStripRecHit2D( lvStereo.first, lvStereo.second,
                         *gdet->stereoDet(),
                         hit.stereoClusterRef());
    
    // return theMatcher->match(&monoHit,&stereoHit,gdet,tkDir,true);
    std::unique_ptr<SiStripMatchedRecHit2D> temp = theMatcher->match(&monoHit,&stereoHit,gdet,tkDir,false);
    if(temp.get() == nullptr) {
      temp = std::unique_ptr<SiStripMatchedRecHit2D>(hit.clone());
    }
    return temp;
}

std::unique_ptr< ProjectedSiStripRecHit2D > TkClonerImpl::operator() ( ProjectedSiStripRecHit2D const &  hit, TrajectoryStateOnSurface const &  tsos  ) const

overridevirtual

Implements TkCloner.

Definition at line 142 of file TkClonerImpl.cc.

References TrackingRecHit::det(), f, TrajectoryStateOnSurface::isValid(), TrajectoryStateOnSurface::localDirection(), StripClusterParameterEstimator::localParameters(), TrackerSingleRecHit::omniCluster(), ProjectedSiStripRecHit2D::originalDet(), GeomDet::position(), LocalError::rotate(), TrackerSingleRecHit::stripCluster(), stripCPE, GeomDet::surface(), GloballyPositioned< T >::toLocal(), PV3DBase< T, PVType, FrameType >::x(), LocalError::xx(), LocalError::xy(), PV3DBase< T, PVType, FrameType >::y(), LocalError::yy(), and PV3DBase< T, PVType, FrameType >::z().

                                                                                                                                                 {
  const SiStripCluster& clust = hit.stripCluster();
  const GeomDetUnit * gdu = reinterpret_cast<const GeomDetUnit *>(hit.originalDet());
  //if (!gdu) std::cout<<"no luck dude"<<std::endl;
  StripClusterParameterEstimator::LocalValues lv = stripCPE->localParameters(clust, *gdu, tsos);

  // project...
  const GeomDet & det = *hit.det();
  const BoundPlane& gluedPlane = det.surface();
  const BoundPlane& hitPlane = gdu->surface();
  LocalVector tkDir = (tsos.isValid() ? tsos.localDirection() : 
               det.surface().toLocal( det.position()-GlobalPoint(0,0,0)));

  auto delta = gluedPlane.localZ( hitPlane.position());
  LocalVector ldir = tkDir;
  LocalPoint lhitPos = gluedPlane.toLocal( hitPlane.toGlobal(lv.first));
  LocalPoint projectedHitPos = lhitPos - ldir * delta/ldir.z();                  

  LocalVector hitXAxis = gluedPlane.toLocal( hitPlane.toGlobal( LocalVector(1.f,0,0)));
  LocalError hitErr = lv.second;
  if (gluedPlane.normalVector().dot( hitPlane.normalVector()) < 0) {
    // the two planes are inverted, and the correlation element must change sign
    hitErr = LocalError( hitErr.xx(), -hitErr.xy(), hitErr.yy());
  }
  LocalError rotatedError = hitErr.rotate( hitXAxis.x(), hitXAxis.y());
  return std::unique_ptr<ProjectedSiStripRecHit2D>{new ProjectedSiStripRecHit2D(projectedHitPos, rotatedError, *hit.det(), *hit.originalDet(), hit.omniCluster())};
}

std::unique_ptr< ProjectedSiStripRecHit2D > TkClonerImpl::project	(	SiStripMatchedRecHit2D const &	hit,
		bool	mono,
		TrajectoryStateOnSurface const &	tsos
	)		const

Definition at line 171 of file TkClonerImpl.cc.

References TrackingRecHit::det(), f, TrajectoryStateOnSurface::isValid(), TrajectoryStateOnSurface::localDirection(), StripClusterParameterEstimator::localParameters(), SiStripMatchedRecHit2D::monoCluster(), SiStripMatchedRecHit2D::monoClusterRef(), GluedGeomDet::monoDet(), GeomDet::position(), LocalError::rotate(), SiStripMatchedRecHit2D::stereoCluster(), SiStripMatchedRecHit2D::stereoClusterRef(), GluedGeomDet::stereoDet(), stripCPE, GeomDet::surface(), GloballyPositioned< T >::toLocal(), PV3DBase< T, PVType, FrameType >::x(), LocalError::xx(), LocalError::xy(), PV3DBase< T, PVType, FrameType >::y(), LocalError::yy(), and PV3DBase< T, PVType, FrameType >::z().

                                                                                                                                                       {
  const GeomDet & det = *hit.det();
  const GluedGeomDet & gdet = static_cast<const GluedGeomDet &> (det);
  const GeomDetUnit * odet = mono ? gdet.monoDet() : gdet.stereoDet();
  const BoundPlane& gluedPlane = det.surface();
  const BoundPlane& hitPlane = odet->surface();


  LocalVector tkDir = (tsos.isValid() ? tsos.localDirection() : 
               det.surface().toLocal( det.position()-GlobalPoint(0,0,0)));

  const SiStripCluster& monoclust   = hit.monoCluster();  
  const SiStripCluster& stereoclust = hit.stereoCluster();

  StripClusterParameterEstimator::LocalValues lv;
  if (tsos.isValid()) 
    lv = mono ?
    stripCPE->localParameters( monoclust, *odet, tsos) :
    stripCPE->localParameters( stereoclust, *odet, gluedToStereo(tsos, &gdet));
  else 
    lv = stripCPE->localParameters( mono ? monoclust : stereoclust, *odet);


  auto delta = gluedPlane.localZ( hitPlane.position());
  LocalVector ldir = tkDir;
  LocalPoint lhitPos = gluedPlane.toLocal( hitPlane.toGlobal(lv.first));
  LocalPoint projectedHitPos = lhitPos - ldir * delta/ldir.z();                  

  LocalVector hitXAxis = gluedPlane.toLocal( hitPlane.toGlobal( LocalVector(1.f,0,0)));
  LocalError hitErr = lv.second;
  if (gluedPlane.normalVector().dot( hitPlane.normalVector()) < 0) {
    // the two planes are inverted, and the correlation element must change sign
    hitErr = LocalError( hitErr.xx(), -hitErr.xy(), hitErr.yy());
  }
  LocalError rotatedError = hitErr.rotate( hitXAxis.x(), hitXAxis.y());
  return std::unique_ptr<ProjectedSiStripRecHit2D>{
    new ProjectedSiStripRecHit2D(projectedHitPos, rotatedError, det, *odet, 
                 mono ? hit.monoClusterRef() : hit.stereoClusterRef() ) };
}

Member Data Documentation

const PixelClusterParameterEstimator* TkClonerImpl::pixelCPE

private

Definition at line 40 of file TkClonerImpl.h.

Referenced by makeShared(), and operator()().

const StripClusterParameterEstimator* TkClonerImpl::stripCPE

private

Definition at line 41 of file TkClonerImpl.h.

Referenced by makeShared(), operator()(), and project().

const SiStripRecHitMatcher* TkClonerImpl::theMatcher

private

Definition at line 42 of file TkClonerImpl.h.

Referenced by operator()().