#################### # [Input parameters # #################### * Reading config file : `config_gw150914_production.ini`. * With sections : ['input', 'Sampler settings', 'Priors', 'Injection', 'Plot']. * I'll be running with the following values: data-H1 : data-L1 : data-V1 : ignore-data-filename : 0 download-data : 1 datalen-download : 4096.0 gw-data-find : 0 gw-data-type-H1 : gw-data-type-L1 : gw-data-type-V1 : tag : CLN channel-H1 : GWOSC channel-L1 : GWOSC channel-V1 : GWOSC config-file : config_gw150914_production.ini run-type : full output : result run-tag : PROD1 screen-output : 0 debug : 0 pesummary : 1 trigtime : 1126259462.423479 detectors : H1,L1 ref-det : H1 sky-frame : equatorial acf-H1 : acf-L1 : acf-V1 : psd-H1 : psd-L1 : psd-V1 : signal-chunksize : 2.0 noise-chunksize : 2.0 window-onsource : 0 window : 1 alpha-window : 0.1 sampling-rate : 16384.0 f-min-bp : 20.0 f-max-bp : 8182.0 bandpassing : 1 noise-averaging-method : mean fft-acf : 1 acf-simple-norm : 1 no-lognorm : 0 truncate : 1 analysis-duration : 0.2 analysis-duration-n : 410 zero-noise : 0 gaussian-noise : gaussian-noise-seed : -1 gaussian-noise-white-sigma : 1e-21 chisquare-computation : 0 non-stationarity-check : 0 maxent-psd : PSD-investigation : 0 onsource-ACF : 0 Dirac-comb : 0 Zeroing-data : 0 likelihood-method : direct-inversion split-inner-products : 0 injection-parameters : None injection-approximant : inject-n-ds-modes : {'s': 0, 'v': 0, 't': 1} inject-area-quantization : 0 inject-braneworld : 0 inject-charge : 0 injection-scaling : 1.0 injection-T : 64.0 template : Kerr single-mode : None n-ds-modes : {'s': 0, 'v': 0, 't': 1} ds-ordering : freq kerr-modes : [(2, 2, 2, 0), (2, 2, 2, 1)] reference-amplitude : 1e-21 spheroidal : 0 qnm-fit : 1 coherent-n : 0 amp-non-prec-sym : 1 max-Kerr-amp-ratio : 0.0 TGR-overtones-ordering : Unordered domega-tgr-modes : None dtau-tgr-modes : None area-quantization : 0 braneworld : 0 tau-AQ : 0 prior-reweight : 0 ParSpec : 0 ParSpec_Dmax_TGR : 2 ParSpec_Dmax_charge : 0 EsGB : 0 charge : 0 gr-time-prior : 1 dist-flat-prior : 0 ds-amp-flat-prior : 0 mf-time-prior : 0.0 analysis-duration-n : 3276 inject-n-ds-modes : {'t': 1} n-ds-modes : {'t': 1} gr-time-prior : 1 dist-flat-prior : 0 ds-amp-flat-prior : 0 mf-time-prior : 67.92493161247017 ########################## # [Setting the start time # ########################## Trigtime in H1 : 1126259462.423479 * Warning: To err on the conservative side, when fixing the start time, pyRing selects the discrete sample immediately after the start time requested in the likelihood. This, added to the data discretisation, implies that the actual start time used in the analysis can be larger than the one requested up to 0.000092 s, i.e. 1.5 * 1/(sampling_rate). If you wish to reduce this error, please increase the sampling rate and the bandpassing frequency accordingly. ################ # [Data reading # ################ ------ - [H1 - ------ * Using GWPY to download data. Fetched 1 URLs from www.gw-openscience.org for [1126257414 .. 1126261510)) Reading data... [Done] * Loaded channel GWOSC starting at 1126257414.0 length 4096.0s. * Bandpassing the raw strain between [20.0, 8182.0] Hz. * Computing the one-sided PSD with the Welch method for comparison with the standard ACF. * No ACF was passed. Estimating ACF. * A Kolmogorov–Smirnov test of whitened data chunks gave 245/20483 normality outliers (at 1.0 % significance). ------ - [L1 - ------ * Using GWPY to download data. Fetched 1 URLs from www.gw-openscience.org for [1126257414 .. 1126261510)) Reading data... [Done] * Loaded channel GWOSC starting at 1126257414.0 length 4096.0s. * Bandpassing the raw strain between [20.0, 8182.0] Hz. * Computing the one-sided PSD with the Welch method for comparison with the standard ACF. * No ACF was passed. Estimating ACF. * A Kolmogorov–Smirnov test of whitened data chunks gave 414/20483 normality outliers (at 1.0 % significance). ######################## # [Model initialisation # ######################## * Using equatorial sky coordinates. ########## # [Priors # ########## I'll be running with the following prior bounds: ------------- - [Extrinsic - ------------- ------------------------- - [Intrinsic: Kerr model - ------------------------- * Running the Kerr model with modes: `[(2, 2, 2, 0), (2, 2, 2, 1)]`. Mf : [20.0, 200.0] af : [0.0, 0.99] ------------------------- - [Intrinsic: Testing-GR - ------------------------- * No TGR parameter was considered. -------------- - [Amplitudes - -------------- A2220 : [0.0, 50.0] A2221 : [0.0, 50.0] phi2220 : [0.0, 6.283185307179586] phi2221 : [0.0, 6.283185307179586] --------- - [Fixed - --------- ra : 1.95 dec : -1.27 psi : 0.82 t : 0.0 cosiota : -1.0 phi : 0.0 ########################## # [Sampler initialisation # ########################## * I'll be running with the following sampler parameters: verbose : 2 poolsize : 128 nthreads : 1 nlive : 4096 maxmcmc : 4096 seed : 1000 resume : 0 periodic_checkpoint_interval : 900 Noise evidence: 297864.014 * Whitened plots using the FD PSD are illustrative, and do not correspond to the full TD treatment applied in the analysis. * Pre-sampling execution time (min): 5.750338443120321 ################## # [Start sampling # ################## * Using CPNest version: `0.11.3`. * The sampling output appears in the `result/Nested_sampler/cpnest.log` file. Signal evidence : 297913.342 ln B : 49.328 Estimated SNR from logB : 11.560 * Sampling execution time (min): 27938419.384084605 ########################### # [Posterior railing tests # ########################### af is railing against the lower prior bound. af is railing against the upper prior bound. A2220 is railing against the lower prior bound. A2221 is railing against the lower prior bound. phi2220 is railing against the lower prior bound. phi2220 is railing against the upper prior bound. phi2221 is railing against the lower prior bound. phi2221 is railing against the upper prior bound. Creating pesummary metafile... * Warning: Information storage in pesumarry failed with error: None. ################### # [SNR computation # ################### Medians: Network matched_filter SNR (TD): 11.395 H1 matched_filter SNR (TD): 8.406 L1 matched_filter SNR (TD): 7.699 Network optimal SNR (TD): 11.362 H1 optimal SNR (TD): 9.251 L1 optimal SNR (TD): 6.593 #################### # [Plots generation # #################### * Skipping sky position plot, since at least one of the sky position parameters was fixed. * Skipping time plot, since time was fixed. * A Kolmogorov–Smirnov test of whitened H1 residuals gave 20829/20829 normality outliers (at 1.0 % significance). * A Kolmogorov–Smirnov test of whitened L1 residuals gave 20829/20829 normality outliers (at 1.0 % significance). No IMR posterior was passed. * Total execution time (min): 5422.679279116789