Changeset 556 for BAORadio/AmasNancay/trunk/etude_onoff.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/etude_onoff.pic (modified) (5 diffs)

BAORadio/AmasNancay/trunk/etude_onoff.pic

@@ -11 +11 @@
 #Plots ON-OFF with all calibrations 
 #####################################
-openppf ${toppath}/${source}/${date}${srclower}/diffOnOff_${date}_${srclower}.ppf
+openppf ${toppath}/${source}/diffOnOff_${date}_${srclower}.ppf
 
 
@@ -34 +34 @@
 settitle "ON-OFF  ${source}/${date} Ch 0 (blue) Ch 1 (red) Mean (black)  All cycles"
 setaxelabels "Freq. (MHz)" "I (Jy)"
+
+#####################################
+#Examine fluctuations sigma for the combination of all cycles
+#####################################
+cp sallcyclev0 specv0
+cp sallcyclev1 specv1
+c++compile reducing        
+c++link reducing.so doreduc
+call doreduc
+
+newwin 1 1
+graphicatt "xylimits=1250,1500,-0.01,0.01"
+ del redMeanv0                                                                      
+ del redMeanv1                                                                      
+ objaoper reducedMeanDiffOnOff row 0 redMeanv0                                      
+ objaoper reducedMeanDiffOnOff row 1 redMeanv1                                      
+ plot2d redMeanv0 (n/32)*250+1250 val n>0 "blue marker=fcircle,7 nsta notit"        
+ plot2d redMeanv1 (n/32)*250+1250 val n>0 "same red marker=fcircle,7 nsta notit"
+settitle "Red Mean ON-OFF  ${source} Ch 0 (blue) Ch 1 (red) ${tmpvec.sum} cycles"
+setaxelabels "Freq. (MHz)" "I (Jy)"
+
+newwin 1 1
+graphicatt "xylimits=1250,1500,0,0.001"
+ del redSigmav0                                                                      
+ del redSigmav1                                                                      
+ objaoper reducedSigmaDiffOnOff row 0 redSigmav0                                      
+ objaoper reducedSigmaDiffOnOff row 1 redSigmav1                                      
+ plot2d redSigmav0 (n/32)*250+1250 val n>0 "blue marker=fcircle,7 nsta notit"        
+ plot2d redSigmav1 (n/32)*250+1250 val n>0 "same red marker=fcircle,7 nsta notit"
+settitle "ON-OFF  ${source} Ch 0 (blue) Ch 1 (red) ${tmpvec.sum} cycles"
+setaxelabels "Freq. (MHz)" "I (Jy)"
+
 
 #####################################
@@ -70 +102 @@
 # Plots RAW no calibration at all
 #####################################
-openppf ${toppath}/${source}/${date}${srclower}/dataRaw_${date}_${srclower}.ppf
+openppf ${toppath}/${source}/dataRaw_${date}_${srclower}.ppf
 newwin 2 2
 graphicatt "xylimits=1250,1500,0.95,1.05"
@@ -141 +173 @@
 # Plots BAO calibration included
 #####################################
-openppf ${toppath}/${source}/${date}${srclower}/dataBAOCalib_${date}_${srclower}.ppf
+openppf ${toppath}/${source}/dataBAOCalib_${date}_${srclower}.ppf
 newwin 2 2
 #graphicatt "xylimits=1250,1500,0.95,1.05"
@@ -214 +246 @@
 # Plots BAO+RT calibration included
 #####################################
-openppf ${toppath}/${source}/${date}${srclower}/dataBAORTCalib_${date}_${srclower}.ppf
+openppf ${toppath}/${source}/dataBAORTCalib_${date}_${srclower}.ppf
 newwin 2 2
 #graphicatt "xylimits=1250,1500,0.95,1.05"