MeasurementByLayerGrouper Class Reference

#include <MeasurementByLayerGrouper.h>

Public Member Functions
	MeasurementByLayerGrouper (const GeometricSearchTracker *search=0)
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

typedef TrajectoryMeasurement MeasurementByLayerGrouper::TM

private

Definition at line 17 of file MeasurementByLayerGrouper.h.

Constructor & Destructor Documentation

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

inlineexplicit

Definition at line 24 of file MeasurementByLayerGrouper.h.

:theGeomSearch(search){};

Member Function Documentation

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()==0){
      LogDebug("MeasurementByLayerGrouper") <<"This hit has no geomdet associated skipping... ";
        return 0;
        }

    //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 0;   
    }

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

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, and command_line::start.

                                                                                                                               {

    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) != 0  //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

const GeometricSearchTracker* MeasurementByLayerGrouper::theGeomSearch

private

Definition at line 18 of file MeasurementByLayerGrouper.h.