PrintoutHelper.h

Go to the documentation of this file.

#ifndef RecoTracker_CkfPattern_PrintoutHelper_h
#define RecoTracker_CkfPattern_PrintoutHelper_h

#include "FWCore/MessageLogger/interface/MessageLogger.h"

#include "TrackingTools/PatternTools/interface/Trajectory.h"
#include "TrackingTools/PatternTools/interface/TrajectoryMeasurement.h"
#include "TrackingTools/PatternTools/interface/bqueue.h"

class TrackerGeometry;

class PrintoutHelper{
 public:
  template< class collection > static std::string dumpCandidates( collection & candidates);
  template< class Candidate> static std::string dumpCandidate( const Candidate & candidate,bool showErrors=false);
  static std::string dumpMeasurements(const std::vector<TrajectoryMeasurement> & v) ;
  static std::string dumpMeasurements(const cmsutils::bqueue<TrajectoryMeasurement> & v);
  static std::string dumpMeasurement(const TrajectoryMeasurement & tm);
  static std::string regressionTest(const TrackerGeometry & tracker,std::vector<Trajectory> & unsmoothedResult);
};

template<class Candidate>
std::string PrintoutHelper::dumpCandidate( const Candidate & traj,bool showErrors ){

  // does not work....
  LogDebug("PrintoutHelperError")<<"switching on error printout"<<(showErrors=true);

  std::stringstream buffer;
  if (!traj.measurements().empty()){
    const TrajectoryMeasurement & last = traj.lastMeasurement();
    
    buffer<<"with: "<<traj.measurements().size()<<" measurements."<< traj.lostHits() << " lost, " << traj.foundHits()<<" found, "
          << traj.trailingFoundHits() << " trailing, " << traj.cccBadHits() << " badCC, "
          << "chi2="<<traj.chiSquared() << ' '
          << int(traj.nLoops())<<" loops\n";
    if (last.updatedState().isValid()) {
      const TrajectoryStateOnSurface & tsos = last.updatedState();
      if (showErrors)
    buffer <<"Last [Updated] state\n : "<<tsos<<"\n";
      else
    buffer <<"Last [Updated] state\n x: "<<tsos.globalPosition()<<"\n p: "<<tsos.globalMomentum()<<"\n";
    } else if(last.forwardPredictedState().isValid()){
      const TrajectoryStateOnSurface & tsos = last.forwardPredictedState();
      if (showErrors)
    buffer <<"Last [fwdPredicted] state\n : "<<tsos<<"\n";
      else
    buffer <<"Last [fwdPredicted] state\n x: "<<tsos.globalPosition()<<"\n p: "<<tsos.globalMomentum()<<"\n";
    } else if (last.predictedState().isValid()){
      const TrajectoryStateOnSurface & tsos = last.predictedState();
      if (showErrors)
    buffer <<"Last [Predicted] state\n : "<<tsos<<"\n";
      else
    buffer <<"Last [Predicted] state\n x: "<<tsos.globalPosition()<<"\n p: "<<tsos.globalMomentum()<<"\n";
    }
    buffer <<" hit is: "<<(last.recHit()->isValid()?"valid":"invalid")<<"\n";
    if (last.recHit()->isValid()) {
      buffer <<"on detId: "<<last.recHit()->geographicalId().rawId()<<"\n";  
      buffer <<"gp: "<<last.recHit()->globalPosition()<<"\n";
    }
  }
  else{
      buffer<<" no measurement. \n";}
  return buffer.str();
}


template< class collection > 
std::string PrintoutHelper::dumpCandidates( collection & candidates) {
  std::stringstream buffer;
  buffer    << "\n____________________________\n";
  unsigned int ic=0;
  for (auto const & traj : candidates) {  
    buffer<<ic++<<"] " << (traj.isValid() ? "valid " : "invalid ");
    buffer<<PrintoutHelper::dumpCandidate(traj);
  }
 buffer << "\n____________________________\n";
  return buffer.str();
}


#endif