| 22 | Avoiding soft and collinear singularities (very mildly): |
| 23 | |
| 24 | {{{ |
| 25 | cuts = all E >= 100 MeV [prt] |
| 26 | and all abs (cos(Theta)) <= 0.99 [prt] |
| 27 | and all M2 >= (1 GeV)^2 [prt, prt] |
| 28 | }}} |
| 29 | |
| 30 | Define title etc. as global variables, that will be used by PLOT: |
| 31 | |
| 32 | {{{ |
| 33 | $description = "A WHIZARD 2.0 Example" |
| 34 | $xlabel = "$\sqrt{s}$/GeV" |
| 35 | $ylabel = "$\sigma(s)$/pb" |
| 36 | }}} |
| 37 | |
| 38 | Allocate one plot |
| 39 | |
| 40 | {{{ |
| 41 | $title = "The Z Lineshape in $e^+e^-\to\mu^+\mu^-$" |
| 42 | plot lineshape_born (88 GeV, 95 GeV) |
| 43 | }}} |
| 44 | |
| 45 | Allocate another plot |
| 46 | |
| 47 | {{{ |
| 48 | $title = "The Z Lineshape in $e^+e^-\to\mu^+\mu^-\gamma$" |
| 49 | plot lineshape_rc (88 GeV, 95 GeV) |
| 50 | }}} |
| 51 | |
| 52 | Scanning in steps over the collider energy, integrating both Born and radiative corrections: |
| 53 | |
| 54 | {{{ |
| 55 | scan sqrts = (88.0 GeV => 90.0 GeV /+ 0.5 GeV, |
| 56 | 90.1 GeV => 91.9 GeV /+ 0.1 GeV, |
| 57 | 92.0 GeV => 95.0 GeV /+ 0.5 GeV) { |
| 58 | beams = e1, E1 |
| 59 | integrate (born) { iterations = 1:2000 } |
| 60 | record lineshape_born (sqrts, integral (born) / 1000) |
| 61 | integrate (rc) { iterations = 5:5000 } |
| 62 | record lineshape_rc (sqrts, integral (rc) / 1000) |
| 63 | } |
| 64 | }}} |
| 65 | |
| 66 | Write out the analysis: |
| 67 | |
| 68 | {{{ |
| 69 | $analysis_filename = "Z-lineshape" |
| 70 | write_analysis |
| 71 | }}} |