ComponentTcad2d.hh

Go to the documentation of this file.

#ifndef G_COMPONENT_TCAD_2D_H
#define G_COMPONENT_TCAD_2D_H
 
#include <array>
#include <memory>
 
#include "Garfield/ComponentTcadBase.hh"
#include "QuadTree.hh"
 
namespace Garfield {

 
class ComponentTcad2d : public ComponentTcadBase<2> {
 public:
  ComponentTcad2d();
  ~ComponentTcad2d() {}

  void SetRangeZ(const double zmin, const double zmax);

  bool GetElement(const size_t i, double& vol, double& dmin, double& dmax,
                  int& type) const;
  bool GetNode(const size_t i, double& x, double& y, double& z) const override;
 
  void ElectricField(const double x, const double y, const double z, double& ex,
                     double& ey, double& ez, double& v, Medium*& m,
                     int& status) override;
 
  void ElectricField(const double x, const double y, const double z, double& ex,
                     double& ey, double& ez, Medium*& m, int& status) override {
    double v = 0.;
    ElectricField(x, y, z, ex, ey, ez, v, m, status);
  }
  using Component::ElectricField;
  Medium* GetMedium(const double x, const double y, const double z) override;
  void DelayedWeightingPotentials(const double x, const double y,
                                  const double z, const std::string& label,
                                  std::vector<double>& dwp) override;
 
  bool GetBoundingBox(double& xmin, double& ymin, double& zmin, double& xmax,
                      double& ymax, double& zmax) override;
  bool GetElementaryCell(double& xmin, double& ymin, double& zmin, double& xmax,
                         double& ymax, double& zmax) override;
 
 private:
  // Bounding box
  bool m_hasRangeZ = false;
 
  // Tetrahedral tree.
  std::unique_ptr<QuadTree> m_tree;
 
  void Reset() override {
    Cleanup();
    m_hasRangeZ = false;
    m_ready = false;
  }
 
  size_t FindElement(const double x, const double y,
                     std::array<double, nMaxVertices>& w) const;
  // Check whether a point is inside a given element and calculate the
  // shape functions if it is.
  bool InElement(const double x, const double y, const Element& element,
                 std::array<double, nMaxVertices>& w) const;
  bool InRectangle(const double x, const double y, const Element& element,
                   std::array<double, nMaxVertices>& w) const;
  bool InTriangle(const double x, const double y, const Element& element,
                  std::array<double, nMaxVertices>& w) const;
  bool OnLine(const double x, const double y, const Element& element,
              std::array<double, nMaxVertices>& w) const;
  bool AtPoint(const double x, const double y, const Element& element,
               std::array<double, nMaxVertices>& w) const;
 
  bool Interpolate(const double x, const double y, const double z,
                   const std::vector<double>& field, double& f) override;
  bool Interpolate(const double x, const double y, const double z,
                   const std::vector<std::array<double, 2> >& field, double& fx,
                   double& fy, double& fz) override;
  void FillTree() override;
};
}  // namespace Garfield
#endif