Cardiff Gravity Exploration Institute

This page is read only. You can view the source, but not change it. Ask your administrator if you think this is wrong.
==== IMRPhenomHM PE Study : Aug 2017 onwards ====

This wiki page is created to note all the motivations for performing specific runs, links to the result pages and codes. The main motivation behind this study is to assess the reliability of IMRPhenomHM waveform to recover physical parameters from injected NR waveforms in zero-noise. We inject NR waveforms in zero-noise so as to suppress all errors from noise contributions to the data. 

=== Choosing the sky position === 

  * 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 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 | |