close
Warning:
Error with navigation contributor "BrowserModule"
- Timestamp:
-
Apr 12, 2010, 10:16:11 AM (15 years ago)
- Author:
-
ohl
- Comment:
-
--
Legend:
- Unmodified
- Added
- Removed
- Modified
-
v2
|
v3
|
|
| 1 | Select the minimal supersymmetric standard model |
1 | 2 | {{{ |
2 | 3 | model = MSSM |
3 | 4 | }}} |
| 5 | and consider a decay cascade of up-squarks into neutralinos and selectrons |
4 | 6 | {{{ |
5 | 7 | process dec_su_q = su1 => u, neu2 |
6 | 8 | process dec_neu_sl2 = neu2 => SE12, e1 |
7 | 9 | }}} |
| 10 | where the up-squark pairs are generated from up quarks |
8 | 11 | {{{ |
9 | 12 | process susybg = u,U => SU1, su1 |
| 13 | }}} |
| 14 | We will compare the results with the full process including all irreducible backgrounds |
| 15 | {{{ |
10 | 16 | process full = u, U => SU1, u, e1, SE12 |
11 | 17 | }}} |
| 18 | Call O'Mega and compile the matrix elements |
12 | 19 | {{{ |
13 | 20 | compile |
14 | 21 | }}} |
| 22 | We use the SLHA convention to choose the SPS1a parameter point |
15 | 23 | {{{ |
16 | 24 | ?slha_read_decays = true |
17 | 25 | read_slha("sps1ap_decays.slha") |
18 | 26 | }}} |
| 27 | Obtain the total widths (simple integration, no singularities, not too many points required) |
19 | 28 | {{{ |
20 | 29 | integrate (dec_su_q, dec_neu_sl2) { iterations = 1:1000 } |
21 | 30 | }}} |
| 31 | We're at the full energy LHC |
22 | 32 | {{{ |
23 | 33 | sqrts = 14000 |
24 | 34 | beams = p, p => lhapdf |
25 | 35 | }}} |
| 36 | Obtain the up-squark production cross section and the cross section for full process |
26 | 37 | {{{ |
27 | 38 | integrate (susybg) { iterations = 5:10000, 2:10000 } |
28 | 39 | integrate (full) |
29 | 40 | }}} |
| 41 | Simulate 10K events |
30 | 42 | {{{ |
31 | 43 | n_events = 10000 |
32 | 44 | }}} |
| 45 | Allocate histograms |
33 | 46 | {{{ |
34 | 47 | $title = "Full process" |
… |
… |
|
53 | 66 | histogram inv_mass1_iso (0,600,20) |
54 | 67 | }}} |
| 68 | Simulate the full process |
55 | 69 | {{{ |
56 | 70 | simulate (full) { |
… |
… |
|
60 | 74 | } |
61 | 75 | }}} |
| 76 | Simulate the cascade decay, including all off-diagonal spin correlations |
62 | 77 | {{{ |
63 | 78 | unstable neu2 (dec_neu_sl2) |
… |
… |
|
70 | 85 | } |
71 | 86 | }}} |
| 87 | Simulate the cascade decay, keeing only the classical (diagonal) spin correlations |
72 | 88 | {{{ |
73 | 89 | ?diagonal_decay = true |
… |
… |
|
80 | 96 | } |
81 | 97 | }}} |
| 98 | Simulate the cascade decay, dropping all spin correlations (i.e. assume isotropic decays) |
82 | 99 | {{{ |
83 | 100 | ?isotropic_decay = true |
… |
… |
|
90 | 107 | } |
91 | 108 | }}} |
| 109 | Write the results to file |
92 | 110 | {{{ |
93 | | write_analysis |
94 | 111 | $analysis_filename = "casc_dec" |
95 | 112 | write_analysis |