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

set source $1
set ncycles $2
set nFreqBin $3
openppf meanDiffOnOffRaw_${source}-${ncycles}Cycles.ppf


##General Mean
del s0 s1
objaoper mean row 0 s0
objaoper mean 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 "<(ON-OFF)[r,c]/Gain[r]>_{r,c} $source (${ncycles} cycles) Ch 0 (blue) Ch 1 (red)"


##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 "<(ON-OFF)[r,c]/Gain[r]>_{r,c} Rebin(${nFreqBin}) $source (${ncycles} cycles) Ch 0 (blue) Ch 1 (red)"

##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 "Sigma[(ON-OFF)[r,c]/Gain[r]] Rebin(${nFreqBin}) $source (${ncycles} cycles) Ch 0 (blue) Ch 1 (red)"