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 };
}  // namespace SurfaceOrientation

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

  ~Surface() override {}

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)
      : ReferenceCountedInConditions(iSurface),
        Base(iSurface),
        theMediumProperties(iSurface.theMediumProperties),
        theBounds(iSurface.theBounds) {}

  Surface(Surface&& iSurface)
      : ReferenceCountedInConditions(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