GeFiCa
Germanium detector Field Calculator
Public Member Functions | Protected Member Functions | List of all members
GeFiCa::Rho Class Reference

1D cylindrical coordinate. More...

#include <Rho.h>

Inheritance diagram for GeFiCa::Rho:
GeFiCa::Grid GeFiCa::Points

Public Member Functions

 Rho (size_t n1=101)
 
void SetupWith (Detector &detector)
 Fix potentials on boundaries based on. More...
 
void SolveAnalytically ()
 
double GetC ()
 
- 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...
 
Gridoperator*= (double scale)
 Potentials at all points are multiplied by. More...
 
Gridoperator+= (Grid &other)
 Potentials of this grid are summed with those of. More...
 
virtual FieldLineGetFieldLineFrom (double c1, double c2, double c3=0, bool positive=true)
 Propogate a field line from (c1,c2,c3). More...
 
- Public Member Functions inherited from GeFiCa::Points
size_t GetN ()
 total number of points More...
 

Protected Member Functions

void OverRelaxAt (size_t idx)
 Over relax potential Vp[. More...
 
 ClassDef (Rho, 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...
 
virtual void CalculateE ()
 Calculate Et, E1, E2, E3 from Vp. 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...
 
DetectorfDetector
 ! Pointer to associated detector object More...
 

Detailed Description

1D cylindrical coordinate.

Examples:
trueCoaxial/checkInitialization.cc, trueCoaxial/compare2analytic.cc, and trueCoaxial/verifyCV.cc.

Definition at line 9 of file Rho.h.

Constructor & Destructor Documentation

§ Rho()

GeFiCa::Rho::Rho ( size_t  n1 = 101)
inline

Definition at line 12 of file Rho.h.

12  :
13  Grid(n1) { fName="rho"; fTitle="1D cylindrical coordinates"; }
Grid(size_t n1=0, size_t n2=0, size_t n3=0)
Default constructor.
Definition: Grid.cc:26

Member Function Documentation

§ ClassDef()

GeFiCa::Rho::ClassDef ( Rho  ,
 
)
protected

§ GetC()

double Rho::GetC ( )
Examples:
trueCoaxial/verifyCV.cc.

Definition at line 53 of file Rho.cc.

54 {
55  Grid::GetC(); // calculate field excluding undepleted region
56 
57  double dV = fDetector->Bias[1]-fDetector->Bias[0]; if (dV<0) dV=-dV;
58  double integral=0;
59  for (size_t i=0; i<GetN(); i++) {
60  integral+=E1[i]*E1[i]*dC1p[i]*C1[i];
61  if (!fIsDepleted[i]) fIsFixed[i]=false; // release undepleted points
62  }
63  double c=integral*2*Pi*epsilon/dV/dV;
64  Info("GetC","%.2f pF/cm",c/pF*cm);
65  return c;
66 }
std::vector< double > Bias
bias on electrodes
Definition: Detector.h:35
static const double Pi
Definition: Units.h:29
std::vector< bool > fIsFixed
true if field values are fixed
Definition: Grid.h:128
static const double cm
centimeter
Definition: Units.h:12
std::vector< double > E1
projection of Et on C1
Definition: Grid.h:17
std::vector< double > dC1p
step length to next point alone C1
Definition: Grid.h:20
size_t GetN()
total number of points
Definition: Grid.h:26
static const double epsilon
permittivity of Ge [C/volt/cm]
Definition: Units.h:27
std::vector< bool > fIsDepleted
true if a grid point is depleted
Definition: Grid.h:129
Detector * fDetector
! Pointer to associated detector object
Definition: Grid.h:131
static const double pF
pico farad
Definition: Units.h:22
double GetC()
Get detector capacitance.
Definition: Grid.cc:198
std::vector< double > C1
the 1st coordinates of the points
Definition: Grid.h:12

§ OverRelaxAt()

void Rho::OverRelaxAt ( size_t  idx)
protectedvirtual

Over relax potential Vp[.

Parameters
idx].

Reimplemented from GeFiCa::Grid.

Definition at line 69 of file Rho.cc.

70 {
71  if (fIsFixed[idx]) return; // no need to calculate on boundaries
72 
73  // calculate Vp[idx] from Vp[idx-1] and Vp[idx+1]
74  double vnew = Src[idx]*dC1p[idx]*dC1m[idx]/2
75  + ((Vp[idx+1]-Vp[idx-1])/C1[idx]/2
76  + Vp[idx+1]/dC1p[idx]+Vp[idx-1]/dC1m[idx])
77  /(1/dC1m[idx]+1/dC1p[idx]);
78  vnew=RelaxationFactor*(vnew-Vp[idx])+Vp[idx]; // over relax
79 
80  // check depletion and update Vp[idx] accordingly
81  double min=Vp[idx-1], max=Vp[idx-1];
82  if (min>Vp[idx+1]) min=Vp[idx+1];
83  if (max<Vp[idx+1]) max=Vp[idx+1];
84  if (vnew<min) {
85  fIsDepleted[idx]=false; Vp[idx]=min;
86  } else if (vnew>max) {
87  fIsDepleted[idx]=false; Vp[idx]=max;
88  } else {
89  fIsDepleted[idx]=true; Vp[idx]=vnew;
90  }
91 
92  // update Vp for impurity-only case even if the point is undepleted
93  if (Vp[0]==Vp[N1-1]) Vp[idx]=vnew;
94 }
std::vector< double > dC1m
step length to previous point alone C1
Definition: Grid.h:21
std::vector< bool > fIsFixed
true if field values are fixed
Definition: Grid.h:128
std::vector< double > Src
-(net impurity concentration)x|Qe|/epsilon
Definition: Grid.h:56
size_t N1
number of points along the 1st coordinate
Definition: Grid.h:53
std::vector< double > dC1p
step length to next point alone C1
Definition: Grid.h:20
std::vector< bool > fIsDepleted
true if a grid point is depleted
Definition: Grid.h:129
std::vector< double > Vp
potential at each point
Definition: Grid.h:15
std::vector< double > C1
the 1st coordinates of the points
Definition: Grid.h:12
double RelaxationFactor
within (0,2), used to speed up convergence
Definition: Grid.h:58

§ SetupWith()

void Rho::SetupWith ( Detector detector)
virtual

Fix potentials on boundaries based on.

Parameters
detectorgeometry. It fills Points data based on
detectorgeometry 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.

Examples:
trueCoaxial/checkInitialization.cc, trueCoaxial/compare2analytic.cc, and trueCoaxial/verifyCV.cc.

Definition at line 6 of file Rho.cc.

7 {
8  Grid::SetupWith(detector); // check number of calls
9 
10  TString type(detector.ClassName());
11  if (type.Contains("TrueCoaxial")==false) {
12  Error("SetupWith", "%s is not expected. "
13  "Please pass in a TrueCoaxial detector.", type.Data());
14  abort();
15  }
16  TrueCoaxial& coaxial = (TrueCoaxial&) detector;
17  coaxial.CheckConfigurations();
18  fDetector = &detector; // for GetC to use fDetector->Bias[]
19 
20  double dR=coaxial.Radius-coaxial.BoreR;
21  for (size_t i=0; i<N1; i++) {
22  dC1p.push_back(dR/(N1-1)); dC1m.push_back(dR/(N1-1));
23  C1.push_back(coaxial.BoreR+i*dC1p[i]);
24  E1.push_back(0); Et.push_back(0);
25  fIsFixed.push_back(false); fIsDepleted.push_back(false);
26  Src.push_back(-coaxial.GetImpurity(C1[i])*Qe/epsilon);
27  }
28  dC1m[0]=0; dC1p[N1-1]=0;
29  // fix 1st and last points
30  fIsFixed[0]=true; fIsFixed[N1-1]=true;
31  // linear interpolation between Bias[0] and Bias[1]
32  double slope = (coaxial.Bias[1]-coaxial.Bias[0])/(N1-1);
33  for (size_t i=0; i<N1; i++) Vp.push_back(coaxial.Bias[0]+slope*i);
34  Vp[N1-1]=coaxial.Bias[1];
35 }
std::vector< double > dC1m
step length to previous point alone C1
Definition: Grid.h:21
std::vector< double > Bias
bias on electrodes
Definition: Detector.h:35
std::vector< bool > fIsFixed
true if field values are fixed
Definition: Grid.h:128
std::vector< double > E1
projection of Et on C1
Definition: Grid.h:17
std::vector< double > Src
-(net impurity concentration)x|Qe|/epsilon
Definition: Grid.h:56
double Radius
radius of the detector
Definition: TrueCoaxial.h:12
Configuration of true coaxial detectors.
Definition: TrueCoaxial.h:9
size_t N1
number of points along the 1st coordinate
Definition: Grid.h:53
virtual void SetupWith(Detector &detector)
Fix potentials on boundaries based on.
Definition: Grid.cc:111
std::vector< double > dC1p
step length to next point alone C1
Definition: Grid.h:20
static const double epsilon
permittivity of Ge [C/volt/cm]
Definition: Units.h:27
std::vector< bool > fIsDepleted
true if a grid point is depleted
Definition: Grid.h:129
Detector * fDetector
! Pointer to associated detector object
Definition: Grid.h:131
double BoreR
radius of the bore
Definition: TrueCoaxial.h:13
double GetImpurity(double height)
Return net impurity concentration at.
Definition: Detector.h:20
void CheckConfigurations()
Definition: TrueCoaxial.cc:10
std::vector< double > Vp
potential at each point
Definition: Grid.h:15
static const double Qe
electron charge in Coulomb [C]
Definition: Units.h:24
std::vector< double > Et
total electric field strength
Definition: Grid.h:16
std::vector< double > C1
the 1st coordinates of the points
Definition: Grid.h:12

§ SolveAnalytically()

void Rho::SolveAnalytically ( )
Examples:
trueCoaxial/compare2analytic.cc.

Definition at line 38 of file Rho.cc.

39 {
40  Grid::SolveAnalytically(); // check if impurity is constant
41  // https://www.wolframalpha.com/input/?i=(V%27(r)r)%27%2Fr%3Da
42  double c1 = (Vp[N1-1]-Vp[0] + Src[0]*(C1[N1-1]*C1[N1-1]-C1[0]*C1[0])/4);
43  c1/=log(C1[N1-1]/C1[0]);
44  double c2 = Vp[N1-1]*log(C1[0]) - Vp[0]*log(C1[N1-1]) +
45  Src[0]*(C1[N1-1]*C1[N1-1]*log(C1[0])-C1[0]*C1[0]*log(C1[N1-1]))/4;
46  c2/=log(C1[0])-log(C1[N1-1]);
47  for (size_t i=0; i<N1; i++)
48  Vp[i] = -Src[0]*C1[i]*C1[i]/4 + c1*log(C1[i]) + c2;
49  CalculateE();
50 }
std::vector< double > Src
-(net impurity concentration)x|Qe|/epsilon
Definition: Grid.h:56
size_t N1
number of points along the 1st coordinate
Definition: Grid.h:53
void SolveAnalytically()
Solve Poisson&#39;s Equation analytically.
Definition: Grid.cc:154
virtual void CalculateE()
Calculate Et, E1, E2, E3 from Vp.
Definition: Grid.cc:551
std::vector< double > Vp
potential at each point
Definition: Grid.h:15
std::vector< double > C1
the 1st coordinates of the points
Definition: Grid.h:12

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