| 1 | model = SM |
|---|
| 2 | |
|---|
| 3 | # top pair production |
|---|
| 4 | process ttbar = gl, gl => t, T |
|---|
| 5 | |
|---|
| 6 | # decay channels |
|---|
| 7 | process tlep = t => b, E1, n1 |
|---|
| 8 | process thad = T => B, U, d |
|---|
| 9 | |
|---|
| 10 | compile |
|---|
| 11 | |
|---|
| 12 | # cuts |
|---|
| 13 | cuts = all Pt > 20 GeV [b:E1:n1:B:U:d] |
|---|
| 14 | and all abs(Eta) < 2. [b:E1:n1:B:U:d] |
|---|
| 15 | |
|---|
| 16 | # integrate decays |
|---|
| 17 | integrate ( tlep, thad ) { iterations = 3:5000,7:10000 } |
|---|
| 18 | |
|---|
| 19 | # beam energy |
|---|
| 20 | sqrts = 1 TeV |
|---|
| 21 | |
|---|
| 22 | # unstable particles and decay chain |
|---|
| 23 | unstable t (tlep) |
|---|
| 24 | unstable T (thad) |
|---|
| 25 | |
|---|
| 26 | integrate ( ttbar ) { iterations = 3:5000,7:10000 } |
|---|
| 27 | |
|---|
| 28 | histogram pt_L ( 0 GeV , 200 GeV , 5 GeV ) |
|---|
| 29 | histogram pt_n ( 0 GeV , 200 GeV , 5 GeV ) |
|---|
| 30 | histogram pt_b ( 0 GeV , 200 GeV , 5 GeV ) |
|---|
| 31 | histogram pt_U ( 0 GeV , 200 GeV , 5 GeV ) |
|---|
| 32 | histogram pt_d ( 0 GeV , 200 GeV , 5 GeV ) |
|---|
| 33 | histogram pt_B ( 0 GeV , 200 GeV , 5 GeV ) |
|---|
| 34 | histogram eta_L ( -5. , 5. , .5 ) |
|---|
| 35 | histogram eta_n ( -5. , 5. , .5 ) |
|---|
| 36 | histogram eta_b ( -5. , 5. , .5 ) |
|---|
| 37 | histogram eta_U ( -5. , 5. , .5 ) |
|---|
| 38 | histogram eta_d ( -5. , 5. , .5 ) |
|---|
| 39 | histogram eta_B ( -5. , 5. , .5 ) |
|---|
| 40 | |
|---|
| 41 | analysis = record pt_L (eval Pt [E1]); |
|---|
| 42 | record pt_n (eval Pt [n1]); |
|---|
| 43 | record pt_b (eval Pt [b]); |
|---|
| 44 | record pt_U (eval Pt [U]); |
|---|
| 45 | record pt_d (eval Pt [d]); |
|---|
| 46 | record pt_B (eval Pt [B]); |
|---|
| 47 | record eta_L (eval Eta [E1]); |
|---|
| 48 | record eta_n (eval Eta [n1]); |
|---|
| 49 | record eta_b (eval Eta [b]); |
|---|
| 50 | record eta_U (eval Eta [U]); |
|---|
| 51 | record eta_d (eval Eta [d]); |
|---|
| 52 | record eta_B (eval Eta [B]) |
|---|
| 53 | |
|---|
| 54 | n_events = 10000 |
|---|
| 55 | |
|---|
| 56 | simulate ( ttbar ) |
|---|
| 57 | |
|---|
| 58 | compile_analysis { $out_file = "foo2.dat" } |
|---|
