saarioka / tinyvdm

README

Lightweight vdM analysis meant to illustrate the process of the vdM framework.

Setup

Will run with a general installation of Python3 found in LXPLUS etc, but needs some additional (common) packages. Depending on your environment, the correct executable might be called either python or python3 (see version with python --version).

Install the needed packages with

python3 -m pip install -U -r requirements.txt

Contents

• minimal_example.ipynb Jupyter Notebook with detailed description about the process. This is a good starting point to run your own test analysis

• minimal_example.html Pre-run version of the previous notebook, featuring printouts. Handy study material if you do not want to set up your environment.

• main.py Script to be called form the command line, has more features than minimal_example

Running

Notebook

After installing the dependencies, start the notebook with jupyter notebook

For SWAN users: In the beginning, start the command line and install the dependencies with pip install --user iminuit mplhep --upgrade. After doing this, run

export PYTHONPATH=$CERNBOX_HOME/.local/lib/python3.9/site-packages:$PYTHONPATH

(substitute Python version number with correct one depending on kernel, here it is the curent default 3.9). Then start the notebook.

CMD

After installing dependencies, run analysis for the preferred example file using main.py, for example

python3 main.py -l pltlumizero /eos/cms/store/group/dpg_bril/comm_bril/vdmdata/2021/original/7525/7525_2110302352_2110310014.hd5

I suggest looking at the code while comparing the printouts to the code statements to get an idea what is being calculated.

If you used the exact command above, look into folder output and compare it to the contents of framework_results_without_beam_current_calibration.

Also try flags -pdf, -cbc, -fit and other values for -l. See help with python3 main.py --help

List of features (CMD)

• Choice of luminometer
• Choice of fit
• Beam current calibration
• Fit pdfs
• Multiple files as input. Give a list of files as unnamed parameter or pipe the output of ls with ls *.hd5 | xargs python3 main.py -l <lumi>