lum_index=0,
n_model='vol_co',
pulse_model='uniform',
pulse_range=[1., 1.],
pulse_mu=1.,
pulse_sigma=0.,
repeat=0.,
si_mu=si,
si_sigma=0.,
z_max=2.5)
pop[si].save()
else:
pop[si] = copy.deepcopy(pop[min(SIS)])
pop[si].frbs.si = np.random.normal(si, 0, n_tot)
pop[si].name = f'si-{si}'
pop[si].save()
pop_obs = {}
if OBSERVE or CREATE:
for si in SIS:
if not CREATE:
pop[si] = unpickle(f'si-{si}')
# Create Survey
perfect = Survey('perfect', gain_pattern='perfect')
# Observe populations
n_model='sfr',
pulse_model='uniform',
pulse_range=[1., 1.],
pulse_mu=1.,
pulse_sigma=0.,
repeat=0.,
si_mu=0.,
si_sigma=0.,
z_max=2.5)
pop[n].save()
else:
pop[n] = copy.deepcopy(pop[35])
n_gen = len(pop[35].frbs.lum_bol)
pop[n].frbs.lum_bol = dis.powerlaw(10**n, 10**n, 0, n_gen)
pop[n].name = f'sc-{n}'
pop[n].save()
pop_obs = {}
if OBSERVE:
for n in (35, 36, 37):
if not CREATE:
pop[n] = unpickle(f'sc-{n}')
# Create Survey
perfect = Survey('perfect-small',
gain_pattern='gaussian',
n_sidelobes=8)
si_sigma=0.,
n_model='vol_co',
days=days,
dm_host_model='normal',
dm_host_mu=0,
dm_host_sigma=0,
dm_igm_index=0,
dm_igm_sigma=0,
dm_mw_model='zero',
emission_range=[10e6, 10e9],
lum_index=0,
pulse_model='uniform',
pulse_range=[1., 1.],
pulse_mu=1.,
pulse_sigma=0.)
population.name = 'test'
# Setup a survey
survey = Survey('perfect', gain_pattern='perfect')
# Observe the FRB population
surv_pop = SurveyPopulation(population, survey)
surv_pop.name = 'lognlogslocal'
surv_pop.save()
else:
surv_pop = unpickle('lognlogslocal')
parms = surv_pop.frbs.fluence
min_p = min(parms)
max_p = max(parms)