24 std::stringstream buffer;
25 buffer<<v.size()<<
" total measurements\n";
26 std::vector<TrajectoryMeasurement>::const_iterator it=v.begin();
27 for (; it!=v.end();++it){
34 std::stringstream buffer;
35 buffer<<v.
size()<<
" total measurements\n";
37 for (; it!=v.
rend();--it){
44 std::stringstream buffer;
46 <<
"layer pointer: "<<tm.
layer()<<
"\n" 49 buffer<<
"updated state: \n" 53 buffer<<
"forward predicted state: \n" 57 buffer<<
"predicted state: \n" 60 else buffer<<
"no valid state\n";
63 <<
"detId: "<<tm.
recHit()->geographicalId().rawId();
64 if (tm.
recHit()->isValid()){
65 buffer<<
"\n hit global x: "<<tm.
recHit()->globalPosition()
66 <<
"\n hit global error: "<<tm.
recHit()->globalPositionError().matrix()
67 <<
"\n hit local x:"<<tm.
recHit()->localPosition()
68 <<
"\n hit local error"<<tm.
recHit()->localPositionError();
69 }
else buffer<<
"\n (-,-,-)";
79 std::vector<Trajectory> & unsmoothedResult
81 std::stringstream buffer;
96 buffer<<
"number of finalTrajectories: " << unsmoothedResult.size()<<std::endl;
97 for (std::vector<Trajectory>::const_iterator it = unsmoothedResult.begin();
98 it != unsmoothedResult.end(); it++) {
99 if (it->lastMeasurement().updatedState().isValid()) {
100 buffer<<
"candidate's n valid and invalid hit, chi2, pt, eta : " 101 << it->foundHits() <<
" , " 102 << it->lostHits() <<
" , " 103 << it->chiSquared() <<
" , " 104 << it->lastMeasurement().updatedState().globalMomentum().perp() <<
" , " 105 << it->lastMeasurement().updatedState().globalMomentum().eta()<<std::endl;
106 }
else if (it->lastMeasurement().predictedState().isValid()) {
107 buffer<<
"candidate's n valid and invalid hit, chi2, pt, eta : " 108 << it->foundHits() <<
" , " 109 << it->lostHits() <<
" , " 110 << it->chiSquared() <<
" , " 111 << it->lastMeasurement().predictedState().globalMomentum().perp() <<
" , " 112 << it->lastMeasurement().predictedState().globalMomentum().eta()<<std::endl;
113 }
else buffer<<
"candidate with invalid last measurement state!"<<std::endl;
115 buffer<<
"=================================================";
const_iterator rend() const
TrajectoryStateOnSurface const & predictedState() const
ConstRecHitPointer const & recHit() const
GlobalPoint globalPosition() const
static std::string dumpMeasurement(const TrajectoryMeasurement &tm)
static std::string dumpMeasurements(const std::vector< TrajectoryMeasurement > &v)
static std::string regressionTest(const TrackerGeometry &tracker, std::vector< Trajectory > &unsmoothedResult)
const DetLayer * layer() const
TrajectoryStateOnSurface const & forwardPredictedState() const
Access to forward predicted state (from fitter or builder)
const_iterator rbegin() const
GlobalVector globalMomentum() const
TrajectoryStateOnSurface const & updatedState() const
TrajectoryStateOnSurface const & backwardPredictedState() const
Access to backward predicted state (from smoother)