source: BAORadio/AmasNancay/v4/plotDiffOnOff.pic@ 654

#$1: source
#$2: nber of cycles
#$3: number of freq. combined for reducing mean/sigma

set source $1
set ncycles $2
set nFreqBin $3
set action $4

set defatt "font=helvetica,bold,20 fixedfontsize"



if ( ${action} == "mean" ) then

openppf meanDiffOnOffRaw_${source}-${ncycles}Cycles.ppf


##General Mean
del s0 s1
objaoper mean row 0 s0
objaoper mean row 1 s1
graphicatt "${defatt} xylimits=1250,1500,-0.01,0.01"
newwin 1 1
plot2d s0 (n/8192)*250+1250 val n>0 "cpts blue nsta notit"
plot2d s1 (n/8192)*250+1250 val n>0 "same cpts red nsta notit"
settitle "$source (${ncycles} cycles) Ch 0 (blue) Ch 1 (red)" ' ' $defatt
setaxelabels "Freq. (MHz)" "<(ON-OFF)[r,c]/Gain[r]>_{r,c}" $defatt



##Reduced Mean
newwin 1 1
del rs0 rs1
objaoper meanred row 0 rs0
objaoper meanred row 1 rs1
plot2d rs0 (n/${nFreqBin})*250+1250 val n>0 "marker=fcircle,5 blue nsta notit"
plot2d rs1 (n/${nFreqBin})*250+1250 val n>0 "same marker=fcircle,5 red nsta notit"
settitle "$source (${ncycles} cycles) Ch 0 (blue) Ch 1 (red)" ' ' $defatt
setaxelabels "Freq. (MHz)" "<(ON-OFF)[r,c]/Gain[r]>_{r,c} Rebin(${nFreqBin})" $defatt


##Reduced Mean
newwin 1 1
del rss0 rss1
objaoper sigmared row 0 rss0
objaoper sigmared row 1 rss1
plot2d rss0 (n/${nFreqBin})*250+1250 val n>0 "marker=fcircle,5 blue nsta notit"
plot2d rss1 (n/${nFreqBin})*250+1250 val n>0 "same marker=fcircle,5 red nsta notit"
settitle "$source (${ncycles} cycles) Ch 0 (blue) Ch 1 (red)" ' ' $defatt 
setaxelabels "Freq. (MHz)" "Sigma[(ON-OFF)[r,c]/Gain[r]] Rebin(${nFreqBin})" $defatt

######
else
######

openppf medianDiffOnOffRaw_${source}-${ncycles}Cycles.ppf


##General Median
del s0 s1
objaoper median row 0 s0
objaoper median row 1 s1
graphicatt "xylimits=1250,1500,-0.01,0.01"
newwin 1 1
plot2d s0 (n/8192)*250+1250 val n>0 "cpts blue nsta notit"
plot2d s1 (n/8192)*250+1250 val n>0 "same cpts red nsta notit"
settitle "$source (${ncycles} cycles) Ch 0 (blue) Ch 1 (red)" ' ' $defatt
setaxelabels "Freq. (MHz)" "Median[(ON-OFF)[r,c]/Gain[r]]"  $defatt


##Reduced Mean of the median
newwin 1 1
del rs0 rs1
objaoper meanmedred row 0 rs0
objaoper meanmedred row 1 rs1
graphicatt "xylimits=1250,1500,-0.01,0.01"
plot2d rs0 (n/${nFreqBin})*250+1250 val n>0 "marker=fcircle,5 blue nsta notit"
plot2d rs1 (n/${nFreqBin})*250+1250 val n>0 "same marker=fcircle,5 red nsta notit"
settitle "$source (${ncycles} cycles) Ch 0 (blue) Ch 1 (red)" ' ' $defatt
setaxelabels "Freq. (MHz)" "Mean Median[(ON-OFF)[r,c]/Gain[r]] Rebin(${nFreqBin})" $defatt


##Reduced Sigma of the median
newwin 1 1
del rss0 rss1
objaoper sigmamedred row 0 rss0
objaoper sigmamedred row 1 rss1
graphicatt "xylimits=1250,1500,0,0.01"
plot2d rss0 (n/${nFreqBin})*250+1250 val n>0 "marker=fcircle,5 blue nsta notit"
plot2d rss1 (n/${nFreqBin})*250+1250 val n>0 "same marker=fcircle,5 red nsta notit"
settitle "$source (${ncycles} cycles) Ch 0 (blue) Ch 1 (red)" ' ' $defatt
setaxelabels "Freq. (MHz)" "Sigma Median[(ON-OFF)[r,c]/Gain[r]] Rebin(${nFreqBin})" $defatt


##Normalized General Median
del ngms0 ngms1
objaoper medianNorm row 0 ngms0
objaoper medianNorm row 1 ngms1
graphicatt "xylimits=1250,1500,-2,2"
newwin 1 1
plot2d ngms0 (n/8192)*250+1250 val n>0 "cpts blue nsta notit"
plot2d ngms1 (n/8192)*250+1250 val n>0 "same cpts red nsta notit"
settitle "$source (${ncycles} cycles) Ch 0 (blue) Ch 1 (red)" ' ' $defatt
setaxelabels "Freq. (MHz)" "Norm Median[(ON-OFF)[r,c]/Gain[r]]"  $defatt

del prod
c++exec TVector<r_4>prod(ngms0.NElts(),BaseArray::RowVector); prod=ngms0&&ngms1;KeepObj(prod);
newwin 1 1
plot2d prod (n/8192)*250+1250 val n>0 "cpts blue nsta notit"
settitle "$source (${ncycles} cycles)" ' ' $defatt
setaxelabels "Freq. (MHz)" "Product Norm Median[(ON-OFF)[r,c]/Gain[r]]"  $defatt
 

endif