def gshift(nu):
T0 = M.T(0, B)
return M.convert(((T0 + e0) / 2.0) + sqrt(((T0 - e0) / 2.0) ** 2 + l * T0 ** 2 * (-0.25 * nu + 1.5)) - e0, 1)
inches_per_pt = 3.4/246.0 # Convert pt to inches
golden_mean = (sqrt(5)-1.0)/2.0 # Aesthetic ratio
fig_width = fig_width_pt*inches_per_pt # width in inches
fig_height =fig_width*golden_mean # height in inches
fig_size = [fig_width,fig_height]
params = {'backend': 'GTKAgg','axes.labelsize': 8,'font.size': 8,'legend.fontsize': 8,'xtick.labelsize': 8,'ytick.labelsize': 8, 'text.usetex': False,'figure.figsize': fig_size}
rcParams.update(params)
"""
# rapidly plot simulation over simulated gatesweep data
try:
delta = double(sys.argv[1])
l = double(sys.argv[2])
shift0 = double(sys.argv[3])
cutoff = int(sys.argv[4])
nrange = arange(double(sys.argv[5]), double(sys.argv[6]), double(sys.argv[7]))
e0 = M.convert(shift0, 2)
B = double(sys.argv[8])
vF = double(sys.argv[9])
l2 = double(sys.argv[10])
except:
print "Usage: ", sys.argv[0], " delta l shift0 cutoff n-low n-high n-step B vF"
sys.exit(1)
width = 0.001
savedir = "./"
# ----------------------------------------------------------------
FitDataPath = "C:/Users/Owner/research/graphene/graphene_hf/results/2009_MagnetoPhonons_121/EXP20090309_121_raman-12_6T-Fit_DoubleLorentz_Intraband.txt"
FitData = loadtxt(FitDataPath)
mask1 = array(
[