diff options
Diffstat (limited to 'src/solid_angle.c')
-rw-r--r-- | src/solid_angle.c | 18 |
1 files changed, 0 insertions, 18 deletions
diff --git a/src/solid_angle.c b/src/solid_angle.c index dbac045..7201533 100644 --- a/src/solid_angle.c +++ b/src/solid_angle.c @@ -24,24 +24,6 @@ static double lookupTable[13][11] = { {1.0041,1.0051,1.0061,1.0070,1.0078,1.0106,1.0120,1.0122,1.0116,1.0097,1.0084} }; -double get_solid_angle_fast(double *pos, double *pmt, double *n, double r) -{ - /* Returns a *very* fast approximation of the solid angle subtended by a - * circular disk of radius r at a position `pmt` with a normal vector `n` from - * a position `pos`. */ - double dir[3]; - double R, cos_theta_pmt; - - SUB(dir,pos,pmt); - - /* Distance to the PMT. */ - R = NORM(dir); - - cos_theta_pmt = fabs(DOT(dir,n)/R); - - return pow(r,2)*cos_theta_pmt/pow(R,2); -} - static double A(double u, double a, double h) { return atan(((a*a-h*h)*u - 2*a*a*h*h)/(2*a*h*sqrt((u-h*h)*(u+a*a)))); |