BH fit example

In [1]:

from bh_molecule.instruments.vis133m import Vis133M
from bh_molecule.physics import BHModel
from bh_molecule.fit import BHFitter
from bh_molecule import use_dark
use_dark()

from bh_molecule.instruments.vis133m import Vis133M from bh_molecule.physics import BHModel from bh_molecule.fit import BHFitter from bh_molecule import use_dark use_dark()

Prep data and model¶

In [2]:

FITS = './data/example.fits'
wavecal_csv = './data/133mVis_wavcal.csv'
s26 = Vis133M(FITS,wavecal_csv)
s26.set_scale(0.0001837)
s26.set_dark(frame=49, channel=43)
s26.set_baseline_zero(True)
s26.set_time_linspace(0, 10)

model = BHModel()

FITS = './data/example.fits' wavecal_csv = './data/133mVis_wavcal.csv' s26 = Vis133M(FITS,wavecal_csv) s26.set_scale(0.0001837) s26.set_dark(frame=49, channel=43) s26.set_baseline_zero(True) s26.set_time_linspace(0, 10) model = BHModel()

Single fit¶

In [3]:

fitr = BHFitter(s26, model, warm_start=True, weight="none")
res = fitr.fit(frame=38, channel=31)
display(res["summary"])
_ = fitr.plot_single(res, title="Example")

fitr = BHFitter(s26, model, warm_start=True, weight="none") res = fitr.fit(frame=38, channel=31) display(res["summary"]) _ = fitr.plot_single(res, title="Example")

	Parameter	Value
0	C	0.262 ± 0.002
1	T_rot	3600 ± 100 K
2	dx	(101 ± 7)e-5 nm
3	w_inst	0.0245 ± 0.0002 nm
4	base	(43 ± 2)e-4
5	I_R7	(25 ± 3)e-5
6	I_R8	(30 ± 3)e-5

Batch fit¶

In [4]:

resb,curves = fitr.batch((25,26,31), range(28,31),return_curves=True)

resb,curves = fitr.batch((25,26,31), range(28,31),return_curves=True)

In [5]:

_ = fitr.plot_grid(curves)

_ = fitr.plot_grid(curves)

Plot Overlay¶

In [6]:

fitr.plot_overlay(curves,frame=31,channels=range(28,31),line_width=1.5)

fitr.plot_overlay(curves,frame=31,channels=range(28,31),line_width=1.5)

Out[6]:

(<Figure size 700x400 with 1 Axes>,
 <Axes: xlabel='wavelength [nm]', ylabel='intensity [arb]'>)