# -*- coding: utf-8 -*-
"""
Created on Fri Aug 27 15:14:48 2021

@author: Joschka
"""
from src import physical_constants as cs
import numpy as np

ge = 1.29
gg = 1.24
Me = 9
Mg = 8

mu_eff = (ge*Me- gg*Mg) * cs.mu_B

detuning = 2 * np.pi * 2.4e6 #Hz
dG = 2.26e-4 #T/cm
z = 1


D_loc = detuning - mu_eff/cs.h_bar * dG * z


print(f"D_loc @ D = 2 pi * 2 MHz = {D_loc/(2*np.pi*1e6)} ")

detuning = 2 * np.pi * 2e6 #Hz
 #G/cm
z = 0.75

D_loc = detuning - mu_eff/cs.h_bar * dG * z
print(f"D_loc @ D = 2 pi * 2 MHz = {D_loc/(2*np.pi*1e6)} ")

detuning = 2 * np.pi * 1e6 #Hz
 #G/cm
z = 0.5

D_loc = detuning - mu_eff/cs.h_bar * dG * z
print(f"D_loc @ D = 2 pi * 2 MHz = {D_loc/(2*np.pi*1e6)} ")


dE = mu_eff/cs.h_bar * dG * 0.5

print(dE/(2*np.pi*1e6))