From fef0e27b103ea85399bb55346e21261243fdecbc Mon Sep 17 00:00:00 2001 From: Matthew Gignac Date: Thu, 1 Oct 2026 07:04:13 -0700 Subject: [PATCH] gsl_ran_gaussian is already a std with width sigma, no need to multipy by sigma again -- which is a bug and makes BS too narrow --- tools/stdhep-tools/src/beam_coords.cc | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/tools/stdhep-tools/src/beam_coords.cc b/tools/stdhep-tools/src/beam_coords.cc index 3961a048e..ff839982c 100644 --- a/tools/stdhep-tools/src/beam_coords.cc +++ b/tools/stdhep-tools/src/beam_coords.cc @@ -147,8 +147,8 @@ int main(int argc, char** argv) nevhep = read_stdhep(&new_event); double shift_x = 0.0, shift_y = 0.0; - if (sigma_x > 0) shift_x = sigma_x * gsl_ran_gaussian(r, sigma_x); - if (sigma_y > 0) shift_y = sigma_y * gsl_ran_gaussian(r, sigma_y); + if (sigma_x > 0) shift_x = gsl_ran_gaussian(r, sigma_x); + if (sigma_y > 0) shift_y = gsl_ran_gaussian(r, sigma_y); double temp_x, temp_y; temp_x = shift_x * cos(theta_z) - shift_y * sin(theta_z);