19
19
#***********************************************************************
20
20
# Tag name for the run (one word) *
21
21
#***********************************************************************
22
tag_1 = run_tag ! name of the run
22
%(run_tag)s = run_tag ! name of the run
23
23
#***********************************************************************
24
24
# Number of LHE events (and their normalization) and the required *
25
25
# (relative) accuracy on the Xsec. *
26
26
# These values are ignored for fixed order runs *
27
27
#***********************************************************************
28
10000 = nevents ! Number of unweighted events requested
29
-1 = req_acc ! Required accuracy (-1=auto determined from nevents)
30
-1 = nevt_job! Max number of events per job in event generation.
28
%(nevents)s = nevents ! Number of unweighted events requested
29
%(req_acc)s = req_acc ! Required accuracy (-1=auto determined from nevents)
30
%(nevt_job)s = nevt_job! Max number of events per job in event generation.
32
32
#***********************************************************************
33
33
# Normalize the weights of LHE events such that they sum or average to *
34
34
# the total cross section *
35
35
#***********************************************************************
36
average = event_norm ! average or sum
36
%(event_norm)s = event_norm ! average or sum
37
37
#***********************************************************************
38
38
# Number of points per itegration channel (ignored for aMC@NLO runs) *
39
39
#***********************************************************************
40
0.01 = req_acc_FO ! Required accuracy (-1=ignored, and use the
40
%(req_acc_fo)s = req_acc_FO ! Required accuracy (-1=ignored, and use the
41
41
! number of points and iter. below)
42
42
# These numbers are ignored except if req_acc_FO is equal to -1
43
5000 = npoints_FO_grid ! number of points to setup grids
44
4 = niters_FO_grid ! number of iter. to setup grids
45
10000 = npoints_FO ! number of points to compute Xsec
46
6 = niters_FO ! number of iter. to compute Xsec
43
%(npoints_fo_grid)s = npoints_FO_grid ! number of points to setup grids
44
%(niters_fo_grid)s = niters_FO_grid ! number of iter. to setup grids
45
%(npoints_fo)s = npoints_FO ! number of points to compute Xsec
46
%(niters_fo)s = niters_FO ! number of iter. to compute Xsec
47
47
#***********************************************************************
48
48
# Random number seed *
49
49
#***********************************************************************
50
0 = iseed ! rnd seed (0=assigned automatically=default))
50
%(iseed)s = iseed ! rnd seed (0=assigned automatically=default))
51
51
#***********************************************************************
52
52
# Collider type and energy *
53
53
#***********************************************************************
54
1 = lpp1 ! beam 1 type (0 = no PDF)
55
1 = lpp2 ! beam 2 type (0 = no PDF)
56
6500 = ebeam1 ! beam 1 energy in GeV
57
6500 = ebeam2 ! beam 2 energy in GeV
54
%(lpp1)s = lpp1 ! beam 1 type (0 = no PDF)
55
%(lpp2)s = lpp2 ! beam 2 type (0 = no PDF)
56
%(ebeam1)s = ebeam1 ! beam 1 energy in GeV
57
%(ebeam2)s = ebeam2 ! beam 2 energy in GeV
58
58
#***********************************************************************
59
59
# PDF choice: this automatically fixes also alpha_s(MZ) and its evol. *
60
60
#***********************************************************************
61
nn23nlo = pdlabel ! PDF set
62
244600 = lhaid ! if pdlabel=lhapdf, this is the lhapdf number
61
%(pdlabel)s = pdlabel ! PDF set
62
%(lhaid)s = lhaid ! if pdlabel=lhapdf, this is the lhapdf number
63
63
#***********************************************************************
64
64
# Include the NLO Monte Carlo subtr. terms for the following parton *
65
65
# shower (HERWIG6 | HERWIGPP | PYTHIA6Q | PYTHIA6PT | PYTHIA8) *
66
66
# WARNING: PYTHIA6PT works only for processes without FSR!!!! *
67
67
#***********************************************************************
68
HERWIG6 = parton_shower
68
%(parton_shower)s = parton_shower
69
69
#***********************************************************************
70
70
# Renormalization and factorization scales *
71
71
# (Default functional form for the non-fixed scales is the sum of *
72
72
# the transverse masses of all final state particles and partons. This *
73
73
# can be changed in SubProcesses/set_scales.f) *
74
74
#***********************************************************************
75
F = fixed_ren_scale ! if .true. use fixed ren scale
76
F = fixed_fac_scale ! if .true. use fixed fac scale
77
91.188 = muR_ref_fixed ! fixed ren reference scale
78
91.188 = muF1_ref_fixed ! fixed fact reference scale for pdf1
79
91.188 = muF2_ref_fixed ! fixed fact reference scale for pdf2
75
%(fixed_ren_scale)s = fixed_ren_scale ! if .true. use fixed ren scale
76
%(fixed_fac_scale)s = fixed_fac_scale ! if .true. use fixed fac scale
77
%(mur_ref_fixed)s = muR_ref_fixed ! fixed ren reference scale
78
%(muf1_ref_fixed)s = muF1_ref_fixed ! fixed fact reference scale for pdf1
79
%(muf2_ref_fixed)s = muF2_ref_fixed ! fixed fact reference scale for pdf2
80
80
#***********************************************************************
81
81
# Renormalization and factorization scales (advanced and NLO options) *
82
82
#***********************************************************************
83
F = fixed_QES_scale ! if .true. use fixed Ellis-Sexton scale
84
91.188 = QES_ref_fixed ! fixed Ellis-Sexton reference scale
85
1 = muR_over_ref ! ratio of current muR over reference muR
86
1 = muF1_over_ref ! ratio of current muF1 over reference muF1
87
1 = muF2_over_ref ! ratio of current muF2 over reference muF2
88
1 = QES_over_ref ! ratio of current QES over reference QES
83
%(fixed_qes_scale)s = fixed_QES_scale ! if .true. use fixed Ellis-Sexton scale
84
%(qes_ref_fixed)s = QES_ref_fixed ! fixed Ellis-Sexton reference scale
85
%(mur_over_ref)s = muR_over_ref ! ratio of current muR over reference muR
86
%(muf1_over_ref)s = muF1_over_ref ! ratio of current muF1 over reference muF1
87
%(muf2_over_ref)s = muF2_over_ref ! ratio of current muF2 over reference muF2
88
%(qes_over_ref)s = QES_over_ref ! ratio of current QES over reference QES
89
89
#***********************************************************************
90
90
# Reweight flags to get scale dependence and PDF uncertainty *
91
91
# For scale dependence: factor rw_scale_up/down around central scale *
92
92
# For PDF uncertainty: use LHAPDF with supported set *
93
93
#***********************************************************************
94
.true. = reweight_scale ! reweight to get scale dependence
95
0.5 = rw_Rscale_down ! lower bound for ren scale variations
96
2.0 = rw_Rscale_up ! upper bound for ren scale variations
97
0.5 = rw_Fscale_down ! lower bound for fact scale variations
98
2.0 = rw_Fscale_up ! upper bound for fact scale variations
99
.false. = reweight_PDF ! reweight to get PDF uncertainty
100
244601 = PDF_set_min ! First of the error PDF sets
101
244700 = PDF_set_max ! Last of the error PDF sets
94
%(reweight_scale)s = reweight_scale ! reweight to get scale dependence
95
%(rw_rscale_down)s = rw_Rscale_down ! lower bound for ren scale variations
96
%(rw_rscale_up)s = rw_Rscale_up ! upper bound for ren scale variations
97
%(rw_fscale_down)s = rw_Fscale_down ! lower bound for fact scale variations
98
%(rw_fscale_up)s = rw_Fscale_up ! upper bound for fact scale variations
99
%(reweight_pdf)s = reweight_PDF ! reweight to get PDF uncertainty
100
%(pdf_set_min)s = PDF_set_min ! First of the error PDF sets
101
%(pdf_set_max)s = PDF_set_max ! Last of the error PDF sets
102
102
#***********************************************************************
103
103
# Merging - WARNING! Applies merging only at the hard-event level. *
104
104
# After showering an MLM-type merging should be applied as well. *
105
105
# See http://amcatnlo.cern.ch/FxFx_merging.htm for more details. *
106
106
#***********************************************************************
107
0 = ickkw ! 0 no merging, 3 FxFx merging, 4 UNLOPS
107
%(ickkw)s = ickkw ! 0 no merging, 3 FxFx merging, 4 UNLOPS
108
108
#***********************************************************************
110
110
#***********************************************************************
111
111
# BW cutoff (M+/-bwcutoff*Gamma) *
112
112
#***********************************************************************
113
%(bwcutoff)s = bwcutoff
114
114
#***********************************************************************
115
115
# Cuts on the jets *
116
116
# Jet clustering is performed by FastJet.
118
118
# considerably softer than the analysis cuts. *
119
119
# (more specific cuts can be specified in SubProcesses/cuts.f) *
120
120
#***********************************************************************
121
1 = jetalgo ! FastJet jet algorithm (1=kT, 0=C/A, -1=anti-kT)
122
0.7 = jetradius ! The radius parameter for the jet algorithm
123
10 = ptj ! Min jet transverse momentum
124
-1 = etaj ! Max jet abs(pseudo-rap) (a value .lt.0 means no cut)
121
%(jetalgo)s = jetalgo ! FastJet jet algorithm (1=kT, 0=C/A, -1=anti-kT)
122
%(jetradius)s = jetradius ! The radius parameter for the jet algorithm
123
%(ptj)s = ptj ! Min jet transverse momentum
124
%(etaj)s = etaj ! Max jet abs(pseudo-rap) (a value .lt.0 means no cut)
125
125
#***********************************************************************
126
126
# Cuts on the charged leptons (e+, e-, mu+, mu-, tau+ and tau-) *
127
127
# (more specific gen cuts can be specified in SubProcesses/cuts.f) *
128
128
#***********************************************************************
129
0 = ptl ! Min lepton transverse momentum
130
-1 = etal ! Max lepton abs(pseudo-rap) (a value .lt.0 means no cut)
131
0 = drll ! Min distance between opposite sign lepton pairs
132
0 = drll_sf ! Min distance between opp. sign same-flavor lepton pairs
133
0 = mll ! Min inv. mass of all opposite sign lepton pairs
134
30 = mll_sf ! Min inv. mass of all opp. sign same-flavor lepton pairs
129
%(ptl)s = ptl ! Min lepton transverse momentum
130
%(etal)s = etal ! Max lepton abs(pseudo-rap) (a value .lt.0 means no cut)
131
%(drll)s = drll ! Min distance between opposite sign lepton pairs
132
%(drll_sf)s = drll_sf ! Min distance between opp. sign same-flavor lepton pairs
133
%(mll)s = mll ! Min inv. mass of all opposite sign lepton pairs
134
%(mll_sf)s = mll_sf ! Min inv. mass of all opp. sign same-flavor lepton pairs
135
135
#***********************************************************************
136
136
# Photon-isolation cuts, according to hep-ph/9801442 *
137
137
# When ptgmin=0, all the other parameters are ignored *
138
138
#***********************************************************************
139
20 = ptgmin ! Min photon transverse momentum
140
-1 = etagamma ! Max photon abs(pseudo-rap)
141
0.4 = R0gamma ! Radius of isolation code
142
1.0 = xn ! n parameter of eq.(3.4) in hep-ph/9801442
143
1.0 = epsgamma ! epsilon_gamma parameter of eq.(3.4) in hep-ph/9801442
144
.true. = isoEM ! isolate photons from EM energy (photons and leptons)
139
%(ptgmin)s = ptgmin ! Min photon transverse momentum
140
%(etagamma)s = etagamma ! Max photon abs(pseudo-rap)
141
%(r0gamma)s = R0gamma ! Radius of isolation code
142
%(xn)s = xn ! n parameter of eq.(3.4) in hep-ph/9801442
143
%(epsgamma)s = epsgamma ! epsilon_gamma parameter of eq.(3.4) in hep-ph/9801442
144
%(isoem)s = isoEM ! isolate photons from EM energy (photons and leptons)
145
145
#***********************************************************************
146
146
# Maximal PDG code for quark to be considered a jet when applying cuts.*
147
147
# At least all massless quarks of the model should be included here. *
148
148
#***********************************************************************
149
%(maxjetflavor)s = maxjetflavor
150
150
#***********************************************************************
151
151
# For aMCfast+APPLGRID use in PDF fitting (http://amcfast.hepforge.org)*
152
152
#***********************************************************************
153
0 = iappl ! aMCfast switch (0=OFF, 1=prepare APPLgrids, 2=fill grids)
153
%(iappl)s = iappl ! aMCfast switch (0=OFF, 1=prepare APPLgrids, 2=fill grids)
154
154
#***********************************************************************