Introduction
This document contains instructions to run a threshold scan of benchmark point two for the CLIC CDR. It is assumed that the user has a valid lxplus account.
Setup
Please prepare your environment, so that it finds all relevant executables and libraries.
export PATH=/afs/cern.ch/sw/lcg/contrib/gcc/4.5.2/x86_64-slc5/bin:${PATH} export LD_LIBRARY_PATH=/afs/cern.ch/sw/lcg/contrib/gcc/4.5.2/x86_64-slc5/lib64:${LD_LIBRARY_PATH} export PATH=/afs/cern.ch/user/j/jfstrube/public/PYTHON/bin:${PATH} export PATH=/afs/cern.ch/user/j/jfstrube/public/Ocaml-3.12.1/bin:${PATH} export LD_LIBRARY_PATH=/afs/cern.ch/user/j/jfstrube/public/Ocaml-3.12.1/lib:${LD_LIBRARY_PATH}
Running the code
The input data
Please copy the directory with the input files from afs
cp -a /afs/cern.ch/user/j/jfstrube/CLIC_Amplitudes/SUSY/ .
The code
The script that does all the work is called compile_thresholdscan.py
. It carries out the following steps
- obtain a list of the subdirectories
- in each of the subdirectories, look for a
.txt
file with the same name as the directory - for each of the given energies, create the Whizard {
.sin
file from the.txt
file - call Whizard on this file
- parse the resulting log file and remember the result
- rename the logfile, so that it ends with the energy at which the cross-section was computed
- move to the next directory
- at the end print the output on the screen and write to file
Call the script with /afs/cern.ch/user/j/jfstrube/public/PYTHON/bin/python compile_thresholdscan.py
Analyzing the output
Conversion to comma-separated values
The first step is to turn the output to a csv file. This is not strictly necessary, but it provides a threshold to involuntarily overwriting the previous output.
/afs/cern.ch/user/j/jfstrube/public/PYTHON/bin/python convert2csv.py
looks for the file results.txt
and creates the file results.csv
Plotting
The csv file can be turned into a plot using /afs/cern.ch/user/j/jfstrube/public/PYTHON/bin/python plot_results.py
. This simple script uses the matplotlib library to create pdf files from results.csv