|
- Timestamp:
-
Apr 12, 2010, 9:28:04 AM (15 years ago)
- Author:
-
ohl
- Comment:
-
Pretty printing
Legend:
- Unmodified
- Added
- Removed
- Modified
-
v2
|
v3
|
|
1 | | Steered by this file |
| 1 | Choose the model |
2 | 2 | {{{ |
3 | 3 | model = Zprime |
4 | | |
| 4 | }}} |
| 5 | Setup aliases for partons (quarks and gluons) and leptons (electrons and positrons) |
| 6 | {{{ |
5 | 7 | alias pr = u:d:s:c:gl:U:D:S:C |
6 | 8 | alias lepton = e1:E1 |
| 9 | }}} |
| 10 | Look at e+e--pairs created by colliding all possible parton pairs |
| 11 | {{{ |
7 | 12 | process zp_drell_yan = pr, pr => e1, E1 |
8 | | |
9 | | # Define model parameters |
| 13 | }}} |
| 14 | Set values for the model parameters (mass and width of the Z') |
| 15 | {{{ |
10 | 16 | mZH = 1500 GeV |
11 | 17 | wZH = 50 GeV |
12 | | |
13 | | # Define reasonable cuts and integrate the cross section(s) |
14 | | # in order to initialize the phase space grids for similation |
| 18 | }}} |
| 19 | Define reasonable cuts |
| 20 | {{{ |
15 | 21 | cuts = all Pt >= 50 GeV [lepton] |
16 | 22 | and all -2.5 <= Eta <= 2.5 [e1:E1] |
17 | 23 | and all M >= 800 GeV [e1,E1] |
18 | | |
19 | | # Define the process scale |
| 24 | }}} |
| 25 | Define the process scale |
| 26 | {{{ |
20 | 27 | scale = 1500 GeV |
21 | | |
22 | | # Beam specification and integration |
| 28 | }}} |
| 29 | Set up the beams (full energy LHC with default lhapdf structure functions |
| 30 | {{{ |
23 | 31 | sqrts = 14 TeV |
24 | | beams = p,p => lhapdf |
| 32 | beams = p,p => lhapdf |
| 33 | }}} |
| 34 | integrate the cross sections in order to initialize the phase space grids for similation |
| 35 | {{{ |
25 | 36 | integrate (zp_drell_yan) { iterations = 6:20000 } |
26 | | |
| 37 | }}} |
| 38 | choose weighted events for smoother histograms |
| 39 | {{{ |
27 | 40 | ?unweighted = false |
28 | | |
29 | | # Define title etc. as global variables, that will be used by PLOT |
| 41 | }}} |
| 42 | define title etc. as global variables, that will be used by PLOT |
| 43 | {{{ |
30 | 44 | $description = |
31 | 45 | "A WHIZARD 2.0 Example. |
… |
… |
|
35 | 49 | distributions." |
36 | 50 | $ylabel = "$N_{\textrm{events}}$" |
37 | | |
38 | | # Allocate plots |
| 51 | }}} |
| 52 | allocate plots |
| 53 | {{{ |
39 | 54 | $title = "$Z'$ Drell-Yan dilepton peak $pp\to \ell\ell$" |
40 | 55 | $xlabel = "$\theta^{\ell}$/GeV" |
… |
… |
|
44 | 59 | $xlabel = "$p_T^{\ell}$/GeV" |
45 | 60 | histogram pt_lepton (100 GeV, 2000 GeV, 100 GeV) |
46 | | |
| 61 | |
47 | 62 | $title = "$Z'$ Drell-Yan: Angular distribution of the $\ell$" |
48 | 63 | $xlabel = "$\theta^{\ell}$/rad" |
49 | 64 | histogram th_lepton (0, 3.1415, 0.157) |
50 | | |
| 65 | }}} |
| 66 | and tell WHIZARD what to store in them |
| 67 | {{{ |
51 | 68 | analysis = record m_lepton (eval M [e1,E1]) and |
52 | 69 | record pt_lepton (eval Pt [e1]) and |
53 | 70 | record th_lepton (eval Theta [e1]) |
54 | | |
| 71 | }}} |
| 72 | simulate 100000 events |
| 73 | {{{ |
55 | 74 | n_events = 100000 |
56 | | |
57 | 75 | simulate (zp_drell_yan) |
58 | | |
| 76 | }}} |
| 77 | write results to file (and run LaTeX to produce plots etc.) |
| 78 | {{{ |
59 | 79 | $analysis_filename = "Zprime" |
60 | 80 | }}} |
| 81 | |
61 | 82 | WHIZARD will produce histograms in the |
62 | 83 | [http://projects.hepforge.org/whizard/trac/attachment/wiki/ExampleZprime/Zprime.pdf PDF output] and as [http://projects.hepforge.org/whizard/trac/attachment/wiki/ExampleZprime/Zprime.dat data file]. The progress can be monitored on the terminal and in the |
|