#################### # [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.4227467 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.1 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)] 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 : 1638 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.422747 * 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 589/40968 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 766/40968 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)]`. Mf : [20.0, 200.0] af : [0.0, 0.99] ------------------------- - [Intrinsic: Testing-GR - ------------------------- * No TGR parameter was considered. -------------- - [Amplitudes - -------------- A2220 : [0.0, 50.0] phi2220 : [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: 148734.574 * 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): 6.932037115097046 ################## # [Start sampling # ################## * Using CPNest version: `0.11.3`. * The sampling output appears in the `result/Nested_sampler/cpnest.log` file. Signal evidence : 148803.934 ln B : 69.360 Estimated SNR from logB : 12.879 * Sampling execution time (min): 27925526.633890197 ########################### # [Posterior railing tests # ########################### Creating pesummary metafile... * Warning: Information storage in pesumarry failed with error: None. ################### # [SNR computation # ################### Medians: Network matched_filter SNR (TD): 12.780 H1 matched_filter SNR (TD): 9.963 L1 matched_filter SNR (TD): 8.013 Network optimal SNR (TD): 12.748 H1 optimal SNR (TD): 10.279 L1 optimal SNR (TD): 7.541 #################### # [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 16507/16507 normality outliers (at 1.0 % significance). * A Kolmogorov–Smirnov test of whitened L1 residuals gave 16507/16507 normality outliers (at 1.0 % significance). No IMR posterior was passed. * Total execution time (min): 106.45562679370245