# 2: all memeber curves with comparison # 3: all memeber curves with Diff keyword = '' #'EJ200SP-1P-N1' #'1X1P N1' # read sample #keyword2 = '1X1P N6' # read sample s = Curve.makeslice(300,600) # x-range (band of wavelength in nm) t = True # transmission or absorption curve; None for absorption newdate = True dateN2 = date(2015,12,19) dateN3 = date(2016,01,14) dateN4 = date(2016,01,11) if newdate: dateN4 = date(2016,03,25) dateCT = date(2015,11,3) # inaccurate allc = curvesKeyword(allcurves, name = keyword) #allc2 = curvesKeyword(allcurves, name = keyword2) #allc += allc2 #allc = curveCreation(date.today(), ['data/EJ_20160328_air_SPN4.csv'], verbose = 1) print 'Drawing curves ...\n' for c in allc: [c.plot,c.plotallcurves,c.plotallcurvesComp,c.plotallcurvesDiff][mode](t=t, ratioyrange=[-0.06,0.06],save = False, plotting = True, s = s)
#allctemp = curveCreation(date(2016,10,21), ['data/EJ_20161021-air_calib.csv'], verbose = 1,skips=2)
#curveCreation(date(2016,10,18),
# ['data/EJ_10-18-2016-air_calib.csv'], verbose = 1,skips=2)
if False:
def avg_gen(c):
# local shorthand
# *** Special: to cancel out effect of higher attenuation on C side ***
nonzero = c.arraynzd()[(0, 1, 3), :]
temp = np.average(c.trsm(nonzero), 0)[s]
return temp
keyword = refkeyword = sys.argv[1] #'EJ200SP-1P-N1'
allcurves = curvesKeyword(allctemp, name=keyword)
refcurves = curvesKeyword(refctemp, name=keyword)
print[c.len() for c in allcurves]
print[c.len() for c in refcurves]
print allcurves
print refcurves
cnum, refnum = 0, 1
s = Curve.makeslice(350, 600)
ratioyrange = [-0.06, 0.06]
c_int = allcurves[cnum]
cref = refcurves[refnum]
fig, (ax1, ax2) = pyt.subplots(2,
1,
sharex=True,
gridspec_kw=shr,
farben = ['blue', 'green', 'red' , 'blueviolet','slategray','darkcyan','magenta']
farben += ['hotpink', 'salmon','springgreen','chocolate', 'crimson']
#pyt.figure(figsize = (8,5))
x= np.linspace(600,200,401)
s = Curve.makeslice(350,550) #slice(None) # x-range (band of wavelength in nm)
t = True
refnum = 0
daysIrr = 7
coldwarm = True
legFont = 13
# plot all curves of material A and day of irradaition D
keyword = '200-1X'
ckeyword = curvesKeyword(curves, char = '-', name = keyword)
allc = curvesAge(ckeyword, day = daysIrr, rincl = False, dayrange = 1)
cref = curvesRef(ckeyword, rred=False)
print cref
#specPlotter(c1, 'EJ260-1X1P on day 24', xlabel='wavelength nm', ylabel='transmission', trsm=1)
batchDes = ['', ' UnIrr.', 'N2 @3Mrad', 'N3 @3Mrad High Rate', 'N4 @3Mrad $N_{2}$', 'CastorT']
if coldwarm:
batchDes = ['', ' UnIrr.']+ [0]*7 +[' UnIrr.', r'23$^{\circ}$C', r'-30$^{\circ}$C']
fig, (ax1, ax2) = pyt.subplots(2,1, sharex=True, gridspec_kw=dict(height_ratios = [3,2]),figsize = (8,8))
tLabel = ' Transmission' if t is not None else ' Absorption'
fig.suptitle('EJ'+ keyword+ tLabel +' Spectrum', y=.95, size=17)
iterfarb = iter(farben)
ax1.plot(x[s],cref[refnum].avgcurve(t=t)[s], label = 'Unirradiated ',
color = iterfarb.next())
## to compare with past refs
## version new: rearranged and merged the curve and ratio plotting procedures
# july 9th: check readCurveSequence for similar capability
import numpy as np
import matplotlib.pyplot as pyt
from datetime import date
import sys
from Curve import (Curve, curvesKeyword, curveCreation, keycGen, specPlotter,
curvesAge, curvesRef, readCurveFile)
execfile('readCurveInput_default.py') # loads default setting
execfile('readCurveInput_260PVTDiff.py')
allc = curvesKeyword(allcurves, name=keyword)
allc2 = curvesKeyword(allcurves, name=keyword2)
refc = curvesKeyword(refcurves, name=refkeyword)
cref = refc[refnum]
fig, (ax1, ax2) = pyt.subplots(2,
1,
sharex=True,
gridspec_kw=shr,
figsize=(8, 8))
## TITLES and labels
fig.suptitle(eval(titleGen), y=.95, size=17)
iterfarb, iterfarb2 = iter(farben), iter(farben2)
# DRAW REFERENCE AND ERROR
import matplotlib.pyplot as pyt
from datetime import date
import sys
from Curve import (Curve, curvesKeyword, curveCreation, keycGen, specPlotter,
curvesAge, curvesRef, readCurveFile)
execfile('readCurveInput_default.py') # loads default setting
allcurves = curveCreation(date(2016, 7, 8), ['data/EJ_20160708_air_test.csv'],
verbose=1,
skips=3,
sepcurve=True)
# separate curves are inputed
## selects only EJ260-1X1P... also only the a side measurements
allcurves = curvesKeyword(allcurves, name='EJ260-1X1P')[:27]
allcurves = [c for c in allcurves if c.view(1) is not 0]
cref = sum(allcurves[4:15])
farbGruppe = [[
'blue', 'darkblue', 'darkslateblue', 'royalblue', 'steelblue', 'darkcyan'
], ['red', 'darkred', 'crimson', 'orangered', 'salmon'
], ['green', 'darkgreen', 'olive', 'darkturquoise', 'greenyellow', 'sage'],
['magenta', 'blueviolet', 'purple', 'darkmagenta', 'darkviolet'],
['slategray', 'dimgrey']]
farben = ['blue'] + farbGruppe[1][:3] + farbGruppe[3][:4]
## Parameters
s = Curve.makeslice(350, 600) ##x-range (band of wavelength in nm)
shr = dict(height_ratios=[2, 3]) #subplot height ratio
figsize = (8, 9) #(8,8)