Skip to content

Changes

Page history
Update Home authored by Haris K's avatar Haris K
Show whitespace changes
Inline Side-by-side
Home.md
View page @ 2ff243cf
...@@ -113,14 +113,32 @@ Investigations for GW190814 result review ...@@ -113,14 +113,32 @@ Investigations for GW190814 result review
| [Haris K](https://git.ligo.org/haris-k), [Nathan K. Johnson-McDaniel](https://git.ligo.org/nathan-johnson-mcdaniel) | Analysis of GW190814 with PE samples as templates. For this analysis data quality flags are not applied. | Done | [Link](https://git.ligo.org/soumen.roy/o3inspiralhom/-/wikis/Distribution-of-on-source-beta) | Soumen | This analysis is done without data quality flags. | :heavy_check_mark: | | [Haris K](https://git.ligo.org/haris-k), [Nathan K. Johnson-McDaniel](https://git.ligo.org/nathan-johnson-mcdaniel) | Analysis of GW190814 with PE samples as templates. For this analysis data quality flags are not applied. | Done | [Link](https://git.ligo.org/soumen.roy/o3inspiralhom/-/wikis/Distribution-of-on-source-beta) | Soumen | This analysis is done without data quality flags. | :heavy_check_mark: |
| [Haris K](https://git.ligo.org/haris-k), [Nathan K. Johnson-McDaniel](https://git.ligo.org/nathan-johnson-mcdaniel) | Injection study for GW190814 event using C00 data | Done | [Link](https://git.ligo.org/soumen.roy/o3inspiralhom/wikis/injection-Study-for-S190814bv-event) | Soumen | This exercise is done on C01 data with the final version of code. We note that the event $`\beta`$ value is on the higher tail of the distribution of $`\beta`$s obtained from the injection campaign. This might come from using the same maxL PE results regardless of the background noise. We suggest investigating this in the future. Further, we observe that the distribution of $\beta$ for h_22+h_33 injections and h_33(residual) injections are consistent. | :heavy_check_mark: | | [Haris K](https://git.ligo.org/haris-k), [Nathan K. Johnson-McDaniel](https://git.ligo.org/nathan-johnson-mcdaniel) | Injection study for GW190814 event using C00 data | Done | [Link](https://git.ligo.org/soumen.roy/o3inspiralhom/wikis/injection-Study-for-S190814bv-event) | Soumen | This exercise is done on C01 data with the final version of code. We note that the event $`\beta`$ value is on the higher tail of the distribution of $`\beta`$s obtained from the injection campaign. This might come from using the same maxL PE results regardless of the background noise. We suggest investigating this in the future. Further, we observe that the distribution of $\beta$ for h_22+h_33 injections and h_33(residual) injections are consistent. | :heavy_check_mark: |
## Updated results:
After the sign-off of the results, we identified a bug in the analysis. The $\beta$ value computed in time-frequency track analysis seems to be highly sensitive to very small changes in the masses and spins. Unfortunately, we didn't notice it earlier (though we knew already that the distribution of $\beta_33$ computed using a subset of PE samples as templates is not very narrow). Soumen and Anand did some investigation to understand the problem. The noise covariance matrix $\sigma(\alpha,\alpha')$ computed from the ensemble of noise vectors $N(\alpha)$ is rank deficient (ill-conditioned). This leads to an unstable inverse matrix (see Supplementary Material:B of https://arxiv.org/pdf/1910.04565.pdf). Replacing the normal inverse function with the pseudo inverse (Penrose Moore pseudo inverse) seems to have solved the problem. The summary of the investigations and fixes are described in [link](https://git.ligo.org/soumen.roy/o3inspiralhom/-/wikis/Issues-with-covaraince-matrix).
To fix the bug the following scripts are added/updated in the pipeline.
| Script | Short description | Status | git hash | Comment | final sign-off |
| ------ | ----------------- | ------ | -------- | ------- | -------------- |
| [inspiralhom_analysis.py](https://git.ligo.org/soumen.roy/o3inspiralhom/blob/master/scripts/inspiralhom_analysis.py) | Stand-alone python script to perform the analysis on an event. | Done | 67846be99ba1a8fa7501be1f659b8f7dd6c0c5c7 | | :X: |
| [read_insphom_results.py](https://git.ligo.org/soumen.roy/o3inspiralhom/-/blob/master/scripts/read_insphom_results.py) | Script for making plots and displaying results. | Done| 67846be99ba1a8fa7501be1f659b8f7dd6c0c5c7 | | :X: |
| [covmat_rank_analysis.py](https://git.ligo.org/soumen.roy/o3inspiralhom/-/blob/master/scripts/covmat_rank_analysis.py) | Stand-alone python script to obtain final results after fixing the covariance matrix inverse | Done | 67846be99ba1a8fa7501be1f659b8f7dd6c0c5c7 | | :X: |
# GW190412
# GW190814
| Script | Short description | Status | git hash | Comment | final sign-off |
| ------ | ----------------- | ------ | -------- | ------- | -------------- |
| [SEOBNRv4HM_test_insphom_GW190814_C01_Fs4096.sh](https://git.ligo.org/soumen.roy/o3inspiralhom/-/blob/master/GW190814/analysis/SEOBNRv4HM/test_insphom_GW190814_C01_Fs4096.sh) | Bash script to produce results (h5 file) for the event GW190814 | Done | c7caff8ea86c90fbbb6ebef932bcfa18c9acac7a | In future, the choice of `track-max-length` can be optimized to get maximum SNR. We still need to check that this is a reviewed PE run. | :heavy_check_mark: |
| [SEOBNRv4HM_plot_script.sh](https://git.ligo.org/soumen.roy/o3inspiralhom/-/blob/19f4c774f84aafa2d9f92ae31345110dd39700ae/GW190814/analysis/SEOBNRv4HM/plot_script.sh) | To produce plots and to display results | Done | c7caff8ea86c90fbbb6ebef932bcfa18c9acac7a | $`\beta=11.18`$, Number of background samples= 4032, Largest background = 6.88, so upper bound on the p-value is $`1/4032 \simeq 2.5\times 10^{-4}`$. | :heavy_check_mark: |
| [SEOBNRv4_highsrate_test_insphom_GW190814_C01_Fs4096.sh](https://git.ligo.org/soumen.roy/o3inspiralhom/-/blob/master/GW190814/analysis/SEOBNRv4_highsrate/test_insphom_GW190814_C01_Fs4096.sh) | Bash script to produce results (h5 file) for the event GW190814 (based on PE results [Prod12](https://ldas-jobs.ligo.caltech.edu/~charlie.hoy/PE/O3/S190814bv/C01/SEOBNRv4_ROM/higher_srate/git_repo/home.html)) | Done | 0c074b0d2a9246fb96440ad856aa849f4a8f95b1 | | :heavy_check_mark: |
| [SEOBNRv4_highsrate_plot_script.sh](https://git.ligo.org/soumen.roy/o3inspiralhom/-/blob/master/GW190814/analysis/SEOBNRv4_highsrate/plot_script.sh) | To produce plots and to display results | Done | 0c074b0d2a9246fb96440ad856aa849f4a8f95b1 | $`\beta=10.83`$, Number of background samples= 4032, Largest background = 7.83, so upper bound on the p-value is $`1/4032 \simeq 2.5\times 10^{-4}`$. | :heavy_check_mark: |
| [IMRPhenomHM_test_insphom_GW190814_C01_Fs4096.sh](https://git.ligo.org/soumen.roy/o3inspiralhom/-/blob/master/GW190814/analysis/IMRPhenomHM/test_insphom_GW190814_C01_Fs4096.sh) | Bash script to produce results (h5 file) for the event GW190814 with IMRPhenomHM PE results | Done | | The PE result used is not the final one. The runs are still in progress. | :x: |
| [IMRPhenomHM_plot_script.sh](https://git.ligo.org/soumen.roy/o3inspiralhom/-/blob/master/GW190814/analysis/IMRPhenomHM/plot_script.sh) | To produce plots and to display results | Done | | $`\beta_{33}=10.64`$, Number of background samples= 4032, Largest background = 6.82. | :x: |
# Final results
| Event | Waveform | on-source $`\beta`$ | p-value | max off-source $`\beta`$ |
| ------ | ------- | ----------------- | ------ | :-----------------------: |
| GW190412 | IMRPhenomHM | 5.94 | 1/1572 | 6.14 |
| GW190412 | SEOBNRv4HM | 5.35 | 4/1572 | 6.0 |
| GW190814 | SEOBNRv4HM | 11.18 | < 1/4032 | 6.6 |
| GW190814 | SEOBNRv4ROM_highsrate | 9.33 | < 1/4032 | 6.4 |
# Link to review telecons minutes # Link to review telecons minutes
......