-
Notifications
You must be signed in to change notification settings - Fork 2
/
Copy pathbatch_sr_vs_mag.pro
159 lines (118 loc) · 7.84 KB
/
batch_sr_vs_mag.pro
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
;pro batch_sr_vs_mag
set3_bin1_h_0529 = obj_new('aodataset', from='20100529_045651', to='20100529_051205')
sr3_bin1_h_0529 = set3_bin1_h_0529->value('irtc.sr_se')
set3_bin1_h_0601 = obj_new('aodataset', from='20100601_034910', to='20100601_035123')
sr3_bin1_h_0601 = set3_bin1_h_0601->value('irtc.sr_se')
set3_bin1_h_0603 = obj_new('aodataset', from='20100603_035522', to='20100603_043715')
sr3_bin1_h_0603 = set3_bin1_h_0603->value('irtc.sr_se')
;set3_bin1_h->merge, set3_bin1_h_bis
set93_bin1_h_0531 = obj_new('aodataset', from='20100531_054115', to='20100531_062602')
sr93_bin1_h_0531 = set93_bin1_h_0531->value('irtc.sr_se')
set93_bin1_j_0531 = obj_new('aodataset', from='20100531_062934', to='20100531_064130')
sr93_bin1_j_0531 = set93_bin1_j_0531->value('irtc.sr_se')
set93_bin2_h_0531 = obj_new('aodataset', from='20100531_071124', to='20100531_072850')
sr93_bin2_h_0531 = set93_bin2_h_0531->value('irtc.sr_se')
set93_bin2_h_0603 = obj_new('aodataset', from='20100603_051742', to='20100603_052149')
sr93_bin2_h_0603 = set93_bin2_h_0603->value('irtc.sr_se')
set63_bin2_h_0531 = obj_new('aodataset', from='20100531_075915', to='20100531_081826')
sr63_bin2_h_0531 = set63_bin2_h_0531->value('irtc.sr_se')
set111_bin2_h_0531 = obj_new('aodataset', from='20100531_085035', to='20100531_091302')
sr111_bin2_h_0531 = set111_bin2_h_0531->value('irtc.sr_se')
set111_bin3_h_0531 = obj_new('aodataset', from='20100531_093542', to='20100531_093859')
sr111_bin3_h_0531 = set111_bin3_h_0531->value('irtc.sr_se')
set122_bin2_h_0603 = obj_new('aodataset', from='20100603_060640', to='20100603_061105')
sr122_bin2_h_0603 = set122_bin2_h_0603->value('irtc.sr_se')
set122_bin3_h_0601 = obj_new('aodataset', from='20100601_065349', to='20100601_065609')
sr122_bin3_h_0601 = set122_bin3_h_0601->value('irtc.sr_se')
set122_bin4_h_0601 = obj_new('aodataset', from='20100601_071105', to='20100601_072414')
sr122_bin4_h_0601 = set122_bin4_h_0601->value('irtc.sr_se')
set84_bin3_h_0601 = obj_new('aodataset', from='20100601_082644', to='20100601_083541')
sr84_bin3_h_0601 = set84_bin3_h_0601->value('irtc.sr_se')
set84_bin2_h_0601 = obj_new('aodataset', from='20100601_085057', to='20100601_090118')
sr84_bin2_h_0601 = set84_bin2_h_0601->value('irtc.sr_se')
set74_bin4_h_0601 = obj_new('aodataset', from='20100601_101108', to='20100601_101332')
sr74_bin4_h_0601 = set74_bin4_h_0601->value('irtc.sr_se')
set151_bin2_h_0603 = obj_new('aodataset', from='20100603_084308', to='20100603_084638')
sr151_bin2_h_0603 = set151_bin2_h_0603->value('irtc.sr_se')
set162_bin1_h_0604 = obj_new('aodataset', from='20100604_033658', to='20100604_034533')
sr162_bin1_h_0604 = set162_bin1_h_0604->value('irtc.sr_se')
set72_bin2_h_0604 = obj_new('aodataset', from='20100604_054156', to='20100604_054348')
sr72_bin2_h_0604 = set72_bin2_h_0604->value('irtc.sr_se')
set162_bin1_h_0619 = obj_new('aodataset', from='20100619_035359', to='20100619_041236')
sr162_bin1_h_0619 = set162_bin1_h_0619->value('irtc.sr_se')
set93_bin1_h_0621 = obj_new('aodataset', from='20100621_035826', to='20100621_040964')
sr93_bin1_h_0621 = set93_bin1_h_0621->value('irtc.sr_se')
set303_bin2_h_0619 = obj_new('aodataset', from='20100619_045408', to='20100619_050620')
sr303_bin2_h_0619 = set303_bin2_h_0619->value('irtc.sr_se')
set122_bin2_h_0621 = obj_new('aodataset', from='20100621_044830', to='20100621_045245')
sr122_bin2_h_0621 = set122_bin2_h_0621->value('irtc.sr_se')
set301_bin2_h_0619 = obj_new('aodataset', from='20100619_080113', to='20100619_081106')
sr301_bin2_h_0619 = set301_bin2_h_0619->value('irtc.sr_se')
set310_bin3_h_0620 = obj_new('aodataset', from='20100620_072834', to='20100620_074345')
sr310_bin3_h_0620 = set310_bin3_h_0620->value('irtc.sr_se')
set84_bin3_h_0619 = obj_new('aodataset', from='20100619_100320', to='20100619_102900')
sr84_bin3_h_0619 = set84_bin3_h_0619->value('irtc.sr_se')
set338_bin3_h_0621 = obj_new('aodataset', from='20100621_060917', to='20100621_062856')
sr338_bin3_h_0621 = set338_bin3_h_0621->value('irtc.sr_se')
set340_bin4_h_0621 = obj_new('aodataset', from='20100621_110843', to='20100621_111105')
sr340_bin4_h_0621 = set340_bin4_h_0621->value('irtc.sr_se')
set323_bin4_h_0620 = obj_new('aodataset', from='20100620_081427', to='20100620_082022')
sr323_bin4_h_0620 = set323_bin4_h_0620->value('irtc.sr_se')
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
mag_simul = findgen(8)+8.5
sr_base = [76., 73, 61, 59, 51, 33, 23, 8.8, 3.2] - 0.75 ; 50% loss due to M2 coating
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
mag3_cat = 4.9
mag93_cat = 7.1
mag162_cat = 8.5
mag63_cat = 9.4
mag303_cat = 10.7
mag111_cat = 11.4
mag122_cat = 11.5
mag151_cat = 12.5
mag301_cat = 12.5
mag310_cat = 13.0
mag72_cat = 13.5
mag84_cat = 13.6
mag74_cat = 14.4
mag338_cat = 14.5
mag340_cat = 16.5
mag323_cat = 17.5
mag_cat = [mag3_cat, mag93_cat, mag162_cat, mag63_cat, mag303_cat, mag111_cat, mag122_cat, mag151_cat, $
mag301_cat, mag310_cat, mag72_cat, mag84_cat, mag74_cat, mag338_cat, mag340_cat, mag323_cat]
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
plot, [4,18], [0,100], /nodata, xra=[4,18], yra=[0,100], charsize=2, symsi=2, thick=2, $
xtitle='R mag', ytitle='SR'
oplot, mag_simul, sr_base, psym=-1, symsi=2, thick=2
oplot, [mag3_cat], [median(sr3_bin1_h_0529)]*100, psym=-4, col=255L, symsi=2, thick=2
oplot, [mag3_cat], [median(sr3_bin1_h_0601)]*100, psym=-4, col=255L, symsi=2, thick=2
oplot, [mag3_cat], [median(sr3_bin1_h_0603)]*100, psym=-4, col=255L, symsi=2, thick=2
oplot, [mag93_cat], [median(sr93_bin1_h_0531)]*100, psym=-4, col=255L, symsi=2, thick=2
oplot, [mag93_cat], [median(sr93_bin2_h_0531)]*100, psym=-2, col=255L, symsi=2, thick=2
oplot, [mag93_cat], [median(sr93_bin2_h_0603)]*100, psym=-2, col=255L, symsi=2, thick=2
oplot, [mag63_cat], [median(sr63_bin2_h_0531)]*100, psym=-2, col=255L, symsi=2, thick=2
oplot, [mag162_cat], [median(sr162_bin1_h_0604)]*100, psym=-4, col=255L, symsi=2, thick=2
oplot, [mag111_cat], [median(sr111_bin2_h_0531)]*100, psym=-2, col=255L, symsi=2, thick=2
oplot, [mag111_cat], [median(sr111_bin3_h_0531)]*100, psym=-1, col=255L, symsi=2, thick=2
oplot, [mag122_cat], [median(sr122_bin2_h_0603)]*100, psym=-2, col=255L, symsi=2, thick=2
oplot, [mag122_cat], [median(sr122_bin3_h_0601)]*100, psym=-1, col=255L, symsi=2, thick=2
oplot, [mag122_cat], [median(sr122_bin4_h_0601)]*100, psym=-6, col=255L, symsi=2, thick=2
oplot, [mag151_cat], [median(sr151_bin2_h_0603)]*100, psym=-2, col=255L, symsi=2, thick=2
oplot, [mag72_cat], [median(sr72_bin2_h_0604)]*100, psym=-2, col=255L, symsi=2, thick=2
oplot, [mag84_cat], [median(sr84_bin3_h_0601)]*100, psym=-1, col=255L, symsi=2, thick=2
oplot, [mag84_cat], [median(sr84_bin2_h_0601)]*100, psym=-2, col=255L, symsi=2, thick=2
oplot, [mag74_cat], [median(sr74_bin4_h_0601)]*100, psym=-6, col=255L, symsi=2, thick=2
oplot, [mag151_cat], [median(sr151_bin2_h_0603)]*100, psym=-2, col=255L, symsi=2, thick=2
oplot, [mag162_cat], [median(sr162_bin1_h_0619)]*100, psym=-4, col=255L, symsi=2, thick=2
oplot, [mag93_cat], [median(sr93_bin1_h_0621)]*100, psym=-4, col=255L, symsi=2, thick=2
oplot, [mag303_cat], [median(sr303_bin2_h_0619)]*100, psym=-2, col=255L, symsi=2, thick=2
oplot, [mag122_cat], [median(sr122_bin2_h_0621)]*100, psym=-2, col=255L, symsi=2, thick=2
oplot, [mag301_cat], [median(sr301_bin2_h_0619)]*100, psym=-2, col=255L, symsi=2, thick=2
oplot, [mag310_cat], [median(sr310_bin3_h_0620)]*100, psym=-1, col=255L, symsi=2, thick=2
oplot, [mag84_cat] , [median(sr84_bin3_h_0619)]*100, psym=-1, col=255L, symsi=2, thick=2
oplot, [mag338_cat], [median(sr338_bin3_h_0621)]*100, psym=-1, col=255L, symsi=2, thick=2
oplot, [mag340_cat], [median(sr340_bin4_h_0621)]*100, psym=-6, col=255L, symsi=2, thick=2
oplot, [mag323_cat], [median(sr323_bin4_h_0620)]*100, psym=-6, col=255L, symsi=2, thick=2
legend,['bin 1','bin 2','bin 3', 'bin 4'], psym=[4,2,1,6], charsi=2, /top, /right
write_jpeg, filepath(root=ao_elabdir(), 'sr_vs_mag.jpeg'), tvrd(true=3), true=3
end