Inside the PU14
directory you can find code we borrowed from the Pileup2014 workshop that takes care of reading in the events and mixing background with signal events.
Inside the include
directory you can find a couple of classes performing background subtraction, grooming jets, and a jet-to-jet matching algorithm.
Follow the installation instructions below to run an example program.
If you are using mac or linux, the steps are relatively straightforward. For windows machines I'm not sure what to do. These are the things you need to install:
- C++ compiler: on mac you could install xcode (found on App Store) to get the g++ compilers
- Latest version of ROOT: follow instructions here https://root.cern.ch/downloading-root
- Fastjet and contrib package: follow the steps in the above about the fastjet installation
- Pythia 8: it can be found here - http://home.thep.lu.se/~torbjorn/Pythia.html
Make sure that the root-config, pythia8-config and fastjet-config executables can be found in the $PATH environment variable. Once the above is done, we can proceed with the compilation of the JetToyHI code:
git clone https://github.com/JetQuenchingTools/JetToyHI.git
cd JetToyHI
echo (FASTJET LOCATION) > .fastjet
echo (PYTHIA8 LOCATION) > .pythia8
cd PU14
echo (FASTJET LOCATION) > .fastjet
echo (PYTHIA8 LOCATION) > .pythia8
./mkmk
make
cd ..
scripts/mkcxx.pl -f -s -1 -r -8 '-IPU14' -l '-LPU14 -lPU14 -lz'
make
./runFromFile -hard samples/PythiaEventsTune14PtHat120.pu14 -pileup samples/ThermalEventsMult12000PtAv0.70.pu14 -nev 10
If you don't have an lxplus account, you can request one here: https://account.cern.ch/account/Externals/
First install fastjet and the contrib package (you only have to do this once)
ssh -Y <username>@lxplus.cern.ch
cd <dirInWhichYouWantToInstall>
source /cvmfs/sft.cern.ch/lcg/external/gcc/4.8.1/x86_64-slc6-gcc48-opt/setup.sh
source /cvmfs/sft.cern.ch/lcg/app/releases/ROOT/6.06.08/x86_64-slc6-gcc48-opt/root/bin/thisroot.sh
#Install fastjet (you just have to do this once)
curl -O http://fastjet.fr/repo/fastjet-3.3.0.tar.gz
tar zxvf fastjet-3.3.0.tar.gz
cd fastjet-3.3.0/
./configure --prefix=$PWD/../fastjet-install
make
make check
make install
FASTJET=$PWD/../fastjet-install
cd ..
export FJ_CONTRIB_VER=1.026
curl -Lo source.tar.gz http://fastjet.hepforge.org/contrib/downloads/fjcontrib-"$FJ_CONTRIB_VER".tar.gz
tar xzf source.tar.gz
cd fjcontrib-"$FJ_CONTRIB_VER"
./configure --fastjet-config=$FASTJET/bin/fastjet-config --prefix=`$FASTJET/bin/fastjet-config --prefix`
make
make install
make fragile-shared #make shared library
make fragile-shared-install
cd ..
Next steps
git clone https://github.com/JetQuenchingTools/JetToyHI.git
cd JetToyHI
. setup.sh #this step you will need to repeat next time you login
echo `$FASTJET/bin/fastjet-config --prefix` > .fastjet
cd PU14
echo `$FASTJET/bin/fastjet-config --prefix` > .fastjet
./mkmk
make
cd ..
scripts/mkcxx.pl -f -s -1 -r -8 '-IPU14' -l '-LPU14 -lPU14 -lz'
make
./runFromFile -hard samples/PythiaEventsTune14PtHat120.pu14 -pileup samples/ThermalEventsMult12000PtAv0.70.pu14 -nev 10
You will have produced a root file with a tree. In this tree properties of jets are stored in std::vector format and are all alligned to the signal jets. So you can for example plot sigJetPt
vs csJetPt
to get the response matrix for constituent-subtracted jets. An example ROOT plotting macro which will draw the jet energy and mass scale can be found here: plot/plotJetEnergyScale.C
.
- If you want to contribute to this code you need to have a github account. Go here to do so: https://github.com/join.
- Fork the original repository. Go to: https://github.com/JetQuenchingTools/JetToyHI and click 'Fork' in the upper right corner.
- Instead of cloning the original repository as shown above, clone your own.
- After committing your changes to your own branch, push them to your own fork. Don't know how to do this, ask your colleages or use google which might bring you here https://services.github.com/on-demand/downloads/github-git-cheat-sheet/
- Do a pull request once you have finished your developements.
Event samples can be found in the jet quenching CERNBOX:
- From lxplus: /eos/project/j/jetquenching/JetWorkshop2017/samples
- Webbrowser: https://cernbox.cern.ch/index.php/s/rNvuBSMHX3hsRhc
- Mount eos on a laptop or local desktop: https://cern.service-now.com/service-portal/article.do?n=KB0003493
You will find samples from various event generators. For underlying event we have: 'thermal' which is independent particle production using a Boltzmann distribution with a fixed multiplicity and mean pT (indicated in the file names). For the hard signal we have PYTHIA8 and JEWEL events with various pT,hat settings.
More details about the available samples can be found here: https://twiki.cern.ch/twiki/bin/view/JetQuenchingTools/PU14Samples