TSiStripMatchedRecHit.cc

Go to the documentation of this file.

#include "RecoTracker/TransientTrackingRecHit/interface/TSiStripMatchedRecHit.h"
#include "DataFormats/TrackerRecHit2D/interface/SiStripMatchedRecHit2D.h"
#include "DataFormats/TrackerRecHit2D/interface/SiStripRecHit2D.h"
#include "RecoLocalTracker/SiStripRecHitConverter/interface/SiStripRecHitMatcher.h"
#include "Geometry/TrackerGeometryBuilder/interface/GluedGeomDet.h"
//#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "RecoTracker/TransientTrackingRecHit/interface/TSiStripRecHit2DLocalPos.h"

#undef RecoTracker_TransientTrackingRecHit_TSiStripMatchedRecHit_RefitProj
#undef RecoTracker_TransientTrackingRecHit_TSiStripMatchedRecHit_RefitLGL
#ifdef RecoTracker_TransientTrackingRecHit_TSiStripMatchedRecHit_RefitLGL 
// Local lo Global lo Local
inline LocalTrajectoryParameters gluedToStereo(const TrajectoryStateOnSurface &tsos, const GluedGeomDet *gdet) {    
    const BoundPlane &stripPlane = gdet->stereoDet()->surface();
    LocalPoint  lp = stripPlane.toLocal(tsos.globalPosition());
    LocalVector ld = stripPlane.toLocal(tsos.globalParameters().momentum());
    return LocalTrajectoryParameters(lp,ld,tsos.charge());
}
#elif defined(RecoTracker_TransientTrackingRecHit_TSiStripMatchedRecHit_RefitProj)
// A la RecHitProjector
inline LocalTrajectoryParameters gluedToStereo(const TrajectoryStateOnSurface &tsos, const GluedGeomDet *gdet) {
    const BoundPlane &stripPlane = gdet->stereoDet()->surface();
    double delta = stripPlane.localZ( tsos.globalPosition());
    LocalVector ld = stripPlane.toLocal(tsos.globalParameters().momentum());
    LocalPoint  lp = stripPlane.toLocal(tsos.globalPosition()) - ld*delta/ld.z();
    return LocalTrajectoryParameters(lp,ld,tsos.charge());
}
#else
// Dummy
inline const LocalTrajectoryParameters & gluedToStereo(const TrajectoryStateOnSurface &tsos, const GluedGeomDet *gdet) {
    return tsos.localParameters();
}
#endif

TSiStripMatchedRecHit::RecHitPointer 
TSiStripMatchedRecHit::clone( const TrajectoryStateOnSurface& ts) const
{
  if (theMatcher != 0) {
    const SiStripMatchedRecHit2D *orig = static_cast<const SiStripMatchedRecHit2D *> (this->hit());
    const GeomDet *det = this->det();
    const GluedGeomDet *gdet = static_cast<const GluedGeomDet *> (det);
    //if ((orig == 0) || (gdet == 0)) return this->clone(); // or just die ?
    LocalVector tkDir = (ts.isValid() ? ts.localDirection() : 
             det->surface().toLocal( det->position()-GlobalPoint(0,0,0)));
    
    if(theCPE != 0){    
      //approximation: the ts parameter on the glued surface are used on the mono
      // and stereo surface to re-evaluate cluster parameter. A further propagation 
      //is slow// and useless (?) in this case.

      const SiStripMatchedRecHit2D* better;

      if(!orig->monoHit()->cluster().isNull()){
    const SiStripCluster& monoclust   = *orig->monoHit()->cluster();  
    const SiStripCluster& stereoclust = *orig->stereoHit()->cluster();

    StripClusterParameterEstimator::LocalValues lvMono = 
      theCPE->localParameters( monoclust, *gdet->monoDet(), ts);
    StripClusterParameterEstimator::LocalValues lvStereo = 
      theCPE->localParameters( stereoclust, *gdet->stereoDet(), gluedToStereo(ts, gdet));
    
    SiStripRecHit2D monoHit = SiStripRecHit2D( lvMono.first, lvMono.second,
                   gdet->monoDet()->geographicalId(),
                   orig->monoHit()->cluster());
    
    SiStripRecHit2D stereoHit = SiStripRecHit2D( lvStereo.first, lvStereo.second,
                     gdet->stereoDet()->geographicalId(),
                     orig->stereoHit()->cluster());
    better =  theMatcher->match(&monoHit,&stereoHit,gdet,tkDir);
      }else{
        const SiStripCluster& monoclust   = *orig->monoHit()->cluster_regional();  
    const SiStripCluster& stereoclust = *orig->stereoHit()->cluster_regional();
    StripClusterParameterEstimator::LocalValues lvMono = 
      theCPE->localParameters( monoclust, *gdet->monoDet(), ts);
    StripClusterParameterEstimator::LocalValues lvStereo = 
      theCPE->localParameters( stereoclust, *gdet->stereoDet(), gluedToStereo(ts, gdet));
    
    SiStripRecHit2D monoHit = SiStripRecHit2D( lvMono.first, lvMono.second,
                           gdet->monoDet()->geographicalId(),
                           orig->monoHit()->cluster_regional());
    
    SiStripRecHit2D stereoHit = SiStripRecHit2D( lvStereo.first, lvStereo.second,
                             gdet->stereoDet()->geographicalId(),
                             orig->stereoHit()->cluster_regional());
    better =  theMatcher->match(&monoHit,&stereoHit,gdet,tkDir);
      }
      
      if (better == 0) {
    //dm::LogWarning("TSiStripMatchedRecHit") << "Refitting of a matched rechit returns NULL";
    return this->clone();
      }

      RecHitPointer result = TSiStripMatchedRecHit::build( gdet, better, theMatcher,theCPE, weight(), getAnnealingFactor()  );
      delete better; //the ownership of the object is passed to the caller of the matcher
      return result;

    }else{
      const SiStripMatchedRecHit2D *better = theMatcher->match(orig,gdet,tkDir);
      if (better == 0) {
    //edm::LogWarning("TSiStripMatchedRecHit") << "Refitting of a matched rechit returns NULL";
    return this->clone();        
      }
      RecHitPointer result = TSiStripMatchedRecHit::build( gdet, better, theMatcher,theCPE, weight(), getAnnealingFactor()  );
      delete better; //the ownership of the object is passed to the caller of the matcher
      return result;
    }
  }
  return this->clone();
   
}



TransientTrackingRecHit::ConstRecHitContainer   
TSiStripMatchedRecHit::transientHits () const {
  ConstRecHitContainer result;

  const GluedGeomDet *gdet = static_cast<const GluedGeomDet *> (this->det());
  const SiStripMatchedRecHit2D *orig = static_cast<const SiStripMatchedRecHit2D *> (this->hit());

  result.push_back(TSiStripRecHit2DLocalPos::build( gdet->monoDet(),orig->monoHit(),theCPE));
  result.push_back(TSiStripRecHit2DLocalPos::build( gdet->stereoDet(),orig->stereoHit(),theCPE));
  return result;
}