~madteam/mg5amcnlo/series2.0

« back to all changes in this revision

Viewing changes to Template/NLO/Cards/run_card.dat

  • Committer: olivier Mattelaer
  • Date: 2015-03-05 00:14:16 UTC
  • mfrom: (258.1.9 2.3)
  • mto: (258.8.1 2.3)
  • mto: This revision was merged to the branch mainline in revision 259.
  • Revision ID: olivier.mattelaer@uclouvain.be-20150305001416-y9mzeykfzwnl9t0j
partial merge

Show diffs side-by-side

added added

removed removed

Lines of Context:
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. 
31
31
                 !  (-1= no split).
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
#***********************************************************************
109
109
#
110
110
#***********************************************************************
111
111
# BW cutoff (M+/-bwcutoff*Gamma)                                       *
112
112
#***********************************************************************
113
 
 15  = bwcutoff
 
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
 
 4 = maxjetflavor
 
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
#***********************************************************************