Plot 88Sr+ ion clock transition frequencies

strontium_freqs.py

# Sr+ frequencies
# AW2019-11-20
# https://www.bipm.org/utils/common/pdf/mep/88Sr+_445THz_2017.pdf

import matplotlib.pyplot as plt
import numpy

f_srs = 444779044095486.50
f_offset = 444779044095000
u_y_srs = 1.5e-15
u_f_srs = u_y_srs*f_srs
print u_f_srs

data = [ [444779044095484.6, 1.5, "Margolis2004", "484.6(1.5)"],
         [444779044095484, 15, "Dube2005", "484(15)"],
         [444779044095485.5, 0.9, "Madej2012", "485.5(9)"],
         [444779044095486.71, 0.24, "Barwood2014", "486.71(24)"],
         [444779044095485.27, 0.75, "Dube2017", "485.27(75)"]
      ]

fig, ax = plt.subplots() #figure()
for (idx, d) in enumerate(data):
    plt.errorbar(idx, d[0]-f_offset, yerr=d[1],fmt='o',capsize=10,label=d[2])
    plt.text(idx,d[0]-f_offset, d[3] )
    
plt.plot([0,5],[f_srs-f_offset, f_srs-f_offset],'g--',label="SRS value, 486.50(67)")
plt.title('88Sr+ clock transition frequency')
ax.axhspan( f_srs-f_offset-u_f_srs, f_srs-f_offset+u_f_srs, alpha=0.1, color='green')
plt.ylim((f_srs-f_offset-5,f_srs-f_offset+5))
plt.grid()
plt.ylabel('f - 444779044095000 / Hz')
plt.legend()
plt.show()