Commit 7f9d40fc authored by Tyson Littenberg's avatar Tyson Littenberg

bug fix in glitch amp proposal (norm on SNR1)

parent a723fef4
......@@ -108,8 +108,11 @@ void draw_glitch_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));
/*
SNR1 == SineGaussianSNR(params)/A
Ay = SNRy / SNR1 = SNRy / (SineGaussianSNR(paramsy)/Ay)
= Ay * SNRy/SineGaussianSNR(paramsy)
*/
params[3] *= SNR/SineGaussianSNR(params, Snf, Tobs);
}
......@@ -175,8 +178,11 @@ 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));
/*
SNR1 == SineGaussianSNR(params)/A
Ay = SNRy / SNR1 = SNRy / (SineGaussianSNR(paramsy)/Ay)
= Ay * SNRy/SineGaussianSNR(paramsy)
*/
params[3] *= SNR/SineGaussianSNR(params, Snf, Tobs);
// FILE *temp=fopen("prior.dat","a");
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment