MeasurementByLayerGrouper Class Reference

#include <MeasurementByLayerGrouper.h>

Public Member Functions
	MeasurementByLayerGrouper (const GeometricSearchTracker *search=nullptr)
std::vector< std::pair< const DetLayer *, std::vector< TM > > >	operator() (const std::vector< TM > &) const

Private Types
typedef TrajectoryMeasurement	TM

Private Member Functions
const DetLayer *	getDetLayer (const TM &tm) const

Private Attributes
const GeometricSearchTracker *	theGeomSearch

Detailed Description

Definition at line 13 of file MeasurementByLayerGrouper.h.

Member Typedef Documentation

TM

typedef TrajectoryMeasurement MeasurementByLayerGrouper::TM

private

Definition at line 15 of file MeasurementByLayerGrouper.h.

Constructor & Destructor Documentation

MeasurementByLayerGrouper()

MeasurementByLayerGrouper::MeasurementByLayerGrouper ( const GeometricSearchTracker *  search = nullptr )

inlineexplicit

Definition at line 21 of file MeasurementByLayerGrouper.h.

: theGeomSearch(search) {}

Member Function Documentation

getDetLayer()

const DetLayer * MeasurementByLayerGrouper::getDetLayer ( const TM &  tm ) const

private

Definition at line 49 of file MeasurementByLayerGrouper.cc.

References Exception, TrajectoryMeasurement::layer(), LogDebug, and TrajectoryMeasurement::recHit().

                                                                         {
  // if the DetLayer is set in the TM...
  if (tm.layer())
    return tm.layer();

  //if it corresponds to an invalid hit with no geomdet associated
  //we can't retrieve the  DetLayer
  //because unfortunately the detlayer is not set in these cases
  //returns 0 for the moment
  //to be revisited

  if (tm.recHit()->det() == nullptr) {
    LogDebug("MeasurementByLayerGrouper") << "This hit has no geomdet associated skipping... ";
    return nullptr;
  }

  //now the cases in which the detid is set

  if (!theGeomSearch) {
    throw cms::Exception("MeasurementByLayerGrouper") << "Impossible to retrieve the det layer because it's not set in "
                                                         "the TM and the pointer to the GeometricSearchTracker is 0 ";
    return nullptr;
  }

  return theGeomSearch->detLayer(tm.recHit()->det()->geographicalId());
}

operator()()

vector< pair< const DetLayer *, vector< TrajectoryMeasurement > > > MeasurementByLayerGrouper::operator()

(

const std::vector< TM > &

vtm

)

const

Definition at line 10 of file MeasurementByLayerGrouper.cc.

References getDetLayer(), watchdog::group, LogTrace, mps_fire::result, isotrackNtupler::start, and trackerHitRTTI::vector.

                                 {
  if (vtm.empty())
    return vector<pair<const DetLayer*, vector<TM> > >();

  vector<pair<const DetLayer*, vector<TM> > > result;
  result.reserve(vtm.size());

  vector<TM>::const_iterator start = vtm.begin();
  //here we assume that the TM on the same detLayer are consecutive (as it should)
  while (start != vtm.end()) {
    vector<TM>::const_iterator ipart = start;
    do {
      ipart++;
    } while (ipart != vtm.end() && getDetLayer(*start) == getDetLayer(*ipart) &&
             getDetLayer(*start) != nullptr  //the returned pointer will be 0 in case
                                             //the measurement contains an invalid hit with no associated detid.
                                             //This kind of hits are at most one per layer.
                                             //this last condition avoids that 2 consecutive measurements of this kind
                                             //are grouped in the same layer.
                                             //it would be useful if invalid hit out of the active area were
                                             //given the detid reserved for the whole layer instead of 0
    );

    vector<TM> group(start, ipart);
    result.push_back(pair<const DetLayer*, vector<TM> >(getDetLayer(*start), group));
    start = ipart;
  }
#ifdef debug_MeasurementByLayerGrouper_
  for (vector<pair<const DetLayer*, vector<TM> > >::const_iterator iter = result.begin(); iter != result.end();
       iter++) {
    LogTrace("MeasurementByLayerGrouper|SiTrackerMultiRecHitUpdator")
        << "DetLayer " << iter->first << " has " << iter->second.size() << " measurements";
  }
#endif

  return result;
}

Member Data Documentation

theGeomSearch

const GeometricSearchTracker* MeasurementByLayerGrouper::theGeomSearch

private

Definition at line 16 of file MeasurementByLayerGrouper.h.