Surface.h

Go to the documentation of this file.

#ifndef Geom_Surface_H
#define Geom_Surface_H

#include "DataFormats/GeometrySurface/interface/ReferenceCounted.h"

#include "DataFormats/GeometrySurface/interface/GloballyPositioned.h"

#include "DataFormats/GeometrySurface/interface/MediumProperties.h"
#include "DataFormats/GeometrySurface/interface/Bounds.h"

#include "FWCore/Utilities/interface/GCC11Compatibility.h"
 
#include "FWCore/Utilities/interface/clone_ptr.h"
#include <algorithm>




namespace SurfaceOrientation {
  enum Side {positiveSide, negativeSide, onSurface};
  enum GlobalFace {outer,inner,zplus,zminus,phiplus,phiminus};
}


//template <class T> class ReferenceCountingPointer;

class Plane;
using TangentPlane = Plane;

class Surface : public ReferenceCountedInConditions
              , public GloballyPositioned<float>
{
public:
  using Side =  SurfaceOrientation::Side;

  using Base = GloballyPositioned<float>;

  virtual ~Surface(){}

protected:
  Surface(){}
  Surface( const PositionType& pos, const RotationType& rot) :
    Base( pos, rot) {}
 
  Surface( const PositionType& pos, const RotationType& rot,
       Bounds* bounds) :
    Base( pos, rot),
    theBounds(bounds)
  {}

  
  Surface( const PositionType& pos, const RotationType& rot,
           MediumProperties mp) :
    Base( pos, rot),
    theMediumProperties(mp)
  {}

  Surface( const PositionType& pos, const RotationType& rot,
           MediumProperties mp,
       Bounds* bounds) :
    Base( pos, rot),
    theMediumProperties(mp),
    theBounds(bounds)
  {}

  

  Surface( const Surface& iSurface ) : 
  Base( iSurface), 
  theMediumProperties(iSurface.theMediumProperties),
  theBounds(iSurface.theBounds)
  {}
  
  Surface(Surface&& iSurface ) : 
  Base(iSurface), 
  theMediumProperties(iSurface.theMediumProperties),
  theBounds(std::move(iSurface.theBounds))
  {}

public:

  virtual Side side( const LocalPoint& p, Scalar tolerance=0) const = 0;
  virtual Side side( const GlobalPoint& p, Scalar tolerance=0) const {
    return side( toLocal(p), tolerance);
  }

  using Base::toGlobal;
  using Base::toLocal;

  GlobalPoint toGlobal( const Point2DBase< Scalar, LocalTag> lp) const {
    return GlobalPoint( rotation().multiplyInverse( lp.basicVector()) +
            position().basicVector());
  }


  const MediumProperties & mediumProperties() const { 
    return theMediumProperties;
  }

  void setMediumProperties( const MediumProperties & mp ) {
    theMediumProperties = mp;
  }

  const Bounds& bounds() const { return *theBounds; }

  // here and not in plane because of PixelBarrelLayer::overlap
  std::pair<float,float> const & phiSpan() const { return bounds().phiSpan(); }
  std::pair<float,float> const & zSpan()   const { return bounds().zSpan(); }
  std::pair<float,float> const & rSpan()   const { return bounds().rSpan(); }
  

  virtual ConstReferenceCountingPointer<TangentPlane> tangentPlane (const GlobalPoint&) const = 0;
  virtual ConstReferenceCountingPointer<TangentPlane> tangentPlane (const LocalPoint&) const = 0;

protected:
  MediumProperties theMediumProperties;
  extstd::clone_ptr<Bounds> theBounds;
};
  

#endif // Geom_Surface_H