Page MenuHomeHEPForge

DrellYanNoMatchingP.sin
No OneTemporary

DrellYanNoMatchingP.sin

# $Id: DrellYanNoMatchingP.sin 2293 2010-04-11 23:57:50Z jr_reuter $
# Simple complete physics example: DY with PYTHIA shower
########################################################################
#
# Copyright (C) 1999-2015 by
# Wolfgang Kilian <kilian@physik.uni-siegen.de>
# Thorsten Ohl <ohl@physik.uni-wuerzburg.de>
# Juergen Reuter <juergen.reuter@desy.de>
# with contributions from
# Sebastian Schmidt
# Christian Speckner <cnspeckn@googlemail.com>
#
# WHIZARD is free software; you can redistribute it and/or modify it
# under the terms of the GNU General Public License as published by
# the Free Software Foundation; either version 2, or (at your option)
# any later version.
#
# WHIZARD is distributed in the hope that it will be useful, but
# WITHOUT ANY WARRANTY; without even the implied warranty of
# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
# GNU General Public License for more details.
#
# You should have received a copy of the GNU General Public License
# along with this program; if not, write to the Free Software
# Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA.
#
########################################################################
model = SM
alias quark = u:d:s:c
alias antiq = U:D:S:C
process qqee = quark:antiq, quark:antiq => e1, E1
ms = 0
mc = 0
compile
sqrts = 1960 GeV
beams = p, pbar => lhapdf
$title = "Pt of Z"
$x_label = "Pt / GeV"
$y_label = "N"
histogram Pt_distribution (0, 250, 1)
ps_isr_primordial_kt_width = 0
ps_isr_tscalefactor = 1
?ps_fsr_active = true
?ps_isr_active = true
?hadronization_active = true
$shower_method = "PYTHIA6"
?ps_isr_pt_ordered = false
?ps_isr_only_onshell_emitted_partons = false
ps_max_n_flavors = 5
ps_mass_cutoff = 1
### Rick Field's CDF Tune D6 using CTEQ6L1
### MSTP(5)=108 : Rick Field's CDF Tune D6 using CTEQ6L1
### MSTP(81)=0 : no multiple interactions
### MSTJ(41)=1 : only QCD branchings
### MSTJ(21)=0 : no decays
### MSTP(68)=0 : hopefully turn off PYTHIA's matching
$ps_PYTHIA_PYGIVE = "MSTP(5)=108;MSTP(81)=0;MSTJ(41)=1;MSTJ(21)=0;MSTP(68)=0"
?mlm_matching = false
cuts = all M > 80 GeV [combine [e1,E1]]
and all M < 100 GeV [combine [e1,E1]]
integrate (qqee) { iterations = 2:20000 }
n_events = 100000
?rebuild_events = true
$sample = "DrellYanNoMatchingP"
sample_format = lhef
analysis = record Pt_distribution (eval Pt [combine[e1,E1]])
checkpoint = 2500
simulate (qqee)
show (results)

File Metadata

Mime Type
text/plain
Expires
Wed, May 14, 11:48 AM (1 h, 55 m)
Storage Engine
blob
Storage Format
Raw Data
Storage Handle
5111520
Default Alt Text
DrellYanNoMatchingP.sin (2 KB)

Event Timeline