Skip to content
Snippets Groups Projects

Draft: Apply detector resolution to offspec simulation (after !2896)

Closed Mikhail Svechnikov requested to merge m.0 into main
36 files
+ 1267
1296
Compare changes
  • Side-by-side
  • Inline
Files
36
+ 3
0
@@ -41,6 +41,9 @@ double Math::Gaussian(double x, double average, double std_dev)
@@ -41,6 +41,9 @@ double Math::Gaussian(double x, double average, double std_dev)
double Math::IntegratedGaussian(double x, double average, double std_dev)
double Math::IntegratedGaussian(double x, double average, double std_dev)
{
{
 
if (std_dev <= std::numeric_limits<double>::min())
 
return Math::Heaviside(x - average);
 
double normalized_x = (x - average) / std_dev;
double normalized_x = (x - average) / std_dev;
static double root2 = std::sqrt(2.0);
static double root2 = std::sqrt(2.0);
return (gsl_sf_erf(normalized_x / root2) + 1.0) / 2.0;
return (gsl_sf_erf(normalized_x / root2) + 1.0) / 2.0;
Loading