# Advanced Calibration demo
# gmoellen (07Oct17)

# BPOLY demo
#  ngc5921
# a pre-B gaincal
default(gaincal);
gaincal(vis='ngc5921.ms',
        caltable='prelim.G',
        field='1331*,1445*',
        spw='0:20~40',
        solint=20.0);

# an ordinary sampled B
default(bandpass);
bandpass(vis='ngc5921.ms',
         caltable='cal.B',
         field='1331*,1445*',
         solint=0.0,
         refant='5',
         gaintable='prelim.G');


# an ordinary sampled B
default(bandpass);
bandpass(vis='ngc5921.ms',
         caltable='cal.Bpoly',
         field='1331*',
         bandtype='BPOLY',
         degamp=25,
         degphase=15,
         refant='5',
         gaintable='prelim.G');

# overlay plot of bandpasses
default(plotcal)
plotsymbol='o';
poln='R';
yaxis='amp';
antenna='3';
plotcal(caltable='cal.B',plotcolor='blue',overplot=F)
plotcal(caltable='cal.Bpoly',plotcolor='red',overplot=T)


# GSPLINE, etc. demo
# NGC4826 BIMA data @ 115 GHz

# solve for a GSPLINE
default(gaincal);
gaincal(vis='n4826_16apr.ms',
        caltable='test.Gspl7200',
        field='1',
        spw='1',
        gaintype='GSPLINE',
        splinetime=7200.0,
        refant='ANT1',
        preavg=300.0);

# solve for an ordinary G
default(gaincal);
gaincal(vis='n4826_16apr.ms',
        caltable='test.G',
        field='1',
        spw='1',
        solint=300.0,
        refant='ANT1');


default(plotcal);
poln='X'
yaxis='phase';
plotsymbol='o';
antenna='ANT3'
plotcal(caltable='test.Gspl7200',yaxis='phase',plotcolor='blue',overplot=F);
plotcal(caltable='test.G',plotcolor='red',overplot=T);



# accum demo
default(accum)
accum(vis='n4826_16apr.ms',
      accumtime=60.0,
      incrtable='test.G',
      caltable='test.Gacc',
      field='1~8',
      calfield='1',
      interp='linear',
      spwmap=[0,1,1,1,1,1]);

default(plotcal);
poln='X'
yaxis='phase';
plotsymbol='o';
antenna='ANT3'
plotcal(caltable='test.Gspl7200',yaxis='phase',plotcolor='blue',overplot=F);
plotcal(caltable='test.G',plotcolor='red',overplot=T);
plotcal(caltable='test.Gacc',plotcolor='green',overplot=T);
      

# smooth demo
default(smoothcal)
smoothcal(vis='n4826_16apr.ms',
          tablein='test.G',
          caltable='test.Gsm',
          smoothtime=600.0)

default(accum)
accum(vis='n4826_16apr.ms',
      accumtime=60.0,
      incrtable='test.Gsm',
      caltable='test.Gsmacc',
      field='1~8',
      calfield='1',
      interp='linear',
      spwmap=[0,1,1,1,1,1]);

default(plotcal);
poln='X'
yaxis='phase';
plotsymbol='o';
antenna='ANT3'
plotcal(caltable='test.Gspl7200',yaxis='phase',plotcolor='blue',overplot=F);
plotcal(caltable='test.G',plotcolor='red',overplot=T);
plotcal(caltable='test.Gsmacc',plotcolor='green',overplot=T);