def __init__(self, name, exp, channel, Delini, interptype="log"):
self.exp = exp
self.channel = channel
self.delini = Delini
self.descr = "No description for this limit"
self.ref = "No reference for this limit"
self.combthr = 1.1 # Standard value used when relevant to combine the upper and lower limit
try:
self.name = name
self.mx_ini, self.lim_ini = uf.LoadData('LimData/' + name + '.dat',
interptype)
# if verbose: print("Loading: ", self.name,self.lim_ini )
except:
self.name = "NotDefined"
self.mx_ini = np.logspace(-3., 1., 30.)
self.lim_ini = np.zeros(np.shape(self.mx_ini))
try:
self.model = model
except:
self.model = "No model for this limit"
self.lim_inifull = {
"V": (self.mx_ini, self.lim_ini),
"AV": (self.mx_ini, self.lim_ini)
}
atlas_monoXlow.ref = "inspirehep.net/record/1635274"
atlas_monoXlow.UpdateLimIni(xi_ATLAS, Lim_ATLAS_Down)
UpdateLimitList(atlas_monoXlow)
atlas_monoXhigh = Limit("", "atlas", "monojet_up", 0.0)
atlas_monoXhigh.descr = """
Limits for ATLAS collaboration, 1711.03301 based on the upper recasted limit from 1807.03817
and our own recast.
"""
atlas_monoXhigh.ref = "inspirehep.net/record/1635274"
atlas_monoXhigh.UpdateLimIni(xi_ATLAS, Lim_ATLAS_Up)
UpdateLimitList(atlas_monoXhigh)
######################### MonoPhoton searches at DELPHI from 1103.0240 #########################
xi_LimLEP_V, LimLEP_V = uf.LoadData('LimData/lep_1103.0240_V.dat', "lin")
xi_LimLEP_AV, LimLEP_AV = uf.LoadData('LimData/lep_1103.0240_AV.dat', "lin")
lep_monoX = Limit("lep_1103.0240", "lep", "monogam", 0.0)
lep_monoX.descr = """
Limits for LEP collaboration, based on the recasting from 1103.0240 for the upper limit.
The lower limit is given by the breakdown of the EFT at the LEP CoM energy.
"""
lep_monoX.ref = "inspirehep.net/record/890992"
lep_monoX.UpdateLimIni(xi_LimLEP_V, LimLEP_V)
lep_monoX.UpdateLimIni(xi_LimLEP_AV, LimLEP_AV, "AV")
UpdateLimitList(lep_monoX)
######################### SN1987 cooling rate constraints #####
if verbose: print("SN1987 cooling ...")