Posted here for reference so that anyone can try it out and poke at the code. These plots, shown in January, depict fits to 54 million positrons above 1.8 GeV. The data is from runs 15923 to 15928 (roughly 5% of the 60 hour dataset).

To reproduce the exact plots I have been showing (after your repository is checked out and set up) do the following:

cd gm2ilratio/fitting/pyfitter
curl -L -O
tar -xvjf clusters_15922-15928.tar.bz2
python -i ./ clusters_*.hist

and all of the plots will pop up.


  • Some results are written out to ./save_fits/, but ONLY if it does not already exist. That means that you must explicitly MOVE, RENAME, OR REMOVE that directory before executing the script. * The plots will not pop up in their own windows if 1) you set matplotlib to use a non-UI backend with something like matplotlib.use('svg'), or 2) you don't have X11 connection for some reason (check your ~/.ssh/config). * If you don't use the -i option then the plots will disappear immediately as the script exits.

Below are some of the results.

N.B. Remember that there are some steps to setting up the environment to run the code. The full list is

  • (if not already done) clone the repository and do git checkout develop (instructions) * (if not already done) build the Blinders library on your system with gm2ilratio/pymodules/Blinders/ (instructions) * set your blinding hash-phrase with export OMEGA_BLINDERS_PHRASE="some text you might remember later" (instructions) * activate the Python environment with source gm2ilratio/py.env (instructions)