SimpleDAFHitCollector Class Reference

#include <SimpleDAFHitCollector.h>

Inheritance diagram for SimpleDAFHitCollector:

Public Member Functions
const MeasurementEstimator *	getEstimator () const
const Propagator *	getPropagator () const
const SiTrackerMultiRecHitUpdator *	getUpdator () const
virtual std::vector < TrajectoryMeasurement >	recHits (const Trajectory &, const MeasurementTrackerEvent *theMTE) const
	SimpleDAFHitCollector (const MeasurementTracker *measurementTracker, const SiTrackerMultiRecHitUpdator *updator, const MeasurementEstimator *est, const Propagator *propagator, bool debug)
virtual	~SimpleDAFHitCollector ()
Public Member Functions inherited from MultiRecHitCollector
const MeasurementTracker *	getMeasurementTracker () const
	MultiRecHitCollector (const MeasurementTracker *meas)

Private Member Functions
void	buildMultiRecHits (const std::vector< TrajectoryMeasurement > &measgroup, std::vector< TrajectoryMeasurement > &result) const

Private Attributes
const bool	debug_
const MeasurementEstimator *	theEstimator
const Propagator *	thePropagator
const SiTrackerMultiRecHitUpdator *	theUpdator

Detailed Description

Definition at line 10 of file SimpleDAFHitCollector.h.

Constructor & Destructor Documentation

SimpleDAFHitCollector::SimpleDAFHitCollector ( const MeasurementTracker *  measurementTracker, const SiTrackerMultiRecHitUpdator *  updator, const MeasurementEstimator *  est, const Propagator *  propagator, bool  debug  )

inlineexplicit

Definition at line 12 of file SimpleDAFHitCollector.h.

                  :MultiRecHitCollector(measurementTracker), theUpdator(updator), theEstimator(est), thePropagator(propagator), debug_(debug){}

virtual SimpleDAFHitCollector::~SimpleDAFHitCollector ( )

inlinevirtual

Definition at line 19 of file SimpleDAFHitCollector.h.

{}

Member Function Documentation

void SimpleDAFHitCollector::buildMultiRecHits ( const std::vector< TrajectoryMeasurement > &  measgroup, std::vector< TrajectoryMeasurement > &  result  ) const

private

const MeasurementEstimator* SimpleDAFHitCollector::getEstimator ( ) const

inline

Definition at line 31 of file SimpleDAFHitCollector.h.

References theEstimator.

{return theEstimator;}

const Propagator* SimpleDAFHitCollector::getPropagator ( ) const

inline

Definition at line 32 of file SimpleDAFHitCollector.h.

References thePropagator.

{return thePropagator;}

const SiTrackerMultiRecHitUpdator* SimpleDAFHitCollector::getUpdator ( ) const

inline

Definition at line 30 of file SimpleDAFHitCollector.h.

References theUpdator.

{return theUpdator;}

vector< TrajectoryMeasurement > SimpleDAFHitCollector::recHits ( const Trajectory &  traj, const MeasurementTrackerEvent *  theMTE  ) const

virtual

Implements MultiRecHitCollector.

Definition at line 16 of file SimpleDAFHitCollector.cc.

References cond::rpcobimon::current, MeasurementDet::geomDet(), tracking::TempMeasurements::hits, i, MeasurementTrackerEvent::idToDet(), TrajectoryStateOnSurface::isValid(), LogTrace, MeasurementDetWithData::mdet(), MeasurementDetWithData::measurements(), Trajectory::measurements(), TrackingRecHit::missing, TrajectoryStateOnSurface::rescaleError(), query::result, and tracking::TempMeasurements::size().

                                                                                                                               {

  LogTrace("MultiRecHitCollector") << " Calling SimpleDAFHitCollector::recHits" << std::endl;

  //WARNING: At the moment the trajectories has the measurements 
  //with reversed sorting after the track smoothing
  const vector<TrajectoryMeasurement>& meas = traj.measurements();

  if (meas.empty()) return vector<TrajectoryMeasurement>();

  //debug
  LogTrace("MultiRecHitCollector") << "Original measurements are:";
  for(vector<TrajectoryMeasurement>::const_iterator itrajmeas = meas.begin(); itrajmeas < meas.end();
      itrajmeas++) {
    if (itrajmeas->recHit()->isValid()){
      LogTrace("MultiRecHitCollector") << "  Valid Hit with DetId " << itrajmeas->recHit()->geographicalId().rawId()
                    << " local position " << itrajmeas->recHit()->hit()->localPosition()
                    << " global position " << itrajmeas->recHit()->hit()->globalPosition() ;
            } else {
              LogTrace("MultiRecHitCollector") << "   Invalid Hit with DetId " << itrajmeas->recHit()->geographicalId().rawId();
            }
  }

  //groups the TrajectoryMeasurements on a sensor by sensor
  //we have to sort the TrajectoryMeasurements in the opposite way in the fitting direction
  vector<TrajectoryMeasurement> result;
  for(vector<TrajectoryMeasurement>::const_reverse_iterator itrajmeas = meas.rbegin(); itrajmeas < meas.rend();
      itrajmeas++) {

      DetId id = itrajmeas->recHit()->geographicalId();
      MeasurementDetWithData measDet = theMTE->idToDet(id);
      tracking::TempMeasurements tmps;
      std::vector<const TrackingRecHit*> hits;

      TrajectoryStateOnSurface current = itrajmeas->updatedState();
      //the error is scaled in order to take more "compatible" hits
      if (current.isValid()) current.rescaleError(10);

      TrajectoryStateOnSurface state = itrajmeas->predictedState();
      if (!state.isValid()){
        LogTrace("MultiRecHitCollector") << "Something wrong! no valid TSOS found in current group ";
        continue;
      }
      //collected hits compatible with the itrajmeas
      if( measDet.measurements(current, *(getEstimator()), tmps)){
        LogTrace("MultiRecHitCollector") << "  Found " << tmps.size() << " compatible measurements";
        for (std::size_t i=0; i!=tmps.size(); ++i){
          DetId idtemps = tmps.hits[i]->geographicalId();

          if( idtemps == id && tmps.hits[i]->hit()->isValid() ) {
            hits.push_back(tmps.hits[i]->hit());
          }
        }

        //I will keep the Invalid hit, IF this is not the first one       
        if (hits.empty()){
          LogTrace("MultiRecHitCollector") << " -> but no valid hits found in current group.";

          if( result.empty() ) continue;

          result.push_back(TrajectoryMeasurement(state,
                                        std::make_shared<InvalidTrackingRecHit>(measDet.mdet().geomDet(), TrackingRecHit::missing)));
        } else {
          //measurements in groups are sorted with increating chi2
          //sort( *hits.begin(), *hits.end(), TrajMeasLessEstim());

          LogTrace("MultiRecHitCollector") << " -> " << hits.size() << " valid hits for this sensor.";

          //building a MultiRecHit out of each sensor group
          result.push_back(TrajectoryMeasurement(state,theUpdator->buildMultiRecHit(hits, state)));
        }
      } else {
          LogTrace("MultiRecHitCollector") << "  No measurements found in current group.";

          if( result.empty() ) continue;

          result.push_back(TrajectoryMeasurement(state,
                                        std::make_shared<InvalidTrackingRecHit>(measDet.mdet().geomDet(), TrackingRecHit::missing)));

      }
      
  }
  LogTrace("MultiRecHitCollector") << " Ending SimpleDAFHitCollector::recHits >> " << result.size();

  //adding a protection against too few hits and invalid hits 
  //(due to failed propagation on the same surface of the original hits)
  if (result.size()>2)
  {
    int hitcounter=0;
    //check if the vector result has more than 3 valid hits
    for (vector<TrajectoryMeasurement>::const_iterator iimeas = result.begin(); iimeas != result.end(); ++iimeas) {
      if(iimeas->recHit()->isValid()) hitcounter++;
    }

    if(hitcounter>2)
      return result;
    else return vector<TrajectoryMeasurement>();
  }

  else{return vector<TrajectoryMeasurement>();}

}

Member Data Documentation

const bool SimpleDAFHitCollector::debug_

private

Definition at line 43 of file SimpleDAFHitCollector.h.

const MeasurementEstimator* SimpleDAFHitCollector::theEstimator

private

Definition at line 40 of file SimpleDAFHitCollector.h.

Referenced by getEstimator().

const Propagator* SimpleDAFHitCollector::thePropagator

private

Definition at line 42 of file SimpleDAFHitCollector.h.

Referenced by getPropagator().

const SiTrackerMultiRecHitUpdator* SimpleDAFHitCollector::theUpdator

private

Definition at line 39 of file SimpleDAFHitCollector.h.

Referenced by getUpdator().