Index
mujpy is here. Documentation on mujpy at ReadTheDocs
Old page: What would be the case for a MuPy (umpy) environment?
- Start from a separate wiki
- Start without gui, edit parameter list, gui later
- The interface should have the following blocks (musrfit-style)
- Data input
- Function assembly
- Parameters
- Plot
- Statistics
- Fourier
- Fit: start from the parameter list of a global fit in a given list of runs
- Fit with different modules: minuit, minuit2, scipy.leastsq, ...
- Standard asymmetry fit
- Tentative Rotating Frame fit (has artifacts)
- Fit globally, across detectors and runs (physical parameters)
- Single detector Minuit fits, along directions {$\alpha=x,y,z$}, according to
{$N_i(t)=N_{0i}e^{-t/\tau_\mu}(1\pm A_iP_\alpha(t))+B_i$}
- Polarization projection functions may be reconstructed in two determinations as
{$ P_{\alpha,\pm}=\pm\frac1{A_i}\left(\frac{N_i-B_i}{N_{0i}}e^{-t/\tau_\mu}-1\right) $}
- Error propagation, with {$s_N=\sqrt(N)$} and {$s_{N0},s_A,s_B$} from the Minuit fit yield the variance on each determination
{$ s_P^2=\left(\frac{e^{-t/\tau_\mu}}{N_0A}\right)^2 \left[N(t)+s_B^2+\frac{N(t)-B}{N_0}s_{N0}^2 \right] + \left(\frac{N-B}{N_0}e^{-t/\tau_\mu}-1\right)^2 s_A^2$}
- Rename the two polarization determinations {$P_\pm$} and their errors {$s_\pm$}; plot polarizations along directions {$\alpha=x,y,z$}, according to
{$P(t)= \frac{s_-^2P_+(t)+s_+^2P_-(t)}{s_+^2+s_-^2}$}
- with error
{$s_P = sqrt2 \frac{s_+s_-}{\sqrt{s_+^2+s_-^2}}$}
- The global fit strategy implies that parameters will belong to one of these categories:
- {$N_{0i},B_i,A_i$}, to be fitted globally across physical parameters (i.e. runs)
- optional polar angles {$\theta,\phi,\cdots$} to relate different polarization projections, according to a specific model (e.g. internal field in a single crystal, muonium etc.), to be also fitted globally across runs
- Run dependent parameters that remain global across detectors
- Global polarization parameters (e.g. temperature independent muon fractions), across both detectors and runs
Index
