set cols ( blue red black green purple orange cyan violet yellow )

set toppath "/sps/baoradio/AmasNancay/JEC"
set source "Abell${1}"
set srclower "abell${1}"
set date   "${2}"


#####################################
#Plots ON-OFF with all calibrations 
#####################################
openppf ${toppath}/${source}/${date}${srclower}/diffOnOff_${date}_${srclower}.ppf


#specONOFF<i>       Matrix 2 x 8192 : i=0...
#specONOFFMean      Matrix 2 x 8192
#specONOFF2ChanMean Vector 8192 x 1


#####################################
#Plots with different Means: per channel and all 
#####################################
del sallcyclev0
objaoper specONOFFMean row 0 sallcyclev0
del sallcyclev1
objaoper specONOFFMean row 1 sallcyclev1

newwin 1 1
graphicatt "xylimits=1250,1500,-0.5,0.5"
plot2d specONOFF2ChanMean (n/8192)*250+1250 val n>0 "black cpts nsta notit"
plot2d sallcyclev0        (n/8192)*250+1250 val n>0 "same blue cpts nsta notit"
plot2d sallcyclev1        (n/8192)*250+1250 val n>0 "same red cpts nsta notit"
settitle "ON-OFF  ${source}/${date} Ch 0 (blue) Ch 1 (red) Mean (black)  All cycles"
setaxelabels "Freq. (MHz)" "I (Jy)"

#####################################
#Plots ON-OFF with all calibrations: chan 0 & 1 separaed and cycle per cycle 
#####################################
newwin 1 1
graphicatt "xylimits=1250,1500,-0.5,0.5"
del sc1v0
objaoper specONOFF1 row 0 sc1v0
plot2d sc1v0 (n/8192)*250+1250 val n>0 "$cols[0] cpts nsta notit"
for ic 2:13
    k = ${ic}-floor(${ic}/${#cols})*${#cols}
    del sc${ic}v0
    objaoper specONOFF${ic} row 0 sc${ic}v0	
    plot2d sc${ic}v0 (n/8192)*250+1250 val n>0 "same $cols[k] cpts nsta notit"
end
settitle "ON-OFF  ${source}/${date} Ch 0 cycle per cycle"
setaxelabels "Freq. (MHz)" "I (Jy)"

newwin 1 1
graphicatt "xylimits=1250,1500,-0.5,0.5"
del sc1v1
objaoper specONOFF1 row 1 sc1v1
plot2d sc1v1 (n/8192)*250+1250 val n>0 "$cols[0] cpts nsta notit"
for ic 2:13
    k = ${ic}-floor(${ic}/${#cols})*${#cols}
    del sc${ic}v1
    objaoper specONOFF${ic} row 1 sc${ic}v1	
    plot2d sc${ic}v1 (n/8192)*250+1250 val n>0 "same $cols[k] cpts nsta notit"
end
settitle "ON-OFF  ${source}/${date} Ch 1 cycle per cycle"
setaxelabels "Freq. (MHz)" "I (Jy)"


#####################################
# Plots RAW no calibration at all
#####################################
openppf ${toppath}/${source}/${date}${srclower}/dataRaw_${date}_${srclower}.ppf
newwin 2 2
graphicatt "xylimits=1250,1500,0.95,1.05"

###         #
### On Ch 0 # On Ch 1 
######################
### Off Ch 0# Off Ch 1
###         #

set mode On	
set first 1
for ic 1:13
  k = ${ic}-floor(${ic}/${#cols})*${#cols}	
    del sm${mode}c${ic}v0
    objaoper specRaw${mode}${ic} row 0 sm${mode}c${ic}v0
    if ( ${first} == 1 ) then
      plot2d sm${mode}c${ic}v0 (n/8192)*250+1250 val n>0 "$cols[k] cpts nsta notit"
    else
      plot2d sm${mode}c${ic}v0 (n/8192)*250+1250 val n>0 "same $cols[k] cpts nsta notit"
    endif
  set first 0
end
settitle "${source}/${date} Raw: On Ch 0"

set first 1
for ic 1:13
  k = ${ic}-floor(${ic}/${#cols})*${#cols}	
    del sm${mode}c${ic}v1
    objaoper specRaw${mode}${ic} row 1 sm${mode}c${ic}v1
    if ( ${first} == 1 ) then
      plot2d sm${mode}c${ic}v1 (n/8192)*250+1250 val n>0 "$cols[k] cpts nsta notit"
    else
      plot2d sm${mode}c${ic}v1 (n/8192)*250+1250 val n>0 "same $cols[k] cpts nsta notit"
    endif
  set first 0
end
settitle "${source}/${date} Raw: On Ch 1"

set mode Off
set first 1
for ic 1:13
  k = ${ic}-floor(${ic}/${#cols})*${#cols}	
    del sm${mode}c${ic}v0
    objaoper specRaw${mode}${ic} row 0 sm${mode}c${ic}v0
    if ( ${first} == 1 ) then
      plot2d sm${mode}c${ic}v0 (n/8192)*250+1250 val n>0 "$cols[k] cpts nsta notit"
    else
      plot2d sm${mode}c${ic}v0 (n/8192)*250+1250 val n>0 "same $cols[k] cpts nsta notit"
    endif
  set first 0
end
settitle "${source}/${date} Raw: Off Ch 0"

set first 1
for ic 1:13
  k = ${ic}-floor(${ic}/${#cols})*${#cols}	
    del sm${mode}c${ic}v1
    objaoper specRaw${mode}${ic} row 1 sm${mode}c${ic}v1
    if ( ${first} == 1 ) then
      plot2d sm${mode}c${ic}v1 (n/8192)*250+1250 val n>0 "$cols[k] cpts nsta notit"
    else
      plot2d sm${mode}c${ic}v1 (n/8192)*250+1250 val n>0 "same $cols[k] cpts nsta notit"
    endif
  set first 0
end
settitle "${source}/${date} Raw: Off Ch 1"

#####################################
# Plots BAO calibration included
#####################################
openppf ${toppath}/${source}/${date}${srclower}/dataBAOCalib_${date}_${srclower}.ppf
newwin 2 2
#graphicatt "xylimits=1250,1500,0.95,1.05"
graphicatt ""


###         #
### On Ch 0 # On Ch 1 
######################
### Off Ch 0# Off Ch 1
###         #

set mode On	
set first 1
for ic 1:13
  k = ${ic}-floor(${ic}/${#cols})*${#cols}	
    del sm${mode}c${ic}v0
    objaoper specBAOCalib${mode}${ic} row 0 sm${mode}c${ic}v0
    if ( ${first} == 1 ) then
      plot2d sm${mode}c${ic}v0 (n/8192)*250+1250 val n>0 "$cols[k] cpts nsta notit"
    else
      plot2d sm${mode}c${ic}v0 (n/8192)*250+1250 val n>0 "same $cols[k] cpts nsta notit"
    endif
  set first 0
end
settitle "${source}/${date} BAOCalib: On Ch 0"

set first 1
for ic 1:13
  k = ${ic}-floor(${ic}/${#cols})*${#cols}	
    del sm${mode}c${ic}v1
    objaoper specBAOCalib${mode}${ic} row 1 sm${mode}c${ic}v1
    if ( ${first} == 1 ) then
      plot2d sm${mode}c${ic}v1 (n/8192)*250+1250 val n>0 "$cols[k] cpts nsta notit"
    else
      plot2d sm${mode}c${ic}v1 (n/8192)*250+1250 val n>0 "same $cols[k] cpts nsta notit"
    endif
  set first 0
end
settitle "${source}/${date} BAOCalib: On Ch 1"

set mode Off
set first 1
for ic 1:13
  k = ${ic}-floor(${ic}/${#cols})*${#cols}	
    del sm${mode}c${ic}v0
    objaoper specBAOCalib${mode}${ic} row 0 sm${mode}c${ic}v0
    if ( ${first} == 1 ) then
      plot2d sm${mode}c${ic}v0 (n/8192)*250+1250 val n>0 "$cols[k] cpts nsta notit"
    else
      plot2d sm${mode}c${ic}v0 (n/8192)*250+1250 val n>0 "same $cols[k] cpts nsta notit"
    endif
  set first 0
end
settitle "${source}/${date} BAOCalib: Off Ch 0"

set first 1
for ic 1:13
  k = ${ic}-floor(${ic}/${#cols})*${#cols}	
    del sm${mode}c${ic}v1
    objaoper specBAOCalib${mode}${ic} row 1 sm${mode}c${ic}v1
    if ( ${first} == 1 ) then
      plot2d sm${mode}c${ic}v1 (n/8192)*250+1250 val n>0 "$cols[k] cpts nsta notit"
    else
      plot2d sm${mode}c${ic}v1 (n/8192)*250+1250 val n>0 "same $cols[k] cpts nsta notit"
    endif
  set first 0
end
settitle "${source}/${date} BAOCalib: Off Ch 1"

#####################################
# Plots BAO+RT calibration included
#####################################
openppf ${toppath}/${source}/${date}${srclower}/dataBAORTCalib_${date}_${srclower}.ppf
newwin 2 2
#graphicatt "xylimits=1250,1500,0.95,1.05"
graphicatt ""


###         #
### On Ch 0 # On Ch 1 
######################
### Off Ch 0# Off Ch 1
###         #

set mode On	
set first 1
for ic 1:13
  k = ${ic}-floor(${ic}/${#cols})*${#cols}	
    del sm${mode}c${ic}v0
    objaoper specBAORTCalib${mode}${ic} row 0 sm${mode}c${ic}v0
    if ( ${first} == 1 ) then
      plot2d sm${mode}c${ic}v0 (n/8192)*250+1250 val n>0 "$cols[k] cpts nsta notit"
    else
      plot2d sm${mode}c${ic}v0 (n/8192)*250+1250 val n>0 "same $cols[k] cpts nsta notit"
    endif
  set first 0
end
settitle "${source}/${date} BAORTCalib: On Ch 0"

set first 1
for ic 1:13
  k = ${ic}-floor(${ic}/${#cols})*${#cols}	
    del sm${mode}c${ic}v1
    objaoper specBAORTCalib${mode}${ic} row 1 sm${mode}c${ic}v1
    if ( ${first} == 1 ) then
      plot2d sm${mode}c${ic}v1 (n/8192)*250+1250 val n>0 "$cols[k] cpts nsta notit"
    else
      plot2d sm${mode}c${ic}v1 (n/8192)*250+1250 val n>0 "same $cols[k] cpts nsta notit"
    endif
  set first 0
end
settitle "${source}/${date} BAORTCalib: On Ch 1"

set mode Off
set first 1
for ic 1:13
  k = ${ic}-floor(${ic}/${#cols})*${#cols}	
    del sm${mode}c${ic}v0
    objaoper specBAORTCalib${mode}${ic} row 0 sm${mode}c${ic}v0
    if ( ${first} == 1 ) then
      plot2d sm${mode}c${ic}v0 (n/8192)*250+1250 val n>0 "$cols[k] cpts nsta notit"
    else
      plot2d sm${mode}c${ic}v0 (n/8192)*250+1250 val n>0 "same $cols[k] cpts nsta notit"
    endif
  set first 0
end
settitle "${source}/${date} BAORTCalib: Off Ch 0"

set first 1
for ic 1:13
  k = ${ic}-floor(${ic}/${#cols})*${#cols}	
    del sm${mode}c${ic}v1
    objaoper specBAORTCalib${mode}${ic} row 1 sm${mode}c${ic}v1
    if ( ${first} == 1 ) then
      plot2d sm${mode}c${ic}v1 (n/8192)*250+1250 val n>0 "$cols[k] cpts nsta notit"
    else
      plot2d sm${mode}c${ic}v1 (n/8192)*250+1250 val n>0 "same $cols[k] cpts nsta notit"
    endif
  set first 0
end
settitle "${source}/${date} BAORTCalib: Off Ch 1"