...
...Link to LikelihoodData I generated...
Using WHIZARD 2 to calculate cross sections
At the time of writing, the Standard Model in MadGraph does not include the vertices mu+ mu- > h or e+ e- > h (See bug report). It is also not a reliable cross section calculator due to the facts that it can only deal with tree-level processes and it does not compute the effects of interference from diagrams that are not explicitly added by the user.
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 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.
| Code Block | ||
|---|---|---|
| ||
# SM_ac allows to change higgs tri-linear self-coupling (fac_gh3)
model = SM_ac
process mumu_nnbbbb = "mu-", "mu+" => numu, numubar, h, h
mH = 125 GeV
wH = 4.07 MeV
fac_gh3 = 1.2
sqrts = 3000 GeV
beams = "mu-", "mu+"
n_events = 10000
iterations = 10:10000, 5:10000, 5:50000
simulate (mumu_hhnn) { $sample = "/path/of/output/data/folder" sample_format = stdhep } |
This SINDARIN script calculates the cross section of double-Higgs production via W-fusion at a high energy muon collider. The cross section is dependent on the Higgs trilinear self-coupling constant, which can be expressed as a ratio to the Standard Model value. In this example the constant is 1.2 times the Standard Model value. Advanced SINDARIN scripts can perform simple analysis and plot histograms.