powheg_atlas_common.py

Go to the documentation of this file.

# Copyright (C) 2002-2021 CERN for the benefit of the ATLAS collaboration
 
"""! Common Standard Model parameters for Powheg.
 
@author James Robinson  <james.robinson@cern.ch>
"""
 

try:
    from EvgenProdTools.offline_dict import parameters
except ImportError:
    parameters = {
        'particles': {
            '11': {'mass': 0.00051, 'width': 0.0},
            '13': {'mass': 0.1057, 'width': 0.0},
            '1': {'mass': 0.32, 'width': 0.0},
            '2': {'mass': 0.32, 'width': 0.0},
            '3': {'mass': 0.5, 'width': 0.0},
            '4': {'mass': 1.55, 'width': 0.0},
            '5': {'mass': 4.95, 'width': 0.0},
            '6': {'mass': 172.5, 'width': 1.32},
            '24': {'mass': 80.399, 'width': 2.085},
            '23': {'mass': 91.1876, 'width': 2.4952},
            '25': {'mass': 125.0, 'width': 0.00407},
        },
        'EW_parameters': {
            'W_to_enu': 0.1082,
            'W_to_leptons': 3 * 0.1082,
            'W_to_hadrons': 1.0 - 3 * 0.1082,
            't_to_Wb': 1.0,
            't_to_Ws': 0.0,
            't_to_Wd': 0.0,
            'Vud': 0.97428,
            'Vus': 0.2253,
            'Vub': 0.00347,
            'Vcd': 0.2252,
            'Vcs': 0.97345,
            'Vcb': 0.041,
            'Vtd': 0.00862,
            'Vts': 0.0403,
            'Vtb': 0.999152,
            'alphaem_0': 0.00729735252,
            'alphaem': 0.00781653039,
            'alphaqcd': 0.1185,
            'G_F': 0.00001166397,
        }
    }
# from ..decorators import Singleton
 
class mass(object):
    """! Particle masses."""
 
    
    e = 0.00051
 
    
    mu = 0.1057
 
    
    tau = 1.777
 
    
    d = 0.32
 
    
    u = 0.32
 
    
    s = 0.5
 
    
    c = 1.55
 
    
    b = 4.95
 
    
    t = 172.5
 
    
    W = 80.399
 
    
    Z = 91.1876
 
    
    H = 125.0
 
class width(object):
    """! Particle widths."""
 
    
    t = 1.32
 
    
    W = 2.085
 
    
    Z = 2.4952
 
    
    H = 0.00407
 
class branching_ratio(object):
    """! Particle branching ratios."""
 
    
    W_to_enu = 0.1082
 
    
    W_to_leptons = 3 * W_to_enu
 
    
    W_to_hadrons = 1.0 - W_to_leptons
 
    
    t_to_Wb = 1.0
 
    
    t_to_Ws = 0.0
 
    
    t_to_Wd = 0.0
 
class CKM(object):
    """! CKM matrix."""
 
    
    Vud = 0.97428
 
    
    Vus = 0.2253
 
    
    Vub = 0.00347
 
    
    Vcd = 0.2252
 
    
    Vcs = 0.97345
 
    
    Vcb = 0.041
 
    
    Vtd = 0.00862
 
    
    Vts = 0.0403
 
    
    Vtb = 0.999152
 
class EW_parameters(object):
    """! EW parameters."""
 
    
    alphaem_0 = 0.00729735252  # 1/137.036
 
    
    alphaem = 0.00781653039  # 1/127.934
 
    
    alphaqcd = 0.1185
 
    
    G_F = 0.00001166397
 
    
    sin2thetaW = 0.23113
 
    
    sin2cabibbo = 0.051
 
class WS_parameters(object):
    """! Wolfenstein parameters."""
    """ would be used i"""
    """ NB: no recalculation is made at the moment if usePDGvalues() is used"""
 
    
    lam = 2.257000e-01
 
    
    A = 8.140000e-01
 
    
    rho = 1.350000e-01
 
    
    eta = 3.490000e-01
    
def usePDGvalues():
    
    mass.e = parameters['particles']['11']['mass']
 
    
    mass.mu = parameters['particles']['13']['mass']
 
    
    mass.tau = parameters['particles']['13']['mass']
 
    
    mass.d = parameters['particles']['1']['mass']
 
    
    mass.u = parameters['particles']['2']['mass']
 
    
    mass.s = parameters['particles']['3']['mass']
 
    
    mass.c = parameters['particles']['4']['mass']
 
    
    mass.b = parameters['particles']['5']['mass']
 
    
    mass.t = parameters['particles']['6']['mass']
 
    
    mass.W = parameters['particles']['24']['mass']
 
    
    mass.Z = parameters['particles']['23']['mass']
 
    
    mass.H = parameters['particles']['25']['mass']
 
    
    width.t = parameters['particles']['6']['width']
 
    
    width.W = parameters['particles']['24']['width']
 
    
    width.Z = parameters['particles']['23']['width']
 
    
    width.H = parameters['particles']['25']['width']
 
    
    branching_ratio.W_to_enu = parameters['EW_parameters']['W_to_enu']
 
    
    branching_ratio.W_to_leptons = parameters['EW_parameters']['W_to_leptons']
 
    
    branching_ratio.W_to_hadrons = parameters['EW_parameters']['W_to_hadrons']
 
    
    branching_ratio.t_to_Wb = parameters['EW_parameters']['t_to_Wb']
 
    
    branching_ratio.t_to_Ws = parameters['EW_parameters']['t_to_Wd']
 
    
    branching_ratio.t_to_Wd = parameters['EW_parameters']['t_to_Ws']
 
    
    CKM.Vud = parameters['EW_parameters']['Vud']
 
    
    CKM.Vus = parameters['EW_parameters']['Vus']
 
    
    CKM.Vub = parameters['EW_parameters']['Vub']
 
    
    CKM.Vcd = parameters['EW_parameters']['Vcd']
 
    
    CKM.Vcs = parameters['EW_parameters']['Vcs']
 
    
    CKM.Vcb = parameters['EW_parameters']['Vcb']
 
    
    CKM.Vtd = parameters['EW_parameters']['Vtd']
 
    
    CKM.Vts = parameters['EW_parameters']['Vts']
 
    
    CKM.Vtb = parameters['EW_parameters']['Vtb']
 
    
    EW_parameters.alphaem_0 = parameters['EW_parameters']['alphaem_0']
 
    
    EW_parameters.alphaem = parameters['EW_parameters']['alphaem']
 
    
    EW_parameters.alphaqcd = parameters['EW_parameters']['alphaqcd']
 
    
    EW_parameters.G_F = parameters['EW_parameters']['G_F']
 
    
    EW_parameters.sin2thetaW = parameters['EW_parameters'][('SIN2THETAW', 'sin2thetaW', 'Sin2ThetaW')]
 
    
    EW_parameters.sin2cabibbo = parameters['EW_parameters']['sin2cabibbo']