-
Notifications
You must be signed in to change notification settings - Fork 0
/
powheg.input
119 lines (99 loc) · 3.87 KB
/
powheg.input
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
!randomseed 352345 ! uncomment to set the random seed to a value of your choice.
! It generates the call RM48IN(352345,0,0) (see the RM48 manual).
! THIS MAY ONLY AFFECTS THE GENERATION OF POWHEG EVENTS!
! If POWHEG is interfaced to a shower MC, refer to the shower MC
! documentation to set its seed.
!Heavy flavour production parameters
numevts 5000 ! number of events to be generated
ih1 1 ! hadron 1
ih2 1 ! hadron 2
#ndns1 131 ! pdf for hadron 1 (hvqpdf numbering)
#ndns2 131 ! pdf for hadron 2
!lhans1 47776000 ! pdf set for hadron 1 (LHA numbering)
!lhans2 47776000 ! pdf set for hadron 2 (LHA numbering)
lhans1 82400 ! pdf set for hadron 1 (LHA numbering)
lhans2 82400 ! pdf set for hadron 2 (LHA numbering)
ebeam1 6500 ! energy of beam 1
ebeam2 6500 ! energy of beam 2
facscfact 1 ! factorization scale factor: mufact=muref*facscfact
renscfact 1 ! renormalization scale factor: muren=muref*renscfact
#fixedscale 1 ! use ref. scale=qmass (default 0, use running scale)
runningscale 0
! Parameters to allow-disallow use of stored data
use-old-grid 1 ! if 1 use old grid if file pwggrids.dat is present (# 1: regenerate)
use-old-ubound 1 ! if 1 use norm of upper bounding function stored in pwgubound.dat, if present; # 1: regenerate
ncall1 30000 ! number of calls for initializing the integration grid
itmx1 1 ! number of iterations for initializing the integration grid
ncall2 50000 ! number of calls for computing the integral and finding upper bound
itmx2 1 ! number of iterations for computing the integral and finding upper bound
foldcsi 1 ! number of folds on x integration
foldy 1 ! number of folds on y integration
foldphi 1 ! number of folds on phi integration
nubound 50000 ! number of bbarra calls to setup norm of upper bounding function
iymax 1 ! <= 10, normalization of upper bounding function in iunorm X iunorm square in y, log(m2qq)
ixmax 1 ! <= 10, normalization of upper bounding function in iunorm X iunorm square in y, log(m2qq)
xupbound 2 ! increase upper bound for radiation generation
manyseeds 0 ! (default 0) allow for the generation of different statistically independent samples (see Docs/Manyseeds.pdf)
parallelstage 0 ! which stage of parallel run 1..4
xgriditeration 0 ! which grid iteration when in parallelstage 1
#evenmaxrat 1
testplots 1
g4 1
mU 2500 ! Mass of the LQ
smear 0
#widthLQ 39.8d0
BWgen 0
LQmasslow 200
LQmasshigh 5000
ptcutlep 0
ptcutjet 0
etacutlep 0
etacutjet 0
mLQcuthi 0
mLQcutlo 0
recombination 1
#bornonly 1
#btildevirt 0
#btildecoll 0
#softmismch 0
#btildeborn 0
#softtest 0
#colltest 0
bornsmartsig 0
! / y_1e y_1m y_1t \ u/d
! | y_2e y_2m y_2t | c/s
! \ y_3e y_3m y_3t / t/b
betaL_1e 0
betaL_2e 0
betaL_3e 0
betaL_1m 0
betaL_2m 0
betaL_3m 0
betaL_1t 0
betaL_2t 0
betaL_3t 0
betaR_1e 0
betaR_2e 0
betaR_3e 0
betaR_1m 0
betaR_2m 0
betaR_3m 0
betaR_1t 0
betaR_2t 0
betaR_3t 0
charge 2 ! Set this to the charge of the desired LQ's absolute charge times 3. Expect 1,2,4 or 5
#------------------------------------------------
# SCALE VARIATION & REWEIGHTING
#------------------------------------------------
storeinfo_rwgt 1 ! store info to allow for reweighting
rwl_file '-'
<initrwgt>
<weightgroup name='First-Weights'>
<weight id='11'> renscfact=1.0 facscfact=1.0 </weight>
<weight id='12'> renscfact=1.0 facscfact=2.0 </weight>
<weight id='21'> renscfact=2.0 facscfact=1.0 </weight>
<weight id='22'> renscfact=2.0 facscfact=2.0 </weight>
<weight id='1H'> renscfact=1.0 facscfact=0.5 </weight>
<weight id='H1'> renscfact=0.5 facscfact=1.0 </weight>
<weight id='HH'> renscfact=0.5 facscfact=0.5 </weight>
</initrwgt>