diff --git a/src/BayesWaveProposal.c b/src/BayesWaveProposal.c
index d07127f443ed529c4ffe3cbfb18fb7a4500d1895..7d535c0d74bdae7ce5da1df704227a6d121704fc 100644
--- a/src/BayesWaveProposal.c
+++ b/src/BayesWaveProposal.c
@@ -174,7 +174,9 @@ void draw_signal_amplitude(double *params, double *Snf, gsl_rng *seed, double To
   }
 
   //SNR defined with Sn(f) but Snf array holdes <n_i^2>
-  params[3] = SNR/sqrt((params[2]/(2.0*RT2PI*params[1]))/(Snf[i]*2.0/Tobs));
+  //params[3] = SNR/sqrt((params[2]/(2.0*RT2PI*params[1]))/(Snf[i]*2.0/Tobs));
+
+  params[3] *= SNR/SineGaussianSNR(params, Snf, Tobs);
   
 //  FILE *temp=fopen("prior.dat","a");
 //  fprintf(temp,"%lg\n",params[3]);