Cu = Stoner.DataFile('/Volumes/data/Projects/Spincurrents/Joe Batley/Measurements/SC004/Transport/Reference Data/Cu spacer resistance/Resistivity/SC004_2_T_Cu_resistivity_vs_T.txt')
Cu.sort('T (K)')
print Cu.column_headers
CuR = interpolate.interp1d(Cu.column('T (K)'),Cu.column(r'$\rho$ ($\Omega$m)'))
# Import Delta R vs T data and group
pattern = re.compile('SC004_(?P<L>\d*)_(?P<Device>\w*)_DeltaRvsT')
folder = DataFolder('/Volumes/data/Projects/Spincurrents/Joe Batley/Measurements/SC004/Transport/DeltaR_vs_Sep/Plot',pattern = pattern)
folder.group('L')
#T = folder[0].column('Sample Temp')
# Walk groups and interpolate
Spinsig = []
Spinsigerr = []
Sep = []
folder.walk_groups(interpSig,group=True,replace_terminal=True)
print Sep
# create a set of Parameters
params = Parameters()
#params.add('Lambda_N', value = 500e-9,min=0)
params.add('Alpha', value = 0.4,min=0,max=1)
T1 = numpy.arange(10,60,5)
T2 = numpy.arange(60,100,10)
T3 = numpy.arange(100,260,20)
T = numpy.concatenate((T1,T2,T3),axis=0)
Lambda = []
avg.apply(err, 1, replace = False, header = 'Err Coef')
return avg
def P_AP_col(folder,keys):
coef = Stoner.DataFile()
for f in folder:
coef.add_column(f.column('Mean Coef'),str(f['IVtemp']))
coef.add_column(f.column('Err Coef'),str(f['IVtemp'])+' Err')
print coef
return coef
pattern = re.compile('_(?P<IVtemp>\d*)K_NLDCIV_300uA_DigFilt10rep_(?P<state>\w*)_')
folder = DataFolder('/Users/Joe/PhD/Measurements/RN0151_4T/NLIVvsT/Both/',pattern = pattern)
folder.group(['state','IVtemp'])
folder.walk_groups(IV_group_Avg,group=True,replace_terminal=True)
folder.walk_groups(P_AP_col,group=True,replace_terminal=True)
print folder['AP']
for f in folder:
for column in f.column_headers:
plt.title(r'NL IV offset vs Temperature for Parallel (red) and Antiparallel (blue) State')
plt.xlabel(r'Temperature (K)')
plt.ylabel(r'$\alpha$ ($\mu$V/A) ')
plt.ticklabel_format(style='plain', scilimits=(3 ,3))
plt.hold(True)
plt.grid(True)
CuRvt.sort("Sample Temp")
Res_Cu = interpolate.interp1d(CuRvt[:, 3], CuRvt[:, 2])
Res_V = interpolate.interp1d(InjV[:, 3], InjV[:, 2])
Res_Py = interpolate.interp1d(InjPy[:, 3], InjPy[:, 2])
Acu = 130e-9 * 150e-9
Asi = 150e-9 * 16e-6
dz = 1000e-9
dx = 425e-9
Seperation = {1: 325e-9, 2: 425e-9, 3: 525e-9, 4: 625e-9, 5: 725e-9, 6: 925e-9, 7: 1125e-9, 8: 1325e-9, 9: 1525e-9}
#################### IMPORTDATA AND WALK GROUPS ####################
pattern = re.compile("_(?P<state>\d*)_(?P<IVtemp>\d*)K_(?P<Inj>\w*)_NLIV_300uA_")
folder = DataFolder(
"/Volumes/data/Projects/Spincurrents/Joe Batley/Measurements/SC004/Transport/SC004_3_T/NLIVvsHvsT_BOTH",
pattern=pattern,
)
Output = Stoner.DataFile()
Output["Sample ID"] = folder[0]["Sample ID"]
folder.group("IVtemp")
# print folder
folder.walk_groups(Py_V, group=True, replace_terminal=True)
Spinsig = []
Spinsig_error = []
beta = []
betaerr = []
Temp = []
heatdata = []
RSP=[]
RSAP=[]
Output = Stoner.DataFile()
Output['Sample ID'] = folder[0]['Sample ID']
folder.group('IVtemp')
print folder
'''
folder.walk_groups(HeatAnalysis,group=True,replace_terminal=True)
Heat = Stoner.DataFile()
Heat.column_headers = ['T (K)','$\Beta$','T$_i$-T$_s$','T$_d$-T$_s$','T$_i$-T$_d$']
Heat['Sample ID'] = Output['Sample ID']
Heat.data = heatdata
Heat.sort('T (K)')
Heat.save('/Volumes/data/Projects/Spincurrents/Joe Batley/Measurements/SC004/Transport/Heat Analysis/' + Heat['Sample ID'] + '_TempdatavsT.txt')
print Heat
p=SP.PlotFile(Heat)
p.setas="x...y"
p.template=SPF.DefaultPlotStyle
title = r'$\Delta$R$_s$/$\rho$ vs temperature'
p.plot(label=p['Sample ID'],figure=1,title=title)
'''
####### CREATE OUTPUT FILE ######
DeltaR = workfile() #Avoid creating temporary python lists, use the final DataFile like object to hold the data
DeltaR.metadata=folder[0].metadata
DeltaR['Sample ID'] = folder[0]['Sample ID']
DeltaR.column_headers=["T","P","Perr","AP","APerr","DR mV","DRerr","Voff","Ptest","APtest"]
#Use the labels attribute to store plot labels that are differnt from column names
DeltaR.labels=[r'T (K)',r'$R_s(P)$ (mV/A)','Perr',r'$R_s(AP)$ (mV/A)','APerr',r'$\Delta R_s$ (mV/A)','DRerr mV',r'$R_s$ offset (mV/A)',"Test Columns"]
####### Group into Temp - walk group and colapse each temp into one file ######
folder.group('IV_Temp')
folder.walk_groups(RsvsH,group=True,replace_terminal=True)
###### Combine different tempertures into one file #######
for a in folder:
DeltaR+=a[0]
DeltaR.sort("T")
DeltaR.template=SPF.JTBPlotStyle
DeltaR.figure() # Creating new figures like this means we don;t reuse windows from run to run
f=plt.gcf()
f.set_size_inches((5.5,3.75),forward=True) # Set for A4 - will make wrapper for this someday
###### Plot ##########
DeltaR.subplot(221)
