diff options
Diffstat (limited to 'utils/plot.py')
| -rw-r--r-- | utils/plot.py | 24 |
1 files changed, 20 insertions, 4 deletions
diff --git a/utils/plot.py b/utils/plot.py index f05b813..4e1a58b 100644 --- a/utils/plot.py +++ b/utils/plot.py @@ -408,6 +408,15 @@ def plot_sens_fixed_angle_pretty(data, outfile, outformat, args): 8: (1.4E-45, 1.00E-34) } + best_limits = { + 3: ('CMB polarization', 1E-43, 'brown'), + 4: ('GRB vacuum birefringence', 1E-38, 'lightgreen'), + 5: ('GRB vacuum birefringence', 1E-34, 'lightgreen'), + 6: ('UHE cosmic ray', 1E-42, 'violet'), + 7: ('GRB vacuum birefringence', 1E-28, 'lightgreen'), + 8: ('gravitational Cherenkov radiation', 1E-46, 'lightblue') + } + show_data = True en = np.log10([1E4, 1E7]) @@ -448,10 +457,6 @@ def plot_sens_fixed_angle_pretty(data, outfile, outformat, args): ax.set_xlim(xlims) ylim = (0.5, 1.5) ax.set_ylim(ylim) - # ax.set_yticks([ylim[0], 1., ylim[1]]) - # ax.set_yticklabels([''] + [yticks[ian]] + [''], fontsize=13) - # ax.yaxis.get_major_ticks()[0].set_visible(False) - # ax.yaxis.get_major_ticks()[2].set_visible(False) ax.set_yticks([1.]) ax.set_yticklabels([yticks[ian]], fontsize=13) ax.yaxis.tick_right() @@ -488,6 +493,8 @@ def plot_sens_fixed_angle_pretty(data, outfile, outformat, args): (bote[dim][1], ylim[0]), bote[dim][0]-bote[dim][1], np.diff(ylim), fill=True, facecolor=col, alpha=alpha, linewidth=0 )) + bc_limit = best_limits[dim] + # ax.axvline(x=np.log10(bc_limit[1]), color=bc_limit[2], alpha=0.7, linewidth=1.5) scales, statistic = data[idim][isrc][ian].T tck, u = interpolate.splprep([scales, statistic], s=0) @@ -532,6 +539,15 @@ def plot_sens_fixed_angle_pretty(data, outfile, outformat, args): fill=False, hatch='\\\\' )) + # fig.text( + # 0.06, 0.5, r'Dimension $d$ of SME coefficient', fontsize=16, rotation='90', + # verticalalignment='center' + # ) + # fig.text( + # 0.96, 0.5, r'New Physics Maximal Mixing Scenario', fontsize=16, + # rotation='270', verticalalignment='center' + # ) + ax.get_xaxis().set_visible(True) ax.set_xlabel(r'${\rm New\:Physics\:Scale}\:[\:{\rm log}_{10} (\Lambda^{-1}\:/\:{\rm GeV}^{-d+4})\: ]$', fontsize=19) |