info.publication = 'Phys. Rev. D 100, 012006 (2019)'
SR = {'obsN' : [7, 1, 5, 7, 6],
'expN' : [8, 2.5, 4.6, 2.8, 5.7],
'bgErr' : [4, 1.3, 1.2, 1, 2.3],
'SR' : ['SRHad-Low', 'SRHad-High', 'SR1Lbb-High', 'SR1Lbb-Medium', 'SR1Lbb-Low']}
DATA = {'expExcl' : ['orig/Expectedlimit0lbb.csv', 'orig/Expectedlimit0lbb.csv', 'orig/Expectedlimit1lbb.csv', 'orig/Expectedlimit1lbb.csv', 'orig/Expectedlimit1lbb.csv'],
'obsExcl' : ['orig/Observedlimit0lbb.csv', 'orig/Observedlimit0lbb.csv', 'orig/Observedlimit1lbb.csv', 'orig/Observedlimit1lbb.csv', 'orig/Observedlimit1lbb.csv'],
'effMap' : ['orig/EffMap_TChiWH_SRHad-Low.txt', 'orig/EffMap_TChiWH_SRHad-High.txt', 'orig/EffMap_TChiWH_SR1Lbb-High.txt', 'orig/EffMap_TChiWH_SR1Lbb-Medium.txt', 'orig/EffMap_TChiWH_SR1Lbb-Low.txt'],
'fig' : ['12a.png', '12a.png', '12b.png', '12b.png', '12b.png']}
for i in range(len(SR['obsN'])):
#+++++++ dataset block ++++++++++++++
dataset = DataSetInput(SR['SR'][i])
dataset.setInfo(dataType = 'efficiencyMap', dataId = SR['SR'][i], observedN = SR['obsN'][i], expectedBG = SR['expN'][i], bgError = SR['bgErr'][i])
#+++++++ next txName block ++++++++++++++
newTx = dataset.addTxName('TChiWH')
newTx.checked = 'No'
newTx.constraint = "[[['W']],[['higgs']]]"
newTx.conditionDescription = None
newTx.condition = None
newTx.source = 'ATLAS'
#+++++++ next mass plane block ++++++++++++++
newPlane = newTx.addMassPlane(2*[[x, y]])
newPlane.figure = 'Fig.12a'
newPlane.figureUrl = 'https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PAPERS/SUSY-2017-01/fig_' + DATA['fig'][i]
newPlane.dataUrl = 'https://www.hepdata.net/record/ins1711261?version=1&table=Upper limit 0lbb'
newPlane.setSources(dataLabels = ['expExclusion', 'obsExclusion', 'efficiencyMap'],
dataFiles = [DATA['expExcl'][i], DATA['obsExcl'][i], DATA['effMap'][i]],
sys.path.append(os.path.abspath(utilsPath))
from smodels_utils.dataPreparation.inputObjects import MetaInfoInput, DataSetInput
from smodels_utils.dataPreparation.databaseCreation import databaseCreator
from smodels_utils.dataPreparation.massPlaneObjects import x, y, z
#+++++++ global info block ++++++++++++++
info = MetaInfoInput('CMS-PAS-SUS-16-022')
info.url = 'http://cms-results.web.cern.ch/cms-results/public-results/preliminary-results/SUS-16-022/'
info.sqrts = 13
info.lumi = 12.9
info.prettyName = '>= 3 leptons'
info.comment = 'https://cds.cern.ch/record/2205165'
#+++++++ dataset block ++++++++++++++
dataset = DataSetInput('data')
dataset.setInfo(dataType='upperLimit', dataId=None)
#+++++++ next txName block ++++++++++++++
T1tttt = dataset.addTxName('T1tttt')
T1tttt.checked = ''
T1tttt.constraint = "[[['t','t']],[['t','t']]]"
T1tttt.conditionDescription = None
T1tttt.condition = None
T1tttt.source = "CMS"
T1tttt.massConstraint = None
T1ttttoff = dataset.addTxName('T1ttttoff')
T1ttttoff.checked = ''
T1ttttoff.constraint = "[[['b','W','b','W']],[['b','W','b','W']]]"
T1ttttoff.conditionDescription = None
T1ttttoff.condition = None
info.arxiv = 'https://arxiv.org/abs/1502.02522'
info.contact ='Andre Lessa <*****@*****.**>; Jan Heisig <*****@*****.**>; SModelS'
info.publication ='https://cds.cern.ch/record/1987723/files/arXiv:1502.02522.pdf'
info.comment ='Search for long-lived charged particles implemented in arXiv:1509.00473. For the topologies with mixed MET-HSCP branches, the MET branch is irrelevant and wildcards are used.'
info.supersedes =''
#++Define list of datasets++
datasetNames = ['c000','c100','c200','c300']
observedNs = [42,7,0,0]
expectedBGs = [44.,5.6,0.56,0.02]
bgErrors = [9.,1.1,0.11,0.004]
obsUpperLimits = ['1.15*fb','0.441*fb','0.16*fb','0.159*fb']
expUpperLimits = ['1.23*fb','0.338*fb','0.16*fb','0.159*fb']
for i,name in enumerate(datasetNames):
#+++++++ dataset block ++++++++++++++
dataset = DataSetInput(name)
dataset.setInfo(dataType = 'efficiencyMap', dataId = name,
observedN=observedNs[i], expectedBG=expectedBGs[i], bgError=bgErrors[i],
upperLimit = obsUpperLimits[i], expectedUpperLimit = expUpperLimits[i])
#+++++++ txnames ++++++++++++++++++++
#+++++++ next txName block ++++++++++++++
HSCPM1 = dataset.addTxName('THSCPM1')
HSCPM1.checked =''
HSCPM1.constraint = "[[],[]]"
HSCPM1.condition =None
HSCPM1.finalState = ['HSCP','HSCP']
HSCPM1.massConstraints = None
HSCPM1.dataUrl = None
HSCPM1.source = 'SModelS'
#+++++++ next mass plane block ++++++++++++++
from smodels_utils.dataPreparation.inputObjects import MetaInfoInput, DataSetInput
from smodels_utils.dataPreparation.databaseCreation import databaseCreator
from smodels_utils.dataPreparation.massPlaneObjects import x, y, z
#+++++++ global info block ++++++++++++++
info = MetaInfoInput('CMS-PAS-EXO-16-036')
info.url = 'http://cms-results.web.cern.ch/cms-results/public-results/preliminary-results/EXO-16-036/index.html'
info.sqrts = 13
info.lumi = 12.9
info.prettyName = 'hscp search'
info.private = False
info.contact = '*****@*****.**'
info.comment = 'Upper limits digitized from Track+TOF png plots. Used conservative gluino bounds (50% gluino-gluon probability) and curve for the cloud model for the squark R-hadron constraints.'
#+++++++ dataset block ++++++++++++++
dataset = DataSetInput('data')
dataset.setInfo(dataType='upperLimit', dataId=None)
#+++++++ txnames ++++++++++++++++++++
#+++++++ next txName block ++++++++++++++
HSCPM1 = dataset.addTxName('THSCPM1b')
HSCPM1.checked = ''
HSCPM1.constraint = "[[],[]]"
HSCPM1.condition = None
HSCPM1.finalState = ['HSCP', 'HSCP']
HSCPM1.massConstraints = None
HSCPM1.dataUrl = 'http://cms-results.web.cern.ch/cms-results/public-results/preliminary-results/EXO-16-036/CMS-PAS-EXO-16-036_Figure_003-b.png'
HSCPM1.source = 'CMS'
#+++++++ next mass plane block ++++++++++++++
plane = HSCPM1.addMassPlane([[x], [x]])
plane.setSources(dataLabels=['upperLimits'],
from smodels_utils.dataPreparation.databaseCreation import databaseCreator
from smodels_utils.dataPreparation.massPlaneObjects import x, y, z
#+++++++ global info block ++++++++++++++
info = MetaInfoInput('ATLAS-SUSY-2016-27')
info.url = "https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PAPERS/SUSY-2016-27/"
info.sqrts = 13
info.lumi = 36.1
info.prettyName = " "
info.private = False
info.arxiv = 'https://arxiv.org/abs/1802.03158'
info.contact = 'ATLAS collaboration'
info.publication = 'Phys. Rev. D 97, 092006 (2018)'
#+++++++ dataset block ++++++++++++++
dataset = DataSetInput('data')
dataset.setInfo(dataType='upperLimit', dataId=None)
#+++++++ next txName block ++++++++++++++
T5Gamma = dataset.addTxName('T5Gamma')
T5Gamma.checked = 'No'
T5Gamma.constraint = "[[['jet','jet'],['photon']],[['jet','jet'],['photon']]]"
T5Gamma.conditionDescription = None
T5Gamma.condition = None
T5Gamma.source = "ATLAS"
#+++++++ next mass plane block ++++++++++++++
T5Gamma_1 = T5Gamma.addMassPlane(2 * [[x, y, 1.]])
T5Gamma_1.figure = 'Fig.8'
T5Gamma_1.figureUrl = 'https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PAPERS/SUSY-2016-27/fig_08.png'
T5Gamma_1.dataUrl = "https://www.hepdata.net/record/ins1654357?version=1&table=Cross section UL 1"
T5Gamma_1.setSources(
from smodels_utils.dataPreparation.massPlaneObjects import x, y, z
#+++++++ global info block ++++++++++++++
info = MetaInfoInput('CMS-SUS-16-035')
info.url = 'http://cms-results.web.cern.ch/cms-results/public-results/publications/SUS-16-035/index.html'
info.sqrts = 13
info.lumi = 35.9
info.prettyName = '2 SS leptons'
info.private = False
info.arxiv = 'https://arxiv.org/abs/1704.07323'
info.contact = '*****@*****.**'
info.publication = 'Eur.Phys.J. C77 (2017) no.9, 578, http://dx.doi.org/10.1140/epjc/s10052-017-5079-z'
info.comment = 'Moriond 2017. Seven SMS interpretations in paper, all implemented.'
#+++++++ dataset block ++++++++++++++
dataset = DataSetInput('data')
dataset.setInfo(dataType='upperLimit', dataId=None)
#+++++++txName block++++++++++++++++++++
T1tttt = dataset.addTxName('T1tttt')
T1tttt.checked = ''
T1tttt.constraint = "[[['t','t']],[['t','t']]]"
T1tttt.condition = None
T1tttt.conditionDescription = None
T1tttt.source = "CMS"
T1tttt.massConstraint = [['dm>=338.0'], ['dm>=338.0']]
T1ttttoff = dataset.addTxName('T1ttttoff')
T1ttttoff.checked = ''
T1ttttoff.constraint = "[[[b,W,b,W]],[[b,W,b,W]]]"
T1ttttoff.condition = None
#+++++++ global info block ++++++++++++++
info = MetaInfoInput('CMS-SUS-16-041')
info.url = 'http://cms-results.web.cern.ch/cms-results/public-results/publications/SUS-16-041/index.html'
info.sqrts = 13
info.lumi = 35.9
info.prettyName = 'Multileptons + jets + Etmiss'
info.private = False
info.arxiv = 'https://arxiv.org/abs/1710.09154'
info.contact = '*****@*****.**'
info.publication = 'Submitted to J. High Energy Phys.'
info.comment = 'Moriond 2017. Omitted the gluino -> q qbar W/Z decay simplified model since the branching ratio information for the gluino to the chargino or neutralino2 was not provided.'
#+++++++ dataset block ++++++++++++++
dataset = DataSetInput('data')
dataset.setInfo(dataType = 'upperLimit', dataId = None)
#+++++++txName block++++++++++++++++++++
T1tttt=dataset.addTxName('T1tttt')
T1tttt.checked=''
T1tttt.constraint="[[['t','t']],[['t','t']]]"
T1tttt.condition=None
T1tttt.conditionDescription = None
T1tttt.source="CMS"
T1tttt.massConstraint=[['dm>=338.0'],['dm>=338.0']]
T1ttttoff=dataset.addTxName('T1ttttoff')
T1ttttoff.checked=''
T1ttttoff.constraint="[[['b','W','b','W']],[['b','W','b','W']]]"
T1ttttoff.condition=None
from smodels_utils.dataPreparation.massPlaneObjects import x, y, z
#+++++++ global info block ++++++++++++++
info = MetaInfoInput('CMS-SUS-16-047')
info.url = 'http://cms-results.web.cern.ch/cms-results/public-results/publications/SUS-16-047/index.html'
info.sqrts = 13
info.lumi = 35.9
info.prettyName = 'Photon + HT'
info.private = False
info.arxiv = 'http://arxiv.org/abs/1707.06193'
info.contact = '*****@*****.**'
info.publication = 'to appear in J. High Energy Phys.'
info.comment = 'Moriond 2017. Several UL maps in the context of gauge-mediated SSB. Implemented T5gg and T6gg; others with mixed BRs are not useable. NB results only for high gluino mass.'
#+++++++ dataset block ++++++++++++++
dataset = DataSetInput('data')
dataset.setInfo(dataType='upperLimit', dataId=None)
#+++++txName block +++++++++++++++++
T5gg = dataset.addTxName('T5gg')
T5gg.checked = ''
T5gg.constraint = "[[['jet','jet'],['photon']],[['jet','jet'],['photon']]]"
T5gg.condition = None
T5gg.conditionDescription = None
T5gg.source = "CMS"
#++++++next mass plane block+++++++++
T5gg_1 = T5gg.addMassPlane(2 * [[x, y, 1.0]])
T5gg_1.figure = 'Fig. 6-c'
#+++++++ global info block ++++++++++++++
info = MetaInfoInput('ATLAS-SUSY-2016-24')
info.url = 'https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PAPERS/SUSY-2016-24/'
info.sqrts = 13
info.lumi = 36.1
info.prettyName = '2+ leptons (e,mu) + jets + Etmiss'
info.private = False
info.arxiv = 'https://arxiv.org/abs/1803.02762'
info.contact = 'ATLAS collaboration'
info.publication ='Eur. Phys. J. C 78 (2018) 995'
#+++++++ dataset block ++++++++++++++
dataset = DataSetInput('data')
dataset.setInfo(dataType = 'upperLimit', dataId = None)
#+++++++ next txName block ++++++++++++++
TChipChimSlepSlep = dataset.addTxName('TChipChimSlepSlep')
TChipChimSlepSlep.checked = 'No'
TChipChimSlepSlep.constraint = "[[['nu'],['e+']],[['nu'],['e-']]]+[[['nu'],['e-']],[['nu'],['e+']]]+[[['nu'],['mu+']],[['nu'],['mu-']]]+[[['nu'],['mu-']],[['nu'],['mu+']]]"
TChipChimSlepSlep.conditionDescription = None
TChipChimSlepSlep.condition = None
TChipChimSlepSlep.source = 'ATLAS'
#+++++++ next mass plane block ++++++++++++++
TChipChimSlepSlep_1 = TChipChimSlepSlep.addMassPlane(2*[[x, 0.5*(x+y), y]])
TChipChimSlepSlep_1.figure = 'Fig.8a'
TChipChimSlepSlep_1.figureUrl = 'https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PAPERS/SUSY-2016-24/fig_08a.png'
TChipChimSlepSlep_1.dataUrl = 'https://www.hepdata.net/record/ins1658902?version=1&table=Table78'
from smodels_utils.dataPreparation.databaseCreation import databaseCreator
from smodels_utils.dataPreparation.massPlaneObjects import x, y, z
#+++++++ global info block ++++++++++++++
info = MetaInfoInput('ATLAS-SUSY-2018-31')
info.url = 'https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PAPERS/SUSY-2018-31/'
info.sqrts = 13
info.lumi = 139.
info.prettyName = 'higgs + b-jets + MET'
info.private = False
info.arxiv = 'arXiv:1908.03122'
info.contact = '*****@*****.**'
info.publication = ''
# +++++++ dataset block ++++++++++++++
SRA_incl = DataSetInput('SRA_incl')
SRA_incl.setInfo(
dataType='efficiencyMap',
dataId="SRA_incl",
observedN=17,
expectedBG=17.1,
bgError=2.8,
# upperLimit = '2.637E-01*fb', expectedUpperLimit = '2.639E-01*fb', ## FIXME this line is nonsense
jsonfile="orig/BkgOnlyA.json")
# +++++++ next txName block ++++++++++++++
T6bbHH = SRA_incl.addTxName("T6bbHH")
T6bbHH.checked = 'NO'
T6bbHH.constraint = "[[['b'],['higgs']],[['b'],['higgs']]]"
T6bbHH.conditionDescription = None
T6bbHH.condition = None
from smodels_utils.dataPreparation.massPlaneObjects import x, y, z
#+++++++ global info block ++++++++++++++
info = MetaInfoInput('CMS-SUS-16-032')
info.url = 'http://cms-results.web.cern.ch/cms-results/public-results/publications/SUS-16-032/index.html'
info.sqrts = 13
info.lumi = 35.9
info.prettyName = 'Sbottom and compressed stop (jets + Etmiss)'
info.private = False
info.arxiv = 'https://arxiv.org/abs/1707.07274'
info.contact = '*****@*****.**'
info.publication = 'Submitted to Phys. Lett. B.'
info.comment = 'Moriond 2017. The name for upper limits in the root files are: hXsec_exp_corr;1 and hXsec_exp_corr;2 with both containing same data (checked by Federico). Here used hXsec_exp_corr;2 for implementation.'
#+++++++ dataset block ++++++++++++++
dataset = DataSetInput('data')
dataset.setInfo(dataType='upperLimit', dataId=None)
#+++++txName block +++++++++++++++++
T2bb = dataset.addTxName('T2bb')
T2bb.checked = ''
T2bb.constraint = "[[['b']],[['b']]]"
T2bb.condition = None
T2bb.conditionDescription = None
T2bb.source = "CMS"
T2bb.massConstraint = None
#++++++next mass plane block+++++++++
T2bb_1 = T2bb.addMassPlane(2 * [[x, y]])
#+++++++ global info block ++++++++++++++
info = MetaInfoInput('ATLAS-SUSY-2015-02')
info.url = 'http://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PAPERS/SUSY-2015-02/'
info.sqrts = 13
info.lumi = 3.2
info.prettyName = 'Top 1l'
info.private = False
info.arxiv = 'https://arxiv.org/abs/1606.03903'
info.source = "ATLAS"
info.publication = 'http://journals.aps.org/prd/abstract/10.1103/PhysRevD.94.052009'
info.comment = 'Very weird exclusion due to best CLs selection, and fluctuation in the Obs data'
#info.supersedes =
#info.supersededBy =
#+++++++ dataset block ++++++++++++++
dataset = DataSetInput("SR2")
dataset.setInfo(dataType='efficiencyMap',
dataId="SR2",
observedN=1,
expectedBG=1.25,
bgError=0.26)
T2tt = dataset.addTxName('T2tt')
T2tt.checked = 'NO'
T2tt.constraint = "[[['t']],[['t']]]"
T2tt.conditionDescription = None
T2tt.condition = None
T2tt.source = "ATLAS"
#+++++++ next mass plane block ++++++++++++++
T2tt_1 = T2tt.addMassPlane([[x, y]] * 2)
#---- new efficiency map -----
#----figure----
from smodels_utils.dataPreparation.databaseCreation import databaseCreator
from smodels_utils.dataPreparation.massPlaneObjects import x, y, z
#+++++++ global info block ++++++++++++++
info = MetaInfoInput('ATLAS-SUSY-2018-31')
info.url = 'https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PAPERS/SUSY-2018-31/'
info.sqrts = 13
info.lumi = 139.
info.prettyName = 'higgs + b-jets + MET'
info.private = False
info.arxiv = 'arXiv:1908.03122'
info.contact = 'ATLAS collaboration'
info.publication = ''
# +++++++ dataset block ++++++++++++++
SRA_L = DataSetInput('SRA_L')
SRA_L.setInfo(
dataType='efficiencyMap',
dataId="SRA_L",
observedN=12,
expectedBG=8.4,
bgError=1.7,
# upperLimit = '2.637E-01*fb', expectedUpperLimit = '2.639E-01*fb', ## FIXME this line is nonsense
jsonfile="orig/BkgOnlyA.json")
# +++++++ next txName block ++++++++++++++
T6bbHH = SRA_L.addTxName("T6bbHH")
T6bbHH.checked = 'NO'
T6bbHH.constraint = "[[['b'],['higgs']],[['b'],['higgs']]]"
T6bbHH.conditionDescription = None
T6bbHH.condition = None
info.lumi = 36.1
info.prettyName = " "
info.private = False
info.arxiv = 'https://arxiv.org/abs/1802.03158'
info.contact = 'ATLAS collaboration'
info.publication = 'Phys. Rev. D 97, 092006 (2018)'
obsN = [0, 0]
expN = [0.5, 0.48]
bgErr = [0.28, 0.275]
SR = ['SRyy-SL', 'SRyy-SH']
for i in range(len(obsN)):
#+++++++ dataset block ++++++++++++++
dataset = DataSetInput(SR[i] + '-q')
dataset.setInfo(dataType = 'efficiencyMap', dataId = SR[i] + '-q', observedN = obsN[i], expectedBG = expN[i], bgError = bgErr[i])
#+++++++ next txName block ++++++++++++++
T6Gamma = dataset.addTxName('T6Gamma')
T6Gamma.checked = 'No'
T6Gamma.constraint = "[[['jet'],['photon']],[['jet'],['photon']]]"
T6Gamma.conditionDescription = None
T6Gamma.condition = None
T6Gamma.source = "ATLAS"
#+++++++ next mass plane block ++++++++++++++
T6Gamma_1 = T6Gamma.addMassPlane(2*[[x, y, 0.]])
T6Gamma_1.figure = 'Fig.9'
T6Gamma_1.figureUrl = 'https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PAPERS/SUSY-2016-27/fig_09.png'
T6Gamma_1.dataUrl = "https://www.hepdata.net/record/ins1654357?version=1&table=Cross section UL 2"
T6Gamma_1.setSources(dataLabels = ['expExclusion', 'obsExclusion', 'efficiencyMap'],
dataFiles = ['orig/HEPData-ins1654357-v1-Exclusion_contour_(expected)_3.csv', 'orig/HEPData-ins1654357-v1-Exclusion_contour_(observed)_4.csv', 'orig/EffMap_T6Gamma_' + SR[i] + '.txt'],
#+++++++ global info block ++++++++++++++
info = MetaInfoInput('ATLAS-SUSY-2016-07')
info.url = "https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PAPERS/SUSY-2016-07/"
info.sqrts = 13
info.lumi = 36.1
info.prettyName = "hadronic jets + MET"
info.private = False
info.arxiv = 'https://arxiv.org/abs/1712.02332'
info.contact = 'ATLAS collaboration'
info.publication = 'Phys. Rev. D 97, 112001 (2018)'
#+++++++ dataset block ++++++++++++++
dataset = DataSetInput('data')
dataset.setInfo(dataType = 'upperLimit', dataId = None)
topos = [ "T1", "T2", "T6WWleft", "T5WW", "T5ZZ", "T5WZh", "T6WZh", "T5WWoff", "T6WWoffleft" ]
# topos = [ "T5WZh", "T6WZh" ]
constraints = { "T2": "[[[q]],[[q]]]", "T1": "[[[q,q]],[[q,q]]]", "T6WWleft": "[[['q'],['W']],[['q'],['W']]]", "T6WWoffleft": "2.23 * [[['jet'],['jet','jet']],[['jet'],['jet','jet']]]", "T5WW": "[[['q','q'],['W']],[['q','q'],['W']]]", "T5ZZ": "[[['q','q'],['Z']],[['q','q'],['Z']]]", "T5WZh": "[[[q,q],[Z]],[[q,q],[W]]]+[[[q,q],[higgs]],[[q,q],[W]]]+[[[q,q],[W]],[[q,q],[W]]]+[[[q,q],[Z]],[[q,q],[Z]]]+[[[q,q],[higgs]],[[q,q],[higgs]]]+[[[q,q],[higgs]],[[q,q],[Z]]]", "T6WZh": "[[[jet],[Z]],[[jet],[W]]]+[[[jet],[higgs]],[[jet],[W]]]+[[[jet],[W]],[[jet],[W]]]+[[[jet],[Z]],[[jet],[Z]]]+[[[jet],[higgs]],[[jet],[higgs]]]+[[[jet],[Z]],[[jet],[higgs]]]", "T5WWoff": "2.23*[[['jet','jet'],['jet','jet']],[['jet','jet'],['jet','jet']]]" }
conditions = { "T5WZh": "Csim([[[q,q],[Z]],[[q,q],[W]]]+[[[q,q],[higgs]],[[q,q],[W]]],2.*([[[q,q],[W]],[[q,q],[W]]]))", "T6WZh": "Csim([[[jet],[Z]],[[jet],[W]]]+[[[jet],[higgs]],[[jet],[W]]],2.*([[[jet],[W]],[[jet],[W]]]))" }
figure = { "T2": "Fig. 76a", "T1": "Fig. 76b", "T6WWleft": ('Fig. 77a', 'Fig. 77b' ), "T5WW": ( 'Fig. 77c', 'Fig. 77d' ), "T5ZZ": ('Fig.79', None) , "T5WZh": ('Fig. 78b', None ), "T6WZh": ('Fig. 78a', None ), "T5WWoff": ( 'Fig. 77c', 'Fig. 77d' ), "T6WWoffleft": ('Fig. 77a', 'Fig. 77b' ) }
baseUrl = 'https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PAPERS/SUSY-2016-07/'
figureUrl = { "T2": 'figaux_076a.png', "T1": 'figaux_076b.png',
"T6WWleft": ( 'figaux_077a.png', 'figaux_077b.png' ),
"T5WW": ( 'figaux_077c.png', 'figaux_077d.png' ),
from smodels_utils.dataPreparation.massPlaneObjects import x, y, z
#+++++++ global info block ++++++++++++++
info = MetaInfoInput('CMS-SUS-16-036')
info.url = 'http://cms-results.web.cern.ch/cms-results/public-results/publications/SUS-16-036/index.html'
info.sqrts = 13
info.lumi = 35.9
info.prettyName = '0L + jets + Etmiss (using MT2)'
info.private = False
info.arxiv = 'https://arxiv.org/abs/1705.04650'
info.contact = '*****@*****.**'
info.publication = 'Eur. Phys. J. C 77 (2017) 710, http://dx.doi.org/10.1140/epjc/s10052-017-5267-x'
info.comment = 'Moriond 2017. The mixed decay topology for stop pair production with stop->t \chi1_0 (50% BR) and stop->b \chi^{\pm} -> b W^{\pm} \chi1_0 (50%), has not been implemented since the branching ratio is fixed and the mass difference between the chargino and neutralino is fixed to 5 GeV. '
#+++++++ dataset block ++++++++++++++
dataset = DataSetInput('data')
dataset.setInfo(dataType='upperLimit', dataId=None)
#+++++txName block +++++++++++++++++
T1bbbb = dataset.addTxName('T1bbbb')
T1bbbb.checked = ''
T1bbbb.constraint = "[[['b','b']],[['b','b']]]"
T1bbbb.condition = None
T1bbbb.conditionDescription = None
T1bbbb.source = "CMS"
#++++++next mass plane block+++++++++
T1bbbb_1 = T1bbbb.addMassPlane(2 * [[x, y]])
T1bbbb_1.figure = 'Fig. 6-a'
#+++++++ global info block ++++++++++++++
info = MetaInfoInput('ATLAS-CONF-2012-166')
info.comment = 'superseding publication contain more Data'
info.sqrts = '8.0'
info.private = False
info.lumi = '13.0'
info.url = 'https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/CONFNOTES/ATLAS-CONF-2012-166/'
info.supersededBy = 'ATLAS-SUSY-2013-15'
info.prettyName = '1 lepton + 4(1 b-)jets + Etmiss'
info.implementedBy = 'MT'
#+++++++ dataset block ++++++++++++++
dataset = DataSetInput('data')
dataset.setInfo(dataType = 'upperLimit', dataId = None)
#+++++++ next txName block ++++++++++++++
T2tt = dataset.addTxName('T2tt')
T2tt.checked ="AL"
T2tt.constraint ="[[['t']],[['t']]]"
T2tt.conditionDescription ="None"
T2tt.condition ="None"
T2tt.source = "ATLAS"
#+++++++ next mass plane block ++++++++++++++
T2tt = T2tt.addMassPlane(2*[[x, y]])
T2tt.figure = 'Fig.(aux) 3'
T2tt.figureUrl = 'https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/CONFNOTES/ATLAS-CONF-2012-166//figaux_03.png'
T2tt.dataUrl = 'Not defined'
T2tt.setSources(dataLabels= ['obsExclusion', 'upperLimits'],
'expExcl': 'orig/HEPData-ins1639856-v4-Table_28.csv',
'obsExcl': 'orig/HEPData-ins1639856-v4-Table_29.csv',
'upLimit': 'orig/HEPData-ins1639856-v4-Table_70.csv'
},
'constraint': [
'[[[b],[W]],[[b],[W]]]',
'[[[b],[L,nu]],[[b],[L,nu]]]+[[[b],[L,nu]],[[b],[q,q]]]+[[[b],[q,q]],[[b],[q,q]]]'
],
'massPlane':
2 * [[x, x - 10, y]]
}
DATA = [T2tt]
#+++++++ dataset block ++++++++++++++
dataset = DataSetInput('data')
dataset.setInfo(dataType='upperLimit', dataId=None)
for TX in DATA:
#+++++++ next txName block ++++++++++++++
newTx = dataset.addTxName(TX['name'][0])
newTx.checked = 'False'
newTx.constraint = TX['constraint'][0]
newTx.conditionDescription = None
newTx.condition = None
newTx.source = 'ATLAS'
#+++++++ next txName block ++++++++++++++
newTxOff1 = dataset.addTxName(TX['name'][1])
newTxOff1.checked = 'False'
newTxOff1.constraint = TX['constraint'][1]
newTxOff1.conditionDescription = None
observedNs = [5, 1, 0, 0] #Extracted from a fit of the data
expectedBGs = [2.63, 0.3377, 0.0258,
0.0045] #Extracted from a log fit of CMS BG
bgErrors = [
0.53, 0.127, 0.0036, 0.001
] #Estimated from the error between a log fit and a linear fit of the CMS BG
#Interpolated SR Upper limits from the PAS (Tables 5 and 6). Note that the SR UL are given by obs_UL*eff (from the table)
obsUpperLimits = ['0.624*fb', '0.306*fb', '0.24*fb', '0.24*fb']
#Interpolated expected SR Upper limits from the PAS (Tables 5 and 6). Note that the SR UL are given by exp_UL*eff (from the table)
expUpperLimits = [
'0.436*fb', '0.242*fb', '0.241*fb', '0.240*fb'
] #I've artificially increase the expUL a bit, so the SR with the largest Mreco is selected
for i, name in enumerate(datasetNames):
#+++++++ dataset block ++++++++++++++
dataset = DataSetInput(name)
dataset.setInfo(dataType='efficiencyMap',
dataId=name,
observedN=observedNs[i],
expectedBG=expectedBGs[i],
bgError=bgErrors[i],
upperLimit=obsUpperLimits[i],
expectedUpperLimit=expUpperLimits[i])
#+++++++ txnames ++++++++++++++++++++
#+++++++ next txName block ++++++++++++++
HSCPM1 = dataset.addTxName('THSCPM1b')
HSCPM1.checked = ''
HSCPM1.constraint = "[[],[]]"
HSCPM1.condition = None
HSCPM1.finalState = ['HSCP', 'HSCP']
'constraint': '[[[b, L, nu]],[[b, jet,jet]]]',
'massConstr': [['dm < 80'], ['dm < 80']],
'massPlanes': [2 * [[x, y]], 2 * [[x, x - y]]]
}
Txnames = [T2tt, T2ttoff, T2bbll]
# Txnames = [T2tt,T2ttoff]
exclusionCurves = {
'obs': ['orig/Table17.csv', 'orig/Table20.csv'],
'exp': ['orig/Table16.csv', 'orig/Table19.csv']
}
for i, sr in enumerate(SR['SR']):
#+++++++ dataset block ++++++++++++++
dataset = DataSetInput(sr)
dataset.setInfo(dataType='efficiencyMap',
dataId=sr,
observedN=SR['obsN'][i],
expectedBG=SR['expN'][i],
bgError=SR['bgErr'][i],
upperLimit=str(SR['NsigUL'][i] / 36.1).strip() + '*fb',
expectedUpperLimit=str(SR['NsigULexp'][i] / 36.1).strip() +
'*fb')
#+++++++ next txName block ++++++++++++++
for TX in Txnames:
#The SRs 'bffN' and 'tN_diag_low' only apply to T2tt and T2bbll
if sr in ['bffN', 'tN_diag_low'] and TX['name'] != 'T2bbll':
continue
newTx = dataset.addTxName(TX['name'])
newTx.checked = 'False'
"T5ZZ": (2 * [[x, y, 1.]], None),
"T5WW": (2 * [[x, 0.5 * (x + y), y]], 2 * [[x, y, 60.]]),
"T5WWoff": (2 * [[x, 0.5 * (x + y), y]], 2 * [[x, y, 60.]]),
"TGQ": ([[x, 0.], [y, 0.]], [[x, 695.], [y, 695.]], [[x, 995.], [y, 995.]])
}
datasets = {}
dsnames = getSignalRegionsEMBaked("orig/%s.embaked" % "T2")
for dsname in dsnames:
#+++++++ dataset block ++++++++++++++
if not dsname in stats:
print("cannot find stats for %s" % dsname)
continue
dataset = DataSetInput(dsname)
dst = stats[dsname]
dataset.setInfo(dataType='efficiencyMap',
dataId=dsname,
observedN=dst["nobs"],
expectedBG=dst["nb"],
bgError=dst["deltanb"])
datasets[dsname] = dataset
rjrs = [1, 2, 3, 5]
for topo in topos:
for dsname in dsnames:
#+++++++ next txName block ++++++++++++++
Tx = datasets[dsname].addTxName(topo)
sys.path.append(os.path.abspath(utilsPath))
from smodels_utils.dataPreparation.inputObjects import MetaInfoInput, DataSetInput
from smodels_utils.dataPreparation.databaseCreation import databaseCreator
from smodels_utils.dataPreparation.massPlaneObjects import x, y, z
#+++++++ global info block ++++++++++++++
info = MetaInfoInput('CMS-PAS-SUS-16-024')
info.url = 'http://cms-results.web.cern.ch/cms-results/public-results/preliminary-results/SUS-16-024/'
info.sqrts = 13
info.lumi = 12.9
info.private = False
info.comment = 'https://cds.cern.ch/record/2205168. Superseded by CMS-SUS-16-039.'
info.supersededBy = 'CMS-SUS-16-039'
#+++++++ dataset block ++++++++++++++
dataset = DataSetInput('data')
dataset.setInfo(dataType='upperLimit', dataId=None)
#+++++++ next txName block ++++++++++++++
TChiChipmStauL = dataset.addTxName('TChiChipmStauL')
TChiChipmStauL.checked = ''
TChiChipmStauL.constraint = "[[['ta'],['ta']],[['nu'],['L']]]"
TChiChipmStauL.conditionDescription = None
TChiChipmStauL.condition = None
TChiChipmStauL.source = "CMS"
#+++++++ next mass plane block ++++++++++++++
TChiChipmStauL_1 = TChiChipmStauL.addMassPlane(2 * [[x, .5 * (x + y), y]])
TChiChipmStauL_1.figure = 'Figure 9'
TChiChipmStauL_1.figureUrl = 'http://cms-results.web.cern.ch/cms-results/public-results/preliminary-results/SUS-16-024/CMS-PAS-SUS-16-024_Figure_009.png'
TChiChipmStauL_1.dataUrl = 'http://cms-results.web.cern.ch/cms-results/public-results/preliminary-results/SUS-16-024/CMS-PAS-SUS-16-024_Figure_009.root'
TChiChipmStauL_1.histoDataUrl = 'http://cms-results.web.cern.ch/cms-results/public-results/preliminary-results/SUS-16-024/CMS-PAS-SUS-16-024_Figure_009.root'
#+++++++ global info block ++++++++++++++
info = MetaInfoInput('ATLAS-SUSY-2015-01')
info.url = "https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PAPERS/SUSY-2015-01/"
info.sqrts = 13
info.lumi = 3.2
info.prettyName = "2b"
info.private = False
info.arxiv = 'https://arxiv.org/abs/1606.08772v2'
info.contact = 'ATLAS collaboration'
info.publication = 'http://link.springer.com/article/10.1140/epjc/s10052-016-4382-4'
#info.comment =
#info.supersedes =
#info.supersededBy =
#+++++++ dataset block ++++++++++++++
dataset = DataSetInput("SRB")
dataset.setInfo(dataType='efficiencyMap',
dataId="SRB",
observedN=6,
expectedBG=12,
bgError=2.5)
T2bb = dataset.addTxName('T2bb')
T2bb.checked = 'NO'
T2bb.constraint = "[[['b']],[['b']]]"
T2bb.conditionDescription = None
T2bb.condition = None
T2bb.source = "ATLAS"
#+++++++ next mass plane block ++++++++++++++
T2bb_1 = T2bb.addMassPlane([[x, y]] * 2)
#---- new efficiency map -----
#----figure----
from smodels_utils.dataPreparation.databaseCreation import databaseCreator
from smodels_utils.dataPreparation.massPlaneObjects import x, y, z
#+++++++ global info block ++++++++++++++
info = MetaInfoInput('ATLAS-SUSY-2016-19')
info.url = "https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PAPERS/SUSY-2016-19/"
info.sqrts = 13
info.lumi = 36.1
info.prettyName = "stops to staus"
info.private = False
info.arxiv = 'https://arxiv.org/abs/1803.10178'
info.contact = 'ATLAS collaboration'
info.publication = 'in process of being published'
#+++++++ dataset block ++++++++++++++
dataset = DataSetInput('data')
dataset.setInfo(dataType='upperLimit', dataId=None)
#+++++++ next txName block ++++++++++++++
T4bnutaubnutau = dataset.addTxName('T4bnutaubnutau')
T4bnutaubnutau.checked = ''
T4bnutaubnutau.constraint = "[[[b,nu],[ta]],[[b,nu],[ta]]]"
T4bnutaubnutau.conditionDescription = None
T4bnutaubnutau.condition = None
T4bnutaubnutau.source = "ATLAS"
#+++++++ next mass plane block ++++++++++++++
T4bnutaubnutau_1 = T4bnutaubnutau.addMassPlane(2 * [[x, y, 0.]])
T4bnutaubnutau_1.figure = 'Fig.7'
T4bnutaubnutau_1.figureUrl = 'https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PAPERS/SUSY-2016-19/fig_07.png'
T4bnutaubnutau_1.setSources(
dataLabels=['expExclusion', 'obsExclusion', 'upperLimits'],
from smodels_utils.dataPreparation.massPlaneObjects import x, y, z
#+++++++ global info block ++++++++++++++
info = MetaInfoInput('CMS-SUS-16-043')
info.url = 'http://cms-results.web.cern.ch/cms-results/public-results/publications/SUS-16-043/index.html'
info.sqrts = 13
info.lumi = 35.9
info.prettyName = 'EWK WH(bb)'
info.private = False
info.arxiv = 'http://arxiv.org/abs/1706.09933'
info.contact = '*****@*****.**'
info.publication = 'JHEP 11 (2017) 029 , http://dx.doi.org/10.1007/JHEP11(2017)029'
info.comment = 'Moriond 2017'
#+++++++ dataset block ++++++++++++++
dataset = DataSetInput('data')
dataset.setInfo(dataType='upperLimit', dataId=None)
#+++++txName block +++++++++++++++++
TChiWH = dataset.addTxName('TChiWH')
TChiWH.checked = ''
TChiWH.constraint = "[[['W']],[['higgs']]]"
TChiWH.condition = None
TChiWH.conditionDescription = None
TChiWH.source = "CMS"
TChiWH.massConstraint = None
#++++++++++mass plane++++++++++++++++++
TChiWH_1 = TChiWH.addMassPlane(2 * [[x, y]])
'observedN': 4,
'expectedBG': 6.9,
'bgError': 1.5,
'index': 5
},
"SR6jt": {
'observedN': 3,
'expectedBG': 4.2,
'bgError': 1.2,
'index': 6
}
}
for SR, data in datasets.items():
#+++++++ dataset block ++++++++++++++
dataset = DataSetInput(SR)
dataset.setInfo(dataType='efficiencyMap',
dataId=SR,
observedN=data['observedN'],
expectedBG=data['expectedBG'],
bgError=data['bgError'])
T1 = dataset.addTxName('T1')
T1.checked = 'NO'
T1.constraint = "[[['jet','jet']],[['jet','jet']]]"
T1.conditionDescription = None
T1.condition = None
T1.source = "SModelS"
#+++++++ next mass plane block ++++++++++++++
T1_1 = T1.addMassPlane([[x, y]] * 2)
T1_1.addSource('efficiencyMap',
'orig/ATLAS_2016_I1458270_eff.dat',
#+++++++ global info block ++++++++++++++
info = MetaInfoInput('ATLAS-SUSY-2015-06')
info.url = 'http://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PAPERS/SUSY-2015-06/'
info.sqrts = 13
info.lumi = 3.2
info.prettyName = '2-6 jets, 0 lep'
info.private = False
info.arxiv = 'https://arxiv.org/abs/1605.03814'
info.contact = 'ATLAS collaboration'
info.publication = 'http://link.springer.com/article/10.1140/epjc/s10052-016-4184-8'
info.comment = 'UL analyses does not have digital data available'
#info.supersedes =
#info.supersededBy =
#+++++++ dataset block ++++++++++++++
dataset = DataSetInput("SR6jm")
dataset.setInfo(dataType='efficiencyMap',
dataId="SR6jm",
observedN=4,
expectedBG=6.9,
bgError=1.5)
T2 = dataset.addTxName('T2')
T2.checked = 'NO'
T2.constraint = "[[['jet']],[['jet']]]"
T2.conditionDescription = None
T2.condition = None
T2.source = "ATLAS"
#+++++++ next mass plane block ++++++++++++++
T2_1 = T2.addMassPlane([[x, y]] * 2)
#---- new efficiency map -----
#----figure----
'''
obsN = [18, 9, 11, 206, 53, 38, 20, 22, 22, 1, 0, 3]
expN = [18.7, 9.3, 8.6, 179., 52.4, 39.3, 20.6, 27.6, 18.9, 7.7, 0.91, 3.64]
bgErr = [2.7, 2.2, 2.1, 26., 7.4, 7.6, 6.5, 4.9, 3.4, 1.2, 0.73, 0.79]
SR = [
'SRAT0', 'SRATW', 'SRATT', 'SRBT0', 'SRBTW', 'SRBTT', 'SRC1', 'SRC2',
'SRC3', 'SRC4', 'SRC5', 'SRE'
]
fig = [
'09c', '09b', '09a', '09c', '09b', '09a', '10a', '10b', '10c', '10d',
'10e', '11b'
]
for i in range(len(obsN)):
#+++++++ dataset block ++++++++++++++
dataset = DataSetInput(SR[i])
dataset.setInfo(dataType='efficiencyMap',
dataId=SR[i],
observedN=obsN[i],
expectedBG=expN[i],
bgError=bgErr[i])
#+++++++ next txName block ++++++++++++++
T2tt = dataset.addTxName('T2tt')
T2tt.checked = 'False'
T2tt.constraint = '[[[t]],[[t]]]'
T2tt.conditionDescription = None
T2tt.condition = None
T2tt.massConstraint = [['dm >= 169.0'], ['dm >= 169.0']]
T2tt.source = 'ATLAS'
#+++++++ next txName block ++++++++++++++
T2ttoff = dataset.addTxName('T2ttoff')
