;
;
;

; look at NIM darks 02sep18

dpath='/usr2/mpenn/02sep20/'

dark140_100=swap_endian(rfits(dpath+'tia140_100.d.i.001'))
dark140_200=swap_endian(rfits(dpath+'tia140_200.d.i.001'))
dark140_300=swap_endian(rfits(dpath+'tia140_300.d.i.001'))
dark150_100=swap_endian(rfits(dpath+'tia150_100.d.i.001'))
dark150_200=swap_endian(rfits(dpath+'tia150_200.d.i.001'))
dark150_300=swap_endian(rfits(dpath+'tia150_300.d.i.001'))
dark160_100=swap_endian(rfits(dpath+'tia160_100.d.i.001'))
dark160_200=swap_endian(rfits(dpath+'tia160_200.d.i.001'))
dark160_300=swap_endian(rfits(dpath+'tia160_300.d.i.001'))
dark180_100=swap_endian(rfits(dpath+'tia180_100.d.i.001'))
dark180_200=swap_endian(rfits(dpath+'tia180_200.d.i.001'))
dark180_300=swap_endian(rfits(dpath+'tia180_300.d.i.001'))


x=[100,200,300]

dmean140=[mean(dark140_100),mean(dark140_200),mean(dark140_300)]
dmean150=[mean(dark150_100),mean(dark150_200),mean(dark150_300)]
dmean160=[mean(dark160_100),mean(dark160_200),mean(dark160_300)]
dmean180=[mean(dark180_100),mean(dark180_200),mean(dark180_300)]

c140=poly_fit(x,dmean140,1)
c150=poly_fit(x,dmean150,1)
c160=poly_fit(x,dmean160,1)
c180=poly_fit(x,dmean180,1)

plot,[140,148,150,160,180],[c140(0),185.6,c150(0),c160(0),c180(0)],psym=2,$
     xtit='TIA VALUE', ytit='Bias value',chars=1.6,xran=[135,185],xstyl=1,$
     yran=[-300,4300],ystyl=1
oplot,[148,150,160],[185.6,c150(0),c160(0)]
oplot,[135,185],[4096,4096],linestyl=2
oplot,[135,185],[0,0],linestyl=2



the_end:
end