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

Timestamp:

Oct 7, 2011, 1:42:23 PM (13 years ago)

Author:

campagne

Message:

generalise the process to download the files for final analysis (jec)

File:

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

BAORadio/AmasNancay/trunk/plotDiffOnOff.pic

@@ -51 +51 @@
 setaxelabels "Freq. (MHz)" "Sigma[(ON-OFF)[r,c]/Gain[r]] Rebin(${nFreqBin})" $defatt
 
+######
 else
+######
 
 openppf medianDiffOnOffRaw_${source}-${ncycles}Cycles.ppf
@@ -93 +95 @@
 
 ##Normalized General Median
-del s0 s1
-objaoper medianNorm row 0 s0
-objaoper medianNorm row 1 s1
+del ngms0 ngms1
+objaoper medianNorm row 0 ngms0
+objaoper medianNorm row 1 ngms1
 graphicatt "xylimits=1250,1500,-2,2"
 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"
+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