• View
• Edit
• History
• Print

Mupy

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