ATLAS Offline Software
Toggle main menu visibility
Loading...
Searching...
No Matches
Generators
PowhegControl
python
processes
powheg
HZj_EW.py
Go to the documentation of this file.
1
# Copyright (C) 2002-2021 CERN for the benefit of the ATLAS collaboration
2
3
from
...
import
Logging
4
from
..powheg_RES
import
PowhegRES
5
6
7
logger = Logging.logging.getLogger(
"PowhegControl"
)
8
9
10
class
HZj_EW
(
PowhegRES
):
11
"""! Default Powheg configuration for associated Higgs boson and Z-boson production plus one jet, with electroweak and quark mass effects.
12
13
Create a configurable object with all applicable Powheg options.
14
@author James Robinson <james.robinson@cern.ch>
15
"""
16
17
def
__init__
(self, base_directory, **kwargs):
18
"""! Constructor: all process options are set here.
19
20
@param base_directory path to PowhegBox code.
21
@param kwargs dictionary of arguments from Generate_tf.
22
"""
23
super(HZj_EW, self).
__init__
(base_directory,
"HZJ_ew"
, **kwargs)
24
25
26
self.
allowed_decay_modes
= [
"z > e+ e-"
,
"z > mu+ mu-"
,
"z > tau+ tau-"
,
"z > ve ve~"
,
"z > vm vm~"
,
27
"z > vt vt~"
,
"z > j j"
,
"z > l l~"
,
"z > v v~"
,
"z > all"
]
28
29
# Add all keywords for this process, overriding defaults if required
30
self.
add_keyword
(
"allrad"
, 1)
31
self.
add_keyword
(
"alphaem"
)
32
self.
add_keyword
(
"bmass"
)
33
self.
add_keyword
(
"bornktmin"
, 0.26)
34
self.
add_keyword
(
"bornonly"
)
35
self.
add_keyword
(
"bornsuppfact"
, 0.0)
36
self.
add_keyword
(
"bornsuppfactV"
, 700.0)
37
self.
add_keyword
(
"bornzerodamp"
)
38
self.
add_keyword
(
"bottomthr"
)
39
self.
add_keyword
(
"bottomthrpdf"
)
40
self.
add_keyword
(
"btildeborn"
)
41
self.
add_keyword
(
"btildecoll"
)
42
self.
add_keyword
(
"btildereal"
)
43
self.
add_keyword
(
"btildevirt"
)
44
self.
add_keyword
(
"btlscalect"
)
45
self.
add_keyword
(
"btlscalereal"
)
46
self.
add_keyword
(
"charmthr"
)
47
self.
add_keyword
(
"charmthrpdf"
)
48
self.
add_keyword
(
"check_bad_st1"
)
49
self.
add_keyword
(
"check_bad_st2"
)
50
self.
add_keyword
(
"chklimseed"
)
51
self.
add_keyword
(
"CKM_diagonal"
)
52
self.
add_keyword
(
"clobberlhe"
)
53
self.
add_keyword
(
"colltest"
)
54
self.
add_keyword
(
"comparison_paper"
)
55
self.
add_keyword
(
"compress_lhe"
)
56
self.
add_keyword
(
"compress_upb"
)
57
self.
add_keyword
(
"compute_rwgt"
)
58
self.
add_keyword
(
"doublefsr"
)
59
self.
add_keyword
(
"enhancereg"
)
60
self.
add_keyword
(
"evenmaxrat"
)
61
self.
add_keyword
(
"ew_renorm_scheme"
)
62
self.
add_keyword
(
"facscfact"
, self.
default_scales
[0])
63
self.
add_keyword
(
"factsc2min"
)
64
self.
add_keyword
(
"fastbtlbound"
)
65
self.
add_keyword
(
"flg_debug"
)
66
self.
add_keyword
(
"foldcsi"
, 5)
67
self.
add_keyword
(
"foldphi"
, 5)
68
self.
add_keyword
(
"foldy"
, 5)
69
self.
add_keyword
(
"for_reweighting"
)
70
self.
add_keyword
(
"frensc2min"
)
71
self.
add_keyword
(
"fullrwgt"
)
72
self.
add_keyword
(
"hdamp"
)
73
self.
add_keyword
(
"hdecaymode"
)
74
self.
add_keyword
(
"hfact"
)
75
self.
add_keyword
(
"hmass"
)
76
self.
add_keyword
(
"icsimax"
)
77
self.
add_keyword
(
"idvecbos"
, hidden=
True
)
78
self.
add_keyword
(
"ih1"
)
79
self.
add_keyword
(
"ih2"
)
80
self.
add_keyword
(
"itmx1"
, 2)
81
self.
add_keyword
(
"itmx1btl"
)
82
self.
add_keyword
(
"itmx1btlbrn"
)
83
self.
add_keyword
(
"itmx1reg"
)
84
self.
add_keyword
(
"itmx1rm"
)
85
self.
add_keyword
(
"itmx2"
, 2)
86
self.
add_keyword
(
"itmx2btl"
)
87
self.
add_keyword
(
"itmx2btlbrn"
)
88
self.
add_keyword
(
"itmx2reg"
)
89
self.
add_keyword
(
"itmx2rm"
)
90
self.
add_keyword
(
"iupperfsr"
)
91
self.
add_keyword
(
"iupperisr"
)
92
self.
add_keyword
(
"iymax"
)
93
self.
add_keyword
(
"kappa_ghb"
)
94
self.
add_keyword
(
"kappa_ght"
)
95
self.
add_keyword
(
"kappa_ghw"
)
96
self.
add_keyword
(
"lambdaHHH"
)
97
self.
add_keyword
(
"lhans1"
, self.
default_PDFs
)
98
self.
add_keyword
(
"lhans2"
, self.
default_PDFs
)
99
self.
add_keyword
(
"LOevents"
)
100
self.
add_keyword
(
"manyseeds"
)
101
self.
add_keyword
(
"max_h_mass"
, 2.0 * list(self.
parameters_by_name
(
"beam_energy"
))[0].value)
102
self.
add_keyword
(
"max_io_bufsize"
)
103
self.
add_keyword
(
"max_z_mass"
, 2.0 * list(self.
parameters_by_name
(
"beam_energy"
))[0].value)
104
self.
add_keyword
(
"maxseeds"
)
105
self.
add_keyword
(
"min_h_mass"
)
106
self.
add_keyword
(
"min_z_mass"
, 10.0)
107
self.
add_keyword
(
"minlo_nnll"
)
108
self.
add_keyword
(
"minlo"
, 1)
109
self.
add_keyword
(
"mint_density_map"
)
110
self.
add_keyword
(
"mintupbratlim"
, 1000)
111
self.
add_keyword
(
"ncall1"
, 120000)
112
self.
add_keyword
(
"ncall1btl"
)
113
self.
add_keyword
(
"ncall1btlbrn"
)
114
self.
add_keyword
(
"ncall1reg"
)
115
self.
add_keyword
(
"ncall1rm"
)
116
self.
add_keyword
(
"ncall2"
, 600000)
117
self.
add_keyword
(
"ncall2btl"
)
118
self.
add_keyword
(
"ncall2btlbrn"
)
119
self.
add_keyword
(
"ncall2reg"
)
120
self.
add_keyword
(
"ncall2rm"
)
121
self.
add_keyword
(
"no_photon_radiation"
)
122
self.
add_keyword
(
"noevents"
)
123
self.
add_keyword
(
"nohad"
)
124
self.
add_keyword
(
"nores"
)
125
self.
add_keyword
(
"novirtual"
)
126
self.
add_keyword
(
"nubound"
, 400000)
127
self.
add_keyword
(
"olpreset"
)
128
self.
add_keyword
(
"olverbose"
)
129
self.
add_keyword
(
"onshellhiggs"
)
130
self.
add_keyword
(
"par_2gsupp"
)
131
self.
add_keyword
(
"par_diexp"
)
132
self.
add_keyword
(
"par_dijexp"
)
133
self.
add_keyword
(
"parallelstage"
)
134
self.
add_keyword
(
"pdfreweight"
)
135
self.
add_keyword
(
"ptsupp"
)
136
self.
add_keyword
(
"ptVhigh"
)
137
self.
add_keyword
(
"ptVlow"
)
138
self.
add_keyword
(
"qed_qcd"
)
139
self.
add_keyword
(
"radregion"
)
140
self.
add_keyword
(
"raisingscales"
)
141
self.
add_keyword
(
"rand1"
)
142
self.
add_keyword
(
"rand2"
)
143
self.
add_keyword
(
"regridfix"
)
144
self.
add_keyword
(
"renscfact"
, self.
default_scales
[1])
145
self.
add_keyword
(
"runningscales"
, 3, description=
"scale at underlying Born level. [0: mh+mw; 1: Ht; 2: sqrt(pt_l1*pt_l2); 3: sqrt((p_H+p_l1+p_l2)^2)]"
)
146
self.
add_keyword
(
"rwl_add"
)
147
self.
add_keyword
(
"rwl_file"
)
148
self.
add_keyword
(
"rwl_format_rwgt"
)
149
self.
add_keyword
(
"rwl_group_events"
)
150
self.
add_keyword
(
"select_EW_virt"
)
151
self.
add_keyword
(
"shortlist_onlyb"
)
152
self.
add_keyword
(
"shortlist"
)
153
self.
add_keyword
(
"smartsig"
)
154
self.
add_keyword
(
"softmismch"
)
155
self.
add_keyword
(
"softonly"
)
156
self.
add_keyword
(
"softtest"
)
157
self.
add_keyword
(
"stage2init"
)
158
self.
add_keyword
(
"storeinfo_rwgt"
)
159
self.
add_keyword
(
"storemintupb"
)
160
self.
add_keyword
(
"sudscalevar"
)
161
self.
add_keyword
(
"testplots"
)
162
self.
add_keyword
(
"testsuda"
)
163
self.
add_keyword
(
"tmass"
)
164
self.
add_keyword
(
"ubexcess_correct"
)
165
self.
add_keyword
(
"ubsigmadetails"
, 0)
# causes crash if enabled
166
self.
add_keyword
(
"use-old-grid"
)
167
self.
add_keyword
(
"use-old-ubound"
)
168
self.
add_keyword
(
"user_reshists_sep"
)
169
self.
add_keyword
(
"vdecaymode"
, self.
allowed_decay_modes
[0], name=
"decay_mode"
, description=
"W boson decay mode. [{}]"
.format(
"; "
.join(self.
allowed_decay_modes
)))
170
self.
add_keyword
(
"virtonly"
)
171
self.
add_keyword
(
"Vstep"
)
172
self.
add_keyword
(
"whichpwhgevent"
)
173
self.
add_keyword
(
"withbtilde"
)
174
self.
add_keyword
(
"withdamp"
)
175
self.
add_keyword
(
"withnegweights"
)
176
self.
add_keyword
(
"withregular"
)
177
self.
add_keyword
(
"withremnants"
)
178
self.
add_keyword
(
"withsubtr"
)
179
self.
add_keyword
(
"xgriditeration"
)
180
self.
add_keyword
(
"xupbound"
, 2)
181
self.
add_keyword
(
"zerowidth"
)
182
183
# Add parameter validation functions
184
self.
validation_functions
.append(
"validate_decays"
)
185
186
def
validate_decays
(self):
187
"""! Validate idvecbos and vdecaymode keywords."""
188
self.
expose
()
# convenience call to simplify syntax
189
self.
check_decay_mode
(self.
decay_mode
, self.
allowed_decay_modes
)
190
# Calculate appropriate decay mode numbers
191
__decay_mode_lookup = {
"e+ e-"
: 1,
"mu+ mu-"
: 2,
"tau+ tau-"
: 3,
192
"ve ve~"
: 4,
"vm vm~"
: 5,
"vt vt~"
: 6,
193
"j j"
: 0,
"l l~"
: 11,
"v v~"
: 12,
"all"
: 10}
194
list(self.
parameters_by_keyword
(
"vdecaymode"
))[0].value = __decay_mode_lookup[self.
decay_mode
.
split
(
"> "
)[1]]
195
python.processes.configurable.Configurable.add_keyword
add_keyword(self, keyword, value=None, name=None, frozen=None, hidden=None, description=None, **kwargs)
Register configurable parameter that is exposed to the user.
Definition
configurable.py:21
python.processes.configurable.Configurable.expose
expose(self)
Add all names to the interface of this object.
Definition
configurable.py:46
python.processes.configurable.Configurable.parameters_by_name
parameters_by_name(self, name)
Retrieve all parameters that use a given name.
Definition
configurable.py:64
python.processes.configurable.Configurable.parameters_by_keyword
parameters_by_keyword(self, keyword)
Retrieve all parameters that use a given keyword.
Definition
configurable.py:57
python.processes.powheg.HZj_EW.HZj_EW
Default Powheg configuration for associated Higgs boson and Z-boson production plus one jet,...
Definition
HZj_EW.py:10
python.processes.powheg.HZj_EW.HZj_EW.decay_mode
decay_mode
Definition
HZj_EW.py:189
python.processes.powheg.HZj_EW.HZj_EW.__init__
__init__(self, base_directory, **kwargs)
Constructor: all process options are set here.
Definition
HZj_EW.py:17
python.processes.powheg.HZj_EW.HZj_EW.allowed_decay_modes
list allowed_decay_modes
List of allowed decay modes.
Definition
HZj_EW.py:26
python.processes.powheg.HZj_EW.HZj_EW.validate_decays
validate_decays(self)
Validate idvecbos and vdecaymode keywords.
Definition
HZj_EW.py:186
python.processes.powheg_RES.PowhegRES
Base class for PowhegBox RES processes.
Definition
powheg_RES.py:7
python.processes.powheg_base.PowhegBase.validation_functions
list validation_functions
List of validation functions to run before preparing runcard.
Definition
powheg_base.py:127
python.processes.powheg_base.PowhegBase.check_decay_mode
check_decay_mode(self, decay_mode, allowed_decay_modes=None)
Check whether a decay mode is allowed an raise an exception if it is not.
Definition
powheg_base.py:231
python.processes.powheg_base.PowhegBase.default_scales
default_scales(self)
Default scale variations for this process.
Definition
powheg_base.py:194
split
std::vector< std::string > split(const std::string &s, const std::string &t=":")
Definition
hcg.cxx:179
Generated on
for ATLAS Offline Software by
1.17.0