TrajectoryToResiduals.cc File Reference

#include "TrajectoryToResiduals.h"
#include "TrackingTools/PatternTools/interface/Trajectory.h"
#include "TrackingTools/TrackFitters/interface/TrajectoryStateCombiner.h"
#include "TrackingTools/TransientTrackingRecHit/interface/HelpertRecHit2DLocalPos.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"

Go to the source code of this file.

Functions
reco::TrackResiduals	trajectoryToResiduals (const Trajectory &trajectory, enum reco::TrackResiduals::ResidualType type)

Function Documentation

reco::TrackResiduals trajectoryToResiduals	(	const Trajectory &	trajectory,
		enum reco::TrackResiduals::ResidualType	type
	)

Definition at line 8 of file TrajectoryToResiduals.cc.

References alongMomentum, TrajectoryMeasurement::backwardPredictedState(), funct::combine(), Trajectory::direction(), TrajectoryMeasurement::forwardPredictedState(), i, TrajectoryStateOnSurface::isValid(), TrajectoryStateOnSurface::localError(), TrajectoryStateOnSurface::localPosition(), Trajectory::measurements(), HelpertRecHit2DLocalPos::parError(), LocalTrajectoryError::positionError(), TrajectoryMeasurement::recHit(), reco::TrackResiduals::setPullXY(), reco::TrackResiduals::setResidualXY(), mathSSE::sqrt(), PV3DBase< T, PVType, FrameType >::x(), vdt::x, reco::TrackResiduals::X_Y_PULLS, reco::TrackResiduals::X_Y_RESIDUALS, LocalError::xx(), detailsBasic3DVector::y, PV3DBase< T, PVType, FrameType >::y(), and LocalError::yy().

Referenced by KfTrackProducerBase::putInEvt().

{
     reco::TrackResiduals residuals(type);
     int i_residual = 0;
     Trajectory::DataContainer::const_iterator i_fwd = 
      trajectory.measurements().begin(); 
     Trajectory::DataContainer::const_reverse_iterator i_bwd = 
      trajectory.measurements().rbegin(); 
     Trajectory::DataContainer::const_iterator i_end = 
      trajectory.measurements().end(); 
     Trajectory::DataContainer::const_reverse_iterator i_rend = 
      trajectory.measurements().rend(); 
     bool forward = trajectory.direction() == alongMomentum;
     for (; forward ? i_fwd != i_end : i_bwd != i_rend; 
      ++i_fwd, ++i_bwd, ++i_residual) {
      const TrajectoryMeasurement *i = forward ? &*i_fwd : &*i_bwd;
      if (!i->recHit()->isValid()||i->recHit()->det()==0) 
           continue;
      TrajectoryStateCombiner combine;
      if (!i->forwardPredictedState().isValid() || !i->backwardPredictedState().isValid())
        {
          edm::LogError("InvalideState")<<"one of the step is invalid";
          continue;
        }

      TrajectoryStateOnSurface combo = combine(i->forwardPredictedState(),
                           i->backwardPredictedState());
      
      if (!combo.isValid()){
        edm::LogError("InvalideState")<<"the combined state is invalid";
        continue;
      }

      LocalPoint combo_localpos = combo.localPosition();
      LocalError combo_localerr = combo.localError().positionError();
      LocalPoint dethit_localpos = i->recHit()->localPosition();     
      LocalError dethit_localerr = i->recHit()->localPositionError();
      HelpertRecHit2DLocalPos helper;
      AlgebraicSymMatrix error_including_alignment = 
           helper.parError(dethit_localerr, *i->recHit()->det());
      switch (type) {
      case reco::TrackResiduals::X_Y_RESIDUALS: 
      {
           double x = (dethit_localpos.x() - combo_localpos.x()) / 
            sqrt(error_including_alignment[0][0]);
           double y = (dethit_localpos.y() - combo_localpos.y()) / 
            sqrt(error_including_alignment[1][1]);
           residuals.setResidualXY(i_residual, x, y);
           break;
      }
      case reco::TrackResiduals::X_Y_PULLS:
      {
           double x = (dethit_localpos.x() - combo_localpos.x()) / 
            sqrt(error_including_alignment[0][0] + combo_localerr.xx());
           double y = (dethit_localpos.y() - combo_localpos.y()) / 
            sqrt(error_including_alignment[1][1] + combo_localerr.yy());
           residuals.setPullXY(i_residual, x, y);
           break;
      }
      default:
           assert(0);
      }
     }
     return residuals;
}