During analysis meeting, few folks were asking about optics codes. I attach a tarball that contains a few scripts to drive the prototype MAUS optics application. In this case I make a plot of matched beta vs momentum.
To reproduce the beta vs p plot, download the attached tarball. Install MAUS version 0.1.1. Then do
source <path to maus version 0.1.1>/env.sh tar -xzf fs2a.tgz cd fs2a python run_optics.py python ana_optics.py
It works in MAUS v0.1.1, probably not in some previous versions (for bureaucratic reasons). The plot looks like:
Just timing it on my quite good desktop PC, it takes 1 min 50 seconds to generate that plot. Probably we can optimise it.
Folks knowledgeable in ROOT can also extract e.g. beta vs z, angular momentum vs z fairly trivially by using the ROOT output.
Beta vs z at 200 MeV/c:
Angular momentum vs z:
There are some further examples at
Really, right now I am treating this like a personal code until we can do the tidy up (i.e. it isn't routinely tested and sometimes it doesn't work just because I forgot to fix this thing or that thing). Peter Lane is kindly applying some effort to tidy it up...