@@ -18,14 +18,14 @@ begin
1818 erange = 10 .^ collect (range (log10 (4 ), log10 (20.0 ), 100 ))
1919 model2 = GaussianLine () + GaussianLine (μ = FitParam (10.0 ))
2020 plot (
21- erange[1 : (end - 1 )],
21+ erange[1 : (end - 1 )],
2222 @time (invokemodel (erange, model)),
2323 marker = :o ,
2424 markersize = 4 ,
2525 label = " XS" ,
2626 )
2727 plot! (
28- erange[1 : (end - 1 )],
28+ erange[1 : (end - 1 )],
2929 @time (invokemodel (erange, model2)),
3030 marker = :o ,
3131 markersize = 4 ,
@@ -40,14 +40,14 @@ model = XS_PowerLaw(a = FitParam(2.0), K = FitParam(1.0))
4040emids = push! (copy (resp. channel_bins_low), last (resp. channel_bins_high))
4141evald = push! (copy (resp. bins_low), last (resp. bins_high))
4242out = resp. matrix * invokemodel (evald, model)
43- plot (emids[1 : (end - 1 )], out)
43+ plot (emids[1 : (end - 1 )], out)
4444
4545
46- sims = @time simulate (model, resp, anc; exposure_time = 1e3 , seed = 42 )
47- plot (sims, xscale = :log10 , yscale = :log10 , xlims = (1.0 , 30.0 ), ylims = (1e -6 , 100.0 ))
46+ sims = @time simulate (model, resp, anc; exposure_time = 1.0e3 , seed = 42 )
47+ plot (sims, xscale = :log10 , yscale = :log10 , xlims = (1.0 , 30.0 ), ylims = (1.0e -6 , 100.0 ))
4848
4949begin
50- sims = @time simulate (model, resp, nothing ; exposure_time = 1e5 , seed = 42 )
50+ sims = @time simulate (model, resp, nothing ; exposure_time = 1.0e5 , seed = 42 )
5151 normalize! (faked)
5252 plot (
5353 faked,
5858 )
5959 plot! (sims)
6060 plot! (
61- emids[1 : (end - 1 )] .+ diff (emids) ./ 2 ,
61+ emids[1 : (end - 1 )] .+ diff (emids) ./ 2 ,
6262 out .* 1.4 ,
6363 xscale = :log10 ,
6464 label = " Hand folded" ,
@@ -71,18 +71,18 @@ result = @time fit(prob, LevenbergMarquadt())
7171plot! (faked, result)
7272
7373
74- sims = @time simulate (model, resp, anc; exposure_time = 1e3 , seed = 42 )
74+ sims = @time simulate (model, resp, anc; exposure_time = 1.0e3 , seed = 42 )
7575plot (
7676 sims,
7777 xscale = :log10 ,
7878 # yscale = :log10,
7979 xlims = (1.0 , 10.0 ),
80- ylims = (1e -6 , 40.0 ),
80+ ylims = (1.0e -6 , 40.0 ),
8181)
8282
8383model = GaussianLine (μ = FitParam (7.0 )) + PowerLaw (a = FitParam (0.2 ))
8484
85- sims = @time simulate (model, resp, nothing ; exposure_time = 1e5 , seed = 42 )
85+ sims = @time simulate (model, resp, nothing ; exposure_time = 1.0e5 , seed = 42 )
8686
8787plot (
8888 sims,
0 commit comments