process = m.group(1)
isToWW = False
isInclusive = False
if 'toWW' in process :
process = process.replace('toWW', '')
isToWW = True
if 'Inclusive' in process:
process = process.replace('Inclusive', '')
isInclusive = True
mass = float(m.group(2))
xSection = 0.
try:
if process == 'VH':
xSection += yrparser13TeV.get(mass)['WH']['sigma']
xSection += yrparser13TeV.get(mass)['ZH']['sigma']
else:
xSection += yrparser13TeV.get(mass)[process]['sigma']
except KeyError:
print 'Higgs mass', mass, 'not found in cross section tables. Interpolating linearly at +- 1 GeV...'
if process=='VH':
xSection += 0.5 * (yrparser13TeV.get(mass-1.)['WH']['sigma'] + xSection + yrparser13TeV.get(mass+1.)['WH']['sigma'])
xSection += 0.5 * (yrparser13TeV.get(mass-1.)['ZH']['sigma'] + yrparser13TeV.get(mass+1.)['ZH']['sigma'])
else:
xSection += 0.5 * (yrparser13TeV.get(mass-1.)[process]['sigma'] + yrparser13TeV.get(mass+1.)[process]['sigma'])
if isToWW :
br = yrparser13TeV.get(mass)['H2B']['WW']
elif isInclusive:
br = 1.
process = m.group(1)
isToWW = False
isInclusive = False
if 'toWW' in process:
process = process.replace('toWW', '')
isToWW = True
if 'Inclusive' in process:
process = process.replace('Inclusive', '')
isInclusive = True
mass = float(m.group(2))
xSection = 0.
try:
if process == 'VH':
xSection += yrparser13TeV.get(mass)['WH']['sigma']
xSection += yrparser13TeV.get(mass)['ZH']['sigma']
else:
xSection += yrparser13TeV.get(mass)[process]['sigma']
except KeyError:
print 'Higgs mass', mass, 'not found in cross section tables. Interpolating linearly at +- 1 GeV...'
if process == 'VH':
xSection += 0.5 * (yrparser13TeV.get(mass - 1.)['WH']['sigma'] +
xSection +
yrparser13TeV.get(mass + 1.)['WH']['sigma'])
xSection += 0.5 * (yrparser13TeV.get(mass - 1.)['ZH']['sigma'] +
yrparser13TeV.get(mass + 1.)['ZH']['sigma'])
else:
xSection += 0.5 * (yrparser13TeV.get(mass - 1.)[process]['sigma'] +
yrparser13TeV.get(mass + 1.)[process]['sigma'])
process = m.group(1)
isToWW = False
isInclusive = False
if "toWW" in process:
process = process.replace("toWW", "")
isToWW = True
if "Inclusive" in process:
process = process.replace("Inclusive", "")
isInclusive = True
mass = float(m.group(2))
xSection = 0.0
try:
if process == "VH":
xSection += yrparser13TeV.get(mass)["WH"]["sigma"]
xSection += yrparser13TeV.get(mass)["ZH"]["sigma"]
else:
xSection += yrparser13TeV.get(mass)[process]["sigma"]
except KeyError:
print "Higgs mass", mass, "not found in cross section tables. Interpolating linearly at +- 1 GeV..."
if process == "VH":
xSection += 0.5 * (
yrparser13TeV.get(mass - 1.0)["WH"]["sigma"] + xSection + yrparser13TeV.get(mass + 1.0)["WH"]["sigma"]
)
xSection += 0.5 * (
yrparser13TeV.get(mass - 1.0)["ZH"]["sigma"] + yrparser13TeV.get(mass + 1.0)["ZH"]["sigma"]
)
else:
xSection += 0.5 * (
yrparser13TeV.get(mass - 1.0)[process]["sigma"] + yrparser13TeV.get(mass + 1.0)[process]["sigma"]