Skip to content
GitLab
Explore
Sign in
Giorgio Calderone
CL_1ES_1927p654
Compare revisions
074db87e4f190f7000529dece5bba79d5bbb6468 to d1975ede42a4e738974e7152d0fef915cec43c6d
Commits on Source (2)
Updated
· 1af8c980
Giorgio Calderone
authored
Aug 12, 2021
1af8c980
Added :floyds job
· d1975ede
Giorgio Calderone
authored
Aug 12, 2021
d1975ede
Show whitespace changes
Inline
Side-by-side
run.jl
View file @
d1975ede
...
...
@@ -124,8 +124,9 @@ function analyze_job(job; Nloop = 6)
viewer
(
res
,
showcomps
=
[
:
qso_cont
,
:
galaxy
,
:
balmer
],
filename
=
"
$
(opt.path)/results_
$(job)
_
$(loop)
.html"
)
viewer
(
res
,
showcomps
=
[
:
qso_cont
,
:
galaxy
,
:
balmer
],
filename
=
"
$
(opt.path)/results_
$(job)
_
$(loop)
_rebin4.html"
,
rebin
=
4
)
else
@info
file
res
=
JLD2
.
load_object
(
file
);
println
(
i
loop
,
" "
,
res
.
fitres
.
gofstat
/
res
.
fitres
.
dof
,
" "
,
res
.
fitres
.
elapsed
/
3600.
)
println
(
loop
,
" "
,
res
.
fitres
.
gofstat
/
res
.
fitres
.
dof
,
" "
,
res
.
fitres
.
elapsed
/
3600.
)
end
OIII_norm
=
Float64
[]
...
...
@@ -200,7 +201,8 @@ epoch_ids = Dict(
:
highres
=>
[
5
,
7
,
13
,
17
,
20
],
:
test
=>
[
7
,
13
],
:
limited
=>
[
1
,
2
,
3
,
13
],
:
all
=>
collect
(
5
:
26
))
:
all
=>
collect
(
5
:
26
),
:
floyds
=>
[
6
,
8
,
9
,
10
,
11
,
12
,
14
,
15
,
16
,
18
,
19
,
21
,
22
,
23
])
# deleteat!(epoch_ids[:all], 4) #insufficient coverage in na_Hg
# Can't use resolution smaller than 150 km / s otherwise some line
...
...
@@ -212,6 +214,7 @@ resolution = Dict()
analyze_job
(
:
test
)
analyze_job
(
:
highres
)
analyze_job
(
:
all
)
analyze_job
(
:
floyds
)
ddd
...
...
@@ -238,38 +241,43 @@ function cont_at(bestfit, λ)
return
λ
*
A
*
B
,
λ
*
(
dd
[
2
]
-
dd
[
1
])
end
# Wrap-up
chosen_job
=
:
all
chosen_iloop
=
5
res
=
JLD2
.
load_object
(
"output/results_
$(chosen_job)
_
$(chosen_iloop)
.dat"
)
tab
=
DataFrame
(
epoch
=
Int
[],
date
=
String
[],
instr
=
String
[],
galaxy
=
Float64
[],
contslope
=
Float64
[],
c3500
=
Float64
[],
c5100
=
Float64
[],
bb_Ha
_norm
=
Float64
[],
bb_Ha
_fwhm
=
Float64
[],
bb_Ha
_voff
=
Float64
[],
br_Ha
_norm
=
Float64
[],
br_Ha
_fwhm
=
Float64
[],
br_Ha
_voff
=
Float64
[],
n
a_
H
a_norm
=
Float64
[],
n
a_
H
a_fwhm
=
Float64
[],
n
a_
H
a_voff
=
Float64
[],
b
b_
H
b_norm
=
Float64
[],
b
b_
H
b_fwhm
=
Float64
[],
b
b_
H
b_voff
=
Float64
[],
br_Hb
_norm
=
Float64
[],
br_Hb
_fwhm
=
Float64
[],
br_Hb
_voff
=
Float64
[],
na_Hb
_norm
=
Float64
[],
na_Hb
_fwhm
=
Float64
[],
na_Hb
_voff
=
Float64
[],
Ha_bb
_norm
=
Float64
[],
Ha_bb
_fwhm
=
Float64
[],
Ha_bb
_voff
=
Float64
[],
Ha_br
_norm
=
Float64
[],
Ha_br
_fwhm
=
Float64
[],
Ha_br
_voff
=
Float64
[],
H
a_
n
a_norm
=
Float64
[],
H
a_
n
a_fwhm
=
Float64
[],
H
a_
n
a_voff
=
Float64
[],
H
b_
b
b_norm
=
Float64
[],
H
b_
b
b_fwhm
=
Float64
[],
H
b_
b
b_voff
=
Float64
[],
Hb_br
_norm
=
Float64
[],
Hb_br
_fwhm
=
Float64
[],
Hb_br
_voff
=
Float64
[],
Hb_na
_norm
=
Float64
[],
Hb_na
_fwhm
=
Float64
[],
Hb_na
_voff
=
Float64
[],
oiii_fwhm
=
Float64
[],
oiii_voff
=
Float64
[])
for
id
in
1
:
length
(
res
.
bestfit
.
models
)
for
id
in
1
:
length
(
res
.
multi
)
push!
(
tab
,
(
id
,
epoch_
filename
s
[
chosen_
epochs
[
id
]][
27
:
end
-
4
],
""
,
res
.
bestfit
[
id
][
:
galaxy
]
.
norm
.
val
,
res
.
bestfit
[
id
][
:
qso_cont
]
.
alpha
.
val
,
cont_at
(
res
.
bestfit
[
id
],
3500
)[
1
],
cont_at
(
res
.
bestfit
[
id
],
5100
)[
1
],
res
.
bestfit
[
id
][
:
bb_Ha
]
.
norm
.
patched
,
res
.
bestfit
[
id
][
:
bb_Ha
]
.
fwhm
.
patched
,
res
.
bestfit
[
id
][
:
bb_Ha
]
.
voff
.
patched
,
res
.
bestfit
[
id
][
:
br_Ha
]
.
norm
.
patched
,
res
.
bestfit
[
id
][
:
br_Ha
]
.
fwhm
.
patched
,
res
.
bestfit
[
id
][
:
br_Ha
]
.
voff
.
patched
,
res
.
bestfit
[
id
][
:
n
a_
H
a
]
.
norm
.
patched
,
res
.
bestfit
[
id
][
:
n
a_
H
a
]
.
fwhm
.
patched
,
res
.
bestfit
[
id
][
:
n
a_
H
a
]
.
voff
.
patched
,
res
.
bestfit
[
id
][
:
b
b_
H
b
]
.
norm
.
patched
,
res
.
bestfit
[
id
][
:
b
b_
H
b
]
.
fwhm
.
patched
,
res
.
bestfit
[
id
][
:
b
b_
H
b
]
.
voff
.
patched
,
res
.
bestfit
[
id
][
:
br_Hb
]
.
norm
.
patched
,
res
.
bestfit
[
id
][
:
br_Hb
]
.
fwhm
.
patched
,
res
.
bestfit
[
id
][
:
br_Hb
]
.
voff
.
patched
,
res
.
bestfit
[
id
][
:
na_Hb
]
.
norm
.
patched
,
res
.
bestfit
[
id
][
:
na_Hb
]
.
fwhm
.
patched
,
res
.
bestfit
[
id
][
:
na_Hb
]
.
voff
.
patched
,
res
.
bestfit
[
id
][
:
OIII_5007
]
.
fwhm
.
patched
,
res
.
bestfit
[
id
][
:
OIII_5007
]
.
voff
.
patched
))
(
epoch_
id
s
[
chosen_
job
][
id
],
string
(
split
(
res
.
source
.
specs
[
id
]
.
label
,
"/"
)[
2
])[
18
:
end
-
4
],
""
,
res
.
multi
[
id
][
:
galaxy
]
.
norm
.
val
,
res
.
multi
[
id
][
:
qso_cont
]
.
alpha
.
val
,
cont_at
(
res
.
multi
[
id
],
3500
)[
1
],
cont_at
(
res
.
multi
[
id
],
5100
)[
1
],
res
.
multi
[
id
][
:
Ha_bb
]
.
norm
.
patched
,
res
.
multi
[
id
][
:
Ha_bb
]
.
fwhm
.
patched
,
res
.
multi
[
id
][
:
Ha_bb
]
.
voff
.
patched
,
res
.
multi
[
id
][
:
Ha_br
]
.
norm
.
patched
,
res
.
multi
[
id
][
:
Ha_br
]
.
fwhm
.
patched
,
res
.
multi
[
id
][
:
Ha_br
]
.
voff
.
patched
,
res
.
multi
[
id
][
:
H
a_
n
a
]
.
norm
.
patched
,
res
.
multi
[
id
][
:
H
a_
n
a
]
.
fwhm
.
patched
,
res
.
multi
[
id
][
:
H
a_
n
a
]
.
voff
.
patched
,
res
.
multi
[
id
][
:
H
b_
b
b
]
.
norm
.
patched
,
res
.
multi
[
id
][
:
H
b_
b
b
]
.
fwhm
.
patched
,
res
.
multi
[
id
][
:
H
b_
b
b
]
.
voff
.
patched
,
res
.
multi
[
id
][
:
Hb_br
]
.
norm
.
patched
,
res
.
multi
[
id
][
:
Hb_br
]
.
fwhm
.
patched
,
res
.
multi
[
id
][
:
Hb_br
]
.
voff
.
patched
,
res
.
multi
[
id
][
:
Hb_na
]
.
norm
.
patched
,
res
.
multi
[
id
][
:
Hb_na
]
.
fwhm
.
patched
,
res
.
multi
[
id
][
:
Hb_na
]
.
voff
.
patched
,
res
.
multi
[
id
][
:
OIII_5007
]
.
fwhm
.
patched
,
res
.
multi
[
id
][
:
OIII_5007
]
.
voff
.
patched
))
end
tab
[
!
,
:
l3500
]
.=
[
3500
.*
Spline1D
(
domain
(
res
.
multi
[
id
])[
:
],
res
.
multi
[
id
]())(
3500.
)
for
id
in
1
:
length
(
chosen_epochs
)]
tab
[
!
,
:
l5100
]
.=
[
5100
.*
Spline1D
(
domain
(
res
.
multi
[
id
])[
:
],
res
.
multi
[
id
]())(
5100.
)
for
id
in
1
:
length
(
chosen_epochs
)]
tab
[
!
,
:
Ha
]
.=
tab
.
br_Ha
_norm
.+
tab
.
bb_Ha
_norm
.+
tab
.
n
a_
H
a_norm
tab
[
!
,
:
Hb
]
.=
tab
.
br_Hb
_norm
.+
tab
.
b
b_
H
b_norm
.+
tab
.
na_Hb
_norm
tab
[
!
,
:
l3500
]
.=
[
3500
.*
Spline1D
(
domain
(
res
.
multi
[
id
])[
:
],
res
.
multi
[
id
]())(
3500.
)
for
id
in
1
:
length
(
res
.
multi
)]
tab
[
!
,
:
l5100
]
.=
[
5100
.*
Spline1D
(
domain
(
res
.
multi
[
id
])[
:
],
res
.
multi
[
id
]())(
5100.
)
for
id
in
1
:
length
(
res
.
multi
)]
tab
[
!
,
:
Ha
]
.=
tab
.
Ha_br
_norm
.+
tab
.
Ha_bb
_norm
.+
tab
.
H
a_
n
a_norm
tab
[
!
,
:
Hb
]
.=
tab
.
Hb_br
_norm
.+
tab
.
H
b_
b
b_norm
.+
tab
.
Hb_na
_norm
dd
=
Date
.
(
tab
.
date
,
Ref
(
"yyyymmdd"
))
tab
[
!
,
:
day
]
.=
getproperty
.
(
dd
-
Date
(
"2018-03-06"
)
.+
Day
(
72
),
:
value
)
...
...
@@ -286,7 +294,7 @@ for i in 1:nrow(tab)
tab
[
i
,
:
instr
]
=
tmp
[
j
,
:
Inst
]
end
f
=
FITS
(
"
$
(opt.path)/1ES_1927p654_results
_
$(job)
.fits"
,
"w"
)
f
=
FITS
(
"
$
(opt.path)/1ES_1927p654_results.fits"
,
"w"
)
write
(
f
,
tab
)
close
(
f
)
...
...
@@ -326,14 +334,14 @@ tab.pt .= getindex.(Ref(cm), tab.instr)
@gp
"set grid"
for
i
in
1
:
nrow
(
tab
)
ss
=
"w l t '
$
(tab[i, :date])' dt
$
(tab[i, :pt]) lw 2"
@gp
:-
domain
(
res
.
m
odel
[
i
])[
:
]
res
.
m
odel
[
i
](
:
qso_cont
)
ss
@gp
:-
domain
(
res
.
m
ulti
[
i
])[
:
]
res
.
m
ulti
[
i
](
:
qso_cont
)
ss
end
save
(
:
prenorm
,
term
=
"png size 800,600"
,
output
=
"
$
(opt.path)/evolution.png"
)
@gp
"set grid"
for
i
in
1
:
length
(
chosen_epochs
)
println
(
res
.
bestfit
[
i
][
:
OIII_5007
])
for
i
in
1
:
length
(
res
.
multi
)
println
(
res
.
multi
[
i
][
:
OIII_5007
])
x
=
domain
(
res
.
multi
[
i
])[
:
]
y
=
res
.
multi
[
i
]()
y0
=
y
.-
res
.
multi
[
i
](
:
OIII_5007
)
...
...