SolidTube.hh

Go to the documentation of this file.

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

 
class SolidTube : public Solid {
 public:
  SolidTube(const double cx, const double cy, const double cz, const double r,
            const double lz);
  SolidTube(const double cx, const double cy, const double cz, const double r,
            const double lz, const double dx, const double dy, const double dz);
  SolidTube(const double cx, const double cy, const double cz, const double ri,
            const double ro, const double lz);
  SolidTube(const double cx, const double cy, const double cz, const double ri,
            const double ro, const double lz, const double dx, const double dy,
            const double dz);
  ~SolidTube() {}
 
  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 IsTube() const override { return true; }

  void SetHalfLength(const double lz);
  void SetRadius(const double r);
  double GetHalfLengthZ() const override { return m_lZ; }
  double GetRadius() const override { return m_rO; }
  double GetOuterRadius() const override { return m_rO; }
  double GetInnerRadius() const override { return m_rI; }

  void SetSectors(const unsigned int n);
  void SetRotation(const double angle) { m_rot = angle; }
  void SetAverageRadius(const bool average) { m_average = average; }
  void SetTopLid(const bool closed) { m_toplid = closed; }
  void SetBottomLid(const bool closed) { m_botlid = closed; }

  unsigned int GetSectors() const { return m_n; }
  double GetRotation() const { return m_rot; }
  bool GetAverage() const { return m_average; }
 
  bool SolidPanels(std::vector<Panel>& panels) override;
  void SetDiscretisationLevel(const double dis) override { m_dis.fill(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_rO;
  double m_rI = 0.;
  double m_lZ;

  double m_rot = 0.;
  unsigned int m_n = 2;
  bool m_average = false;
  double m_rpO;
  double m_rpI = 0.;
  double m_riO;
  double m_riI = 0.;
  std::vector<double> m_xpO;
  std::vector<double> m_ypO;
  std::vector<double> m_xpI;
  std::vector<double> m_ypI;

  bool m_toplid = true;
  bool m_botlid = true;

  std::array<double, 3> m_dis{{-1., -1., -1.}};
 
  void UpdatePolygon();
};
}  // namespace Garfield
 
#endif