/data/refman/pasoursint/CMSSW_5_3_0/src/TrackPropagation/NavGeometry/src/LinearSearchNavSurfaceImpl.cc

Go to the documentation of this file.

#include "TrackPropagation/NavGeometry/interface/LinearSearchNavSurfaceImpl.h"
#include "DataFormats/GeometrySurface/interface/Bounds.h"

const NavVolume* LinearSearchNavSurfaceImpl::nextVolume( const NavSurface::LocalPoint& point, 
                                                         SurfaceOrientation::Side side) const
{
    if (side == SurfaceOrientation::positiveSide) {
        return findVolume( point, thePosVolumes);
    }
    else if (side == SurfaceOrientation::negativeSide) {
        return findVolume( point, theNegVolumes);
    }
    else return 0;  // should not be called with SurfaceOrientation::onSurface
}

const NavVolume* LinearSearchNavSurfaceImpl::findVolume( const NavSurface::LocalPoint& point,
                                                         const VolumeContainer& vols) const
{
// simple linear search for volume who's bounds contain the point

//MM:return 0 if no volume was defined on this side! 
  if (vols.size()==0) return 0;

    for (VolumeContainer::const_iterator i=vols.begin(); i!=vols.end(); i++) {
        if (i->second->inside(point)) return i->first;
    }
    return 0; // if point outside of all bounds on this side
}

const Bounds* LinearSearchNavSurfaceImpl::bounds( const NavVolume* vol)
{
    for (VolumeContainer::const_iterator i=thePosVolumes.begin(); i!=thePosVolumes.end(); i++) {
        if (i->first == vol) return i->second;
    }
    for (VolumeContainer::const_iterator i=theNegVolumes.begin(); i!=theNegVolumes.end(); i++) {
        if (i->first == vol) return i->second;
    }
    return 0; // if volume not found
}

void LinearSearchNavSurfaceImpl::addVolume( const NavVolume* vol, const Bounds* bounds, 
                                            SurfaceOrientation::Side side)
{
    if (side == SurfaceOrientation::positiveSide) {
        thePosVolumes.push_back( VolumeAndBounds( vol, bounds->clone()));
    }
    else {
        theNegVolumes.push_back( VolumeAndBounds( vol, bounds->clone()));
    }
}