whizard is hosted by Hepforge, IPPP Durham

Changes between Version 2 and Version 3 of ExampleZprime


Ignore:
Timestamp:
Apr 12, 2010, 9:28:04 AM (14 years ago)
Author:
ohl
Comment:

Pretty printing

Legend:

Unmodified
Added
Removed
Modified
  • ExampleZprime

    v2 v3  
    1 Steered by this file
     1Choose the model
    22{{{
    33model = Zprime
    4 
     4}}}
     5Setup aliases for partons (quarks and gluons) and leptons (electrons and positrons)
     6{{{
    57alias pr = u:d:s:c:gl:U:D:S:C
    68alias lepton = e1:E1
     9}}}
     10Look at e+e--pairs created by colliding all possible parton pairs
     11{{{
    712process zp_drell_yan = pr, pr => e1, E1
    8 
    9 # Define model parameters
     13}}}
     14Set values for the model parameters (mass and width of the Z')
     15{{{
    1016mZH = 1500 GeV
    1117wZH = 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}}}
     19Define reasonable cuts
     20{{{
    1521cuts = all Pt >= 50 GeV [lepton]
    1622       and all -2.5 <= Eta <= 2.5 [e1:E1]
    1723       and all M >= 800 GeV [e1,E1]
    18 
    19 # Define the process scale
     24}}}
     25Define the process scale
     26{{{
    2027scale = 1500 GeV
    21 
    22 # Beam specification and integration
     28}}}
     29Set up the beams (full energy LHC with default lhapdf structure functions
     30{{{
    2331sqrts = 14 TeV
    24 beams = p,p => lhapdf
     32beams = p,p => lhapdf
     33}}}
     34integrate the cross sections in order to initialize the phase space grids for similation
     35{{{
    2536integrate (zp_drell_yan) { iterations = 6:20000 }
    26 
     37}}}
     38choose weighted events for smoother histograms
     39{{{
    2740?unweighted = false
    28 
    29 # Define title etc. as global variables, that will be used by PLOT
     41}}}
     42define title etc. as global variables, that will be used by PLOT
     43{{{
    3044$description =
    3145  "A WHIZARD 2.0 Example.
     
    3549   distributions."
    3650$ylabel = "$N_{\textrm{events}}$"
    37 
    38 # Allocate plots
     51}}}
     52allocate plots
     53{{{
    3954$title = "$Z'$ Drell-Yan dilepton peak $pp\to \ell\ell$"
    4055$xlabel = "$\theta^{\ell}$/GeV"
     
    4459$xlabel = "$p_T^{\ell}$/GeV"
    4560histogram pt_lepton (100 GeV, 2000 GeV, 100 GeV)
    46  
     61
    4762$title = "$Z'$ Drell-Yan: Angular distribution of the $\ell$"
    4863$xlabel = "$\theta^{\ell}$/rad"
    4964histogram th_lepton (0, 3.1415, 0.157)
    50 
     65}}}
     66and tell WHIZARD what to store in them
     67{{{
    5168analysis = record m_lepton (eval M [e1,E1]) and
    5269           record pt_lepton (eval Pt [e1]) and
    5370           record th_lepton (eval Theta [e1])
    54 
     71}}}
     72simulate 100000 events
     73{{{
    5574n_events = 100000
    56 
    5775simulate (zp_drell_yan)
    58 
     76}}}
     77write results to file (and run LaTeX to produce plots etc.)
     78{{{
    5979$analysis_filename = "Zprime"
    6080}}}
     81
    6182WHIZARD will produce histograms in the
    6283[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