Changeset 611 for BAORadio/AmasNancay/trunk/etude_mergeAna.pic

Timestamp:

Nov 28, 2011, 10:04:36 AM (13 years ago)

Author:

torrento

Message:

Add pic analysis files and emptydir function for proc_specmfib

File:

• BAORadio/AmasNancay/trunk/etude_mergeAna.pic (modified) (12 diffs)

BAORadio/AmasNancay/trunk/etude_mergeAna.pic

@@ +1 @@
+
+set path "/sps/baoradio/AmasNancay/AST"
 
 set source  $1
@@ -9 +11 @@
 
 set fcalib 1410
-set f1420low 1420.2
-set f1420high 1420.6
-
-
+
+# Abell 85, 500 cycles
+if ( ${source} == "abell85" ) then
+  set f1420low 1420.2
+  set f1420high 1420.6
+  set packVal "1 10 25 50 100"
+endif
+
+# Abell 1205, 600 cycles
+if ( ${source} == "abell1205" ) then
+  set f1420low 1420.1
+  set f1420high 1420.6
+  set packVal "1 10 25 50 100 120"
+endif
+
+# Abell 2440, 210 cycles
+if ( ${source} == "abell2440" ) then
+  set f1420low 1420.4
+  set f1420high 1420.7
+  set packVal "1 10 25 42"
+endif
+
+echo "f1420low=${f1420low}, f1420high=${f1420high}"
+echo "Packing values="${packVal}
 
 nbins1420 = floor((${f1420high}-${f1420low})/250*8192)
-
 echo $nbins1420
-
 
 clearscript diffOnOff
@@ -26 +46 @@
 clearscript sigmaRaw1400a1420
 clearscript evolAt1410a1415OffOn
+
 ############################################################################
 defscript evolAtCalibFreq
@@ -90 +111 @@
 
 
-graphicatt "xylimits=0,500,-0.01,0.01"
+graphicatt "xylimits=0,${ncycles},-0.02,0.02"
 
 newwin 1 2
@@ -159 +180 @@
 line2vec vecCh1All $linCh1All
 
-#packing values only valid for 500 cycles...
-set packVal "1 10 25 50 100"
 set sigmaCh0 ""
 set errsigCh0 ""
 set sigmaCh1 ""
 set errsigCh1 ""
-
 
 foreach ipack ( $packVal )
@@ -236 +254 @@
 line2vec vecCh1All $linCh1All
 
-#packing values only valid for 500 cycles...
-set packVal "1 10 25 50 100"
 set sigmaCh0 ""
 set errsigCh0 ""
@@ -329 +345 @@
 
 
-graphicatt "xylimits=0,500,-0.01,0.01"
+graphicatt "xylimits=0,${ncycles},-0.01,0.01"
 
 newwin 1 2
@@ -394 +410 @@
 
 
-graphicatt "xylimits=0,500,-0.1,0.1"
+graphicatt "xylimits=0,${ncycles},-0.1,0.1"
 
 newwin 1 2
@@ -487 +503 @@
 setaxelabels "Freq. (MHz)" "I (a.u)" $axedefatt
 
-
-
 graphicatt "xylimits=1250,1500,-0.01,0.01"
 del snocal0 snocal1
@@ -499 +513 @@
 setaxelabels "Freq. (MHz)" "I (a.u)" $axedefatt
 
-del sruncal0 sruncal1
-objaoper meanPerRunCalib row 0 sruncal0
-objaoper meanPerRunCalib row 1 sruncal1
-newwin 1 1
-plot2d sruncal0 (n/8192)*250+1250 val n>0 "blue cpts notit nsta"
-plot2d sruncal1 (n/8192)*250+1250 val n>0 "red same cpts notit nsta"
-settitle "ON-OFF Calib Run  ${source} Ch 0 (blue) Ch 1 (red) ${ncycles}cycles" ' ' $defatt
-setaxelabels "Freq. (MHz)" "I (a.u)" $axedefatt
-
-newwin 1 1
-del scyclecal0 scyclecal1
-objaoper meanPerCycleCalib row 0 scyclecal0
-objaoper meanPerCycleCalib row 1 scyclecal1
-plot2d scyclecal0 (n/8192)*250+1250 val n>0 "blue cpts notit nsta"
-plot2d scyclecal1 (n/8192)*250+1250 val n>0 "red same cpts notit nsta"
-settitle "ON-OFF Calib Cycle  ${source} Ch 0 (blue) Ch 1 (red) ${ncycles}cycles" ' ' $defatt
-setaxelabels "Freq. (MHz)" "I (a.u)" $axedefatt
+#del sruncal0 sruncal1
+#objaoper meanPerRunCalib row 0 sruncal0
+#objaoper meanPerRunCalib row 1 sruncal1
+#newwin 1 1
+#plot2d sruncal0 (n/8192)*250+1250 val n>0 "blue cpts notit nsta"
+#plot2d sruncal1 (n/8192)*250+1250 val n>0 "red same cpts notit nsta"
+#settitle "ON-OFF Calib Run  ${source} Ch 0 (blue) Ch 1 (red) ${ncycles}cycles" ' ' $defatt
+#setaxelabels "Freq. (MHz)" "I (a.u)" $axedefatt
+
+#newwin 1 1
+#del scyclecal0 scyclecal1
+#objaoper meanPerCycleCalib row 0 scyclecal0
+#objaoper meanPerCycleCalib row 1 scyclecal1
+#plot2d scyclecal0 (n/8192)*250+1250 val n>0 "blue cpts notit nsta"
+#plot2d scyclecal1 (n/8192)*250+1250 val n>0 "red same cpts notit nsta"
+#settitle "ON-OFF Calib Cycle  ${source} Ch 0 (blue) Ch 1 (red) ${ncycles}cycles" ' ' $defatt
+#setaxelabels "Freq. (MHz)" "I (a.u)" $axedefatt
 
 
@@ -526 +540 @@
 #use minuscule for source label
 
-
-
-
-openppf onoffsurvey_${source}-${ncycles}Cycles.ppf
+openppf ${path}/onoffsurvey_${source}-${ncycles}Cycles.ppf
 print onoffevol
+
+#openppf ${path}/rawOnOffDiff_${source}-${ncycles}Cycles.ppf
+#ls
 
 set defatt "font=helvetica,bold,20 fixedfonsize"
@@ -546 +560 @@
 #evolAt1420SideFreq
 #IntAt1420
+#sigmaRaw1420Side
+#sigmaRaw1400a1420
+#evolAt1400a1420OffOn