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 https://cdcvs.fnal.gov/redmine/attachments/download/52138/clusters_15922-15928.tar.bz2 tar -xvjf clusters_15922-15928.tar.bz2 python -i ./data_fit_test.py 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
-ioption 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/build.sh(instructions) * set your blinding hash-phrase with
export OMEGA_BLINDERS_PHRASE="some text you might remember later"(instructions) * activate the Python environment with