from measurement.measurement import Measurement
from measurement.measurement import UpperLimitData
import magellan.config as magconf
import hepdata.hepdata as hepdata
limit_vs_mA_GammaA_path = os.path.join(magconf.paths['measurements'], 'ATLAS_EXOT_2016_10', './limit_vs_mA_Gamma_A.dat')
m = Measurement(
data = UpperLimitData.from_datafile(limit_vs_mA_GammaA_path, independent_variables=['mA','GammaA_div_mA']),
description = "ATLAS Higgs A -> Zh, sqrt(s) = 13 TeV, lumi = 36 fb-1",
tag = "ATLAS_EXOT_2016_10",
label = "ATLAS-EXOT-2016-10",
process = r"$A \rightarrow Zh$",
factor = 1.0,
observables = {
"xsec_sushi_ggh_A_NNLO_x_br_A_Zh":
{"upper_limit" : {"function" : None,
"arguments" : ["mA" , "GammaA_div_mA", "br_h_bb"]}}
},
CL = 0.95,
sqrt_s = 13.0,
luminosity = 36.0,
url = "https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PAPERS/EXOT-2016-10/",
comment = ""
)
def limit_vs_mA_GammaA_interpolation_function(vars, limit_type='limit_exp'):
interpolation_function = m.data['interpolation'][limit_type]['function']
if isinstance(vars, pd.DataFrame):
'ATLAS_HIGG_2016_12', './ggF_limit.yaml')
m = Measurement(
data=UpperLimitData.from_hepdatafile(hepdata_limit,
independent_variables=['mass']),
description="ATLAS Higgs H/A -> tautau, sqrt(s) = 13 TeV, lumi = 36 fb-1",
tag="ATLAS_HIGG_2016_12",
label="ATLAS-HIGG-2016-12",
process=r"$A/H \rightarrow \tau\tau$",
factor=1.0,
observables={
"xsec_sushi_ggh_A_NNLO_x_br_A_tautau": {
"upper_limit": {
"function": None,
"arguments": ["mA", "GammaA_div_mA"]
}
},
"xsec_sushi_ggh_H_NNLO_x_br_H_tautau": {
"upper_limit": {
"function": None,
"arguments": ["mH", "GammaH_div_mH"]
}
},
},
CL=0.95,
sqrt_s=13.0,
luminosity=36.0,
url="https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PAPERS/HIGG-2016-12/",
comment="")
from measurement.measurement import UpperLimitData
hepdata_limit = os.path.join(magconf.paths['measurements'],
'ATLAS_EXOT_2016_31', './scalar.yaml')
m = Measurement(
data=UpperLimitData.from_hepdatafile(hepdata_limit,
independent_variables=['mH']),
description=
"ATLAS Higgs H -> hh -> bbbb, sqrt(s) = 13 TeV, lumi = 24.3 fb-1",
tag="ATLAS_EXOT_2016_31",
label="ATLAS-EXOT-2016-31",
process=r"$H \rightarrow hh \rightarrow bbbb$",
factor=1.0,
observables={
"xsec_sushi_ggh_H_NNLO_x_br_H_hh_x_br_hh_bbbb": {
"upper_limit": {
"function": None,
"arguments": ["mH", "GammaH_div_mH"]
}
}
},
CL=0.95,
sqrt_s=13.0,
luminosity=24.3,
url="https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PAPERS/EXOT-2016-31/",
comment="")
def limit_vs_mass_func(vars, limit_type='limit_exp'):
interpolation_function = m.data['interpolation'][limit_type]['function']
