Author: Stephane Willocq

Sherpa 2.1.1 is installed on titan. For information about Sherpa, see http://www.sherpa-mc.de/.

Run instructions below are for the worker nodes (with SLC6 OS).

Define environment variables to point to the Sherpa and Root installation directories for convenience

• export SHERPAROOT=/home/willocq/bin/sherpa/SHERPA-MC-2.1.1
• export ROOTSYS=/home/willocq/root/root_v5.34.21/root

Add the location of the Sherpa and Root binaries to your PATH environment variable

• export PATH=$SHERPAROOT/bin:$ROOTSYS/bin:$PATH

Make sure that the Sherpa, Root and hepmc libraries can be found by adding them to the LD_LIBRARY_PATH environment variable

• export LD_LIBRARY_PATH=$SHERPAROOT/lib/SHERPA-MC:/home/willocq/root/root_v5.34.21/root/lib:/home/willocq/bin/hepmc/x86_64-slc6-gcc45-opt/lib:$LD_LIBRARY_PATH

For a quicker setup of all the required environment variables do as follows:

• source /home/willocq/bin/sherpa/setup_sherpa.sh

Copy a set of run cards to your own directory. For example use the dijet process by copying the following cards: $SHERPAROOT/Examples/Jets_at_HadronColliders/LHC_Jets_MEPS/Run.dat. This particular process takes several days to compute matrix elements though! So it may be better to try a different example like the one in $SHERPAROOT/Examples/Models/SM_ZPrime/, (the latter requires executing Sherpa then ./makelibs before running again with Sherpa)

To run, type

• Sherpa

T0 specify the type of output, see https://sherpa.hepforge.org/doc/SHERPA-MC-2.1.1.html#Event-output-formats. Unfortunately, the software to produce Delphes output from Sherpa is incompatible with the latest Delphes releases (3.X.Y).

xAOD Analysis

• The most up-to-date and detailed tutorial for xAOD analysis is available at https://twiki.cern.ch/twiki/bin/viewauth/AtlasComputing/SoftwareTutorialxAODAnalysisInROOT.
• First time setup
  • setupATLAS
  • cd to your working release directory
  • rcSetup Base,2.3.16
  • rc find_packages
  • rc compile
• Setup for every session
  • setupATLAS
  • cd to your working release directory
  • rcSetup
• To add a package
  • rc checkout_pkg PackageName
• To start clean and recompile a whole working area
  • cd to your working release directory
  • rc clean
  • rc find_packages
  • rc compile
• To run
  • cd Run/     (from your working release directory)
  • root -l -b -q $ROOTCOREDIR/scripts/load_packages.C ‘ATestRun.cxx(“outputDir”)’

Truth Smearing Functions

Smearing functions for physics upgrade studies are available in package TruthToRecoFunctions, see https://svnweb.cern.ch/trac/atlasgroups/browser/PhysicsAnalysis/EuropeanStrategy/TruthToRecoFunctions/trunk.

 

To get started with research on ATLAS, here are some pointers to learn some of the more technical aspects of the data analysis.

1. Computer account

Our compute cluster is part of a larger system named NET3. Instructions on how to request and set up an account are available at https://websites.umass.edu/willocq/2018/09/05/request-account-on-net3/.

General information about the cluster is available at https://twiki.cern.ch/twiki/bin/view/Main/UMassNET3. In particular, you will see instructions to set up your account for ATLAS analysis.

2. UNIX

The titan computer cluster runs with the UNIX operating system. If you are not familiar with it, check http://www.ee.surrey.ac.uk/Teaching/Unix/.

3. C++

To help firm up your knowledge about programming languages, see the book in room LGRT 1036: “C++ How To Program”. Chapters 1 through 11 are most relevant. It’s a good idea to do some of the exercises (unless you already are very familiar with the material).

Otherwise, a good alternate option is the following online tutorial: http://www.cplusplus.com/doc/tutorial/.

This is a fair amount of material so you may want to do this in a number of steps in between ROOT tutorial sessions.

4. ROOT tutorial

The main data analysis tool is ROOT. For general information, see http://root.cern.ch/. You can download the user guide and binaries from that web site. It may be good for you to try running on your laptop if you wish.

A good primer is provided on the ROOT web site, see https://root.cern.ch/root/htmldoc/guides/primer/ROOTPrimer.html.

To run root on titan you will need to type the following every time you login:

lsetup root

I suggest you create root and tutorial subdirectories and work inside of the tutorial directory as you follow each of the tutorials above in turn.

cd
mkdir -p root/tutorial
cd root/tutorial

To run just enter:

root

5. Other tools

If you plan to use your laptop you will need to install:

• text editor (or your choice, I like emacs)
• remote connection program like PuTTy (if you run Windows) to run ssh, otherwise macs already have ssh
• X11 graphics program to view graphics output from titan, e.g. Cygwin for Windows: http://x.cygwin.com/
• root (from root.cern.ch)

The above will give you plenty to do for a couple of weeks but will be valuable to establish a foundation for your research.

Stephane W.

From titan, the following steps need to be done for every new terminal session:

• cd root/2013    (0r choose your main directory; all directories mentioned below are relative to this directory)
• source script.sh

Files needed for analysis:

•  Directory ContactInteractionAnalysis contains the main analysis code class implemented in ContactInteractionAnalysis/Analyzer.cxx and ContactInteractionAnalysis/ContactInteractionAnalysis/Analysis.h.
• input.txt lists the input data sample for analysis.
• Input ntuple files can be listed at /titan_rw/atlas/common/D3PD/mc12_8TeV/slim_p1328 for Monte Carlo simulation and at /titan_rw/atlas/common/D3PD/data12_8TeV/slim_p1345 for real data.

To compile the code:

• cd ContactInteractionAnalysis/cmt
• make

To run the analysis:

• cd ContactInteractionAnalysis
• ./bin/runClassTest -f input.txt -o outputFile.root –isMC      (for MC, note the double dash in front of isMC)
• ./bin/runClassTest -f input.txt -o outputFile.root      (for data)

Relevant variables in ntuple files:

• mu_muid_pt          muon transverse momentum
• mu_muid_eta        muon pseudorapidity
• mu_muid_charge  muon electric charge
• mu_muid_id_qoverp   associated ID track q/p
• mu_muid_id_theta      associated ID track polar angle
• mu_muid_id_phi         associated ID tracks
• mu_muid_me_qoverp   associated MS track q/p
• mu_muid_me_theta      associated MS track polar angle
• mu_muid_me_phi         associated MS tracks

Histograms to produce:

• mu_muid_pt
  • for all muons in the event and those with opposite q/p for ID and MS tracks
  • for all muons passing the selection criteria (goodSmearedMuons for MC) and those with opposite q/p for ID and MS tracks
  • the goal is to plot the rate of charge flips between the ID and MS tracks as a function of muon pt for both data and MC simulation
• mu_muid_eta
  • same categories as for pt listed above

• Code for the Z_KK to ZH decay analysis is on titan in /home/willocq/Snowmass2013/root/. Create a work directory named Z_KK and copy the contents of the original directory into it:
  • cd
  • mkdir -p root/Z_KK
  • cd ~/root/Z_KK
  • cp /home/willocq/Snowmass2013/root/*.C .
• Create a link to an include file required for compilation
  • mkdir classes
  • cd ~/root/Z_KK/classes
  • ln -s /scratch/willocq/Snowmass2013/Delphes-3.0.9/classes/SortableObject.h
• Run the analysis code
  • cd ~/root/Z_KK
  • root -b -q runMyExample1.C
• A series of histograms is produced and stored in an output file with a name specified in runMyExample1.C, it’ll be of the form output_filename.root, where filename is the name of a given input file. The input files are stored in /scratch/willocq/Snowmass2013/Delphes-3.0.9/.
• The histograms can be plotted with the help of the myPlots.C macro; for example start up root and give it the histogram file name as input, as follows
  • root output_filename.root
• At the root prompt you can then execute the plotting macro, which actually runs a number of functions defined in myPlots.C
  • .x runMyPlots.C

The content of ntuples produced by the Delphes fast simulation program is described in https://cp3.irmp.ucl.ac.be/projects/delphes/wiki/WorkBook/RootTreeDescription.

Sample analysis programs are available on titan in directory /home/willocq/Snowmass2013/root, see the files myExample1.C and myExample2.C. The first is a rather simple example of analysis code to produce some simple histograms. The second example uses a more complex program structure as well as Delphes analysis classes like a histogram manager.

Root files for analysis are located in /scratch/willocq/Snowmass2013/Delphes-3.0.9.

• On titan, copy the following file ~willocq/root/example/simple_plot_truth.C to your own subdirectory. After setting up root with localSetupROOT, the simple plotting macro can be run with

root simple_plot_truth.C

or you can just start root and at the root prompt type:

.x simple_plot_truth.C

• Once this is working, extend the macro to create additional canvases like c1 to plot the following variables stored in the truth ntuple file that is opened in the macro:
  • jetAntiKt4TruthJets_pt          transverse momentum of truth jets
  • jetAntiKt4TruthJets_eta        pseudorapidity  of truth jets
  • el_pt                                           transverse momentum of truth electrons
  • el_eta                                         pseudorapidity of truth electrons
  • mu_muid_pt                             transverse momentum of truth muons
  • mu_muid_eta                           pseudorapidity of truth muons
• Do the above plots with a logarithmic vertical scale
• Produce a second set of those above plots with a minimum pt cut of 25000 MeV for  jets and 5000 MeV for both electrons and muons, i.e. in the example there is a cut on the barcode, simply replace that with a cut on the variable mu_muid_pt>5000 for muons, and similarly for electrons and jets
• Produce a plot of the pseudorapidity of electrons and muons (separately) showing an overlay of negatively and postively charged particles (histograms with different colors).