/data/refman/pasoursint/CMSSW_6_1_2_SLHC2/src/RecoTracker/TrackProducer/src/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(), 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(), x, reco::TrackResiduals::X_Y_PULLS, reco::TrackResiduals::X_Y_RESIDUALS, LocalError::xx(), PV3DBase< T, PVType, FrameType >::y(), detailsBasic3DVector::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;
}