Garfield::ComponentParallelPlate Class Reference

Component for parallel-plate geometries. More...

#include <ComponentParallelPlate.hh>

Inheritance diagram for Garfield::ComponentParallelPlate:

Public Member Functions
	ComponentParallelPlate ()
	Constructor.
	~ComponentParallelPlate ()
	Destructor.
void	Setup (const int N, std::vector< double > eps, std::vector< double > d, const double V, std::vector< int > sigmaIndex={})
void	ElectricField (const double x, const double y, const double z, double &ex, double &ey, double &ez, Medium *&m, int &status) 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
	Calculate the drift field [V/cm] and potential [V] at (x, y, z).
double	WeightingPotential (const double x, const double y, const double z, const std::string &label) override
bool	GetVoltageRange (double &vmin, double &vmax) override
	Calculate the voltage range [V].
void	AddPixel (double x, double z, double lx, double lz, const std::string &label, bool fromAnode=true)
void	AddStrip (double z, double lz, const std::string &label, bool fromAnode=true)
	Add strip electrode.
void	AddPlane (const std::string &label, bool fromAnode=true)
void	SetMedium (Medium *medium)
	Setting the medium.
void	SetWeightingPotentialGrid (const double xmin, const double xmax, const double xsteps, const double ymin, const double ymax, const double ysteps, const double zmin, const double zmax, const double zsteps, const std::string &label)
void	SetWeightingPotentialGrids (const double xmin, const double xmax, const double xsteps, const double ymin, const double ymax, const double ysteps, const double zmin, const double zmax, const double zsteps)
void	LoadWeightingPotentialGrid (const std::string &label)
Medium *	GetMedium (const double x, const double y, const double z) override
	Get the medium at a given location (x, y, z).
bool	GetBoundingBox (double &xmin, double &ymin, double &zmin, double &xmax, double &ymax, double &zmax) override
	Get the bounding box coordinates.
bool	getLayer (const double y, int &m, double &epsM)
int	NumberOfLayers ()
void	SetIntegrationPrecision (const double eps)
void	SetIntegrationUpperbound (const double p)
void	DisablePotentialCalculationOutsideGasGap ()
std::array< double, 3 >	ElectricField (const double x, const double y, const double z)
	Calculate the drift field [V/cm] at (x, y, z).
Public Member Functions inherited from Garfield::Component
	Component ()=delete
	Default constructor.
	Component (const std::string &name)
	Constructor.
virtual	~Component ()
	Destructor.
virtual void	SetGeometry (Geometry *geo)
	Define the geometry.
virtual void	Clear ()
	Reset.
std::array< double, 3 >	ElectricField (const double x, const double y, const double z)
	Calculate the drift field [V/cm] at (x, y, z).
virtual double	ElectricPotential (const double x, const double y, const double z)
	Calculate the (drift) electrostatic potential [V] at (x, y, z).
virtual void	WeightingField (const double x, const double y, const double z, double &wx, double &wy, double &wz, const std::string &label)
virtual void	DelayedWeightingField (const double x, const double y, const double z, const double t, double &wx, double &wy, double &wz, const std::string &label)
virtual double	DelayedWeightingPotential (const double x, const double y, const double z, const double t, const std::string &label)
virtual void	MagneticField (const double x, const double y, const double z, double &bx, double &by, double &bz, int &status)
void	SetMagneticField (const double bx, const double by, const double bz)
	Set a constant magnetic field.
virtual bool	IsReady ()
	Ready for use?
virtual bool	GetElementaryCell (double &xmin, double &ymin, double &zmin, double &xmax, double &ymax, double &zmax)
	Get the coordinates of the elementary cell.
double	IntegrateFluxCircle (const double xc, const double yc, const double r, const unsigned int nI=50)
double	IntegrateFluxSphere (const double xc, const double yc, const double zc, const double r, const unsigned int nI=20)
double	IntegrateFluxParallelogram (const double x0, const double y0, const double z0, const double dx1, const double dy1, const double dz1, const double dx2, const double dy2, const double dz2, const unsigned int nU=20, const unsigned int nV=20)
double	IntegrateWeightingFluxParallelogram (const std::string &label, const double x0, const double y0, const double z0, const double dx1, const double dy1, const double dz1, const double dx2, const double dy2, const double dz2, const unsigned int nU=20, const unsigned int nV=20)
double	IntegrateFluxLine (const double x0, const double y0, const double z0, const double x1, const double y1, const double z1, const double xp, const double yp, const double zp, const unsigned int nI, const int isign=0)
virtual bool	CrossedWire (const double x0, const double y0, const double z0, const double x1, const double y1, const double z1, double &xc, double &yc, double &zc, const bool centre, double &rc)
virtual bool	InTrapRadius (const double q0, const double x0, const double y0, const double z0, double &xw, double &yw, double &rw)
virtual bool	CrossedPlane (const double x0, const double y0, const double z0, const double x1, const double y1, const double z1, double &xc, double &yc, double &zc)
void	EnablePeriodicityX (const bool on=true)
	Enable simple periodicity in the direction.
void	EnablePeriodicityY (const bool on=true)
	Enable simple periodicity in the direction.
void	EnablePeriodicityZ (const bool on=true)
	Enable simple periodicity in the direction.
void	IsPeriodic (bool &perx, bool &pery, bool &perz)
	Return periodicity flags.
void	EnableMirrorPeriodicityX (const bool on=true)
	Enable mirror periodicity in the direction.
void	EnableMirrorPeriodicityY (const bool on=true)
	Enable mirror periodicity in the direction.
void	EnableMirrorPeriodicityZ (const bool on=true)
	Enable mirror periodicity in the direction.
void	IsMirrorPeriodic (bool &perx, bool &pery, bool &perz)
	Return mirror periodicity flags.
void	EnableAxialPeriodicityX (const bool on=true)
	Enable axial periodicity in the direction.
void	EnableAxialPeriodicityY (const bool on=true)
	Enable axial periodicity in the direction.
void	EnableAxialPeriodicityZ (const bool on=true)
	Enable axial periodicity in the direction.
void	IsAxiallyPeriodic (bool &perx, bool &pery, bool &perz)
	Return axial periodicity flags.
void	EnableRotationSymmetryX (const bool on=true)
	Enable rotation symmetry around the axis.
void	EnableRotationSymmetryY (const bool on=true)
	Enable rotation symmetry around the axis.
void	EnableRotationSymmetryZ (const bool on=true)
	Enable rotation symmetry around the axis.
void	IsRotationSymmetric (bool &rotx, bool &roty, bool &rotz)
	Return rotation symmetry flags.
void	EnableDebugging ()
	Switch on debugging messages.
void	DisableDebugging ()
	Switch off debugging messages.
virtual bool	HasMagneticField () const
	Does the component have a non-zero magnetic field?
virtual bool	HasTownsendMap () const
	Does the component have maps of the Townsend coefficient?
virtual bool	HasAttachmentMap () const
	Does the component have attachment maps?
virtual bool	HasVelocityMap () const
	Does the component have velocity maps?
virtual bool	ElectronAttachment (const double, const double, const double, double &eta)
	Get the electron attachment coefficient.
virtual bool	HoleAttachment (const double, const double, const double, double &eta)
	Get the hole attachment coefficient.
virtual bool	ElectronTownsend (const double, const double, const double, double &alpha)
	Get the electron Townsend coefficient.
virtual bool	HoleTownsend (const double, const double, const double, double &alpha)
	Get the hole Townsend coefficient.
virtual bool	ElectronVelocity (const double, const double, const double, double &vx, double &vy, double &vz)
	Get the electron drift velocity.
virtual bool	HoleVelocity (const double, const double, const double, double &vx, double &vy, double &vz)
	Get the hole drift velocity.
virtual bool	GetElectronLifetime (const double, const double, const double, double &etau)
virtual bool	GetHoleLifetime (const double, const double, const double, double &htau)
virtual double	StepSizeHint ()

Additional Inherited Members
Protected Attributes inherited from Garfield::Component
std::string	m_className = "Component"
	Class name.
Geometry *	m_geometry = nullptr
	Pointer to the geometry.
std::array< double, 3 >	m_b0 = {{0., 0., 0.}}
	Constant magnetic field.
bool	m_ready = false
	Ready for use?
bool	m_debug = false
	Switch on/off debugging messages.
std::array< bool, 3 >	m_periodic = {{false, false, false}}
	Simple periodicity in x, y, z.
std::array< bool, 3 >	m_mirrorPeriodic = {{false, false, false}}
	Mirror periodicity in x, y, z.
std::array< bool, 3 >	m_axiallyPeriodic = {{false, false, false}}
	Axial periodicity in x, y, z.
std::array< bool, 3 >	m_rotationSymmetric = {{false, false, false}}
	Rotation symmetry around x-axis, y-axis, z-axis.

Detailed Description

Component for parallel-plate geometries.

Definition at line 17 of file ComponentParallelPlate.hh.

Constructor & Destructor Documentation

ComponentParallelPlate()

Garfield::ComponentParallelPlate::ComponentParallelPlate

(

)

Constructor.

Definition at line 15 of file ComponentParallelPlate.cc.

: Component("ParallelPlate") {}

~ComponentParallelPlate()

Garfield::ComponentParallelPlate::~ComponentParallelPlate ( )

inline

Destructor.

Definition at line 22 of file ComponentParallelPlate.hh.

{}

Member Function Documentation

AddPixel()

void Garfield::ComponentParallelPlate::AddPixel	(	double	x,
		double	z,
		double	lx,
		double	lz,
		const std::string &	label,
		bool	fromAnode = true )

Add a pixel electrode.

Parameters

x,z	position of the center of the electrode in the xz-plane.
lx	width in the along .
lz	width in the along .
label	give name using a string.
fromAnode	is is the electrode is the andode and if it is the cathode.

Definition at line 290 of file ComponentParallelPlate.cc.

                                                  {
 
  // Here I switch conventions back with the y-axis the direction of drift.
 
  const auto it = std::find(m_readout.cbegin(), m_readout.cend(), label);
  if (it != m_readout.end() && m_readout.size() > 0) {
    std::cerr << m_className << "::AddPixel:\n"
              << "Note that the label " << label << " is already in use.\n";
  }
  Electrode pixel;
  pixel.label = label;
  pixel.ind = structureelectrode::Pixel;
  pixel.xpos = z;
  pixel.ypos = x;
  pixel.lx = lz_input;
  pixel.ly = lx_input;
 
  pixel.formAnode = anode;
 
  m_readout.push_back(label);
  m_readout_p.push_back(std::move(pixel));
  std::cout << m_className << "::AddPixel: Added pixel electrode.\n";
}

AddPlane()

void Garfield::ComponentParallelPlate::AddPlane	(	const std::string &	label,
		bool	fromAnode = true )

Add plane electrode, if you want to read the signal from the cathode set the second argument to false.

Definition at line 337 of file ComponentParallelPlate.cc.

                                                                        {
  const auto it = std::find(m_readout.cbegin(), m_readout.cend(), label);
  if (it != m_readout.end() && m_readout.size() > 0) {
    std::cerr << m_className << "::AddPlane:\n"
              << "Note that the label " << label << " is already in use.\n";
  }
  Electrode plate;
  plate.label = label;
  plate.ind = structureelectrode::Plane;
 
  plate.formAnode = anode;
 
  m_readout.push_back(label);
  m_readout_p.push_back(std::move(plate));
 
  std::cout << m_className << "::AddPlane: Added plane electrode.\n";
}

AddStrip()

void Garfield::ComponentParallelPlate::AddStrip	(	double	z,
		double	lz,
		const std::string &	label,
		bool	fromAnode = true )

Add strip electrode.

Definition at line 316 of file ComponentParallelPlate.cc.

                                                                          {
  const auto it = std::find(m_readout.cbegin(), m_readout.cend(), label);
  if (it != m_readout.end() && m_readout.size() > 0) {
    std::cerr << m_className << "::AddStrip:\n"
              << "Note that the label " << label << " is already in use.\n";
  }
  Electrode strip;
  strip.label = label;
  strip.ind = structureelectrode::Strip;
  strip.xpos = z;
  strip.lx = lz_input;
 
  strip.formAnode = anode;
 
  m_readout.push_back(label);
  m_readout_p.push_back(std::move(strip));
 
  std::cout << m_className << "::AddStrip: Added strip electrode.\n";
}

DisablePotentialCalculationOutsideGasGap()

void Garfield::ComponentParallelPlate::DisablePotentialCalculationOutsideGasGap ( )

inline

Definition at line 137 of file ComponentParallelPlate.hh.

                                                  {
    m_getPotentialInPlate = false;
  }

ElectricField() [1/3]

std::array< double, 3 > Garfield::Component::ElectricField	(	const double	x,
		const double	y,
		const double	z )

Calculate the drift field [V/cm] at (x, y, z).

Definition at line 55 of file Component.cc.

                                                    {
  double ex = 0., ey = 0., ez = 0.;
  Medium* medium = nullptr;
  int status = 0;
  ElectricField(x, y, z, ex, ey, ez, medium, status);
  std::array<double, 3> efield = {ex, ey, ez};
  return efield;
}

ElectricField() [2/3]

void Garfield::ComponentParallelPlate::ElectricField ( const double x, const double y, const double z, double & ex, double & ey, double & ez, double & v, Medium *& m, int & status )

overridevirtual

Calculate the drift field [V/cm] and potential [V] at (x, y, z).

Implements Garfield::Component.

Definition at line 185 of file ComponentParallelPlate.cc.

                                                                    {
  // Here I switch conventions back with the y-axis the direction of drift.
 
  ex = ey = ez = v = 0;
 
  int im = -1;
  double epsM = -1;
  if (!getLayer(y, im, epsM)) {
    if (m_debug)
      std::cout << m_className << "::ElectricField: Not inside geometry.\n";
    status = -6;
    return;
  }
 
  ey = constEFieldLayer(im);
 
  v = -m_V - (y - m_z[im - 1]) * constEFieldLayer(im);
  for (int i = 1; i <= im - 1; i++) {
    v -= (m_z[i] - m_z[i - 1]) * constEFieldLayer(i);
  }
 
  m = m_geometry ? m_geometry->GetMedium(x, y, z) : m_medium;
 
  if (!m) {
    if (m_debug) {
      std::cout << m_className << "::ElectricField: No medium at (" << x << ", "
                << y << ", " << z << ").\n";
    }
    status = -6;
    return;
  }
 
  if (epsM == 1) {
    status = 0;
  } else {
    status = -5;
  }
}

ElectricField() [3/3]

void Garfield::ComponentParallelPlate::ElectricField ( const double x, const double y, const double z, double & ex, double & ey, double & ez, Medium *& m, int & status )

overridevirtual

Calculate the drift field at given point.

Parameters

x,y,z	coordinates [cm].
ex,ey,ez	components of the electric field [V/cm].
m	pointer to the medium at this location.
status	status flag

Status flags:

        0: Inside an active medium
      > 0: Inside a wire of type X
-4 ... -1: On the side of a plane where no wires are
       -5: Inside the mesh but not in an active medium
       -6: Outside the mesh
      -10: Unknown potential type (should not occur)
    other: Other cases (should not occur)

Implements Garfield::Component.

Definition at line 148 of file ComponentParallelPlate.cc.

                                                        {
  // Here I switch conventions back with the y-axis the direction of drift.
 
  ex = ey = ez = 0;
 
  int im = -1;
  double epsM = -1;
  if (!getLayer(y, im, epsM)) {
    if (m_debug)
      std::cout << m_className << "::ElectricField: Not inside geometry.\n";
    status = -6;
    return;
  }
 
  ey = constEFieldLayer(im);
 
  m = m_geometry ? m_geometry->GetMedium(x, y, z) : m_medium;
 
  if (!m) {
    if (m_debug) {
      std::cout << m_className << "::ElectricField: No medium at (" << x << ", "
                << y << ", " << z << ").\n";
    }
    status = -6;
    return;
  }
 
  if (epsM == 1) {
    status = 0;
  } else {
    status = -5;
  }
}

GetBoundingBox()

bool Garfield::ComponentParallelPlate::GetBoundingBox ( double & xmin, double & ymin, double & zmin, double & xmax, double & ymax, double & zmax )

overridevirtual

Get the bounding box coordinates.

Reimplemented from Garfield::Component.

Definition at line 81 of file ComponentParallelPlate.cc.

                                                        {
  // If a geometry is present, try to get the bounding box from there.
 
  // Here I switch conventions back with the y-axis the direction of drift.
  if (m_geometry) {
    if (m_geometry->GetBoundingBox(x0, y0, z0, x1, y1, z1)) return true;
  }
  z0 = -std::numeric_limits<double>::infinity();
  x0 = -std::numeric_limits<double>::infinity();
  z1 = +std::numeric_limits<double>::infinity();
  x1 = +std::numeric_limits<double>::infinity();
 
  y0 = 0.;
  y1 = m_z.back();
  return true;
}

getLayer()

bool Garfield::ComponentParallelPlate::getLayer ( const double y, int & m, double & epsM )

inline

Definition at line 115 of file ComponentParallelPlate.hh.

                                                      {
 
    int mholer = -1;
 
    for (int i = 1; i < m_N; i++) {
      if (y <= m_z[i]) {
        mholer = i;
        break;
      }
    }
    if (mholer == -1) return false;
 
    m = mholer;
    epsM = m_epsHolder[m - 1];
    return true;
  }

Referenced by Garfield::AvalancheGrid::AsignLayerIndex(), ElectricField(), and ElectricField().

GetMedium()

Medium * Garfield::ComponentParallelPlate::GetMedium ( const double x, const double y, const double z )

overridevirtual

Get the medium at a given location (x, y, z).

Reimplemented from Garfield::Component.

Definition at line 355 of file ComponentParallelPlate.cc.

                                                          {
  if (m_geometry) {
    return m_geometry->GetMedium(x, y, z);
  } else if (m_medium) {
    return m_medium;
  }
  return nullptr;
}

GetVoltageRange()

bool Garfield::ComponentParallelPlate::GetVoltageRange ( double & vmin, double & vmax )

overridevirtual

Calculate the voltage range [V].

Implements Garfield::Component.

Definition at line 227 of file ComponentParallelPlate.cc.

                                                                       {
  if (m_V == 0) return false;
 
  if (m_V < 0) {
    vmin = m_V;
    vmax = 0;
  } else {
    vmin = 0;
    vmax = m_V;
  }
  return true;
}

LoadWeightingPotentialGrid()

void Garfield::ComponentParallelPlate::LoadWeightingPotentialGrid ( const std::string & label )

inline

This will load a previously calculated grid of time-dependant weighting potential values.

Definition at line 91 of file ComponentParallelPlate.hh.

                                                          {
    for (auto &electrode : m_readout_p) {
      if (electrode.label == label) {
        if (electrode.grid.LoadWeightingField(label + "map", "xyz", true)) {
          std::cerr
              << m_className
              << "::LoadWeightingPotentialGrid: Weighting potential set for "
              << label << ".\n";
          electrode.m_usegrid = true;
          return;
        }
      }
    }
    std::cerr << m_className
              << "::LoadWeightingPotentialGrid: Could not find file for "
              << label << ".\n";
  }

Referenced by SetWeightingPotentialGrid().

NumberOfLayers()

int Garfield::ComponentParallelPlate::NumberOfLayers ( )

inline

Definition at line 131 of file ComponentParallelPlate.hh.

{ return m_N - 1; }

Referenced by Garfield::AvalancheGrid::AsignLayerIndex().

SetIntegrationPrecision()

void Garfield::ComponentParallelPlate::SetIntegrationPrecision ( const double eps )

inline

Definition at line 133 of file ComponentParallelPlate.hh.

{ m_precision = eps; }

SetIntegrationUpperbound()

void Garfield::ComponentParallelPlate::SetIntegrationUpperbound ( const double p )

inline

Definition at line 135 of file ComponentParallelPlate.hh.

{ m_upperBoundIntegration = p; }

SetMedium()

void Garfield::ComponentParallelPlate::SetMedium ( Medium * medium )

inline

Setting the medium.

Definition at line 66 of file ComponentParallelPlate.hh.

{ m_medium = medium; }

Setup()

void Garfield::ComponentParallelPlate::Setup	(	const int	N,
		std::vector< double >	eps,
		std::vector< double >	d,
		const double	V,
		std::vector< int >	sigmaIndex = {} )

Define the geometry.

Parameters

N	amount of layers in the geometry, this includes the gas gaps .
d	thickness of the layers starting from the bottom to the top lauer along .
eps	relative permittivities of the layers starting from the bottom to the top lauer along . Here, the gas gaps having a value of 1.
sigmaIndex	Indices of the resistive layers (optional).
V	applied potential difference between the parallel plates.

Definition at line 17 of file ComponentParallelPlate.cc.

                                                              {
 
  // Here I switch conventions with the z-axis the direction of drift.
  std::vector<double> placeHolder(N + 1, 0);
 
  const int Nholder1 = eps.size();
  const int Nholder2 = d.size();
  if (N != Nholder1 || N != Nholder2) {
    std::cout << m_className
              << "::Inconsistency between the number of layers, permittivities "
                 "and thicknesses given.\n";
    return;
  } else if (N < 2) {
 
    std::cout << m_className
              << "::Setup:: Number of layers must be larger then 1.\n";
    return;
  }
 
  if (m_debug) std::cout << m_className << "::Setup:: Loading parameters.\n";
  m_epsHolder = eps;
  m_eps = placeHolder;
 
  m_dHolder = d;
  m_d = placeHolder;
  m_N = N + 1;
  m_V = V;
 
  if (sigmaIndex.size() == 0) {
    for (int i = 0; i < N; i++) {
      if (eps[i] != 1) sigmaIndex.push_back(i + 1);
    }
  }
 
  m_sigmaIndex = sigmaIndex;
 
  std::vector<double> m_zHolder(N + 1);
  m_zHolder[0] = 0;
  for (int i = 1; i <= N; i++) {
    m_zHolder[i] = m_zHolder[i - 1] + m_dHolder[i - 1];
 
    if (m_debug)
      std::cout << m_className << "Setup:: layer " << i
                << ":: z = " << m_zHolder[i]
                << ", epsr = " << m_epsHolder[i - 1] << ".\n";
  }
  m_z = m_zHolder;
 
  if (m_debug) std::cout << m_className << "Setup:: Constructing matrices.\n";
  constructGeometryMatrices(m_N);
 
  if (m_debug)
    std::cout << m_className
              << "Setup:: Computing weighting potential functions.\n";
  setHIntegrand();
  setwpStripIntegrand();
  setwpPixelIntegrand();
 
  std::cout << m_className << "Setup:: Geometry with N = " << N
            << " layers set.\n";
}

SetWeightingPotentialGrid()

void Garfield::ComponentParallelPlate::SetWeightingPotentialGrid	(	const double	xmin,
		const double	xmax,
		const double	xsteps,
		const double	ymin,
		const double	ymax,
		const double	ysteps,
		const double	zmin,
		const double	zmax,
		const double	zsteps,
		const std::string &	label )

Calculate time-dependent weighting potential on a grid.

Parameters

xmin,ymin,zmin	minimum value of the interval in the -, - and -direction.
xmax,ymax,zmax	maximum value of the interval in the -, - and -direction.
xsteps,ysteps,zsteps	mumber of grid nodes in the -, - and -direction.
label	give name using a string.

Definition at line 595 of file ComponentParallelPlate.cc.

                            {
 
  for (auto &electrode : m_readout_p) {
    if (electrode.label == label) {
      if (electrode.m_usegrid) {
        std::cerr << m_className
                  << "::SetWeightingPotentialGrid: Overwriting grid.\n";
      }
 
      if (electrode.grid.SetMesh(xsteps, ysteps, zsteps, xmin, xmax, ymin, ymax,
                                 zmin, zmax)) {
        std::cerr << m_className << "::SetWeightingPotentialGrid: Mesh set for "
                  << label << " (for " << xsteps *ysteps *zsteps
                  << " points).\n";
      }
 
      electrode.grid.SaveWeightingField(this, label, label + "map", "xyz");
 
      LoadWeightingPotentialGrid(label);
    }
  }
}

Referenced by SetWeightingPotentialGrids().

SetWeightingPotentialGrids()

void Garfield::ComponentParallelPlate::SetWeightingPotentialGrids	(	const double	xmin,
		const double	xmax,
		const double	xsteps,
		const double	ymin,
		const double	ymax,
		const double	ysteps,
		const double	zmin,
		const double	zmax,
		const double	zsteps )

This will calculate all electrodes time-dependent weighting potential on the specified grid.

Definition at line 622 of file ComponentParallelPlate.cc.

                                                               {
 
  for (auto &electrode : m_readout_p) {
    SetWeightingPotentialGrid(xmin, xmax, xsteps, ymin, ymax, ysteps, zmin,
                              zmax, zsteps, electrode.label);
  }
}

WeightingPotential()

double Garfield::ComponentParallelPlate::WeightingPotential ( const double x, const double y, const double z, const std::string & label )

overridevirtual

Calculate the weighting potential at a given point.

Parameters

x,y,z	coordinates [cm].
label	name of the electrode.

Returns

weighting potential [dimensionless].

Reimplemented from Garfield::Component.

Definition at line 240 of file ComponentParallelPlate.cc.

                                                                          {
 
  // Here I switch conventions back with the y-axis the direction of drift.
 
  double ret = 0.;
 
  for (auto &electrode : m_readout_p) {
    if (electrode.label == label) {
      double yin = y;
      if (!electrode.formAnode) yin = m_z.back() - y;
      if (!electrode.m_usegrid) {
        ret += IntegratePromptPotential(electrode, z, x, yin);
      } else {
        ret += FindWeightingPotentialInGrid(electrode, z, x, yin);
      }
    }
  }
  return ret;
}

---

The documentation for this class was generated from the following files:

• ComponentParallelPlate.hh
• ComponentParallelPlate.cc