00001 # in advance switch of Process semihard QCD 2 -> 2 Process
00002 pysubs msub 96 0
00003
00004 # Selecting jet + gamma*/Z process (number 13):
00005 pysubs msel 13
00006
00007 # randomnumber
00008 pydatr mrpy 1 97823417
00009
00010
00011 pydat3 mdme 1 174 0
00012 pydat3 mdme 1 175 0
00013 pydat3 mdme 1 176 0
00014 pydat3 mdme 1 177 0
00015 pydat3 mdme 1 178 0
00016 pydat3 mdme 1 179 0
00017 pydat3 mdme 1 180 0
00018 pydat3 mdme 1 181 0
00019 # electrons
00020 pydat3 mdme 1 182 1
00021 pydat3 mdme 1 183 0
00022 # muons
00023 pydat3 mdme 1 184 1
00024 pydat3 mdme 1 185 0
00025 #taus
00026 pydat3 mdme 1 186 0
00027 pydat3 mdme 1 187 0
00028 pysubs ckin 1 60.
00029
00030 # set pdfsets
00031 pypars mstp 51 7
00032 pypars mstp 52 1
00033
00034 # needed for CTEQ6L
00035 #pypars mstp 51 10042
00036 #pypars mstp 52 2
00037 #pypars mstp 53 10042
00038 #pypars mstp 54 2
00039 #pypars mstp 55 0042
00040 #pypars mstp 56 2
00041
00042 pydat2.pmas[6-1][1-1] = 175.;
00043 pydat2.pmas[24-1][1-1] =80.42;
00044 pydat2.pmas[24-1][2-1] =2.124;
00045 pydat2.pmas[23-1][1-1] =91.19;
00046 pydat2.pmas[23-1][2-1] =2.495;
00047
00048 # fix junk output for documentary particles
00049 pypars mstp 128 1
00050 # error handling switch
00051 pydat1 mstu 21 1
00052
00053 # ATLAS default for MI and also for the associated
00054 pypars mstp 81 21
00055
00056 # Turn off Multiple Interactions (MI)
00057 pypars mstp 81 0
00058 # Turn off Final state Radiation (FSR)
00059 pypars mstp 71 0
00060 # Turn off Initial state Radiation (ISR)
00061 pypars mstp 61 1
00062
00063 # select peterson for charm fragmentation
00064 pydat1 mstj 11 3
00065 # Make K_S, Lambda stable
00066 pydat1 mstj 22 2
00067 #c hadronization
00068 pydat1 parj 54 -0.07
00069 #b hadronization
00070 pydat1 parj 55 -0.006
00071
00072 pypars parp 89 1800.
00073
00074 pypars mstp 81 1
00075 pypars parp 67 3.
00076 pypars parp 82 2.2
00077 pypars parp 83 0.5
00078 pypars parp 84 0.3
00079 pypars parp 85 0.33
00080 pypars parp 86 0.66
00081 pypars parp 90 0.24
00082
00083 #mutiple interaction model
00084 pypars mstp 82 4
00085 pypars mstp 86 2
00086