-
Notifications
You must be signed in to change notification settings - Fork 1
/
analysis_wip.pro
82 lines (68 loc) · 3.01 KB
/
analysis_wip.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
pro analysis_wip
; ################################################################
; Produce the maps for each P1-P4 in >2keV, 2-4 keV and 4-6 keV
for pid=1,4 do make_map_nov14o4,pid=pid,erang=2
for pid=1,4 do make_map_nov14o4,pid=pid,erang=[2,4]
for pid=1,4 do make_map_nov14o4,pid=pid,erang=[4,6]
; ################################################################
; Figure 1
plot_livetime_orbit4_hsi
plot_maps_nov14o4
; ################################################################
; Figure 2
; need the GOES/SXI and SDO/AIA maps not included
plot_aiasxi_maps_nov14o4
; ################################################################
; ################################################################
; ################################################################
; Spectra and livetime - analysis and some playing about with the grades/pile-up
; Firstly make sure the *.dat with the regions of interest co-ords are saved out
make_roi_nov14
; Then make the grade histograms per fov and region
; need IDL >8.0 for the plots to work
make_gradehist_nov14o4,pid=1,fid='A'
make_gradehist_nov14o4,pid=1,fid='B'
make_gradehist_nov14o4,pid=2,fid='A'
make_gradehist_nov14o4,pid=2,fid='B'
; Then make the spectra per fov and regions
; again need IDL >8.0 for the plots to work
make_specs_nov14o4_de,pid=1,fid='A'
make_specs_nov14o4_de,pid=1,fid='B'
make_specs_nov14o4_de,pid=2,fid='A'
make_specs_nov14o4_de,pid=2,fid='B'
; Then test to see how much of the spectrum is piled-up
specs_pileup_test,pid=1,fid='A'
specs_pileup_test,pid=1,fid='B'
specs_pileup_test,pid=2,fid='A'
specs_pileup_test,pid=2,fid='B'
; Make the srm for ospex
regnm=['D1','D2','L1','L2','L3']
for i=0, n_elements(regnm)-1 do make_ns_srm,pid=1,fid='A',regnm=regnm[i]
for i=0, n_elements(regnm)-1 do make_ns_srm,pid=1,fid='B',regnm=regnm[i]
make_ns_srm,pid=2,fid='A',regnm='D1'
make_ns_srm,pid=2,fid='A',regnm='L1'
make_ns_srm,pid=2,fid='B',regnm='D1'
make_ns_srm,pid=2,fid='B',regnm='L1'
; Fit an isothermal model to the spectrum
for i=0, 4 do fitvth_spec_nov14o4,pid=1,fid='A',regid=i
for i=0, 4 do fitvth_spec_nov14o4,pid=1,fid='B',regid=i
fitvth_spec_nov14o4,pid=2,fid='A',regid=0
fitvth_spec_nov14o4,pid=2,fid='A',regid=2
fitvth_spec_nov14o4,pid=2,fid='B',regid=0
fitvth_spec_nov14o4,pid=2,fid='B',regid=2
; ################################################################
; ################################################################
; ################################################################
; Figure 3
plot_fitted_specs_p1
; ################################################################
; Figure 4
plot_fitted_specs_p2
; ################################################################
; ################################################################
; ################################################################
; Default work out the limits
for i=0, 4 do mc_2ndvth_uplims,regid=i
for i=0, 4 do mc_2ndvth_uplims,regid=i,lime=6.0
stop
end