/afs/cern.ch/work/a/aaltunda/public/www/CMSSW_6_2_7/src/MagneticField/VolumeGeometry/src/MagVolume6Faces.cc

Go to the documentation of this file.

#include "MagneticField/VolumeGeometry/interface/MagVolume6Faces.h"

MagVolume6Faces::MagVolume6Faces( const PositionType& pos,
                                  const RotationType& rot, 
                                  DDSolidShape shape,
                                  const std::vector<VolumeSide>& faces,
                                  const MagneticFieldProvider<float> * mfp,
                                  double sf)
  : MagVolume(pos,rot,shape,mfp,sf),  volumeNo(0), copyno(0), theFaces(faces)
{}

bool MagVolume6Faces::inside( const GlobalPoint& gp, double tolerance) const 
{

  // check if the point is on the correct side of all delimiting surfaces
  for (std::vector<VolumeSide>::const_iterator i=theFaces.begin(); i!=theFaces.end(); ++i) {
    Surface::Side side = i->surface().side( gp, tolerance);
    if ( side != i->surfaceSide() && side != SurfaceOrientation::onSurface) return false;
  }
  return true;
}