NavVolume6Faces.h

Go to the documentation of this file.

#ifndef NavVolume6Faces_H
#define NavVolume6Faces_H

#include "TrackPropagation/NavGeometry/interface/NavVolume.h"
#include "TrackPropagation/NavGeometry/interface/NavVolumeSide.h"

#include <vector>

class Bounds;
class Plane;
class TrajectoryStateOnSurface;

class NavVolume6Faces final : public NavVolume {
public:

    NavVolume6Faces( const PositionType& pos, const RotationType& rot, 
             DDSolidShape shape, const std::vector<NavVolumeSide>& faces,
             const MagneticFieldProvider<float> * mfp);

    explicit NavVolume6Faces( const MagVolume& magvol, const bool isIron=false);

    Container nextSurface( const NavVolume::LocalPoint& pos, 
                   const NavVolume::LocalVector& mom,
                   double charge, PropagationDirection propDir=alongMomentum) const override;

    Container nextSurface( const NavVolume::LocalPoint& pos, 
                   const NavVolume::LocalVector& mom,
                   double charge, PropagationDirection propDir,
                   ConstReferenceCountingPointer<Surface> NotThisSurfaceP) const override;

    VolumeCrossReturnType crossToNextVolume( const TrajectoryStateOnSurface& currentState,
                             const Propagator& prop) const override;

    const std::vector<VolumeSide>& faces() const override { return theFaces; }

    bool isIron() const override { return isThisIron; }

    using MagVolume::inside;
    bool inside( const GlobalPoint& gp, double tolerance) const override; 

private:

    std::vector<VolumeSide> theFaces;
    Container theNavSurfaces;
    bool isThisIron; 

    void computeBounds(const std::vector<NavVolumeSide>& faces);
    Bounds* computeBounds( int index, const std::vector<const Plane*>& bpc);
    Bounds* computeBounds( int index, const std::vector<NavVolumeSide>& faces);

};

/* bool NavVolume6Faces::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;
}
*/

#endif