1 | model = SM |
---|
2 | |
---|
3 | alias n = n1:n2:n3 |
---|
4 | alias N = N1:N2:N3 |
---|
5 | alias l = e1:e2 |
---|
6 | alias L = E1:E2 |
---|
7 | alias q = u:d:s:c:b |
---|
8 | alias Q = U:D:S:C:B |
---|
9 | |
---|
10 | # Set the input parameters |
---|
11 | # Would-be 115 GeV Higgs at LEP |
---|
12 | mH = 125 GeV |
---|
13 | wH = 4.30 MeV |
---|
14 | # Running b mass |
---|
15 | me = 0 |
---|
16 | mmu = 0 |
---|
17 | ms = 0 |
---|
18 | mc = 0 |
---|
19 | mb = 0 |
---|
20 | |
---|
21 | process qqh = e1, E1 => q, Q, h |
---|
22 | process hdecayl1 = h => l, L, n, N{$restrictions = "2+3~Z && 4+5~Z"} |
---|
23 | |
---|
24 | sqrts = 250 |
---|
25 | |
---|
26 | compile |
---|
27 | ?vis_history = false |
---|
28 | seed = 10000 |
---|
29 | |
---|
30 | ?ps_fsr_active = true |
---|
31 | $shower_method = "PYTHIA6" |
---|
32 | ?hadronization_active = true |
---|
33 | |
---|
34 | ps_max_n_flavors = 5 |
---|
35 | ps_mass_cutoff = 1 |
---|
36 | |
---|
37 | integrate(qqh){ beams = e1, E1 => isr isr_mass = 0.000511} |
---|
38 | $circe2_file = "cepc250.circe" |
---|
39 | $circe2_design = "CEPC" |
---|
40 | integrate(hdecayl1) |
---|
41 | |
---|
42 | unstable h(hdecayl1) |
---|
43 | |
---|
44 | n_events = 1 |
---|
45 | ?write_raw = false |
---|
46 | |
---|
47 | sample_format = lhef, debug, stdhep |
---|
48 | simulate (qqh) |
---|
49 | |
---|
50 | show ( results ) |
---|