Changeset 552 for BAORadio/AmasNancay/trunk/plotDiffOnOff.pic

Timestamp:

Oct 5, 2011, 3:13:36 PM (13 years ago)

Author:

campagne

Message:

improvment (jec)

File:

• BAORadio/AmasNancay/trunk/plotDiffOnOff.pic (modified) (4 diffs)

BAORadio/AmasNancay/trunk/plotDiffOnOff.pic

@@ -6 +6 @@
 set ncycles $2
 set nFreqBin $3
+set action $4
+
+set defatt "font=helvetica,bold,20"
+
+
+
+if ( ${action} == "mean" ) then
+
 openppf meanDiffOnOffRaw_${source}-${ncycles}Cycles.ppf
 
@@ -13 +21 @@
 objaoper mean row 0 s0
 objaoper mean row 1 s1
-graphicatt "xylimits=1250,1500,-0.01,0.01"
+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 "<(ON-OFF)[r,c]/Gain[r]>_{r,c} $source (${ncycles} cycles) Ch 0 (blue) Ch 1 (red)"
+settitle "$source (${ncycles} cycles) Ch 0 (blue) Ch 1 (red)" ' ' $defatt
+setaxelabels "Freq. (MHz)" "<(ON-OFF)[r,c]/Gain[r]>_{r,c}" $defatt
+
 
 
@@ -27 +37 @@
 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)"
+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
@@ -36 +48 @@
 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)"
+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
+ 
+
+endif
+