GeFiCa::RTheta Class Reference

2D spherical coordinates. More...

#include <RTheta.h>

Inheritance diagram for GeFiCa::RTheta:

Collaboration diagram for GeFiCa::RTheta:

Public Member Functions
	RTheta (size_t n1=101, size_t n2=181)
void	SetupWith (Detector &detector)
	Fix potentials on boundaries based on. More...
Public Member Functions inherited from GeFiCa::Grid
	Grid (size_t n1=0, size_t n2=0, size_t n3=0)
	Default constructor. More...
virtual	~Grid ()
void	SuccessiveOverRelax ()
	Successively over-relax potentials on grid points. More...
void	SolveAnalytically ()
	Solve Poisson's Equation analytically. More...
double	GetV (double c1, double c2=0, double c3=0) const
	Get potential at (c1,c2,c3) by interpolation. More...
double	GetE (double c1, double c2=0, double c3=0) const
double	GetE1 (double c1, double c2=0, double c3=0) const
double	GetE2 (double c1, double c2=0, double c3=0) const
double	GetE3 (double c1, double c2=0, double c3=0) const
double	GetC ()
	Get detector capacitance. More...
TTree *	GetTree (bool createNew=false)
	Create &/or return a TTree with field data. More...
bool	IsDepleted ()
	Check if every grid point is depleted. More...
Grid &	operator*= (double scale)
	Potentials at all points are multiplied by. More...
Grid &	operator+= (Grid &other)
	Potentials of this grid are summed with those of. More...
virtual FieldLine *	GetFieldLineFrom (double c1, double c2, bool positive=true)
	Propogate a field line from (c1,c2,c3). More...
virtual FieldLine *	GetFieldLineFrom (double c1, double c2, double c3, bool positive=true)
Public Member Functions inherited from GeFiCa::Points
size_t	GetN ()
	total number of points More...

Protected Member Functions
virtual void	OverRelaxAt (size_t idx)
	Over relax potential Vp[. More...
void	CalculateE ()
	Calculate Et, E1, E2, E3 from Vp. More...
	ClassDef (RTheta, 1)
Protected Member Functions inherited from GeFiCa::Grid
size_t	GetIdxOfPointToTheRightOf (double c1, size_t begin, size_t end) const
	Get index of point near. More...
size_t	GetIdxOfPointToTheRightOf (double c1, double c2, size_t begin, size_t end) const
size_t	GetIdxOfPointToTheRightOf (double c1, double c2, double c3, size_t begin, size_t end) const
virtual double	GetData (const std::vector< double > &data, double c1, double c2, double c3) const
	Interpolate grid data at (c1,c2,c3). More...
size_t	GetIdxOfMaxV ()
	Get index of the grid point with max potential. More...
size_t	GetIdxOfMinV ()
	Get index of the grid point with min potential. More...
double	twopoint (double dataset[2], double tarlocationset, double pointxset[2]) const
	Calculate interpolation value between two point. More...
double	threepoint (double dataset[3], double tarlocationset[2], double pointxset[3], double pointyset[3]) const
	Calculate interpolation value between three point (triangle) More...
double	fourpoint (double dataset[4], double tarlocationset[2], double pointxset[4], double pointyset[4]) const
	Calculate interpolation value between four point (rectangle) More...
	ClassDef (Grid, 1)

Additional Inherited Members
Public Attributes inherited from GeFiCa::Grid
size_t	N1
	number of points along the 1st coordinate More...
size_t	N2
	number of points along the 2nd coordinate More...
size_t	N3
	number of points along the 3rd coordinate More...
std::vector< double >	Src
	-(net impurity concentration)x|Qe|/epsilon More...
size_t	MaxIterations
	maximal iterations of SOR to be performed More...
double	RelaxationFactor
	within (0,2), used to speed up convergence More...
double	Tolerance
	SOR stops when error<Tolerance. More...
size_t	Iterations
	number of iterations of SOR performed More...
Public Attributes inherited from GeFiCa::Points
std::vector< double >	C1
	the 1st coordinates of the points More...
std::vector< double >	C2
	the 2nd coordinates of the points More...
std::vector< double >	C3
	the 3rd coordinates of the points More...
std::vector< double >	Vp
	potential at each point More...
std::vector< double >	Et
	total electric field strength More...
std::vector< double >	E1
	projection of Et on C1 More...
std::vector< double >	E2
	projection of Et on C2 More...
std::vector< double >	E3
	projection of Et on C3 More...
std::vector< double >	dC1p
	step length to next point alone C1 More...
std::vector< double >	dC1m
	step length to previous point alone C1 More...
std::vector< double >	dC2p
	step length to next point along C2 More...
std::vector< double >	dC2m
	step length to previous point along C2 More...
std::vector< double >	dC3p
	step length to next point alone C3 More...
std::vector< double >	dC3m
	step length to previous point alone C3 More...
Protected Attributes inherited from GeFiCa::Grid
std::vector< bool >	fIsFixed
	true if field values are fixed More...
std::vector< bool >	fIsDepleted
	true if a grid point is depleted More...
TTree *	fTree
	! ROOT tree to visualize fields More...
Detector *	fDetector
	! Pointer to associated detector object More...

Detailed Description

2D spherical coordinates.

Definition at line 9 of file RTheta.h.

Constructor & Destructor Documentation

RTheta()

GeFiCa::RTheta::RTheta ( size_t  n1 = 101, size_t  n2 = 181  )

inline

Definition at line 12 of file RTheta.h.

                                           : Grid(n1, n2) {
         fName="rt"; fTitle="2D spherical coordinate"; }

Member Function Documentation

CalculateE()

void RTheta::CalculateE ( )

protectedvirtual

Calculate Et, E1, E2, E3 from Vp.

Reimplemented from GeFiCa::Grid.

Definition at line 63 of file RTheta.cc.

{
   double N=N1*N2;
   for (size_t i=1; i<N-1; i++) {
      if(i<N1-1)
         E2[i]=(C2[i+N1]*Vp[i+N1]-C2[i]*Vp[i])/(dC2p[i])/C2[i]/C1[i];
      else if(i>N-N1-1)
         E2[i]=(C2[i]*Vp[i]-C2[i-N1]*Vp[i-N1])/(dC2m[i])/C2[i]/C1[i];
      else 
         E2[i]=(C2[i+N1]*Vp[i+N1]-C2[i-N1]*Vp[i-N1])/(dC2p[i]+dC2m[i])/C2[i]/C1[i];
      if(i%N1==0)
         E1[i]=(C1[i+1]*Vp[i+1]-C1[i]*Vp[i])/(dC1p[i])/C1[i];
      else if(i%N1==N1-1)
         E1[i]=(C1[i]*Vp[i]-C1[i-1]*Vp[i-1])/(dC1m[i])/C1[i];
      else E1[i]=(C1[i+1]*Vp[i+1]-C1[i-1]*Vp[i-1])/(dC1p[i]+dC1m[i])/C1[i];
      Et[i]=sqrt(E1[i]*E1[i]+E2[i]*E2[i]);
   }
}

ClassDef()

GeFiCa::RTheta::ClassDef ( RTheta  , 1    )

protected

OverRelaxAt()

void RTheta::OverRelaxAt ( size_t  idx )

protectedvirtual

Over relax potential Vp[.

Parameters

idx].

Reimplemented from GeFiCa::Grid.

Definition at line 38 of file RTheta.cc.

{
   if (fIsFixed[idx]) return; // no need to calculate on boundaries

   double vnew=Src[idx]*dC1m[idx]*dC1p[idx]/2+(1/C1[idx]*(Vp[idx+1]-Vp[idx-1])
         +Vp[idx+1]/dC1p[idx]+Vp[idx-1]/dC1m[idx])/(1/dC1m[idx]+1/dC1p[idx]);
   vnew=RelaxationFactor*(vnew-Vp[idx])+Vp[idx];

   // check depletion and update Vp[idx] accordingly
   double min=Vp[idx-1], max=Vp[idx-1];
   if (min>Vp[idx+1]) min=Vp[idx+1];
   if (max<Vp[idx+1]) max=Vp[idx+1];
   if (vnew<min) {
      fIsDepleted[idx]=false; Vp[idx]=min;
   } else if (vnew>max) {
      fIsDepleted[idx]=false; Vp[idx]=max;
   } else {
      fIsDepleted[idx]=true; Vp[idx]=vnew;
   }

   // update Vp for impurity-only case even if the point is undepleted
   if (Vp[0]==Vp[N1-1]) Vp[idx]=vnew;
}

SetupWith()

void RTheta::SetupWith ( Detector &  detector )

virtual

Fix potentials on boundaries based on.

Parameters

detector	geometry. It fills Points data based on
detector	geometry and N1, N2 and/or N3, and raises the flag fIsFixed for points on/outside electrodes. It has to be called before SuccessiveOverRelax().

Reimplemented from GeFiCa::Grid.

Definition at line 6 of file RTheta.cc.

{
   Grid::SetupWith(detector); // check number of calls

   TString type(detector.ClassName());
   if (type.Contains("Hemispherical")==false) {
      Error("SetupWith", "%s is not expected. "
            "Please pass in a Hemispherical detector.", type.Data());
      abort();
   }
   Hemispherical& hemi = (Hemispherical&) detector;
   hemi.CheckConfigurations();
   fDetector = &detector; // for GetC to use fDetector->Bias[]

   double dR=hemi.Height-hemi.PointContactR;
   for (size_t i=0; i<N1; i++) {
      dC1p.push_back(dR/(N1-1)); dC1m.push_back(dR/(N1-1));
      C1.push_back(hemi.PointContactR+i*dC1p[i]);
      E1.push_back(0); Et.push_back(0);
      fIsFixed.push_back(false); fIsDepleted.push_back(false);
      Src.push_back(-hemi.GetImpurity(C1[i])*Qe/epsilon);
   }
   dC1m[0]=0; dC1p[N1-1]=0;
   // fix 1st and last points
   fIsFixed[0]=true; fIsFixed[N1-1]=true;
   // linear interpolation between Bias[0] and Bias[1]
   double slope = (hemi.Bias[1]-hemi.Bias[0])/(N1-1);
   for (size_t i=0; i<N1; i++) Vp.push_back(hemi.Bias[0]+slope*i);
   Vp[N1-1]=hemi.Bias[1];
}

---

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

• RTheta.h
• RTheta.cc