asic='r18'

;---------------------------------------------------
loadct,39
basedir='~/nustar/asic_screening/data/probe_station/'
reportdir='~/nustar/asic_screening/reports/probe_station/'+asic+'/'
   ; Make the directory, if necessary
file_mkdir, reportdir

       ; Input files for all of the tests
sfx_ndis='_a_ndisctest_log.txt'
sfx_tadc='_a_testadc_log.txt'
sfx_wpre='_a_pr_osc.txt'
sfx_wsha='_a_sh_osc.txt'
sfx_wcpm='_a_cp_osc.txt'
sfx_off='_a_pixeldc_log.txt'

asicdir=strcompress(basedir+asic+'/',/remove)
csize =1.0
csize1=1.0

;-----------------   NDISC test  ------------------
filename=strcompress(asicdir+asic+sfx_ndis,/remove)


outps1 =strcompress(reportdir+asic+'_ndisc.ps',/remove)
outpdf1 =strcompress(reportdir+asic+'_ndisc.pdf',/remove)
print,filename
.run setps
.run parse_ndisc
.run fit_erf
.run plot_erf_results
.run endps

;-----------------   Threshold ADC  ---------------
filename=strcompress(asicdir+asic+sfx_tadc,/remove)
print,filename
outps1 =strcompress(reportdir+asic+'_adc.ps',/remove)
outpdf1 =strcompress(reportdir+asic+'_adc.pdf',/remove)
.run setps
.run parse_adc_lin
.run plot_adc_lin
.run endps

;-----------------   Waveform plot  ---------------
filename=strcompress(asicdir+asic+sfx_wpre,/remove)
print,filename
outps1 =strcompress(reportdir+asic+'_preamp.ps',/remove)
outpdf1 =strcompress(reportdir+asic+'_preamp.pdf',/remove)
.run setps

yrange = [-.1, 0.3]
.run parse_preamp.pro
.run endps







filename=strcompress(asicdir+asic+sfx_wsha,/remove)
print,filename
yrange = [-.3, 0.8]
outps1 =strcompress(reportdir+asic+'_waveform.ps',/remove)
outpdf1 =strcompress(reportdir+asic+'_waveform.pdf',/remove)

.run setps
.run parse_waveform.pro
.run endps








filename=strcompress(asicdir+asic+sfx_wcpm,/remove)
print,filename
yrange = [-.1, 0.2]
outps1 =strcompress(reportdir+asic+'_cpmode.ps',/remove)
outpdf1 =strcompress(reportdir+asic+'_cpmode.pdf',/remove)
.run setps
.run parse_cpmode.pro
.run endps






filename=strcompress(asicdir+asic+sfx_off,/remove)
outps1 =strcompress(reportdir+asic+'_pixeldc.ps',/remove)
outpdf1 =strcompress(reportdir+asic+'_pixeldc.pdf',/remove)
print,filename
.run setps
.run parse_offset
.run fit_offset
.run endps