SimpleDAFHitCollector Class Reference

#include <SimpleDAFHitCollector.h>

Inheritance diagram for SimpleDAFHitCollector:

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

Private Member Functions
std::unique_ptr< TrackingRecHit >	clone (BaseTrackerRecHit const &hit2D) const
std::unique_ptr< TrackingRecHit >	rightdimension (TrackingRecHit const &hit) const

Private Attributes
const bool	debug_
const MeasurementEstimator *	theEstimator
TkClonerImpl	theHitCloner
const Propagator *	thePropagator
const TrackerTopology *	theTopology
const SiTrackerMultiRecHitUpdator *	theUpdator

Detailed Description

Definition at line 17 of file SimpleDAFHitCollector.h.

Constructor & Destructor Documentation

SimpleDAFHitCollector()

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

inlineexplicit

Definition at line 19 of file SimpleDAFHitCollector.h.

References runTheMatrix::const, SiTrackerMultiRecHitUpdator::getBuilder(), theHitCloner, and theUpdator.

      : MultiRecHitCollector(measurementTracker),
        theTopology(trackerTopology),
        theUpdator(updator),
        theEstimator(est),
        thePropagator(propagator),
        debug_(debug) {
    theHitCloner = static_cast<TkTransientTrackingRecHitBuilder const*>(theUpdator->getBuilder())->cloner();
  }

~SimpleDAFHitCollector()

SimpleDAFHitCollector::~SimpleDAFHitCollector ( )

inlineoverride

Definition at line 34 of file SimpleDAFHitCollector.h.

{}

Member Function Documentation

clone()

std::unique_ptr<TrackingRecHit> SimpleDAFHitCollector::clone ( BaseTrackerRecHit const &  hit2D ) const

inlineprivate

Definition at line 75 of file SimpleDAFHitCollector.h.

References TrackingRecHit::clone(), TrackingRecHit::det(), TrackingRecHit::detUnit(), makeMuonMisalignmentScenario::endcap, BaseTrackerRecHit::firstClusterRef(), BaseTrackerRecHit::localPosition(), BaseTrackerRecHit::localPositionError(), SiStripPI::max, and LocalError::xx().

Referenced by rightdimension().

                                                                            {
    auto const& detU = *hit2D.detUnit();
    //Use 2D SiStripRecHit in endcap
    bool endcap = detU.type().isEndcap();
    if (endcap)
      return std::unique_ptr<TrackingRecHit>{hit2D.clone()};
    return std::unique_ptr<TrackingRecHit>{
        new SiStripRecHit1D(hit2D.localPosition(),
                            LocalError(hit2D.localPositionError().xx(), 0.f, std::numeric_limits<float>::max()),
                            *hit2D.det(),
                            hit2D.firstClusterRef())};
  }

Debug()

void SimpleDAFHitCollector::Debug

(

const std::vector< TrajectoryMeasurement >

TM

)

const

Definition at line 156 of file SimpleDAFHitCollector.cc.

References MuonSubdetId::CSC, DetId::det(), MuonSubdetId::DT, LogTrace, DetId::Muon, PixelSubdetector::PixelBarrel, PixelSubdetector::PixelEndcap, MuonSubdetId::RPC, DetId::subdetId(), StripSubdetector::TEC, StripSubdetector::TIB, StripSubdetector::TID, StripSubdetector::TOB, and DetId::Tracker.

                                                                                 {
#ifdef EDM_ML_DEBUG
  for (vector<TrajectoryMeasurement>::const_iterator itrajmeas = TM.begin(); itrajmeas < TM.end(); itrajmeas++) {
    if (itrajmeas->recHit()->isValid()) {
      LogTrace("MultiRecHitCollector") << "  Valid Hit with DetId " << itrajmeas->recHit()->geographicalId().rawId()
                                       << " and dim:"
                                       << itrajmeas->recHit()->dimension()
                                       //<< " type " << typeid(itrajmeas->recHit()).name()
                                       << " local position " << itrajmeas->recHit()->hit()->localPosition()
                                       << " global position " << itrajmeas->recHit()->hit()->globalPosition()
                                       << " and r " << itrajmeas->recHit()->hit()->globalPosition().perp();

      DetId hitId = itrajmeas->recHit()->geographicalId();

      if (hitId.det() == DetId::Tracker) {
        switch (hitId.subdetId()) {
          case StripSubdetector::TIB:
            LogTrace("MultiRecHitCollector") << "  I am TIB " << theTopology->tibLayer(hitId);
            break;
          case StripSubdetector::TOB:
            LogTrace("MultiRecHitCollector") << "  I am TOB " << theTopology->tobLayer(hitId);
            break;
          case StripSubdetector::TEC:
            LogTrace("MultiRecHitCollector") << "  I am TEC " << theTopology->tecWheel(hitId);
            break;
          case StripSubdetector::TID:
            LogTrace("MultiRecHitCollector") << "  I am TID " << theTopology->tidWheel(hitId);
            break;
          case PixelSubdetector::PixelBarrel:
            LogTrace("MultiRecHitCollector") << "  I am PixBar " << theTopology->pxbLayer(hitId);
            break;
          case PixelSubdetector::PixelEndcap:
            LogTrace("MultiRecHitCollector") << "  I am PixFwd " << theTopology->pxfDisk(hitId);
            break;
          default:
            LogTrace("MultiRecHitCollector") << "  UNKNOWN TRACKER HIT TYPE ";
            break;
        }
      } else if (hitId.det() == DetId::Muon) {
        if (hitId.subdetId() == MuonSubdetId::DT)
          LogTrace("MultiRecHitCollector") << "  I am DT " << DTWireId(hitId);
        else if (hitId.subdetId() == MuonSubdetId::CSC)
          LogTrace("MultiRecHitCollector") << "  I am CSC " << CSCDetId(hitId);
        else if (hitId.subdetId() == MuonSubdetId::RPC)
          LogTrace("MultiRecHitCollector") << "  I am RPC " << RPCDetId(hitId);
        else
          LogTrace("MultiRecHitCollector") << "  UNKNOWN MUON HIT TYPE ";
      } else
        LogTrace("MultiRecHitCollector") << "  UNKNOWN HIT TYPE ";

      LogTrace("MultiRecHitCollector") << "  TSOS predicted_fwd " << itrajmeas->predictedState().localPosition();
      LogTrace("MultiRecHitCollector") << "  TSOS predicted_bwd "
                                       << itrajmeas->backwardPredictedState().localPosition();
      LogTrace("MultiRecHitCollector") << "  TSOS smoothtsos " << itrajmeas->updatedState().localPosition();
    } else {
      LogTrace("MultiRecHitCollector") << "  Invalid Hit with DetId " << itrajmeas->recHit()->geographicalId().rawId();
    }
    LogTrace("MultiRecHitCollector") << "\n";
  }
#endif
}

getEstimator()

const MeasurementEstimator* SimpleDAFHitCollector::getEstimator ( ) const

inline

Definition at line 46 of file SimpleDAFHitCollector.h.

References theEstimator.

{ return theEstimator; }

getPropagator()

const Propagator* SimpleDAFHitCollector::getPropagator ( ) const

inline

Definition at line 47 of file SimpleDAFHitCollector.h.

References thePropagator.

{ return thePropagator; }

getUpdator()

const SiTrackerMultiRecHitUpdator* SimpleDAFHitCollector::getUpdator ( ) const

inline

Definition at line 45 of file SimpleDAFHitCollector.h.

References theUpdator.

{ return theUpdator; }

recHits()

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

overridevirtual

Implements MultiRecHitCollector.

Definition at line 28 of file SimpleDAFHitCollector.cc.

References heavyFlavorValidationHarvestingSequence_cff::combiner, Debug, MeasurementDet::geomDet(), hfClusterShapes_cfi::hits, tracking::TempMeasurements::hits, mps_fire::i, TrajectoryStateOnSurface::isValid(), LogTrace, MeasurementDetWithData::mdet(), MeasurementDetWithData::measurements(), Trajectory::measurements(), TrackingRecHit::missing, eostools::move(), mps_fire::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();
  unsigned int hitcounter = 1;

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

  LogTrace("MultiRecHitCollector") << "  Original measurements are:";
  Debug(meas);

  //groups hits on a sensor by sensor with same Id of previous TM
  //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++, hitcounter++) {
    DetId id = itrajmeas->recHit()->geographicalId();
    MeasurementDetWithData measDet = theMTE->idToDet(id);
    tracking::TempMeasurements tmps;

    std::vector<const TrackingRecHit*> hits;
    std::vector<std::unique_ptr<const TrackingRecHit>> hitsOwner;

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

    TrajectoryStateOnSurface predtsos_fwd = itrajmeas->predictedState();
    TrajectoryStateOnSurface predtsos_bwd = itrajmeas->backwardPredictedState();
    if (!predtsos_fwd.isValid() || !predtsos_bwd.isValid()) {
      LogTrace("MultiRecHitCollector") << "Something wrong! no valid TSOS found in current group ";
      continue;
    }

    TrajectoryStateCombiner combiner;
    TrajectoryStateOnSurface combtsos;
    if (hitcounter == meas.size())
      combtsos = predtsos_fwd;
    else if (hitcounter == 1)
      combtsos = predtsos_bwd;
    else
      combtsos = combiner(predtsos_bwd, predtsos_fwd);

    //collected hits compatible with the itrajmeas
    if (measDet.measurements(smoothtsos, *(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()) {
          LogTrace("MultiRecHitCollector") << "  This is valid with position " << tmps.hits[i]->hit()->localPosition()
                                           << " and error " << tmps.hits[i]->hit()->localPositionError();

          TransientTrackingRecHit::RecHitPointer transient = theUpdator->getBuilder()->build(tmps.hits[i]->hit());
          TrackingRecHit::ConstRecHitPointer preciseHit = theHitCloner.makeShared(transient, combtsos);
          auto righthit = rightdimension(*preciseHit);
          hitsOwner.push_back(std::move(righthit));
          hits.push_back(hitsOwner.back().get());
        }
      }

      //the error was scaled, now is scaled back (even if this is not tightly necessary)
      //        if (smoothtsos.isValid()) smoothtsos.rescaleError(0.1);

      //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(
            predtsos_fwd, std::make_shared<InvalidTrackingRecHit>(measDet.mdet().geomDet(), TrackingRecHit::missing)));
      } else {
        //measurements in groups are sorted with increating chi2
        //sort( *hits.begin(), *hits.end(), TrajMeasLessEstim());
        if (!itrajmeas->recHit()->isValid())
          LogTrace("MultiRecHitCollector")
              << "  -> " << hits.size() << " valid hits for this sensor. (IT WAS INVALID!!!)";
        else
          LogTrace("MultiRecHitCollector") << "  -> " << hits.size() << " valid hits for this sensor.";

        //building a MultiRecHit out of each sensor group
        result.push_back(TrajectoryMeasurement(predtsos_fwd, theUpdator->buildMultiRecHit(hits, combtsos, measDet)));
      }
    } else {
      LogTrace("MultiRecHitCollector") << "  No measurements found in current group.";
      //the error was scaled, now is scaled back (even if this is not tightly necessary)
      //          if (smoothtsos.isValid()) smoothtsos.rescaleError(0.1);

      if (result.empty())
        continue;

      result.push_back(TrajectoryMeasurement(
          predtsos_fwd, std::make_shared<InvalidTrackingRecHit>(measDet.mdet().geomDet(), TrackingRecHit::missing)));
    }
  }
  LogTrace("MultiRecHitCollector") << " Ending SimpleDAFHitCollector::recHits >> " << result.size();

  //LogTrace("MultiRecHitCollector") << "  New measurements are:";
  //Debug(result);

  //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>();
  }
}

rightdimension()

std::unique_ptr<TrackingRecHit> SimpleDAFHitCollector::rightdimension ( TrackingRecHit const &  hit ) const

inlineprivate

Definition at line 55 of file SimpleDAFHitCollector.h.

References clone(), runTheMatrix::const, and LogDebug.

                                                                                {
    if (!hit.isValid() || (hit.dimension() != 2)) {
      return std::unique_ptr<TrackingRecHit>{hit.clone()};
    }
    auto const& thit = static_cast<BaseTrackerRecHit const&>(hit);
    auto const& clus = thit.firstClusterRef();
    if (clus.isPixel())
      return std::unique_ptr<TrackingRecHit>{hit.clone()};
    else if (clus.isPhase2())
      return std::unique_ptr<TrackingRecHit>{hit.clone()};
    else if (thit.isMatched()) {
      LogDebug("MultiRecHitCollector") << " SiStripMatchedRecHit2D to check!!!";
      return std::unique_ptr<TrackingRecHit>{hit.clone()};
    } else if (thit.isProjected()) {
      edm::LogError("MultiRecHitCollector") << " ProjectedSiStripRecHit2D should not be present at this stage!!!";
      return std::unique_ptr<TrackingRecHit>{hit.clone()};
    } else
      return clone(thit);
  }

Member Data Documentation

debug_

const bool SimpleDAFHitCollector::debug_

private

Definition at line 95 of file SimpleDAFHitCollector.h.

theEstimator

const MeasurementEstimator* SimpleDAFHitCollector::theEstimator

private

Definition at line 91 of file SimpleDAFHitCollector.h.

Referenced by getEstimator().

theHitCloner

TkClonerImpl SimpleDAFHitCollector::theHitCloner

private

Definition at line 94 of file SimpleDAFHitCollector.h.

Referenced by SimpleDAFHitCollector().

thePropagator

const Propagator* SimpleDAFHitCollector::thePropagator

private

Definition at line 93 of file SimpleDAFHitCollector.h.

Referenced by getPropagator().

theTopology

const TrackerTopology* SimpleDAFHitCollector::theTopology

private

Definition at line 89 of file SimpleDAFHitCollector.h.

theUpdator

const SiTrackerMultiRecHitUpdator* SimpleDAFHitCollector::theUpdator

private

Definition at line 90 of file SimpleDAFHitCollector.h.

Referenced by getUpdator(), and SimpleDAFHitCollector().