001 path enable AllPath
002 path list
003 module input WGRAD
004 #
005 mod talk WGRAD
006 debug set f
007 usePythia set t
008 cardsMenu
009 W_Charge set 1
010 widthTreatment set 1
011 DecayLepton set 1
012 QED_RadType set 2
013 W_Mass set 80.4
014 W_Width set 2.1
015 CME set 1800.0
016 Order set 2
017 Scale set 0
018 fsrCollinearCut set 0
019 PhotonEnergyCut set 0.01
020 PhotonCollinearCut set 0.001
021 Lepton_ptmin set 25.0
022 Lepton_ymax set 5.0
023 Neutrino_ptmin set 25.0
024 # Lepton_ptmin set 25.0
025 # Lepton_ymax set 1.1
026 # Neutrino_ptmin set 25.0
027 mT_low set 0.0
028 mT_high set 200.0
029 itmaxGrid set 5
030 ncallGrid set 5000
031 itmaxPrecision set 10
032 ncallPrecision set 10000
033 # itmaxGrid set 5
034 # ncallGrid set 10000
035 # itmaxPrecision set 10
036 # ncallPrecision set 250000
037 exit
038 pythiaInterfaceMenu
039 #######################################################
040 # If not set, Pythia defaults will be substituted #
041 #######################################################
042 # Fragmentation and decay : 0 - off, 1 - on (default) :
043 mstp(111) set 1
044 # 0 - no showers, 1 - QCD only, 2 - QCD+QED (default) :
045 mstj(41) set 1
046 # Random number seed :
047 mrpy(1) set 280273
048 # multiple interactions : 0 - off, 1 - on (default) :
049 mstp(81) set 0
050 # intrinsic kT : 0 - none, 1 - Gaussian (default), 2 - exponential :
051 mstp(91) set 1
052 # kT Gaussian width (GeV) :
053 parp(91) set 4.0
054 # kT exponential width (GeV) :
055 parp(92) set 1.63
056 # Upper bound on intrinsic kT (GeV) :
057 parp(93) set 20.0
058 exit
059 exit
060 #
061 mod talk HepRootManager
062 # Name of file containing output histograms :
063 histfile set wgrad.root
064 exit
065 #
066 ev begin -nev 1000
067 # show timer
068 exit
Send problems or questions to cdfcode@fnal.gov