source: BAORadio/AmasNancay/v3/etude_mergeAna.pic @ 564

Last change on this file since 564 was 563, checked in by campagne, 13 years ago

improve mean ON-OFF with Calib BAO (jec)

File size: 2.9 KB
Line 
1
2set source  $1
3set ncycles $2
4
5openppf onoffsurvey_${source}-${ncycles}Cycles.ppf
6print onoffevol
7
8set defatt "font=helvetica,bold,20 fixedfonsize"
9set axedefatt  "font=helvetica,bold,30 fixedfonsize grid"
10graphicatt "$defatt"
11setaxesatt "$axedefatt"
12
13newwin 1 2
14n/pl onoffevol.onoffRaw0%cycle ! ! "blue cpts notit nsta"
15settitle "Raw ON-OFF  ${source} Ch 0 (blue) Ch 1 (red) ${ncycles}cycles " ' ' $defatt
16n/pl onoffevol.onoffRaw1%cycle ! ! "red cpts notit nsta"
17setaxelabels "cycle" "I[1346MHz/6.25MHz] (a.u)" $axedefatt
18
19newwin 1 2
20n/proj hRaw0 onoffevol.onoffRaw0 ! ! ! "blue cpts notit"
21settitle "Raw ON-OFF  ${source} Ch 0 (blue) Ch 1 (red) ${ncycles}cycles " ' ' $defatt
22n/proj hRaw1 onoffevol.onoffRaw1 ! ! ! "red cpts notit"
23
24
25
26       
27newwin 1 2
28n/pl onoffevol.onoffRun0%cycle ! ! "blue cpts notit nsta"
29settitle "ON-OFF Calib Run ${source} Ch 0 (blue) Ch 1 (red) ${ncycles}cycles" ' ' $defatt
30n/pl onoffevol.onoffRun1%cycle ! ! "red cpts notit nsta"
31setaxelabels "cycle" "I[1346MHz/6.25MHz] (a.u)" $axedefatt
32
33newwin 1 2
34n/proj hRun0 onoffevol.onoffRun0 ! ! ! "blue cpts notit"
35settitle "ON-OFF Calib Run ${source} Ch 0 (blue) Ch 1 (red) ${ncycles}cycles" ' ' $defatt
36n/proj hRun1 onoffevol.onoffRun1 ! ! ! "red cpts notit"
37
38
39newwin 1 2
40n/pl onoffevol.onoffCycle0%cycle ! ! "blue cpts notit nsta"
41settitle "ON-OFF Calib Cycle ${source} Ch 0 (blue) Ch 1 (red) ${ncycles}cycles" ' ' $defatt
42n/pl onoffevol.onoffCycle1%cycle ! ! "red cpts notit nsta"
43setaxelabels "cycle" "I[1346MHz/6.25MHz] (a.u)" $axedefatt
44
45newwin 1 2
46n/proj hCycle0 onoffevol.onoffCycle0 ! ! ! "blue cpts notit"
47settitle "ON-OFF Calib Cycle ${source} Ch 0 (blue) Ch 1 (red) ${ncycles}cycles" ' ' $defatt
48n/proj hCycle1 onoffevol.onoffCycle1 ! ! ! "red cpts notit"
49
50
51
52graphicatt "xylimits=1250,1500,-0.01,0.01"
53del snocal0 snocal1
54objaoper meanNoCalib row 0 snocal0
55objaoper meanNoCalib row 1 snocal1
56newwin 1 1
57plot2d snocal0 (n/8192)*250+1250 val n>0 "blue cpts notit nsta"
58plot2d snocal1 (n/8192)*250+1250 val n>0 "red same cpts notit nsta"
59settitle "Raw ON-OFF  ${source} Ch 0 (blue) Ch 1 (red) ${ncycles}cycles" ' ' $defatt
60setaxelabels "Freq. (MHz)" "I (a.u)" $axedefatt
61
62
63
64del sruncal0 sruncal1
65objaoper meanPerRunCalib row 0 sruncal0
66objaoper meanPerRunCalib row 1 sruncal1
67newwin 1 1
68plot2d sruncal0 (n/8192)*250+1250 val n>0 "blue cpts notit nsta"
69plot2d sruncal1 (n/8192)*250+1250 val n>0 "red same cpts notit nsta"
70settitle "ON-OFF Calib Run  ${source} Ch 0 (blue) Ch 1 (red) ${ncycles}cycles" ' ' $defatt
71setaxelabels "Freq. (MHz)" "I (a.u)" $axedefatt
72
73newwin 1 1
74del scyclecal0 scyclecal1
75objaoper meanPerCycleCalib row 0 scyclecal0
76objaoper meanPerCycleCalib row 1 scyclecal1
77plot2d scyclecal0 (n/8192)*250+1250 val n>0 "blue cpts notit nsta"
78plot2d scyclecal1 (n/8192)*250+1250 val n>0 "red same cpts notit nsta"
79settitle "ON-OFF Calib Cycle  ${source} Ch 0 (blue) Ch 1 (red) ${ncycles}cycles" ' ' $defatt
80setaxelabels "Freq. (MHz)" "I (a.u)" $axedefatt
Note: See TracBrowser for help on using the repository browser.