SolidSphere.hh

Go to the documentation of this file.

#ifndef G_SOLID_SPHERE_H
#define G_SOLID_SPHERE_H
 
#include <mutex>
 
#include "Garfield/Solid.hh"
 
namespace Garfield {

 
class SolidSphere : public Solid {
 public:
  SolidSphere(const double cx, const double cy, const double cz,
              const double r);
  SolidSphere(const double cx, const double cy, const double cz,
              const double rmin, const double rmax);
  ~SolidSphere() {}
 
  bool IsInside(const double x, const double y, const double z,
                const bool tesselated) const override;
  bool GetBoundingBox(double& xmin, double& ymin, double& zmin, double& xmax,
                      double& ymax, double& zmax) const override;
  bool IsSphere() const override { return true; }

  void SetRadius(const double r);
  void SetRadii(const double rmin, const double rmax);
 
  double GetRadius() const override { return m_rMax; }
  double GetInnerRadius() const override { return m_rMin; }
  double GetOuterRadius() const override { return m_rMax; }

  void SetMeridians(const unsigned int n);
 
  bool SolidPanels(std::vector<Panel>& panels) override;
  void SetDiscretisationLevel(const double dis) override { m_dis = dis; }
  double GetDiscretisationLevel(const Panel& panel) override;
 
  void Cut(const double x0, const double y0, const double z0, const double xn,
           const double yn, const double zn,
           std::vector<Panel>& panels) override;
 
 private:
  std::mutex m_mutex;

  double m_rMin = 0., m_rMax = 1.;

  unsigned int m_n = 10;

  double m_dis = -1.;

  std::vector<Panel> m_panelsO;
  std::vector<Panel> m_panelsI;
 
  void UpdatePanels();
  void MakePanels(const int vol, const double r, const bool out,
                  std::vector<Panel>& panels) const;
};
}  // namespace Garfield
 
#endif