| Table of Contents |
|---|
Overview
This page will guide you through the process of installing and using some of the software we use to generate, fully simulate and reconstruct events in the MCDRCal01 detector for a Higgs Factory Muon Collider. We assume a Linux /UNIX environment.
We start with the MadGraph5 matrix element generator, a simple program that allows you to define a process or group of processes, calculate their cross sections, generate events and seamlessly interface with Pythia for Monte Carlo shower simulation (hadronization).
...
Note: If you already have the .hep or .stdhep files you want, or if you just want to use the particle gun in SLIC, you do not need to do the steps in this section; go on to Installing ILCSoft Programs with ilcinstall. If you have fully simulated .slcio events, go on to Analysis.
Download and Install MadGraph5 and Pythia
Download MadGraph 5 from the MadGraph 5 downloads page after making a free account and unpack the file somewhere.
| No Format |
|---|
$ tar xzvf MadGraph5_v1_5_11.tar.gz
|
Next, use MadGraph to download and install Pythia:
| No Format |
|---|
$ cd MadGraph5_v1_5_11
$ ./bin/mg5
mg5> install pythia
|
Note: If you have Python 3 installed on your system you may need to run MadGraph and associated programs with Python 2, eg.
| No Format |
|---|
$ python2 ./bin/mg5
|
Generate Events:
For this tutorial we will be simulating the process of two muons colliding with a center of mass energy of 125 GeV, creating a Z* boson which decays into two light quarks. This data will be used to generate detector-specific information that is used by the PandoraPFA photon reconstruction algorithm. A copy of this 'likelihood data' will be provided for the mcdrcal01 detector. For a more complete introduction to MadGraph, simply type mg5> tutorial or refer to the MadGraph 5 Wiki.
| No Format |
|---|
mg5> define uds = u d s u~ d~ s~ # define the multiparticle uds as any one of the u, d and s quarks or their antiparticles mg5> generate mu+ mu- > Z > uds uds # calculate the diagrams for the process # use `add process mu+ mu- > Z > b b~' to add more processes. mg5> output my-mumu-z-uds # create a folder in the current directory. mg5> open index.html # optional; view the generated diagrams and other information mg5> exit |
Now we can generate events and simulate their hadronization with Pythia.
| No Format |
|---|
> launch $> cd my-mumu-z-uds $ ./bin/madevent MGME5> launch MGME5> 22 # select to run MadEvent and Pythia MGME5> (select 1 to only calculate cross sections) > 2 # edit run card. eg. Change beam type (0 for muons or electrons, ie. no PDF), # set energy (62.5GeV each), edit cuts, number of events MGME5> 0 # none of the cuts apply to uds quarks, so in this example we will ignore them. > 1 # (optional) edit parameter card to change physics parameters such as Higgs mass and width. > 0 # Launches event generator and then runs Pythia. MGME5>> exit $ cd my-mumu-z-uds/ $ gzip -d Events/run_01/tag_1_pythia_events.hep.gz |
Now we have a .hep file which can be read into SLIC for full detector simulation. If you want to open the events in JAS3, change the extension to `.stdhep'
Installing ILCSoft programs with ilcinstall
...
ilcinstall is installed with the command
| No Format |
|---|
$ svn co https://svnsrv.desy.de/public/ilctools/ilcinstall/trunk/ ilcinstall
|
This should create a folder called ilcinstall in the directory from which you ran the command. The version used for this tutorial was v01-17-0204, rev 28202914. For more help installing and using ilcinstall, see the ilcinstall documentation.
...
Java_versioncan be found by running$ java -version. Here we use 1.7.0_40.Java_pathis the directory containingbin/java.java, which can be found with$ which javailcsoft_install_prefixwill be the name of the directory in which it installs programs, eg./home/$USER/ilcsoft/.
...
| Code Block | ||||
|---|---|---|---|---|
| ||||
/physics/select FTFP_BERT
/lcio/path /path/to/eventdata
/lcio/filename Z_uds_125GeV_iso
/lcio/fileExists delete
/lcdd/url /path/to/detectors/mcdrcal01/mcdrcal01.lcdd
/run/initialize
/vis/open HepRepFile
/vis/heprep/setFileDir /path/to/detectors/mcdrcal01/
/vis/heprep/setFileName mcdrcal01
/vis/heprep/setCullInvisibles true
/vis/scene/create
/vis/sceneHandler/attach
/vis/viewer/set/culling global
/vis/enable
/vis/viewer/flush
/lcdd/dumpGDML /path/to/detectors/mcdrcal01/mcdrcal01.gdml
/generator/select stdhep
/generator/filename /path/to/mumu-Z-uds-125GeV.hep
/run/beamOn 5000
|
(...Running on mcdrd or OSG?...)
| No Format |
|---|
<mcdrd.fnal.gov> nohup nice -m ~/path/to/slic/script
|
...
| Code Block | ||||||||
|---|---|---|---|---|---|---|---|---|
| ||||||||
<\!-\- Standalone photon clustering \-->
<algorithm type = "PhotonReconstruction">
<algorithm type = "ConeClustering" description = "PhotonClusterFormation">
<ClusterSeedStrategy>0</ClusterSeedStrategy>
<ShouldUseTrackSeed>false</ShouldUseTrackSeed>
<ShouldUseOnlyECalHits>true</ShouldUseOnlyECalHits>
<ConeApproachMaxSeparation>250.</ConeApproachMaxSeparation>
</algorithm>
<ClusterListName>PhotonClusters</ClusterListName>
<ShouldMakePdfHistograms>true</ShouldMakePdfHistograms>
<NEnergyBins>9</NEnergyBins>
<EnergyBinLowerEdges>0 0.5 1 1.5 2.5 5 10 20 50 </EnergyBinLowerEdges>
<HistogramFile>MyPath/MyLikelihoodData9EBins.xml</HistogramFile>
</algorithm>
|
...