|
|
Gravelamps code is developed by Michael Wright. The review of Gravelamps is undertaken by Suvodip Mukherjee. The review chair is Martin Hendry.
|
|
|
|
|
|
**Code Review**
|
|
|
|
|
|
A google document containing the minutes of the Code Review meetings is located [**here**](https://docs.google.com/document/d/1MtFdW1W6pQlDVMhevHDy-J6CEilzCwN51kF-I-6w6BY/edit)**.**
|
|
|
|
|
|
Materials for the Code Review can be found in the `review_materials` folder of this repository located [here](https://git.ligo.org/michael.wright/Gravelamps/-/tree/main/review_materials). Inside of this are scripts containing tests for each of the functions within `gravelamps.inference` and `gravelamps.lensing` the two submodules of Gravelamps' code.
|
|
|
|
|
|
**Reviewer's comments:** The Gravelamps code jointly infers GW source parameters and the parameters associated with lens within a Bayesian setup. This version of the code uses only the distortion in the GW waveform due to lensing to infer the parameters related to lens and source. It does not explore the time delay between multiple images. So, this code will be useful to infer the properties of the GW sources which are already confirmed lensed events.
|
|
|
|
|
|
I have checked the code and in this process, we have obtained results from the following modules: [utils_testing](https://git.ligo.org/michael.wright/Gravelamps/-/blob/main/review_materials/utils_testing/utils_testing.py), [lens_generation_testing](https://git.ligo.org/michael.wright/Gravelamps/-/tree/main/review_materials/lens_generation_testing), [helpers_testing](https://git.ligo.org/michael.wright/Gravelamps/-/tree/main/review_materials/helpers_testing), [file_generators_testing ](https://git.ligo.org/michael.wright/Gravelamps/-/tree/main/review_materials/file_generators_testing), [gravelamps_lensing_testing](https://git.ligo.org/michael.wright/Gravelamps/-/tree/main/review_materials/gravelamps_lensing_testing)
|
|
|
|
|
|
The input values match with the output results from these codes.
|
|
|
|
|
|
Next, we reviewed the Bayesian setup and joint inference of the parameters for three different sets of injected parameters such as luminosity distance, and inclination angle. As the analysis is computationally costly, we have used three sets of runs to validate the parameter estimation of the code.
|
|
|
|
|
|
The contour plots agree to corresponding injected values for all the parameters except for the parameter \\phi. Though the parameter \\phi is usually poorly constrained and it will not affect the lensing inference significantly, I expect this is not going to be a problem for the O3 analysis. However, for O4 analysis, it will be important to study the reason for this discrepancy. One further clarification, in the contour plots, the lines corresponding to the injected value for chirp mass and mass ratio are not present. However, the posterior distribution on these two parameters matches the injected values (check the input file).
|
|
|
|
|
|
The results from all these analyses are present in the tables below.
|
|
|
|
|
|
Examples of generated amplification factors can be found at the following locations:
|
|
|
| Lens Type |
|
|
|
|-----------|
|
|
|
| [Point Mass Lens - Wave Optics](https://git.ligo.org/michael.wright/Gravelamps/-/blob/main/review_materials/lens_generation_testing/pointlens-wave_plot.pdf) |
|
|
|
| [Point Mass Lens - Geometric Optics](https://git.ligo.org/michael.wright/Gravelamps/-/blob/main/review_materials/lens_generation_testing/pointlens-geo_plot.pdf) |
|
|
|
| [SIngular Isothermal Sphere - Wave Optics](https://git.ligo.org/michael.wright/Gravelamps/-/blob/main/review_materials/lens_generation_testing/sislens-wave_plot.pdf) |
|
|
|
| [Singular Isothermal Sphere - Geometric Optics](https://git.ligo.org/michael.wright/Gravelamps/-/blob/main/review_materials/lens_generation_testing/sislens-geo_plot.pdf) |
|
|
|
| [Navarro, Frenk, White - Wave Optics](https://git.ligo.org/michael.wright/Gravelamps/-/blob/main/review_materials/lens_generation_testing/nfwlens-wave_plot.pdf) |
|
|
|
| [Navarro, Frenk, White - Geometric Optics](https://git.ligo.org/michael.wright/Gravelamps/-/blob/main/review_materials/lens_generation_testing/nfwlens-geo_plot.pdf) |
|
|
|
|
|
|
Results of example test runs can be found in the following locations:
|
|
|
|
|
|
<table>
|
|
|
<tr>
|
|
|
<th>Run Type</th>
|
|
|
<th>Run Directory</th>
|
|
|
<th>Corner Plot Link</th>
|
|
|
</tr>
|
|
|
<tr>
|
|
|
<td>GW150914-like Injection</td>
|
|
|
<td>
|
|
|
|
|
|
[using_reviewed_lens_settings/pointlens_injection](https://ligo.gravity.cf.ac.uk/\~michael.wright/review_runs/using_revied_lens_settings/pointlens_injection/)
|
|
|
</td>
|
|
|
<td>
|
|
|
|
|
|
[injection_review.png](https://ligo.gravity.cf.ac.uk/\~michael.wright/injection_review.png)
|
|
|
</td>
|
|
|
</tr>
|
|
|
<tr>
|
|
|
<td>
|
|
|
|
|
|
GW150914-like Injection<br>with double distance
|
|
|
</td>
|
|
|
<td>
|
|
|
|
|
|
[higher_distance](https://ligo.gravity.cf.ac.uk/\~michael.wright/review_runs/higher_distance/)
|
|
|
</td>
|
|
|
<td>
|
|
|
|
|
|
[higher_distance_corner.png](https://ligo.gravity.cf.ac.uk/\~michael.wright/higher_distance_corner.png)
|
|
|
</td>
|
|
|
</tr>
|
|
|
<tr>
|
|
|
<td>
|
|
|
|
|
|
GW150914-like Injection<br>with higher inclination angle
|
|
|
</td>
|
|
|
<td>
|
|
|
|
|
|
[inclination_angle_change](https://ligo.gravity.cf.ac.uk/\~michael.wright/review_runs/inclination_angle_change/)
|
|
|
</td>
|
|
|
<td>
|
|
|
|
|
|
[high_inclination_corner.png](https://ligo.gravity.cf.ac.uk/\~michael.wright/high_inclination.corner.png)
|
|
|
</td>
|
|
|
</tr>
|
|
|
<tr>
|
|
|
<td>GW170809 Real Event Run</td>
|
|
|
<td>
|
|
|
|
|
|
[using_reviewed_lens_settings/pointlens_realevent](https://ligo.gravity.cf.ac.uk/\~michael.wright/review_runs/using_revied_lens_settings/pointlens_realevent/)
|
|
|
</td>
|
|
|
<td>
|
|
|
|
|
|
[real_event.png](https://ligo.gravity.cf.ac.uk/\~michael.wright/real_event.png)
|
|
|
</td>
|
|
|
</tr>
|
|
|
<tr>
|
|
|
<td>
|
|
|
|
|
|
GW150914-like Injection<br>with 10^6 Solar Mass lens<br>using direct injection (SIS)
|
|
|
</td>
|
|
|
<td>
|
|
|
|
|
|
[using_direct_implementation/sislens_injection_direct](https://ligo.gravity.cf.ac.uk/\~michael.wright/using_direct_implementation/sislens_injection_direct/)
|
|
|
</td>
|
|
|
<td>
|
|
|
|
|
|
[sislens_direct_corner.png](https://ligo.gravity.cf.ac.uk/\~michael.wright/sislens_direct_corner.png)
|
|
|
</td>
|
|
|
</tr>
|
|
|
<tr>
|
|
|
<td>
|
|
|
|
|
|
GW150914-like Injection
|
|
|
|
|
|
with 10^6 Solar Mass lens using direct injection (NFW)
|
|
|
</td>
|
|
|
<td>
|
|
|
|
|
|
[using_direct_implementation/nfwlens_injection_direct](https://ligo.gravity.cf.ac.uk/\~michael.wright/using_direct_implementation/nfwlens_injection_direct/)
|
|
|
</td>
|
|
|
<td>
|
|
|
|
|
|
[nfw_direct_corner.png](https://ligo.gravity.cf.ac.uk/\~michael.wright/nfw_direct_corner.png)
|
|
|
</td>
|
|
|
</tr>
|
|
|
</table>
|
|
|
|
|
|
**Results Review** |
|
|
\ No newline at end of file |