Collection of scripts for performing Muon Collider studies. These scripts assume you're inside the singularity image used in the Fermilab Muon Collider tutorial which does some path mapping for you on the Snowmass cluster. In particular, you should run something along the lines of:
apptainer build k4toroid.sif docker://madbaron/k4test-ubuntu:latest
apptainer run --no-home -B /ospool/uc-shared/project/futurecolliders/data:/data -B /scratch/$USER -B /home/$USER /scratch/trholmes/mucol/k4toroid.sif
You can update this command to use your own location instead of mine (although I believe it will work if you just run my .sif file and skip the first line).
and then from inside the image, run:
source /opt/spack/share/spack/setup-env.sh
spack env activate -V k4prod
source /opt/spack/opt/spack/linux-ubuntu22.04-x86_64/gcc-11.3.0/mucoll-stack-2023-07-30-ysejlaccel4azxh3bxzsdb7asuzxbfof/setup.sh
export MARLIN_DLL=$MARLIN_DLL:/opt/MyBIBUtils/lib/libMyBIBUtils.so
(FYI, putting these two in one shell script and sourcing it will not work; it won't execute the second command until you exit the singularity. You can also update to a newer release, which you can find here, but I ran into trouble running with 1.7 with my current setup.)
For a very simple version of reading an slcio file, see makeMuonPlots_simple.py.
For a slightly more advanced version that handles making multiple histograms more elegantly, see makeMuonPlots.py. It may be helpful to look at these side-by-side if you're trying to learn what's going on here.
Both of these scripts will put plots in a plots/ directory - make one if you don't have it already or it will crash.
To work with key4hep files, look at makeMuonPlots_k4hep.py. This is just a skeleton file; it doesn't show you as much as even the makeMuonPlots_simple.py file when it comes to e.g. plotting, but it does give a simple recipe for how to open and loop over these files with the key4hep syntax. Be sure to read the comments at the top of this file for specific instructions on how to run the docker image for this.
For the SLCIO files, to get a list of the colletions you can access, run anajob <filename> on one of your files.
Use the COLLECTION NAME to access them.
To understand what kinds of functions you can use on the particles in these collections, look at the COLLECTION TYPE and look up its functions here.
To understand in more depth how reconstruction works, you can look at the code here. To understand what options were passed to this code, you'll need to look at the steering files that were used to run it. Fede's are here, but you'd have to know which one he used. For more generic ones, have a look at the ones used in the tutorial.