Файл: Mplwp ramsey fringes monochromatic.svg

#! / usr / bin / python # - * - кодировка: utf8 - * -import matplotlib.pyplot as pltimport matplotlib as mplimport numpy as npfrom math import *code_website = 'http://commons.wikimedia.org/wiki/User:Geek3/mplwp'try: import mplwpexcept ImportError, er: print 'ImportError:', er print 'You need to download mplwp.py from', code_website exit(1)name = 'mplwp_ramsey_fringes_monochromatic.svg'fig = mplwp.fig_standard(mpl)xlim = -5,5; fig.gca().set_xlim(xlim)ylim = -0.05,1.15; fig.gca().set_ylim(ylim)mplwp.mark_axeszero(fig.gca())# free evolution time vs. interaction pulse time ratioT_over_tau = 8def flip(delta): # return spin flip probability of atom beam in Ramsey device # delta is the detuning in terms of the Rabi frequency Omega0 Omega_rel = sqrt(1.0 + delta**2) phase_pulse = 0.25 * pi * Omega_rel phase_free = 0.25 * pi * delta * T_over_tau Pflip = (sin(phase_pulse) * (Omega_rel * cos(phase_pulse) * cos(phase_free) - delta * sin(phase_pulse) * sin(phase_free)) * 2. / Omega_rel**2)**2 return Pflipx = np.linspace(xlim[0], xlim[1], 5001)y = [flip(xx) for xx in x]mpl.rc('legend', borderaxespad=1.0)plt.plot(x, y, label='$P_{flip}(\Delta\omega/\Omega_0)$')plt.legend(loc='upper right').get_frame().set_alpha(0.9) plt.savefig(name)mplwp.postprocess(name)