def fig_prs_ratio():
kinds = ['conv_nue', 'pr_nue', 'conv_numu', 'pr_numu']
cos_ths = [0.25, 0.85]
labels = [r'ALLM97', r'BB']
for kind in kinds:
plt.figure()
plt.title(titling[kind])
for cos_th in cos_ths:
clabel = r'$\cos \theta_z = {}$'.format(cos_th)
allm97 = plots.fn(cos_th)(cos_th,
kind,
prpl='ice_allm97_step_1',
label=clabel)
bb = plots.fn(cos_th)(cos_th,
kind,
prpl='ice_bb_step_1',
linestyle='--',
color=allm97.get_color())
plt.plot(bb.get_data()[0],
allm97.get_data()[1] / bb.get_data()[1],
':',
color=allm97.get_color())
plt.axvline(np.nan, color='k', linestyle='-', label=labels[0])
plt.axvline(np.nan, color='k', linestyle='--', label=labels[1])
plt.axvline(np.nan, color='k', linestyle=':', label='ALLM97/BB')
plt.legend()
plt.tight_layout(0.3)
save('fig/prs_ratio_{}.eps'.format(kind))
def fig_pmodels():
kinds = ['conv_nue', 'pr_nue', 'conv_numu', 'pr_numu']
pmodels = [
(pm.HillasGaisser2012, 'H3a', 'H3a'),
# (pm.GlobalSplineFitBeta, None, 'GSF spl'),
(pm.GaisserStanevTilav, '4-gen', 'GST 4-gen'),
# (pm.PolyGonato, False, 'poly-gonato'),
# (pm.GaisserHonda, None, 'GH'),
(pm.ZatsepinSokolskaya, 'default', 'ZS')
]
cos_ths = [0.25, 0.85]
for kind in kinds:
plt.figure()
plt.title(titling[kind])
for idx, pmodel in enumerate(pmodels):
for cos_th in cos_ths:
clabel = r'$\cos \theta_z = {}$'.format(
cos_th) if idx == 0 else None
plots.fn(cos_th)(cos_th,
kind,
pmodel=pmodel[:-1],
label=clabel,
linestyle=linestyles[idx])
plt.axvline(np.nan,
color='k',
linestyle=linestyles[idx],
label=pmodel[-1])
plt.gca().set_prop_cycle(None)
plt.legend()
plt.tight_layout(0.3)
save('fig/pmodels_{}.eps'.format(kind, cos_th))
def fig_extsv():
kinds = ['conv_nue', 'pr_nue', 'conv_numu', 'pr_numu']
cos_ths = [0.25, 0.85]
ens = np.logspace(2, 9, 100)
useexts = [False, True]
labels = ['This work', 'GJKvS']
for kind in kinds:
plt.figure()
plt.title(titling[kind])
for idx, useext in enumerate(useexts):
for cos_th in cos_ths:
if useext:
emu = extsv.minimum_muon_energy(extsv.overburden(cos_th))
plt.plot(ens,
exthp.passrates(kind)(ens, emu, cos_th),
linestyles[idx])
else:
clabel = r'$\cos \theta_z = {}$'.format(
cos_th) if idx == 0 else None
plots.fn(cos_th)(cos_th, kind, label=clabel)
plt.axvline(np.nan,
color='k',
linestyle=linestyles[idx],
label=labels[idx])
plt.gca().set_prop_cycle(None)
plt.legend()
plt.tight_layout(0.3)
save('fig/extsv_{}.eps'.format(kind))
def fig_density():
kinds = ['conv_nue', 'pr_nue', 'conv_numu', 'pr_numu']
dmodels = [('CORSIKA', ('SouthPole', 'June'), 'MSIS-90-E SP/Jul'),
('CORSIKA', ('SouthPole', 'December'), 'MSIS-90-E SP/Dec'),
('MSIS00', ('Karlsruhe', 'July'), 'NRLMSISE-00 KR/Jul'),
('MSIS00', ('Karlsruhe', 'December'), 'NRLMSISE-00 KR/Dec')]
cos_ths = [0.25, 0.85]
for kind in kinds:
plt.figure()
plt.title(titling[kind])
for idx, dmodel in enumerate(dmodels):
for cos_th in cos_ths:
clabel = r'$\cos \theta_z = {}$'.format(
cos_th) if idx == 0 else None
plots.fn(cos_th)(cos_th,
kind,
density=dmodel[:2],
label=clabel,
linestyle=linestyles[idx])
plt.axvline(np.nan,
color='k',
linestyle=linestyles[idx],
label='{}'.format(dmodel[2]))
plt.gca().set_prop_cycle(None)
plt.legend()
plt.tight_layout(0.3)
save('fig/dmodels_{}.eps'.format(kind, cos_th))
def fig_hadrs():
kinds = ['conv_nue', 'pr_nue', 'conv_numu', 'pr_numu']
hadrs_prompt = ['SIBYLL2.3c', 'SIBYLL2.3', 'DPMJET-III']
hadrs_conv = ['SIBYLL2.3c', 'SIBYLL2.3', 'QGSJET-II-04', 'EPOS-LHC']
cos_ths = [0.25, 0.85]
for kind in kinds:
plt.figure()
plt.title(titling[kind])
hadrs = hadrs_conv if kind.split('_')[0] == 'conv' else hadrs_prompt
for idx, hadr in enumerate(hadrs):
for cos_th in cos_ths:
clabel = r'$\cos \theta_z = {}$'.format(
cos_th) if idx == 0 else None
plots.fn(cos_th)(cos_th,
kind,
hadr=hadr,
label=clabel,
linestyle=linestyles[idx])
plt.axvline(np.nan,
color='k',
linestyle=linestyles[idx],
label=hadr)
plt.gca().set_prop_cycle(None)
plt.legend()
plt.tight_layout(0.3)
save('fig/hadrs_{}.eps'.format(kind))
def fig_medium():
kinds = ['conv_nue', 'pr_nue', 'conv_numu', 'pr_numu']
cos_ths = [0.25, 0.85]
prpls = [('ice_allm97_step_1', 1.95 * Units.km),
('water_allm97_step_1', 1.95 * Units.km),
('ice_allm97_step_1', 3.5 * Units.km),
('water_allm97_step_1', 3.5 * Units.km)]
labels = [
r'Ice $1.95$ km', r'Water $1.95$ km', r'Ice $3.5$ km',
r'Water $3.5$ km'
]
for kind in kinds:
plt.figure()
plt.title(titling[kind])
for idx, prpl in enumerate(prpls):
for cos_th in cos_ths:
clabel = r'$\cos \theta_z = {}$'.format(
cos_th) if idx == 0 else None
plots.fn(cos_th)(cos_th,
kind,
prpl=prpl[0],
depth=prpl[1],
label=clabel,
linestyle=linestyles[idx])
plt.axvline(np.nan,
color='k',
linestyle=linestyles[idx],
label=labels[idx])
plt.gca().set_prop_cycle(None)
plt.legend()
plt.tight_layout(0.3)
save('fig/medium_{}.eps'.format(kind))
def fig_pls():
kinds = ['conv_nue', 'pr_nue', 'conv_numu', 'pr_numu']
cos_ths = [0.25, 0.85]
prpls = [
'ice_allm97_step_1', 'ice_allm97_step_0.75',
'ice_allm97_sigmoid_0.75_0.25'
]
labels = [
r'$\Theta(E_\mu^{\rm f} - 1\,{\rm TeV})$',
r'$\Theta(E_\mu^{\rm f} - 0.75\,{\rm TeV})$',
r'$\Phi\left(\frac{E_\mu^{\rm f} - 0.75\,{\rm TeV}}{0.25\,{\rm TeV}}\right)$'
]
for kind in kinds:
plt.figure()
plt.title(titling[kind])
for idx, prpl in enumerate(prpls):
for cos_th in cos_ths:
clabel = r'$\cos \theta_z = {}$'.format(
cos_th) if idx == 0 else None
plots.fn(cos_th)(cos_th,
kind,
prpl=prpl,
label=clabel,
linestyle=linestyles[idx])
plt.axvline(np.nan,
color='k',
linestyle=linestyles[idx],
label=labels[idx])
plt.gca().set_prop_cycle(None)
plt.legend()
plt.tight_layout(0.3)
save('fig/pls_{}.eps'.format(kind))
def fig_prs():
kinds = ['conv_nue', 'pr_nue', 'conv_numu', 'pr_numu']
cos_ths = [0.25, 0.85]
prpls = ['ice_allm97_step_1', 'ice_bb_step_1']
labels = [r'ALLM97', r'BB']
for kind in kinds:
plt.figure()
plt.title(titling[kind])
for idx, prpl in enumerate(prpls):
for cos_th in cos_ths:
clabel = r'$\cos \theta_z = {}$'.format(
cos_th) if idx == 0 else None
plots.fn(cos_th)(cos_th,
kind,
prpl=prpl,
label=clabel,
linestyle=linestyles[idx])
plt.axvline(np.nan,
color='k',
linestyle=linestyles[idx],
label=labels[idx])
plt.gca().set_prop_cycle(None)
plt.legend()
plt.tight_layout(0.3)
save('fig/prs_{}.eps'.format(kind))
