Fixed incorrect LSF width used in convolution

cschreib · Aug 4, 2023 · fc7e38e · fc7e38e
1 parent e67b18c
commit fc7e38e
Show file tree

Hide file tree

Showing 2 changed files with 1 addition and 2 deletions.
diff --git a/src/fast++-fitter.cpp b/src/fast++-fitter.cpp
@@ -544,7 +544,6 @@ void fitter_t::fit_galaxies(const model_t& model, uint_t i0, uint_t i1) {
 
             wsp.wflux[iw] = input.flux.safe(is,il)*wsp.weight[iw];
             wsp.wmodel[iw] = model.flux.safe[il]*wsp.weight[iw];
-
         };
 
         for (uint_t iw : range(iflx, nscale)) {

diff --git a/src/fast++-read_input.cpp b/src/fast++-read_input.cpp
@@ -1816,7 +1816,7 @@ bool read_lsf_spec(const options_t& opts, input_state_t& state) {
         vec1d l = rgen(l0-dll/2, l0+dll/2, npt);
 
         f.wl = l;
-        f.tr = 0.5*(erf((l0+dl/2-l)/(sqrt(2)*s)) - erf((l0-dl/2-l)/(sqrt(2)*s)));
+        f.tr = 0.5*(erf((l0+dl/2-l)/(sqrt(2)*lsf)) - erf((l0-dl/2-l)/(sqrt(2)*lsf)));
 
         double ttot = integrate(f.wl, f.tr);
         if (!is_finite(ttot) || ttot == 0) {