Using Rivet in ALICE
Rivet has been integrated into the ALICE build system and is maintained there. Thus, it is also deployed to cvmfs, from where it can be used by Grid workernodes, interactive lxplus sessions, or other machines using software from cvmfs (e.g. CERN Linux desktops). It can also be installed locally using aliBuild. However, if you just get started and want to gain experience with Rivet, we recommend to use the cvmfs installation on lxplus, as this is maintained centrally and known to work (if it is not, please report!).
Using Rivet on lxplus
In order to use Rivet on lxplus, you just have to logon and use alienv to enter a shell with the proper environment loaded (same procedure as for AliPhysics):
alienv enter Rivet
You can then check which packages and versions are loaded using:
Now you have all the rivet tools (rivet, rivet-build, ...) available. To list the available analyses (those distributed with the used version of Rivet), you can use:
rivet --list-analyses ALICE_
where the second command only lists analyses for ALICE publications. You can display the details of one of these analyses (here ALICE_2010_S8625980) using:
rivet --show-analysis ALICE_2010_S8625980
Assuming you have some output from a Monte Carlo generator in a HepMC file (e.g. /eos/project/a/alipwgmm/rivet/hepmc/pythia6_pp7000.hepmc or see here for instructions on running MC generators), you can run an analysis (e.g. ALICE_2010_S8625980) on this input using:
rivet -a ALICE_2010_S8625980 /eos/project/a/alipwgmm/rivet/hepmc/pythia6_pp7000.hepmc
where -a is used to specify the analysis you want to run (you can use any analysis that was distributed with Rivet). You can also list more than one analysis.
Beware the absence (!) of extension of your analysis name specified at the rivet command, i.e. no "rivet -a ALICE_2010_S8625980.cc ...". This would make your analysis not recognisable by Rivet, and it will complain so.
The command line above runs the analysis on the given input (~1000 generated events, here) and stores the produced output, typically histograms, in the file Rivet.yoda. Based on this yoda file, you can produce plots by running:
This command puts the plots in the subdirectory rivet-plots/. The directory also contains an index.html, which you can open in a browser to look and navigate through the plots. If you have run Rivet on an lxplus node you probably must copy this directory to your local PC and then point your browser to the directory, i.e.
scp -r lxplus.cern.ch:<myRivetDir>/rivet-plots /tmp/myRivetPlots
open /tmp/myRivetPlots/index.html (on Mac OS X, otherwise point browser to file:///tmp/myRivetPlots/index.html)
If on lxplus you use a directory under /tmp you have to copy from the exact host you are working on instead of lxplus.cern.ch, e.g. lxplus012.cern.ch. You can get actual name of the hosting machine with
We recommend not to use direclty the directory /tmp to avoid clashes with other user. Please use /tmp/
Running Rivet together with a generator
You can use Rivet to run on the output of an arbitrary generator for which an interface to HepMC exists - e.g. those which have been integrated into the ALICE build system and are, thus, available on cvmfs or for local installation. There are two common options to run the generator and Rivet:
- intermediate file:
- run generator writing to hepmc file
- run Rivet reading from hepmc file (see example above)
- FIFO (first-in-first-out pipe, on Linux)
- create FIFO (NB: not possible on afs, need to use /tmp on lxplus)
- run generator writing to FIFO (process won't finish)
- run Rivet reading from FIFO
The latter option avoids writing a file to disk and is beneficial if the generator is fast. You can run a generator directly or use
aligenmc after you have set it up with:
alienv enter AliGenerators
Then, you would have to do something like:
aligenmc -g pythia8 -E 7000 -N 10000 -o /tmp/$USER/fifo_$UID.hepmc &> /tmp/$USER/aligenmc_$UID.log &
rivet -a ALICE_2010_S8625980 /tmp/$USER/fifo_$UID.hepmc
aligenmc is the script developed within ALICE that does the running of a chosen generator, with a dedicated configuration file, and writing the output to the specified file/FIFO. For further instructions on how to run different generators, see:
- MC generator documentation
aligenmc is not yet updated to generate heavy-ion collisions with the Pythia8/Angantyr model. In case you want to generate Pb-Pb collisions with Pythia8/Angantyr model you can
run-pythia -n 100 -e 2760 -o /tmp/$USER/fifo_$UID.hepmc -c "Beams:idA 1000822080" -c "Beams:idB 1000822080" -c "HeavyIon:SigFitNGen = 0" -c "HeavyIon:SigFitDefPar = 12.32,1.63,0.19,0.0,0.0,0.0,0.0,0.0"
Beware that this line is valid only for Pb-Pb 2.76 TeV collisions. Should you need to change the collision energy (the the beams), make sure you run only the initial part
run-pythia -n 100 -e 2760 -o /tmp/$USER/fifo_$UID.hepmc -c "Beams:idA 1000822080" -c "Beams:idB 1000822080"