A Partial-Wave/Amplitude Analysis Software Framework

 

The PyPWA Project
Thomas Jefferson National Accelerator Facility
Newport News, VA

MCMC (Markov Chain Monte Carlo) Tutorial (code by Peter Pauli, Glasgow, UK)

You need to run first a normal Minuit fitting to get the seeds (initial values) for the MCMC.

Perhaps you will need to include the library "corner" in your enviroment (conda install -c conda-forge corner)

The two notebook example files are:

MCMCexample_minuit.ipynt

MCMCexample.ipynt

and two associated amplitudes (a la JPAC)

AmplitudeJPACfit.py

AmplitudeJPACfitAngles.py

The difference between those amplitudes is that the first treat amplitudes as two real numbers (components)

and the second as magnitude and angle.