asicbase='H26'


;---------------------------------------------------
loadct,39
asicdir='~/nustar/asic_screening/data/temperature/'
basereportdir='~/nustar/asic_screening/reports/temperature/'

sfx_pdc='PIXELDC'
sfx_intlin='INTLIN'
sfx_ndisc='DISC'


csize =1.0
csize1=1.0




temp = 'T-15C'
asic = asicbase+temp
datadir=strcompress(asicdir+asic+'/',/remove)
reportdir=strcompress(basereportdir+asic+'/', /remove)


;-----------------   Pixel DC Test  ------------------

; Make the directory, if necessary
file_mkdir, reportdir

filename=strcompress(datadir+sfx_pdc+'_'+asic+'.TXT',/remove)
outps1 =strcompress(reportdir+asic+'_pixeldc.ps',/remove)
outpdf1 =strcompress(reportdir+asic+'_pixeldc.pdf',/remove)
print,filename

 .run setps
.run parse_pixeldc
.run fit_pixeldc
.run endps

;-----------------   NDISC Test  ------------------
filename=strcompress(datadir+sfx_ndisc+'_'+asic+'.TXT',/remove)
outps1 =strcompress(reportdir+asic+'_ndisc.ps',/remove)
outpdf1 =strcompress(reportdir+asic+'_ndisc.pdf',/remove)

.run setps
print,filename
.run parse_ndisc
.run fit_erf
.run plot_erf_results
.run endps




;-----------------   Integral Linearity Test  ------------------

filename=strcompress(datadir+sfx_intlin+'_'+asic+'.TXT',/remove)
print,filename
outps1 =strcompress(reportdir+asic+'_intlin.ps',/remove)
outpdf1 =strcompress(reportdir+asic+'_intlin.pdf',/remove)
.run setps
.run parse_intlin
.run plot_intlin
.run endps








   ; ----------- Done with -15 C. Do 0C -------------------



temp = 'T0C'
asic = asicbase+temp
datadir=strcompress(asicdir+asic+'/',/remove)
reportdir=strcompress(basereportdir+asic+'/', /remove)


;-----------------   Pixel DC Test  ------------------

; Make the directory, if necessary
file_mkdir, reportdir

filename=strcompress(datadir+sfx_pdc+'_'+asic+'.TXT',/remove)
outps1 =strcompress(reportdir+asic+'_pixeldc.ps',/remove)
outpdf1 =strcompress(reportdir+asic+'_pixeldc.pdf',/remove)
print,filename

 .run setps
.run parse_pixeldc
.run fit_pixeldc
.run endps

;-----------------   NDISC Test  ------------------
filename=strcompress(datadir+sfx_ndisc+'_'+asic+'.TXT',/remove)
outps1 =strcompress(reportdir+asic+'_ndisc.ps',/remove)
outpdf1 =strcompress(reportdir+asic+'_ndisc.pdf',/remove)

.run setps
print,filename
.run parse_ndisc
.run fit_erf
.run plot_erf_results
.run endps




;-----------------   Integral Linearity Test  ------------------

filename=strcompress(datadir+sfx_intlin+'_'+asic+'.TXT',/remove)
print,filename
outps1 =strcompress(reportdir+asic+'_intlin.ps',/remove)
outpdf1 =strcompress(reportdir+asic+'_intlin.pdf',/remove)
.run setps
.run parse_intlin
.run plot_intlin
.run endps





   ; ----------- Done with -0 C. Do 40C -------------------



temp = 'T40C'
asic = asicbase+temp
datadir=strcompress(asicdir+asic+'/',/remove)
reportdir=strcompress(basereportdir+asic+'/', /remove)


;-----------------   Pixel DC Test  ------------------

; Make the directory, if necessary
file_mkdir, reportdir

filename=strcompress(datadir+sfx_pdc+'_'+asic+'.TXT',/remove)
outps1 =strcompress(reportdir+asic+'_pixeldc.ps',/remove)
outpdf1 =strcompress(reportdir+asic+'_pixeldc.pdf',/remove)
print,filename

 .run setps
.run parse_pixeldc
.run fit_pixeldc
.run endps

;-----------------   NDISC Test  ------------------
filename=strcompress(datadir+sfx_ndisc+'_'+asic+'.TXT',/remove)
outps1 =strcompress(reportdir+asic+'_ndisc.ps',/remove)
outpdf1 =strcompress(reportdir+asic+'_ndisc.pdf',/remove)

.run setps
print,filename
.run parse_ndisc
.run fit_erf
.run plot_erf_results
.run endps




;-----------------   Integral Linearity Test  ------------------

filename=strcompress(datadir+sfx_intlin+'_'+asic+'.TXT',/remove)
print,filename
outps1 =strcompress(reportdir+asic+'_intlin.ps',/remove)
outpdf1 =strcompress(reportdir+asic+'_intlin.pdf',/remove)
.run setps
.run parse_intlin
.run plot_intlin
.run endps