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( [