Differences

This shows you the differences between two versions of the page.

--- imrphenomhm_pe_study [2017/09/15 11:47] – chinmay.kalaghatgi@LIGO.ORG
+++ imrphenomhm_pe_study [2017/09/18 14:11] (current) – [PE Runs] chinmay.kalaghatgi@LIGO.ORG
@@ Line 7: / Line 7: @@
   * The total response of a detector network (2 detector) is taken to be : sqrt( F<sup>2</sup><sub>+, H</sub> + F<sup>2</sup><sub>x, H</sub> + F<sup>2</sup><sub>+, L</sub> + F<sup>2</sup><sub>x, L</sub>). If a third detector is added, then we include F<sub>+/x</sub> for that detector to compute the total response.
-  * The response function depends on the ra, dec values of the sky-position along with the polarisation of the waveform and gps trigger time. So, we fix the polarisation to 1.75 radians (~ 100 deg) and the gps time to 1187128000. Constraining the values of ra, dec between (0, 2 pi) and (-pi/2, pi/2) respectively, the plot of response function for HL and HLV. The blue point is the chosen sky-position for setting up the PE runs.
+  * The response function depends on the ra, dec values of the sky-position along with the polarisation of the waveform and gps trigger time. So, we fix the polarisation to 1.75 radians (~ 100 deg) and the gps time to 1187128000. Constraining the values of ra, dec between (0, 2 pi) and (-pi/2, pi/2) respectively, the plot of response function for HL and HLV. The blue point is the chosen sky-position for setting up the PE runs. The blue dot has the values __(ra, dec) = (0.1843687374, -0.68820842)__ and we will be using these values for creating NR frame files.
   * {{:hl_resp.png?500|}}  HL detector network response
   * {{:hlv_resp.png?500|}} HLV detector network response
+=== Creating frame files ===
+  * First, create a h5 file from the psi4 data for the different simulations. The script to create the h5 file is at : __/home/c1549390/work/PhenomHM_pe/src/nrutils-make-nrh5-BAM.py__ for the BAM waveforms, and  __/home/c1549390/work/PhenomHM_pe/src/nrutils-make-nrh5-GaTech.py__ for the Georgia tech waveforms. All modes upto and including l=4 are stored in the h5 file.
+  * As of 18/09/2017; I have created h5 files for q1, q2, q4, q8 and q18 non-spinning NR waveforms. q1 and q2 are Gerogia Tech waveforms while q4, q8 and q18 are BAM.
+  * I wanted to keep the f_low (as calculated from the NR waveform) of the waveform for all the systems to be the same, so I adjust the mass of the injection to adjust the f_low ~ 20Hz for all runs.
+  * Created HL and HVL (inclination series) frame files for the q4 and q8 non-spinning systems. To ensure f_low is same, q4 is injected at 80 Msun and q8 at 115 Msun.
+  * The two detector network injection (H-L) frames can be found at : __/home/c1549390/work/PhenomHM_pe/frames/HL_Frames__
+  * The three detector network injection (H-L-V) frames can be found at: __/home/c1549390/work/PhenomHM_pe/frames/HLV_Frames__
+=== PE Runs ===
+  * In some previous runs, we had noticed that the improvement of PhenomHM vs PhenomD to constrain the distance was increased by ~ 50 % (i.e., the posteriors of PhenomHM were narrower by 50% compared to PhenomD ) for the two detector network and for the HLV configuration, the same value was ~ 40%.
+  * One caveat between those two runs was that for the HL run, the source was placed at a sky-position where the response was the maximum, but for the HLV run the source was at a sky-position where the response was about half of that.
+  * I wanted to check the effect of changing the inclination of the source on the recovery and to check for any biases on the posteriors and to check whether including a third detector would improve the results.
+  * With that in mind, I have started a set of runs for the BAM q4 & q8 non-spinning waveforms for two and three detector configuration, at three different inclinations for each. So it's a total of 24 runs.
+  * I've used the aLIGOzeroDetHighPower and aVIRGO PSD's for the runs.
+  * The base directory for all the run directories is at: __/home/c1549390/work/PhenomHM_pe/runs/18_09_2017/__
+^  Simulation h5 file  ^  Detector configuration  ^  Inclination  ^  Waveform  ^ Date started  ^  Run Directory : basedir + ^  Status  ^  Results page  ^
+|  q4a0_T_112_448_fix_scale_Lmax_4_Allm.h5 | HL  |  0 | IMRPhenomD | 18/09/17 | BAM_q4a0a0_SNR_25_mode_lMax_4_inc_0.0_M80_HL_PhenomD | Running | |
+|  q4a0_T_112_448_fix_scale_Lmax_4_Allm.h5 | HL  |  60 | IMRPhenomD | 18/09/17 | BAM_q4a0a0_SNR_25_mode_lMax_4_inc_60.0_M80_HL_PhenomD | Running | |
+|  q4a0_T_112_448_fix_scale_Lmax_4_Allm.h5 | HL  |  90 | IMRPhenomD | 18/09/17 | BAM_q4a0a0_SNR_25_mode_lMax_4_inc_90.0_M80_HL_PhenomD | Running | |
+|  q4a0_T_112_448_fix_scale_Lmax_4_Allm.h5 | HL  |  0 | IMRPhenomHM | 18/09/17 | BAM_q4a0a0_SNR_25_mode_lMax_4_inc_0.0_M80_HL_PhenomHM | Running | |
+|  q4a0_T_112_448_fix_scale_Lmax_4_Allm.h5 | HL  |  60 | IMRPhenomHM | 18/09/17 | BAM_q4a0a0_SNR_25_mode_lMax_4_inc_60.0_M80_HL_PhenomHM | Running | |
+|  q4a0_T_112_448_fix_scale_Lmax_4_Allm.h5 | HL  |  90 | IMRPhenomHM | 18/09/17 | BAM_q4a0a0_SNR_25_mode_lMax_4_inc_90.0_M80_HL_PhenomHM | Running | |
+|  q4a0_T_112_448_fix_scale_Lmax_4_Allm.h5 | HLV  |  0 | IMRPhenomD | 18/09/17 | BAM_q4a0a0_SNR_25_mode_lMax_4_inc_0.0_M80_HLV_PhenomD | Running | |
+|  q4a0_T_112_448_fix_scale_Lmax_4_Allm.h5 | HLV  |  60 | IMRPhenomD | 18/09/17 | BAM_q4a0a0_SNR_25_mode_lMax_4_inc_60.0_M80_HLV_PhenomD | Running | |
+|  q4a0_T_112_448_fix_scale_Lmax_4_Allm.h5 | HLV  |  90 | IMRPhenomD | 18/09/17 | BAM_q4a0a0_SNR_25_mode_lMax_4_inc_90.0_M80_HLV_PhenomD | Running | |
+|  q4a0_T_112_448_fix_scale_Lmax_4_Allm.h5 | HLV  |  0 | IMRPhenomHM | 18/09/17 | BAM_q4a0a0_SNR_25_mode_lMax_4_inc_0.0_M80_HLV_PhenomHM | Running | |
+|  q4a0_T_112_448_fix_scale_Lmax_4_Allm.h5 | HLV  |  60 | IMRPhenomHM | 18/09/17 | BAM_q4a0a0_SNR_25_mode_lMax_4_inc_60.0_M80_HLV_PhenomHM | Running | |
+|  q4a0_T_112_448_fix_scale_Lmax_4_Allm.h5 | HLV  |  90 | IMRPhenomHM | 18/09/17 | BAM_q4a0a0_SNR_25_mode_lMax_4_inc_90.0_M80_HLV_PhenomHM | Running | |
+|  q8a0a0_T_112_588n768_fix_scale_Lmax_4_Allm.h5 | HL  |  0 | IMRPhenomD | 18/09/17 | BAM_q8a0a0_SNR_25_mode_lMax_4_inc_0.0_M115_HL_PhenomD | Running | |
+|  q8a0a0_T_112_588n768_fix_scale_Lmax_4_Allm.h5 | HL  |  60 | IMRPhenomD | 18/09/17 | BAM_q8a0a0_SNR_25_mode_lMax_4_inc_60.0_M115_HL_PhenomD | Running | |
+|  q8a0a0_T_112_588n768_fix_scale_Lmax_4_Allm.h5 | HL  |  90 | IMRPhenomD | 18/09/17 | BAM_q8a0a0_SNR_25_mode_lMax_4_inc_90.0_M115_HL_PhenomD | Running | |
+|  q8a0a0_T_112_588n768_fix_scale_Lmax_4_Allm.h5 | HL  |  0 | IMRPhenomHM | 18/09/17 | BAM_q8a0a0_SNR_25_mode_lMax_4_inc_0.0_M115_HL_PhenomHM | Running | |
+|  q8a0a0_T_112_588n768_fix_scale_Lmax_4_Allm.h5 | HL  |  60 | IMRPhenomHM | 18/09/17 | BAM_q8a0a0_SNR_25_mode_lMax_4_inc_60.0_M115_HL_PhenomHM | Running | |
+|  q8a0a0_T_112_588n768_fix_scale_Lmax_4_Allm.h5 | HL  |  90 | IMRPhenomHM | 18/09/17 | BAM_q8a0a0_SNR_25_mode_lMax_4_inc_90.0_M115_HL_PhenomHM | Running | |
+|  q8a0a0_T_112_588n768_fix_scale_Lmax_4_Allm.h5 | HLV  |  0 | IMRPhenomD | 18/09/17 | BAM_q8a0a0_SNR_25_mode_lMax_4_inc_0.0_M115_HLV_PhenomD | Running | |
+|  q8a0a0_T_112_588n768_fix_scale_Lmax_4_Allm.h5 | HLV  |  60 | IMRPhenomD | 18/09/17 | BAM_q8a0a0_SNR_25_mode_lMax_4_inc_60.0_M115_HLV_PhenomD | Running | |
+|  q8a0a0_T_112_588n768_fix_scale_Lmax_4_Allm.h5 | HLV  |  90 | IMRPhenomD | 18/09/17 | BAM_q8a0a0_SNR_25_mode_lMax_4_inc_90.0_M115_HLV_PhenomD | Running | |
+|  q8a0a0_T_112_588n768_fix_scale_Lmax_4_Allm.h5 | HLV  |  0 | IMRPhenomHM | 18/09/17 | BAM_q8a0a0_SNR_25_mode_lMax_4_inc_0.0_M115_HLV_PhenomHM | Running | |
+|  q8a0a0_T_112_588n768_fix_scale_Lmax_4_Allm.h5 | HLV  |  60 | IMRPhenomHM | 18/09/17 | BAM_q8a0a0_SNR_25_mode_lMax_4_inc_60.0_M115_HLV_PhenomHM | Running | |
+|  q8a0a0_T_112_588n768_fix_scale_Lmax_4_Allm.h5 | HLV  |  90 | IMRPhenomHM | 18/09/17 | BAM_q8a0a0_SNR_25_mode_lMax_4_inc_90.0_M115_HLV_PhenomHM | Running | |