...
| No Format |
|---|
$ cd MadGraph5_v1_5_11
$ ./bin/mg5
mg5> install pythia-pgs
|
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 |
|---|
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, eg., `add process mu+ mu- > Z > b b~' to add moreanother 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
|
Important note: If you want to do mu+ mu- > h > x x events to load this interaction vertex, be sure to execute `mg5> import model sm-full' before launching. The mu+ mu- > h vertex is disabled in the default sm model for performance reasons.
Now we can generate events and simulate their hadronization with Pythia.
...
| Code Block | ||||
|---|---|---|---|---|
| ||||
/physics/select FTFP_BERT
/physics/setPDGFile /path/to/slic/v00-00/slic/v03-01-04/data/particle.tbl
/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
|
...
(...Installation of slicPandora and PandoraSettings...)
| 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>
|
...
WHIZARD 2 is a good tool for more precise cross section calculations because it computes every possible channel between the given initial and final states. One necessary restriction of this is that unlike MadGraph, you do not have control over the intermediate channels, such as whether your mu+ mu- > b b~ process is mediated by a Z or h boson. In addition, stdhep files generated by WHIZARD do not retain little Monte Carlo information about the event. Refer to the website for instructions on installing and using WHIZARD 2. We provide here a sample SINDARIN input file for a muon collider at 3TeV.
...