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>
"""
 

from EvgenProdTools.offline_dict import parameters
# 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']