Skip to content

Commit

Permalink
Small update so Github isnt annoying about image file
Browse files Browse the repository at this point in the history
  • Loading branch information
vvelan committed Sep 20, 2024
1 parent 423c58a commit 98b3609
Showing 1 changed file with 16 additions and 11 deletions.
27 changes: 16 additions & 11 deletions examples/plot_sapphire_NR.py
Original file line number Diff line number Diff line change
Expand Up @@ -9,8 +9,8 @@
mpl.rcParams.update({'font.size': 20})
mpl.rcParams.update({'axes.linewidth': 2})
fig, ax = plt.subplots(1, 1, figsize=(11, 11))
xmin = 1e1; xmax = 1e6;
ymin = 1e-48; ymax = 1e-28;
xmin = 1e-2; xmax = 1e7;
ymin = 1e-48; ymax = 1e-20;

# Existing limits
m_limit, x_limit = np.loadtxt(curves_dir + 'SENSEI_Migdal_NR.txt').transpose()
Expand All @@ -23,15 +23,20 @@
ax.plot(m_limit*1e3, x_limit, '--', lw=1.5, label='LZ')

# Simulation
m_limit, x_limit = np.loadtxt('sapphire/results_sapphire_oi_phonon_massless_001_days/HeRALD_FC_1d_1device_1fold_100mus.txt').transpose()
ax.plot(m_limit*1e3, x_limit, 'b--', lw=4, label='1 day, multiphonon')
m_limit, x_limit = np.loadtxt('sapphire/results_sapphire_oi_phonon_massless_010_days/HeRALD_FC_10d_1device_1fold_100mus.txt').transpose()
ax.plot(m_limit*1e3, x_limit, 'r--', lw=4, label='10 days, multiphonon')
# Also massive
m_limit, x_limit = np.loadtxt('sapphire/results_sapphire_oi_scan_NR_001days/HeRALD_FC_1d_1device_1fold_100mus.txt').transpose()
ax.plot(m_limit*1e3, x_limit, 'b:', lw=4, alpha=0.5, label='1 day, NR')
m_limit, x_limit = np.loadtxt('sapphire/results_sapphire_oi_scan_NR_010days/HeRALD_FC_10d_1device_1fold_100mus.txt').transpose()
ax.plot(m_limit*1e3, x_limit, 'r:', lw=4, alpha=0.5, label='1 day, NR')
m_limit, x_limit = np.loadtxt('results/sapphire_oi_NR_1h_low_bkgd/HeRALD_FC_0d_1device_1fold_100mus.txt').transpose()
ax.plot(m_limit*1e3, x_limit, 'k--', lw=3, label='NR Low Bkgd')

m_limit, x_limit = np.loadtxt('results/sapphire_oi_phonon_massive_1h_low_bkgd/HeRALD_FC_0d_1device_1fold_100mus.txt').transpose()
ax.plot(m_limit*1e3, x_limit, 'g--', lw=3, label='Multiphonon Low Bkgd')

m_limit, x_limit = np.loadtxt('results/sapphire_oi_NR_1h/HeRALD_FC_0d_1device_1fold_100mus.txt').transpose()
ax.plot(m_limit*1e3, x_limit, 'k-', lw=3, label='NR')

m_limit, x_limit = np.loadtxt('results/sapphire_oi_phonon_massive_1h/HeRALD_FC_0d_1device_1fold_100mus.txt').transpose()
ax.plot(m_limit*1e3, x_limit, 'g-', lw=3, label='Multiphonon')

m_limit, x_limit = np.loadtxt('results/limit.txt').transpose()
ax.plot(m_limit*1e3, x_limit, 'r-', lw=3, label='Test')

ax.set_xscale('log')
ax.set_xlim([xmin, xmax])
Expand Down

0 comments on commit 98b3609

Please sign in to comment.