Skip to content
Snippets Groups Projects
Commit 88c52d2a authored by Colm Talbot's avatar Colm Talbot
Browse files

BUGFIX: add case of just time and calibration marginalzation

parent ca4b03f1
No related branches found
No related tags found
1 merge request!1234Time calibration
Pipeline #506076 passed
......@@ -767,6 +767,8 @@ class GravitationalWaveTransient(Likelihood):
log_l_tc_array = self.phase_marginalized_likelihood(
d_inner_h=d_inner_h_tc_array,
h_inner_h=h_inner_h)
elif self.calibration_marginalization:
log_l_tc_array = np.real(d_inner_h_tc_array) - h_inner_h[:, np.newaxis] / 2
else:
log_l_tc_array = np.real(d_inner_h_tc_array) - h_inner_h / 2
return logsumexp(log_l_tc_array, b=time_prior_array, axis=-1)
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment