1 | model = PSSSM |
---|
2 | |
---|
3 | process dec_lqino_tb = lqino1 => u,e1,neu1 { $restrictions = "3+4~'se1-'" } |
---|
4 | process sig1 = u,g => SE11, lqino1 |
---|
5 | |
---|
6 | compile |
---|
7 | |
---|
8 | ?slha_read_decays = true |
---|
9 | read_slha("sps1aNMSSMed.in") |
---|
10 | |
---|
11 | !if integrate is commented out, the histogram gets faulty and event generation takes longer by a factor of 100?! |
---|
12 | !integrate (dec_lqino_tb) |
---|
13 | |
---|
14 | mlqino_1 = 600 |
---|
15 | wlqino_1 = 5.87 |
---|
16 | yuk_nd = 0.312 |
---|
17 | yuk_d1 = 0.312 |
---|
18 | yuk_d2 = 0 |
---|
19 | yuk_dc = 0.312 |
---|
20 | |
---|
21 | sqrts = 14000 |
---|
22 | beams = p, p => lhapdf |
---|
23 | |
---|
24 | integrate (sig1) { iterations = 5:50000, 2:50000 } |
---|
25 | |
---|
26 | n_events = 10000 |
---|
27 | |
---|
28 | histogram inv_mass1_tb (0,1000,10) |
---|
29 | |
---|
30 | unstable lqino1 (dec_lqino_tb) |
---|
31 | |
---|
32 | simulate (sig1) { |
---|
33 | |
---|
34 | $sample = "test" |
---|
35 | sample_format = lhef |
---|
36 | checkpoint = 1000 |
---|
37 | |
---|
38 | analysis = |
---|
39 | record inv_mass1_tb (eval M [combine[u,e1]]) |
---|
40 | } |
---|
41 | |
---|
42 | write_analysis |
---|
43 | compile_analysis { $out_file = "sigtehest.dat"} |
---|
44 | |
---|
45 | |
---|
46 | |
---|
47 | |
---|
48 | |
---|
49 | |
---|