#
# 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 46 passband spectrum

d.plotr(46,1,,1)
xlbl := "Barycentric frequency (MHz)"
ylbl := "Intensity (CNTS)"
file := "figure1.ps"
toplbl := "Fig 1: TSPECTEST_02 Scan 46 U9179o"
plotit(1412,1427,-0.2,3,vel=F,temp=T,factor=factor,xlbl=xlbl, ylbl=ylbl,
       toplbl=toplbl,file=file)

print readline("continue:")

# plot the start of linear X poln of scan 46 reduced spectrum against
# frequency to show noisy behaviour
 
d.plotc(46,,1)
xlbl := "Barycentric frequency (MHz)"
ylbl := "Intensity (TA*)"
file := "figure2.ps"
toplbl := "Fig 2: TSPECTEST_02 Scan 46 U9179o"
plotit(1412.5,1413.6,-2.2,2.2,vel=F,temp=T,factor=factor, 
      xlbl=xlbl, ylbl=ylbl, toplbl=toplbl,file=file)

print readline("continue:")

# scan 61 reduced spectrum to show one of the "worst" spectra

d.plotc(61,,1)
xlbl := "Barycentric frequency (MHz)"
ylbl := "Intensity (TA*)"
file := "figure3.ps"
toplbl := "Fig 3: TSPECTEST_02 Scan 61 U9179o"
plotit(1413.5,1425.5,-0.5,2,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

   d.nfit(1)
   d.nregion(1700,3500,6000,7500)

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

d.plotc(46)

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

print readline("continue:")

d.sclear()
for (i in [46,51,53,55,57,59,61]) {
   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_02 Scans 46-61 U9179o"
plotit(1413.5,1425.5,-0.5,2,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 radio velocity (km/s)"
ylbl := "Intensity (Jy)"
file := "figure5.ps"
toplbl := "Fig 5: TSPECTEST_02 Scans 46-61 U9179o"
plotit(100,500,-0.2,1,vel=T,temp=F,factor=factor,xlbl=xlbl, ylbl=ylbl,
       toplbl=toplbl,file=file)