#
# glish script to plot a bunch of results from TSPECTEST_02
# to allow final plots to be re-done as necessary
#
# R. Prestage 9th September 2002
#
# requires "plotit.g" from this directory to be already loaded.
#
# call args: plotit (xmin, xmax, ymin, ymax, vel=T, temp=T, factor=1.91,
#      xlbl = "X", ylbl = "Y", toplbl = "Plot", file ="plot.ps")
# 

# start up the plotter

# set factor for all plots

factor := 1.91

# scan 33 passband spectrum

d.plotr(33,1,,1)
xlbl := "Barycentric frequency (MHz)"
ylbl := "Intensity (CNTS)"
file := "figure1.ps"
toplbl := "Fig 1: TSPECTEST_03 Scan 33 SL1720+6236"
plotit(1375,1430,-0.5,4,vel=F,temp=T,factor=factor,xlbl=xlbl, ylbl=ylbl,
       toplbl=toplbl,file=file)

print readline("continue:")

# plot the start of linear Y poln of scan 33 reduced spectrum against
# frequency to show noisy behaviour

d.plotc(33,,2)
xlbl := "Barycentric frequency (MHz)"
ylbl := "Intensity (TA*)"
file := "figure2.ps"
toplbl := "Fig 2: TSPECTEST_02 Scan 33 SL1720+623"
plotit(1378,1381,-0.4,0,vel=F,temp=T,factor=factor, 
      xlbl=xlbl, ylbl=ylbl, toplbl=toplbl,file=file)

print readline("continue:")

# plot reduced scan 35 y poln to show a bad baseline

d.plotc(35,,2)
xlbl := "Barycentric frequency (MHz)"
ylbl := "Intensity (TA*)"
file := "figure3.ps"
toplbl := "Fig 3: TSPECTEST_03 Scan 35 SL1720+6236"
plotit(1375,1430,-0.5,0.1,vel=F,temp=T,factor=factor,xlbl=xlbl, ylbl=ylbl,
       toplbl=toplbl,file=file)

print readline("continue:")


# Now make the average of the seven good scans 46 through 61


# above nregion didn't seem to work - need to set it manually!

d.plotc(37)

   d.nfit(1)
   d.nregion(232, 989, 1677, 3588, 5430, 6876)

#print " "
#print "set two baselines regions now...."
#print " "
#d.setregion()

print "look at the plot, and change scale to pixels...."

d.plotc(37,,1)

print readline("continue:")

d.sclear()
for (i in [37,39,41,43,45,47,49,51]) {
   print "processing scan", i
   d.getc(i);d.baseline();d.accum()
}
d.ave();d.avgFeeds();d.show()

# plot the whole thing

xlbl := "Barycentric frequency (MHz)"
ylbl := "Intensity (TA*)"
file := "figure4.ps"
toplbl := "Fig 4: TSPECTEST_03 Scans 37-51 SL1720+6236"
plotit(1385,1418,-0.02,0.02,vel=F,temp=T,factor=factor,xlbl=xlbl, ylbl=ylbl,
       toplbl=toplbl,file=file)

print readline("continue:")

# and a zoomed in region for comparison with Rick's results

xlbl := "Barycentric frequency (MHz)"
ylbl := "Intensity (Jy)"
file := "figure5.ps"
toplbl := "Fig 5: TSPECTEST_03 Scans 37-51 SL1720+6236"
plotit(1402.6,1404.4,-0.013,0.043,vel=F,temp=F,factor=factor,xlbl=xlbl, 
        ylbl=ylbl,
       toplbl=toplbl,file=file)