For Robin's work on NMSSM H1 -> 2a1 -> 4tau

Monte Carlo stuff is in pythia8 and Delphes3 folders. You'll need Boost both Pythia and Delphes.

Sorry in advance if this is a bit rubbish - it's constantly in a state of flux, and no-one but me uses it. (Standard weak HEP excuse)

Also see runningOnBlueCrystal.md

• BOOST (for Pythia event generation and Delphes)
• ROOT (for Delphes)
• Pythia8 (for event generation)
• Delphes 3.1.X, X>2 (for converting HepMC to ROOT)

• See below

• Download tarball, extract
• Follow README.HepMC

• mkdir root
• cd root
• Get latest 5-34-XX version from git, see website (e.g. git clone http://root.cern.ch/git/root.git root-v5-34, git checkout -b v5-34-19 v5-34-19)
• mkdir install
• mkdir config
• cd root-v5-34
• ./configure --prefix=$HOME/root/install --etcdir=/gpfs/cluster/phys/ra12451/root/config/
• make -j2
• make install
• add source $HOME/root/root-v5-34/bin/thisroot.sh to ~/.bashrc or ~/.bash_profile

** Requires version 3.1.2 or greater** (basically need track class to have impact parameter variables Dxy,Xd, Yd, etc - check in DelphesClasses.h)

• Download tarball, extract
• cd into folder, run make
• Edit and run setupDelphesLinks.sh to create symlinks to files in your Delphes folder

Signal:

• Make signal MC using CalcHEP. Outputs as LHE file.
• Hadronise it using Pythia:

make mainLHEhadronise
./mainLHEhadronise.exe # to show usage
./mainLHEhadronise.exe <LHEfilename & path> <output HEPMC filename>

This outputs a HEPMC file.

QCD[b/c/q-g scatter]:

• Make it all in Pythia using one executable

make mainQCD
./mainQCD.exe # to show usage

Once you have a HEPMC file, you now run it through Delphes to produce a ROOT file:

• Run ./DelphesHepMC for relevant options.
• Generally, you'll want to use DelphesCards\delphes_card_cms_bare.tcl. If you're using a version of Delphes older than 3.1.2, use delphes_card_CMS_bare_OLD.tcl.
• Use Scripts/submitDelphesLots.sh or Scripts/submitDelphesSingle.sh for use on PBS batch system, or Scripts/runDelphesLots.sh for lots of local jobs.

• basicScript: template script for making a new program (see below)
• mainAnalysis: does lots of plots, like pT, track eta Vs phi, soft track distributions. Used for QCDb rejection studies
• massPlots: makes plots of invariant masses of mu+tk, and calculates & plots correlation coefficients for backgrounds
• cutFlow: prints out cutflow of various cuts stages in signal selection
• IP: plots impact parameters of things. Not really kept up to date. OUT OF DATE

• Look at Scripts/createNewScript.sh - you'll need to edit the path where Delphes is installed.
• Run createNewScript.sh <myscriptName>. Make sure it puts the new scripts in your Delphes/examples folder
• Go to your Delphes directory, run ./configure to pick up new script
• Look at Scripts/createMakefile.sh, modify paths as necessary
• Run Scripts/createMakefile.sh to modify the Delphes makefile to add support for C++11 and Boost.
• Now run make
• To run your program, do ./runMyCoolNewProgram in the main Delphes installation folder.

You can now edit your new analysis code. To re-make it, just run makeScripts. This will re-make ALL the analysis programs. See makeMain and makeMass for examples on how to write a script to re-make just one program.

• There are old Pythia programs, mainQCDb, mainQCDScatter, etc. These are are all old, and buggy. Do not use them. Only use mainQCD (plus, it does everything with command-line options, so win).

• cuts.h is loats of central functions that test muons, tracks, etc. Not the best, but is centralised.

• Delphes distance units for Xd, Yd, Dxy, etc are mm. You can see this in ParticlePropagator.cc

• By default Dxy is NOT stored properly. Add the line:

candidate->Dxy = dxy*1.0E3;

after

candidate->Zd = zd*1.0E3;

http://www.boost.org/doc/libs/1_55_0/more/getting_started/unix-variants.html

1. Download boost_1_55_0.tar.bz2.

2. Extract tar --bzip2 -xf /path/to/boost_1_55_0.tar.bz2 (say to ~/boost_1_55_0).

3. Now need to build the separate Boost libraries

   • The only Boost libraries that must be built separately are:

   • Boost.Chrono

   • Boost.Context

   • Boost.Filesystem

   • Boost.GraphParallel

   • Boost.IOStreams

   • Boost.Locale

   • Boost.MPI

   • Boost.ProgramOptions

   • Boost.Python (see the Boost.Python build documentation before building and installing it)

   • Boost.Regex

   • Boost.Serialization

   • Boost.Signals

   • Boost.System

   • Boost.Thread

   • Boost.Timer

   • Boost.Wave

4. Create installation folder for those libraries: mkdir ~/boost_1_55_0_install

5. In the original boost folder, do ./bootsrap.sh --prefix=${HOME}/boost_1_55_0_install

6. Then run ./b2 install. This should install libraries to ~/boost_1_55_0_install/lib

7. In the Delphes Makefile, change the following:

   • Add -I $(HOME)/boost_1_55_0 -I $(HOME)/boost_1_55_0_install/include to CCXFLAGS (line 18)
   • Add -L/panfs/panasas01/phys/ra12451/boost_1_55_0_install/lib -lboost_program_options to DELPHES_LIBS (line 19)
   • Also adding the lib folder to LD_LIBRARY_PATH in your ~/.bashrc:

      export LD_LIBRARY_PATH=$LD_LIBRARY_PATH:/panfs/panasas01/phys/ra12451/boost_1_55_0_install/lib/
     

8. Note that if you use another boost library (like Filesystem) you'll need to add -lboost_XXX to DELPHES_LIBS

9. For Boost Program Options, you can now follow the tutorial.

LD_LIBRARY_PATH=/cm/shared/apps/gcc/4.7.0/lib:/cm/shared/apps/gcc/4.7.0/lib64:/cm/shared/languages/Python-2.7.6/lib:/cm/shared/apps/torque/4.2.4.1/lib:/cm/shared/apps/moab/7.2.2/lib:/cm/shared/tools/git-1.8.4.2/lib:/panfs/panasas01/phys/ra12451/boost_1_55_0_install/lib/:/panfs/panasas01/phys/ra12451/root/root/lib/

First time setup:

svn co -N svn+ssh://svn.cern.ch/reps/tdr2 myDir
# svn co -N svn+ssh://raggleto@svn.cern.ch/reps/tdr2 myDir
cd myDir
svn update utils
svn update -N notes
svn update notes/AN-14-007
eval `notes/tdr runtime -sh` # (for bash. use -csh for tcsh.)
cd notes/AN-14-007/trunk/
tdr --style=an b AN-14-007

To update/remake:

cd myDir
svn update notes/AN-14-007
eval `notes/tdr runtime -sh` # (for bash. use -csh for tcsh.)
cd notes/AN-14-007/trunk/
tdr --style=an b AN-14-007

See https://twiki.cern.ch/twiki/bin/view/CMS/Internal/TdrProcessing