Generating Higgs To 4 Muons at NIKHEF

An exercise to simulated Higgs production events at the LHC, where the Higgs boson decays into 2 Z bosons that each decay into 2 muons.

$H\rightarrow ZZ^{*}\rightarrow \mu ^{+}\mu ^{-}\mu ^{+}\mu ^{-}$

The exercise is ment as a starting point for the 'monkey-see monkey-do' technique. In this example we will use AtlFast for the detector simulation and reconstruction. We will produce an AOD that contains the MC truth and reconstructed AtlFast objects. Since the AOD is in pool format we will also transform the AOD into an Ntuple that allows a simple analysis program to be constructed in Root.

Note: We assume you have the CMT and Athena set-up at NIKHEF in ordnung Starting with CMT and Athena at NIKHEF

1) Setting up the ATLAS environment at NIKHEF

a) Setting up the general ATLAS environment at NIKHEF (first time only)

First set up the ATLAS environment at NIKHEF. Follow the instructions on: ATLAS setup at NIKHEF.

b) Setting up the Package required to produce Ntuples from the AOD (first time only)

To produce Ntuples from an AOD you'll need to add an additional package created at NIKHEF.

1. Go to the directory where you want to install the package: cd /project/atlas/users/<your_login_name>
2. Check out the TTBarAnalysis package from the NIKHEF/ATLAS CVS repository: cvs -d /project/atlas/cvs co TTBarAnalysis
3. Go to the cmt directory: cd TTBarAnalysis/cmt
4. Execute cmt config
5. Execute source setup.csh
6. Build the library: gmake

You can also get a more detailed set of instructions from Installing the AOD->Ntuple (TTBarabalysis) package.

Once this is set-up you can produce TopNtuples from an AOD

c) Setting up all required packages (every time)

On every login you should now make sure the shell knows where to get the various programs, which means both the ATLAS general and the Ntuple Make program. You can do this by simply sourcing a script similar to init1002.csh. Simply source it in every window where you want to do the generation: source init1002.csh

2) Generating Higgs events

a) Download the scripts

Go to your favorite area and create a running directory and download the code. At NIKHEF a logical place would be again your project disk:

cd /project/atlas/users/<your_login_name>

cvs -d /project/atlas/cvs co Higgs4MuonAnalysis

cd Higgs4MuonAnalysis

Let's have a look at what files are in the package.

Athena requires steering files telling it what to do. These files are called joboptions files and since this exercise is made up of 2 steps we need 2 joboptions files: 1) From Pythia -> AOD: jobOptions_Pythia_To_Atlfast_To_AOD_BASIC.py 2) From AOD -> TopNtuple: jobOptions_AOD_to_Ntuple_BASIC.py

Then there is the main running script and some strange file Athena requires for various tasks:

3) ShipOff_Pythia.py The script that generates events 4) PDGTABLE.MeV A steering file required for MC production in Athena (not to be edited)

Finished!

Finally, chaning the Number of events to 1000 and the number of jobs to 10, you will produce 10,000 events with $H\rightarrow ZZ^{*}\rightarrow \mu ^{+}\mu ^{-}\mu ^{+}\mu ^{-}$.