ATLAS Offline Software
W_SMEFT.py
Go to the documentation of this file.
1 # Copyright (C) 2002-2021 CERN for the benefit of the ATLAS collaboration
2 
3 from AthenaCommon import Logging
4 from ..powheg_V2 import PowhegV2
5 import itertools
6 
7 
8 logger = Logging.logging.getLogger("PowhegControl")
9 
10 
12  """!
13  Powheg configuration for W boson production in Standard Model Effective Field Theory
14 
15  This process is an implementation of the W+/- cross sections pp -> W -> l v by SMEFT operators
16 
17  It includes 196 SMEFT coefficients that the user can set. For symmetry, the authors included coefficients
18  for operators involving top-quarks. These should not be used as they are only relevant when initial-state
19  tops are considered, which is never the case.
20 
21  Publication (please cite if using this process): https://arxiv.org/abs/1804.07407 https://arxiv.org/abs/1703.04751
22  Brief manual (highly recommended reading if using this process): accessible on LXPlus etc. at $POWHEGPATH/POWHEG-BOX-V2/W_smeft/Docs/manual-BOX-W_smeft.pdf
23 
24  @author Timoth\'ee Theveneaux-Pelzer <tpelzer@cern.ch>
25  """
26 
27  def __init__(self, base_directory, **kwargs):
28  """!
29  Constructor: all process options are set here.
30 
31  @param base_directory: path to PowhegBox code.
32  @param kwargs dictionary of arguments from Generate_tf.
33  """
34  super(W_SMEFT, self).__init__(base_directory, "W_smeft", **kwargs)
35 
36  # Add parameter validation functions
37  self.validation_functions.append("validate_decays")
38 
39 
40  self.allowed_decay_modes = ["w- > e- ve~", "w- > mu- vm~", "w- > tau- vt~", "w+ > e+ ve", "w+ > mu+ vm", "w+ > tau+ vt"]
41 
42  # Add all keywords for this process, overriding defaults if required
43  self.add_keyword("alphaem")
44  self.add_keyword("bmass_lhe")
45  self.add_keyword("bornktmin")
46  self.add_keyword("bornonly")
47  self.add_keyword("bornsuppfact")
48  self.add_keyword("bornzerodamp")
49  self.add_keyword("bottomthr")
50  self.add_keyword("bottomthrpdf")
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  # clobberlhe: needed to make reweighting work --- otherwise Powheg crashes because it
58  # cannot overwrite the LHE file with events before reweighting.
59  # The value given to the parameter here doesn't matter, it just needs to be present in
60  # the powheg.input file. The reweighter postprocessor will replace the value with 1
61  # anyway.
62  self.add_keyword("clobberlhe")
63  self.add_keyword("cmass_lhe")
64  self.add_keyword("colltest")
65  self.add_keyword("compute_rwgt")
66  self.add_keyword("dim6", 1)
67  self.add_keyword("facscfact", self.default_scales[0])
68  self.add_keyword("fcnc", 0)
69  self.add_keyword("flg_debug")
70  self.add_keyword("foldcsi", 5)
71  self.add_keyword("foldphi", 2)
72  self.add_keyword("foldy", 2)
73  self.add_keyword("fullrwgt")
74  self.add_keyword("fullrwgtmode")
75  self.add_keyword("hdamp")
76  self.add_keyword("hfact")
77  self.add_keyword("icsimax")
78  self.add_keyword("idvecbos", hidden=True)
79  self.add_keyword("ih1")
80  self.add_keyword("ih2")
81  self.add_keyword("itmx1", 5)
82  self.add_keyword("itmx2", 2)
83  self.add_keyword("iupperfsr")
84  self.add_keyword("iupperisr")
85  self.add_keyword("iymax")
86  self.add_keyword("LambdaNP", 1000.0)
87  self.add_keyword("lhans1", self.default_PDFs)
88  self.add_keyword("lhans2", self.default_PDFs)
89  self.add_keyword("lhrwgt_descr")
90  self.add_keyword("lhrwgt_group_combine")
91  self.add_keyword("lhrwgt_group_name")
92  self.add_keyword("lhrwgt_id")
93  self.add_keyword("manyseeds")
94  self.add_keyword("mass_high", 2.0 * list(self.parameters_by_name("beam_energy"))[0].value)
95  self.add_keyword("mass_low", 60.0)
96  self.add_keyword("masswindow_high")
97  self.add_keyword("masswindow_low")
98  self.add_keyword("ncall1", 40000)
99  self.add_keyword("ncall2", 70000)
100  self.add_keyword("nubound", 10000)
101  self.add_keyword("par_2gsupp")
102  self.add_keyword("par_diexp")
103  self.add_keyword("par_dijexp")
104  self.add_keyword("parallelstage")
105  self.add_keyword("pdfreweight")
106  self.add_keyword("ptsqmin")
107  self.add_keyword("ptsupp")
108  self.add_keyword("radregion")
109  self.add_keyword("rand1")
110  self.add_keyword("rand2")
111  self.add_keyword("renscfact", self.default_scales[1])
112  self.add_keyword("running_width")
113  self.add_keyword("runningscale", 1) # if 1, central scale equal to Z/gamma virtuality. If 0, central scale equal to M_Z
114  self.add_keyword("rwl_add")
115  self.add_keyword("rwl_file")
116  self.add_keyword("rwl_format_rwgt")
117  self.add_keyword("rwl_group_events")
118  self.add_keyword("smartsig")
119  self.add_keyword("softtest")
120  self.add_keyword("sthw2")
121  self.add_keyword("storeinfo_rwgt")
122  self.add_keyword("testplots")
123  self.add_keyword("testsuda")
124  self.add_keyword("ubsigmadetails")
125  self.add_keyword("use-old-grid")
126  self.add_keyword("use-old-ubound")
127  self.add_keyword("vdecaymode", self.allowed_decay_modes[0], name="decay_mode")
128  self.add_keyword("withdamp")
129  self.add_keyword("withnegweights")
130  self.add_keyword("withsubtr")
131  self.add_keyword("Wmass")
132  self.add_keyword("Wwidth")
133  self.add_keyword("xupbound")
134  self.add_keyword("Zmass")
135  self.add_keyword("Zwidth")
136  # SMEFT operator coefficients
137  for coefficient_name in ["ReGEw", "ImGEw"]:
138  self.add_keyword(coefficient_name)
139  for coefficient in ["ReGUw_{flavours}", "ImGUw_{flavours}", "ReLeQu3_{flavours}", "ImLeQu3_{flavours}", "ReLeQu_{flavours}", "ImLeQu_{flavours}"]:
140  for flavours in ["{}{}".format(*x) for x in itertools.product(["u", "c", "t"], repeat=2)]:
141  coefficient_name = coefficient.format(flavours=flavours)
142  self.add_keyword(coefficient_name)
143  for coefficient in ["ReGDw_{flavours}", "ImGDw_{flavours}", "ReLedQ_{flavours}", "ImLedQ_{flavours}"]:
144  for flavours in ["{}{}".format(*x) for x in itertools.product(["d", "s", "b"], repeat=2)]:
145  coefficient_name = coefficient.format(flavours=flavours)
146  self.add_keyword(coefficient_name)
147  for coefficient in ["QphiU_{flavours}", "QLu_{flavours}"]:
148  for flavours in ["uu", "uc", "ut", "cc", "ct", "tt"]:
149  coefficient_name = coefficient.format(flavours=flavours)
150  self.add_keyword(coefficient_name)
151  for coefficient in ["QphiD_{flavours}", "QLd_{flavours}"]:
152  for flavours in ["dd", "ds", "db", "ss", "sb", "bb"]:
153  coefficient_name = coefficient.format(flavours=flavours)
154  self.add_keyword(coefficient_name)
155  for coefficient in ["ReXi_{flavours}", "ImXi_{flavours}"]:
156  for flavours in ["{}{}".format(*x) for x in itertools.product(["u", "c", "t"], ["d", "s", "b"])]:
157  coefficient_name = coefficient.format(flavours=flavours)
158  self.add_keyword(coefficient_name)
159 
160 
161 
162  def validate_decays(self):
163  """! Validate vdecaymode keyword."""
164  self.expose() # convenience call to simplify syntax
165  self.check_decay_mode(self.decay_mode, self.allowed_decay_modes)
166  # Calculate appropriate decay mode numbers
167  list(self.parameters_by_keyword("idvecbos"))[0].value = 24 * [-1, +1][self.decay_mode.startswith("w+")]
168  __decay_mode_lookup = {"e- ve~": 1, "mu- vm~": 2, "tau- vt~": 3, "e+ ve": 1, "mu+ vm": 2, "tau+ vt": 3}
169  list(self.parameters_by_keyword("vdecaymode"))[0].value = __decay_mode_lookup[self.decay_mode.split("> ")[1]]
python.processes.configurable.Configurable.expose
def expose(self)
Add all names to the interface of this object.
Definition: configurable.py:46
vtune_athena.format
format
Definition: vtune_athena.py:14
python.processes.configurable.Configurable.parameters_by_name
def parameters_by_name(self, name)
Retrieve all parameters that use a given name.
Definition: configurable.py:64
python.processes.powheg_base.PowhegBase.default_scales
def default_scales(self)
Default scale variations for this process.
Definition: powheg_base.py:246
python.processes.powheg.W_SMEFT.W_SMEFT.validate_decays
def validate_decays(self)
Validate vdecaymode keyword.
Definition: W_SMEFT.py:162
python.processes.configurable.Configurable.add_keyword
def 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.powheg.W_SMEFT.W_SMEFT.allowed_decay_modes
allowed_decay_modes
List of allowed decay modes.
Definition: W_SMEFT.py:40
dumpHVPathFromNtuple.append
bool append
Definition: dumpHVPathFromNtuple.py:91
python.processes.configurable.Configurable.parameters_by_keyword
def parameters_by_keyword(self, keyword)
Retrieve all parameters that use a given keyword.
Definition: configurable.py:57
python.processes.powheg_base.PowhegBase.validation_functions
validation_functions
List of validation functions to run before preparing runcard.
Definition: powheg_base.py:179
python.processes.powheg_base.PowhegBase.check_decay_mode
def 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:283
python.processes.powheg.W_SMEFT.W_SMEFT.__init__
def __init__(self, base_directory, **kwargs)
Constructor: all process options are set here.
Definition: W_SMEFT.py:27
histSizes.list
def list(name, path='/')
Definition: histSizes.py:38
python.processes.powheg_base.PowhegBase.default_PDFs
def default_PDFs(self)
Default PDFs for this process.
Definition: powheg_base.py:241
python.processes.powheg.W_SMEFT.W_SMEFT
Powheg configuration for W boson production in Standard Model Effective Field Theory.
Definition: W_SMEFT.py:11
python.processes.powheg_V2.PowhegV2
Base class for PowhegBox V2 processes.
Definition: powheg_V2.py:6
Trk::split
@ split
Definition: LayerMaterialProperties.h:38