Changes between Initial Version and Version 1 of ExampleWendpoint

Timestamp:

Apr 12, 2010, 10:21:45 AM (15 years ago)

Author:

ohl

Comment:

--

ExampleWendpoint

@@ +1 @@
+{{{
+model = SM
+}}}
+Set up the process
+{{{
+alias parton = u:U:d:D:g
+alias jet = parton
+}}}
+We could have added strange and charm jets:
+{{{
+# alias jet = parton:s:S:c:C
+}}}
+{{{
+alias lepton = e1:e2
+alias neutrino = n1:N1:n2:N2
+}}}
+{{{
+process enj = parton, parton => lepton, neutrino, jet
+compile
+}}}
+LHC design energy
+{{{
+sqrts = 14 TeV
+beams = p, p => lhapdf { $lhapdf_file = "cteq5l.LHgrid" }
+}}}
+Define reasonable cuts and integrate the cross section(s) in order to initialize the phase space grids for similation
+{{{
+cuts = all Pt >= 10 GeV [jet:lepton]
+integrate (enj) { iterations = 5:20000 }
+}}}
+Define title etc. as global variables, that will be used by PLOT
+{{{
+$description =
+  "A WHIZARD 2.0 Example."
+$ylabel = "$N_{\textrm{events}}$"
+}}}
+Allocate plots
+{{{
+$title = "$W$ Endpoint (``Jacobian Peak'') in $pp\to \ell\bar\nu j$"
+$xlabel = "$p_T^\ell$/GeV"
+histogram pt_lepton (0 GeV, 80 GeV, 2 GeV)
+}}}
+{{{
+$title = "Jet-$p_T$ in $pp\to \ell\bar\nu j$"
+$xlabel = "$p_T^j$/GeV"
+histogram pt_jet (0 GeV, 80 GeV, 2 GeV)
+}}}
+{{{
+$title = "Lepton Energy in $pp\to \ell\bar\nu j$"
+$xlabel = "$E^\ell$/GeV"
+histogram e_lepton (0 GeV, 160 GeV, 4 GeV)
+}}}
+{{{
+$title = "Jet Energy in $pp\to \ell\bar\nu j$"
+$xlabel = "$E^j$/GeV"
+histogram e_jet (0 GeV, 80 GeV, 2 GeV)
+}}}
+Record transverse momenta and energies
+{{{
+analysis = record pt_lepton (eval Pt [extract index 1 [sort by Pt [lepton]]])
+       and record pt_jet (eval Pt [extract index 1 [sort by Pt [jet]]])
+       and record e_lepton (eval E [extract index 1 [sort by Pt [lepton]]])
+       and record e_jet (eval E [extract index 1 [sort by Pt [jet]]])
+}}}
+{{{
+simulate (enj) { n_events = 1000 }
+}}}
+{{{
+$analysis_filename = "W-endpoint"
+write_analysis
+}}}