blase
blase copied to clipboard
An attribute (not nn.Param) for Equivalent Width and FWHM
What we want
An attribute for Equivalent Width. It will have the same shape as lam_centers
, amp
, sigma
, and gamma
, but it will not be an nn.Parameter
. Just a derived parameter.
How to compute it
Two options:
-
Analytically This turns out to be impossible for Voigt profiles because the Lorenzian long tails integrate to infinity. However, we could simply provide a number of FWHMs out to which we integrate. This choice has the benefit of operating on a table of derived properties $\sigma$, $A$, $\gamma$.
-
Numerically Simply
torch.sum()
the n_lines axis offlux2D
, including the size of the pixels in angstroms. Probably np.trapz or equivalent would work better (or scipy quad).