Event Rates with 3G detectors

[1]:

import numpy as np
from ler.rates import LeR

[2]:

ler = LeR(batch_size=50000, ifos=['CE', 'ET'], verbose=False)

[3]:

ler.unlensed_cbc_statistics(size=1000000, save_batch=False, resume=True, output_jsonfile='unlensed_3G_1M.json');

unlensed params will be store in unlensed_3G_1M.json
chosen batch size = 50000 with total size = 1000000
There will be 20 batche(s)
resuming from ./ler_data/unlensed_3G_1M.json
saving all unlensed_params in ./ler_data/unlensed_3G_1M.json...

[4]:

ler.unlensed_rate();

getting unlensed_params from json file unlensed_3G_1M.json...
total unlensed rate (yr^-1) (with step function): 91843.65042655931
number of simulated unlensed detectable events: 887182
number of all simulated unlensed events: 1000000
storing detectable unlensed params in unlensed_param_detectable.json

[5]:

ler.unlensed_rate(snr_threshold=6);

getting unlensed_params from json file unlensed_3G_1M.json...
total unlensed rate (yr^-1) (with step function): 94258.63255779142
number of simulated unlensed detectable events: 910510
number of all simulated unlensed events: 1000000
storing detectable unlensed params in unlensed_param_detectable.json

[7]:

ler.lensed_cbc_statistics(size=1000000, save_batch=False, resume=True, output_jsonfile='lensed_3G_1M.json');

lensed params will be store in lensed_3G_1M.json
chosen batch size = 50000 with total size = 1000000
There will be 20 batche(s)
resuming from ./ler_data/lensed_3G_1M.json
saving all lensed_params in ./ler_data/lensed_3G_1M.json...

[8]:

ler.lensed_rate();

getting lensed_params from json file lensed_3G_1M.json...
total lensed rate (yr^-1) (with step function): 152.3718678085319
number of simulated lensed detectable events: 613339
number of simulated all lensed events: 1000000
storing detectable lensed params in lensed_param_detectable.json

[9]:

ler.lensed_rate(snr_threshold=[6,6], num_img=[1,1]);

getting lensed_params from json file lensed_3G_1M.json...
total lensed rate (yr^-1) (with step function): 167.8888117226168
number of simulated lensed detectable events: 675799
number of simulated all lensed events: 1000000
storing detectable lensed params in lensed_param_detectable.json

[8]:

ler.json_file_names

[8]:

{'ler_param': './ler_params.json',
 'unlensed_param': 'unlensed_3G_1M.json',
 'unlensed_param_detectable': './unlensed_param_detectable.json',
 'lensed_param': 'lensed_3G_1M.json',
 'lensed_param_detectable': './lensed_param_detectable.json'}

[11]:

ler.rate_ratio();

unlensed_rate: 91843.65042655931
lensed_rate: 152.3718678085319
ratio: 602.7598909660187

With ANN snr (IMRphenomXPHM)

(under development)

[5]:

from ler.rates import LeR
ler = LeR(
    batch_size=50000,
    verbose=False,
    spin_zero=False,
    spin_precession=True,
    # ifos=['CE', 'ET'], # This will be implemented in the next version
    snr_type='ann',
    waveform_approximant='IMRPhenomXPHM'
    )

[10]:

ler.lensed_cbc_statistics(size=100000, save_batch=False);

lensed params will be store in ./lensed_param.json
chosen batch size = 50000 with total size = 100000
There will be 2 batche(s)
Batch no. 1
sampling lensed params...
solving lens equations...

100%|███████████████████████████████████████████████████████| 50000/50000 [00:14<00:00, 3389.90it/s]

calculating snrs...
Batch no. 2
sampling lensed params...
solving lens equations...

100%|███████████████████████████████████████████████████████| 50000/50000 [00:14<00:00, 3358.70it/s]

calculating snrs...

[2]:

ler.lensed_rate();

getting lensed_params from json file ./lensed_param.json...
total lensed rate (yr^-1) (with step function): 1.488096286349158
number of simulated lensed detectable events: 599
number of simulated all lensed events: 100000
storing detectable lensed params in ./lensed_param_detectable.json

[3]:

ler.lensed_rate(snr_recalculation=True,);

getting lensed_params from json file ./lensed_param.json...
calculating snrs...

100%|████████████████████████████████████████████████████████████| 840/840 [00:04<00:00, 177.26it/s]
100%|████████████████████████████████████████████████████████████| 840/840 [00:04<00:00, 201.31it/s]
100%|████████████████████████████████████████████████████████████| 729/729 [00:03<00:00, 193.07it/s]
100%|████████████████████████████████████████████████████████████| 570/570 [00:03<00:00, 160.27it/s]

total lensed rate (yr^-1) (with step function): 1.1800429649680302
number of simulated lensed detectable events: 475
number of simulated all lensed events: 100000
storing detectable lensed params in ./lensed_param_detectable.json