ATLAS Offline Software
Toggle main menu visibility
Loading...
Searching...
No Matches
Generators
PowhegControl
python
processes
powheg
bblvlv_modified.py
Go to the documentation of this file.
1
# Copyright (C) 2002-2024 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
bblvlv_modified
(
PowhegRES
):
11
"""! Default Powheg configuration for top pair and Wt production including non-resonant and interference effects.
12
This is the modified version that contains all dilepton decay channels in one sample.
13
Create a configurable object with all applicable Powheg options.
14
15
@author Andrej Saibel <Andrej.Saibel@cern.ch>
16
"""
17
18
def
__init__
(self, base_directory, **kwargs):
19
"""! Constructor: all process options are set here.
20
21
@param base_directory: path to PowhegBox code.
22
@param kwargs dictionary of arguments from Generate_tf.
23
"""
24
super(bblvlv_modified, self).
__init__
(base_directory,
"b_bbar_4l_modified"
, **kwargs)
25
26
# make sure some warnings end up being displayed as such
27
self.
warning_output
= [
"inverseMappingFSR: warning"
,
"inverseMappingISR: warning"
]
28
29
# Add parameter validation functions
30
self.
validation_functions
.append(
"validate_decays"
)
31
32
33
self.
allowed_decay_modes
= [
"b mu+ vmu b~ e- ve~"
,
"b e+ ve b~ mu- vmu~"
,
"b emu+ vemu b~ emu- vemu~"
,\
34
"b tau+ vtau b~ e- ve~"
,
"b e+ ve b~ tau- vtau~"
, \
35
"b mu+ vmu b~ tau- vtau~"
,
"b tau+ vtau b~ mu- vmu~"
,
"b l+ vl b~ l- vl~"
]
36
37
# Add all keywords for this process, overriding defaults if required
38
self.
add_keyword
(
"allrad"
, 1)
39
self.
add_keyword
(
"alpha"
)
40
self.
add_keyword
(
"bmass"
)
41
self.
add_keyword
(
"bornktmin"
)
42
self.
add_keyword
(
"bornonly"
)
43
self.
add_keyword
(
"bornsuppfact"
)
44
self.
add_keyword
(
"bornzerodamp"
)
45
self.
add_keyword
(
"bottomthr"
)
46
self.
add_keyword
(
"bottomthrpdf"
)
47
self.
add_keyword
(
"btildeborn"
)
48
self.
add_keyword
(
"btildecoll"
)
49
self.
add_keyword
(
"btildereal"
)
50
self.
add_keyword
(
"btildevirt"
)
51
self.
add_keyword
(
"btlscalect"
)
52
self.
add_keyword
(
"btlscalereal"
)
53
self.
add_keyword
(
"charmthr"
)
54
self.
add_keyword
(
"charmthrpdf"
)
55
self.
add_keyword
(
"check_bad_st1"
)
56
self.
add_keyword
(
"check_bad_st2"
)
57
self.
add_keyword
(
"chklimseed"
)
58
self.
add_keyword
(
"clobberlhe"
)
59
self.
add_keyword
(
"colltest"
)
60
self.
add_keyword
(
"complexGFermi"
)
61
self.
add_keyword
(
"compress_lhe"
)
62
self.
add_keyword
(
"compress_upb"
)
63
self.
add_keyword
(
"compute_rwgt"
)
64
self.
add_keyword
(
"doublefsr"
)
65
self.
add_keyword
(
"enhancereg"
)
66
self.
add_keyword
(
"evenmaxrat"
)
67
self.
add_keyword
(
"ewscheme"
)
68
self.
add_keyword
(
"facscfact"
, self.
default_scales
[0])
69
self.
add_keyword
(
"fastbtlbound"
)
70
self.
add_keyword
(
"fixedscale"
)
71
self.
add_keyword
(
"flg_debug"
)
72
self.
add_keyword
(
"foldcsi"
, 2)
73
self.
add_keyword
(
"foldphi"
, 5)
74
self.
add_keyword
(
"foldy"
, 5)
75
self.
add_keyword
(
"for_reweighting"
)
76
self.
add_keyword
(
"fullrwgt"
)
77
self.
add_keyword
(
"hdamp"
, 172.5)
78
self.
add_keyword
(
"hfact"
)
79
self.
add_keyword
(
"hmass"
)
80
self.
add_keyword
(
"hwidth"
)
81
self.
add_keyword
(
"icsimax"
)
82
self.
add_keyword
(
"ih1"
)
83
self.
add_keyword
(
"ih2"
)
84
self.
add_keyword
(
"itmx1"
)
85
self.
add_keyword
(
"itmx1btl"
)
86
self.
add_keyword
(
"itmx1btlbrn"
)
87
self.
add_keyword
(
"itmx1rm"
)
88
self.
add_keyword
(
"itmx2"
, 6)
89
self.
add_keyword
(
"itmx2btl"
)
90
self.
add_keyword
(
"itmx2btlbrn"
)
91
self.
add_keyword
(
"itmx2rm"
)
92
self.
add_keyword
(
"iupperfsr"
)
93
self.
add_keyword
(
"iupperisr"
)
94
self.
add_keyword
(
"iymax"
)
95
self.
add_keyword
(
"lhans1"
, self.
default_PDFs
)
96
self.
add_keyword
(
"lhans2"
, self.
default_PDFs
)
97
self.
add_keyword
(
"lhrwgt_descr"
)
98
self.
add_keyword
(
"lhrwgt_group_combine"
)
99
self.
add_keyword
(
"lhrwgt_group_name"
)
100
self.
add_keyword
(
"lhrwgt_id"
)
101
self.
add_keyword
(
"LOevents"
)
102
self.
add_keyword
(
"manyseeds"
)
103
self.
add_keyword
(
"max_io_bufsize"
)
104
self.
add_keyword
(
"maxseeds"
)
105
self.
add_keyword
(
"minlo"
)
106
self.
add_keyword
(
"mint_density_map"
)
107
self.
add_keyword
(
"mintupbratlim"
)
108
self.
add_keyword
(
"MSbarscheme"
)
109
self.
add_keyword
(
"ncall1"
, 120000)
110
self.
add_keyword
(
"ncall1btl"
)
111
self.
add_keyword
(
"ncall1btlbrn"
)
112
self.
add_keyword
(
"ncall1rm"
)
113
self.
add_keyword
(
"ncall2"
, 180000)
114
self.
add_keyword
(
"ncall2btl"
)
115
self.
add_keyword
(
"ncall2btlbrn"
)
116
self.
add_keyword
(
"ncall2rm"
)
117
self.
add_keyword
(
"ncallfrominput"
)
118
self.
add_keyword
(
"noevents"
)
119
self.
add_keyword
(
"nores"
)
120
self.
add_keyword
(
"novirtual"
)
121
self.
add_keyword
(
"nubound"
, 100000)
122
self.
add_keyword
(
"olpreset"
)
123
self.
add_keyword
(
"olverbose"
)
124
self.
add_keyword
(
"openloops-stability"
)
125
self.
add_keyword
(
"openloopsreal"
)
126
self.
add_keyword
(
"openloopsvirtual"
)
127
self.
add_keyword
(
"par_2gsupp"
)
128
self.
add_keyword
(
"par_diexp"
)
129
self.
add_keyword
(
"par_dijexp"
)
130
self.
add_keyword
(
"parallelstage"
)
131
self.
add_keyword
(
"pdfreweight"
)
132
self.
add_keyword
(
"ptsqmin"
)
133
self.
add_keyword
(
"ptsupp"
)
134
self.
add_keyword
(
"radregion"
)
135
self.
add_keyword
(
"rand1"
)
136
self.
add_keyword
(
"rand2"
)
137
self.
add_keyword
(
"regridfix"
)
138
self.
add_keyword
(
"renscfact"
, self.
default_scales
[1])
139
self.
add_keyword
(
"rwl_add"
)
140
self.
add_keyword
(
"rwl_file"
)
141
self.
add_keyword
(
"rwl_format_rwgt"
)
142
self.
add_keyword
(
"rwl_group_events"
)
143
self.
add_keyword
(
"smartsig"
)
144
self.
add_keyword
(
"softmismch"
)
145
self.
add_keyword
(
"softonly"
)
146
self.
add_keyword
(
"softtest"
)
147
self.
add_keyword
(
"stage2init"
)
148
self.
add_keyword
(
"storeinfo_rwgt"
)
149
self.
add_keyword
(
"storemintupb"
)
150
self.
add_keyword
(
"testplots"
)
151
self.
add_keyword
(
"testsuda"
)
152
self.
add_keyword
(
"tmass_phsp"
)
153
self.
add_keyword
(
"tmass"
)
154
self.
add_keyword
(
"twidth"
)
155
self.
add_keyword
(
"twidth_phsp"
)
156
self.
add_keyword
(
"ubexcess_correct"
)
157
self.
add_keyword
(
"ubsigmadetails"
, -1)
# disable cross-section output to avoid Fortran crash
158
self.
add_keyword
(
"use-old-grid"
)
159
self.
add_keyword
(
"use-old-ubound"
)
160
self.
add_keyword
(
"user_reshists_sep"
)
161
self.
add_keyword
(
"verytinypars"
)
162
self.
add_keyword
(
"virtonly"
)
163
self.
add_keyword
(
"whichpwhgevent"
)
164
self.
add_keyword
(
"withbtilde"
)
165
self.
add_keyword
(
"withdamp"
, 1)
166
self.
add_keyword
(
"withnegweights"
)
167
self.
add_keyword
(
"withremnants"
)
168
self.
add_keyword
(
"withsubtr"
)
169
self.
add_keyword
(
"wmass"
)
170
self.
add_keyword
(
"wwidth"
)
171
self.
add_keyword
(
"xgriditeration"
)
172
self.
add_keyword
(
"xupbound"
, 2)
173
self.
add_keyword
(
"zerowidth"
)
174
self.
add_keyword
(
"zmass"
)
175
self.
add_keyword
(
"zwidth"
)
176
self.
add_keyword
(
"channel"
, self.
allowed_decay_modes
[7], name=
"decay_mode"
)
177
178
def
validate_decays
(self):
179
"""! Validate decay_mode keyword."""
180
self.
expose
()
# convenience call to simplify syntax
181
if
self.
decay_mode
not
in
self.
allowed_decay_modes
:
182
logger.warning(
"Decay mode {} not recognised!"
.format(self.
decay_mode
))
183
raise
ValueError(
"Decay mode {} not recognised!"
.format(self.
decay_mode
))
184
# Calculate appropriate decay mode numbers
185
186
__decay_mode_lookup = {
"b mu+ vmu b~ e- ve~"
: 0,
"b e+ ve b~ mu- vmu~"
: 1,
"b emu+ vemu b~ emu- vemu~"
: 2,\
187
"b tau+ vtau b~ e- ve~"
: 3,
"b e+ ve b~ tau- vtau~"
: 4, \
188
"b mu+ vmu b~ tau- vtau~"
: 5,
"b tau+ vtau b~ mu- vmu~"
: 6,
"b l+ vl b~ l- vl~"
: 7}
189
self.
parameters_by_keyword
(
"channel"
)[0].value = __decay_mode_lookup[self.
decay_mode
]
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_keyword
parameters_by_keyword(self, keyword)
Retrieve all parameters that use a given keyword.
Definition
configurable.py:57
python.processes.powheg.bblvlv_modified.bblvlv_modified
Default Powheg configuration for top pair and Wt production including non-resonant and interference e...
Definition
bblvlv_modified.py:10
python.processes.powheg.bblvlv_modified.bblvlv_modified.allowed_decay_modes
list allowed_decay_modes
List of allowed decay modes.
Definition
bblvlv_modified.py:33
python.processes.powheg.bblvlv_modified.bblvlv_modified.__init__
__init__(self, base_directory, **kwargs)
Constructor: all process options are set here.
Definition
bblvlv_modified.py:18
python.processes.powheg.bblvlv_modified.bblvlv_modified.validate_decays
validate_decays(self)
Validate decay_mode keyword.
Definition
bblvlv_modified.py:178
python.processes.powheg.bblvlv_modified.bblvlv_modified.decay_mode
decay_mode
Definition
bblvlv_modified.py:182
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.default_scales
default_scales(self)
Default scale variations for this process.
Definition
powheg_base.py:194
Generated on
for ATLAS Offline Software by
1.17.0