if args.smodelsPath:
sys.path.append(os.path.abspath(args.smodelsPath))
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
databaseCreator.ncpus = 1
#+++++++ global info block ++++++++++++++
info = MetaInfoInput('ATLAS-SUSY-2016-14')
info.url = 'http://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PAPERS/SUSY-2016-14/'
info.sqrts = 13
info.lumi = 36.1
info.prettyName = '2 same-sign or 3 leptons + jets + MET'
info.private = False
info.arxiv = 'https://arxiv.org/abs/1706.03731'
info.contact = '*****@*****.**'
info.publication = 'https://link.springer.com/article/10.1007/JHEP09(2017)084'
info.comment = 'Moriond 2017. Omitted RPV SUSY and long cascade topologies.'
#+++++++ 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
from utilsPath import utilsPath
utilsPath = databaseRoot + utilsPath
if args.smodelsPath:
sys.path.append(os.path.abspath(args.smodelsPath))
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')
if args.smodelsPath:
sys.path.append(os.path.abspath(args.smodelsPath))
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
databaseCreator.ncpus = 1
#+++++++ global info block ++++++++++++++
info = MetaInfoInput('CMS-SUS-15-008')
info.url = 'http://cms-results.web.cern.ch/cms-results/public-results/publications/SUS-15-008/index.html'
info.sqrts = 13
info.lumi = 2.3
info.prettyName = 'same-sign dilepton'
info.private = False
info.arxiv = 'http://arxiv.org/abs/1605.03171'
info.publication = 'http://dx.doi.org/10.1140/epjc/s10052-016-4261-z'
info.comment = 'Only T1tttt can be implemneted (other topologies have only up to 2 mass planes). Superseded by CMS-SUS-16-035.'
info.supersededBy = 'CMS-SUS-16-035'
#+++++++ dataset block ++++++++++++++
dataset = DataSetInput('data')
dataset.setInfo(dataType='upperLimit', dataId=None)
#+++++++ next txName block ++++++++++++++
T1tttt = dataset.addTxName('T1tttt')
T1tttt.constraint = "[[['t','t']],[['t','t']]]"
T1tttt.conditionDescription = None
T1tttt.condition = None
utilsPath = databaseRoot + utilsPath
if args.smodelsPath:
sys.path.append(os.path.abspath(args.smodelsPath))
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('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)
from utilsPath import utilsPath
utilsPath = databaseRoot + utilsPath
if args.smodelsPath:
sys.path.append(os.path.abspath(args.smodelsPath))
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('ATLAS-SUSY-2015-09')
info.url = 'https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PAPERS/SUSY-2015-09/'
info.sqrts = 13
info.lumi = 3.2
info.prettyName = 'jets + 2 SS or >=3 leptons'
info.private = False
info.arxiv = 'http://arxiv.org/abs/1602.09058'
info.contact = ''
info.publication = ''
info.comment = ''
info.supersededBy = ''
#+++++++ 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']]]"
from utilsPath import utilsPath
utilsPath = databaseRoot + utilsPath
if args.smodelsPath:
sys.path.append(os.path.abspath(args.smodelsPath))
sys.path.append(os.path.abspath(utilsPath))
from smodels_utils.dataPreparation.inputObjects import MetaInfoInput, DataSetInput, getSignalRegionsEMBaked, getStatsEMBaked
from smodels_utils.dataPreparation.databaseCreation import databaseCreator
from smodels_utils.dataPreparation.massPlaneObjects import x, y, z
#+++++++ 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)'
stats = getStatsEMBaked()
topos = [
"T1", "T2", "T6WWleft", "T5WW", "T5ZZ", "T5WZh", "T6WZh", "T5WWoff",
"T6WWoffleft", "TGQ"
]
# topos = [ "TGQ" ]
#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]]]", "T6WZh": "[[[q],[Z]],[[q],[W]]]+[[[q],[higgs]],[[q],[W]]]", "T5WWoff": "2.23*[[['jet','jet'],['jet','jet']],[['jet','jet'],['jet','jet']]]", "TGQ": "[[[q]],[[q,q]]]" }
constraints = {
"T2": "[[[q]],[[q]]]",
from utilsPath import utilsPath
utilsPath = databaseRoot + utilsPath
if args.smodelsPath:
sys.path.append(os.path.abspath(args.smodelsPath))
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-014')
info.url = 'http://cms-results.web.cern.ch/cms-results/public-results/preliminary-results/SUS-16-014/'
info.sqrts = 13
info.lumi = 12.9
info.prettyName = 'jets + Etmiss, HT'
info.private = False
info.comment = 'Only CDS entry:https://cds.cern.ch/record/2205158. Superseded by CMS-SUS-16-033.'
info.supersededBy = 'CMS-SUS-16-033'
#+++++++ 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"
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('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 ++++++++++++++
from utilsPath import utilsPath
utilsPath = databaseRoot + utilsPath
if args.smodelsPath:
sys.path.append(os.path.abspath(args.smodelsPath))
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-SUS-17-001')
info.url = 'http://cms-results.web.cern.ch/cms-results/public-results/publications/SUS-17-001/index.html'
info.sqrts = 13
info.lumi = 35.9
info.prettyName = 'Stop search in dilepton + jets + Etmiss final state'
info.private = False
info.arxiv = 'http://arxiv.org/abs/1711.00752'
info.contact = ''
info.publication = 'to appear in Phys. Rev. D '
info.comment = 'Moriond 2017. A long cascade decay not implemented: stop -> b chargino, chargino-> slepton v, slepton -> l lsp with chargino fixed halfway between the stop and lsp and 3 mass planes provided for the intermediate slepton at x = 0.05,0.5 and 0.95 the mass of the chargino. '
#+++++++ dataset block ++++++++++++++
dataset = DataSetInput('data')
dataset.setInfo(dataType='upperLimit', dataId=None)
#+++++txName block +++++++++++++++++
T2tt = dataset.addTxName('T2tt')
T2tt.checked = ''
T2tt.constraint = "[[['t']],[['t']]]"
from utilsPath import utilsPath
utilsPath = databaseRoot + utilsPath
if args.smodelsPath:
sys.path.append(os.path.abspath(args.smodelsPath))
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-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']]]"
if args.smodelsPath:
sys.path.append(os.path.abspath(args.smodelsPath))
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
databaseCreator.ncpus = 1
#+++++++ global info block ++++++++++++++
info = MetaInfoInput('ATLAS-SUSY-2016-17')
info.url = 'http://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PAPERS/SUSY-2016-17/'
info.sqrts = 13
info.lumi = 36.1
info.prettyName = '2 opposite sign leptons + MET'
info.private = False
info.arxiv = 'https://arxiv.org/abs/1708.03247'
info.contact = '*****@*****.**'
info.publication = 'https://link.springer.com/article/10.1140/epjc/s10052-017-5445-x'
info.comment = ''
#+++++++ dataset block ++++++++++++++
dataset = DataSetInput('data')
dataset.setInfo(dataType = 'upperLimit', dataId = None)
#+++++++txName block++++++++++++++++++++
T2tt=dataset.addTxName('T2tt')
T2tt.checked=''
T2tt.constraint="[[['t']],[['t']]]"
from utilsPath import utilsPath
utilsPath = databaseRoot + utilsPath
if args.smodelsPath:
sys.path.append(os.path.abspath(args.smodelsPath))
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-17-004')
info.url = 'http://cms-results.web.cern.ch/cms-results/public-results/preliminary-results/SUS-17-004/index.html'
info.sqrts = 13
info.lumi = 35.9
info.prettyName = 'Multilepton EWK searches'
info.private = False
info.arxiv = ''
info.contact = ''
info.publication = ''
info.comment = 'TChiHH, TChiZZ mass planes not provided. TChiWH/Z not implemented as BR(Chi20->H Chi10) = BR(Chi20->Chi10 Z)=0.5 provided for a single mass plane.'
#+++++++ dataset block ++++++++++++++
dataset = DataSetInput('data')
dataset.setInfo(dataType='upperLimit', dataId=None)
#+++++txName block +++++++++++++++++
TChiWZ = dataset.addTxName('TChiWZ')
TChiWZ.checked = ''
TChiWZ.constraint = "[[['W']],[['Z']]]"
if args.smodelsPath:
sys.path.append(os.path.abspath(args.smodelsPath))
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
databaseCreator.ncpus = 1
#+++++++ global info block ++++++++++++++
info = MetaInfoInput('CMS-PAS-SUS-16-015')
info.url = 'http://cms-results.web.cern.ch/cms-results/public-results/preliminary-results/SUS-16-015/'
info.sqrts = 13
info.lumi = 12.9
info.prettyName = 'jets + Etmiss, MT2'
info.private = False
info.comment = 'https://cds.cern.ch/record/2205162. Superseded by CMS-SUS-16-036.'
info.supersededBy = 'CMS-SUS-16-036'
#+++++++ 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"
if args.smodelsPath:
sys.path.append(os.path.abspath(args.smodelsPath))
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('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
if args.smodelsPath:
sys.path.append(os.path.abspath(args.smodelsPath))
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('ATLAS-CONF-2013-007')
info.comment = 'A technical problem has been found in the pseudo-experiments used to make Table 4 in the conference note of March 1. (Partly) superseded by ATLAS-SUSY-2013-09, fastlim maps contain more topologies for ATLAS-CONF-2013-007. Topologies included here are, T1, T1bbbb, T1bbbt, T1bbqq, T1bbtt, T1btbt, T1btqq, T1bttt, T1qqtt, T1tttt, T2, T2bb, T2bt, T2tt, T5bbbb, T5bbbt, T5btbt, T5tbtb, T5tbtt, T5tttt, TGQ, TGQbbq, TGQbtq, TGQQtt.'
info.sqrts = '8.0'
info.private = False
info.lumi = '20.7'
info.url = 'https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/CONFNOTES/ATLAS-CONF-2013-007/'
info.prettyName = '2 SS leptons + 0-3 b-jets + Etmiss'
info.implementedBy = 'MT'
#+++++++ dataset block ++++++++++++++
dataset = DataSetInput('data')
dataset.setInfo(dataType='upperLimit', dataId=None)
#+++++++ next txName block ++++++++++++++
T6ttWW = dataset.addTxName('T6ttWW')
T6ttWW.checked = 'VM'
T6ttWW.constraint = "[[['t+'],['W-']],[['t-'],['W+']]]"
T6ttWW.conditionDescription = "None"
T6ttWW.condition = "None"
T6ttWW.massConstraint = None
T6ttWW.source = 'ATLAS'
T6ttWWoff = dataset.addTxName('T6ttWWoff')
sys.path.append(os.path.abspath(args.smodelsPath))
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.datasetCreation import DatasetsFromLatex
from smodels_utils.dataPreparation.massPlaneObjects import x, y, z
DataSetInput.ntoys = 50000
#+++++++ global info block ++++++++++++++
info = MetaInfoInput('CMS-PAS-SUS-16-052-agg')
info.url = 'http://cms-results.web.cern.ch/cms-results/public-results/preliminary-results/SUS-16-052/index.html'
info.sqrts = '13.0*TeV'
info.lumi = 35.9
info.prettyName = 'soft lepton, <= 2 jets'
info.private = False
# info.arxiv = 'https://arxiv.org/abs/1704.07781'
info.comment = 'https://cds.cern.ch/record/2273394, http://inspirehep.net/record/1609006'
info.implementedBy = 'WW'
info.contact = 'SModelS'
# info.combinableWith = "CMS-PAS-SUS-16-050" ## FIXME can be combined with a lot, since it vetoes > 2 jets!
max_datasets = None
aggregate = None
## all 44 signal regions. agreement 97%. around 3.3s per point.
## corr > .6: 23 agg regions, agreement 86%. around 2s per point on commodore.
# aggregate = [[1], [2, 3, 4, 18, 19, 21, 22], [5], [6, 7, 8], [9], [10, 11, 12], [13], [14, 15, 16], [17, 20], [23], [24, 25, 26], [27], [28, 29, 30], [31], [32, 33, 34], [35], [36, 37, 38], [39], [40], [41], [42], [43], [44]]
from utilsPath import utilsPath
utilsPath = databaseRoot + utilsPath
if args.smodelsPath:
sys.path.append(os.path.abspath(args.smodelsPath))
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('ATLAS-SUSY-2016-16')
info.url = 'https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PAPERS/SUSY-2016-16/'
info.sqrts = 13
info.lumi = 36.1
info.prettyName = 'stop to lsp + MET'
info.private = False
info.arxiv = 'https://arxiv.org/abs/1711.11520'
info.contact = 'ATLAS collaboration'
info.publication = 'JHEP 06 (2018) 108'
SR = {
'obsN': [8, 50, 19, 115, 34, 68, 70],
'expN': [3.8, 36.3, 18.3, 115, 30.3, 71, 60.5],
'bgErr': [1, 6.6, 2.2, 31, 5.9, 16, 6.1],
'SR': [
'tN_high', 'tN_med', 'tN_diag_high', 'tN_diag_med', 'tN_diag_low',
'bWN', 'bffN'
],
'NsigUL': [10.0, 31.0, 11.0, 58.0, 17.0, 31.0, 28.0],
'NsigULexp': [5.8, 19.0, 11.0, 58.0, 19.0, 33.0, 21.0]
from utilsPath import utilsPath
utilsPath = databaseRoot + utilsPath
if args.smodelsPath:
sys.path.append(os.path.abspath(args.smodelsPath))
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-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
from utilsPath import utilsPath
utilsPath = databaseRoot + utilsPath
if args.smodelsPath:
sys.path.append(os.path.abspath(args.smodelsPath))
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-SUS-16-034')
info.url = 'http://cms-results.web.cern.ch/cms-results/public-results/publications/SUS-16-034/index.html'
info.sqrts = 13
info.lumi = 35.9
info.prettyName = '2 OSSF leptons'
info.private = False
info.arxiv = 'https://arxiv.org/abs/1709.08908'
info.contact = '*****@*****.**'
info.publication = 'Submitted to J. High Energy Phys. '
info.comment = 'Moriond 2017. Fig. 10 with long cascade decay via chargino and slepton not implemented because of assumtions on masses and BRs.'
#+++++++ dataset block ++++++++++++++
dataset = DataSetInput('data')
dataset.setInfo(dataType='upperLimit', dataId=None)
#+++++txName block +++++++++++++++++
TChiWZ = dataset.addTxName('TChiWZ')
TChiWZ.checked = ''
TChiWZ.constraint = "[[['W']],[['Z']]]"
TChiWZ.condition = None
if args.smodelsPath:
sys.path.append(os.path.abspath(args.smodelsPath))
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-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']]]"
sys.path.append(os.path.abspath(utilsPath))
from smodels_utils.dataPreparation.inputObjects import MetaInfoInput,DataSetInput
from smodels_utils.dataPreparation.commandlineArgs import setEnv
setEnv ( args )
from smodels_utils.dataPreparation.databaseCreation import databaseCreator
from smodels_utils.dataPreparation.datasetCreation import DatasetsFromLatex,createAggregationList
from smodels_utils.dataPreparation.massPlaneObjects import x, y, z
# DataSetInput.ntoys = 50
#+++++++ global info block ++++++++++++++
info = MetaInfoInput('CMS-SUS-16-050-agg')
info.url = 'http://cms-results.web.cern.ch/cms-results/public-results/publications/SUS-16-050/index.html'
info.sqrts = '13.0*TeV'
info.lumi = 35.9
info.prettyName = 'hadronic top tagging'
info.private = False
info.arxiv = 'https://arxiv.org/abs/1719.11188'
info.comment = ''
info.implementedBy = 'HARG'
info.contact = 'SModelS'
info.supersedes = "CMS-PAS-SUS-16-050"
## we dont yet perform combinations
# info.combinableWith = "CMS-PAS-SUS-16-052, CMS-PAS-SUS-16-024, CMS-PAS-SUS-16-022"
max_datasets = None
## all 84 signal regions take around 16s on commodore??
## corr>.4: 60 aggregate regions, agreements 98/97/97. Takes around 30s(?) per point on wnouc.
from utilsPath import utilsPath
utilsPath = databaseRoot + utilsPath
if args.smodelsPath:
sys.path.append(os.path.abspath(args.smodelsPath))
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-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']
from utilsPath import utilsPath
utilsPath = databaseRoot + utilsPath
if args.smodelsPath:
sys.path.append(os.path.abspath(args.smodelsPath))
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('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
sys.path.append(os.path.abspath(args.smodelsPath))
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
databaseCreator.ncpus = 1
#+++++++ global info block ++++++++++++++
info = MetaInfoInput('CMS-PAS-SUS-16-019')
info.url = 'http://cms-results.web.cern.ch/cms-results/public-results/preliminary-results/SUS-16-019/index.html'
info.sqrts = 13
info.lumi = 12.9
info.prettyName = 'jets + 1 lepton'
info.private = False
info.comment = 'https://cds.cern.ch/record/2204932. Superseded by CMS-SUS-16-037.'
info.supersededBy = 'CMS-SUS-16-037'
#+++++++ 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
if args.smodelsPath:
sys.path.append(os.path.abspath(args.smodelsPath))
sys.path.append(os.path.abspath(utilsPath))
from smodels_utils.dataPreparation import dataHandlerObjects
dataHandlerObjects.allowTrimming = True
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-SUS-16-033')
info.url = 'http://cms-results.web.cern.ch/cms-results/public-results/publications/SUS-16-033/index.html'
info.sqrts = 13
info.lumi = 35.9
info.prettyName = '0L + jets + MET (using MHT)'
info.private = False
info.arxiv = '1704.07781'
info.contact = '*****@*****.**'
info.publication = 'Phys. Rev D 96 (2017) 032003, http://dx.doi.org/10.1103/PhysRevD.96.032003'
info.comment = 'Moriond 2017. Not implemented: T5qqqqVV (Fig.12d) because decays to W and Z are summed over and T1tbtb because of fixed chargino mass (not specified). Added expected ULs (WW, june 2019).'
#+++++++ dataset block ++++++++++++++
dataset = DataSetInput('data')
dataset.setInfo(dataType='upperLimit', dataId=None)
#+++++txName block +++++++++++++++++
T1 = dataset.addTxName('T1')
T1.checked = ''
T1.constraint = "[[['jet','jet']],[['jet','jet']]]"
T1.condition = None
from utilsPath import utilsPath
utilsPath = databaseRoot + utilsPath
if args.smodelsPath:
sys.path.append(os.path.abspath(args.smodelsPath))
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('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
from utilsPath import utilsPath
utilsPath = databaseRoot + utilsPath
if args.smodelsPath:
sys.path.append(os.path.abspath(args.smodelsPath))
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-SUS-16-051')
info.url = 'http://cms-results.web.cern.ch/cms-results/public-results/publications/SUS-16-051/index.html'
info.sqrts = 13
info.lumi = 35.9
info.prettyName = '1L stop'
info.private = False
info.arxiv = 'http://arxiv.org/abs/1706.04402'
info.contact = '*****@*****.**'
info.publication = 'JHEP 10 (2017) 019, http://dx.doi.org/10.1007/JHEP10(2017)019'
info.comment = 'Moriond 2017. Fig. 7 with asymmetric decay of stop to b+chargino and top+neutralino not implemented because both, BRs and chargino mass are fixed.'
#+++++++ dataset block ++++++++++++++
dataset = DataSetInput('data')
dataset.setInfo(dataType='upperLimit', dataId=None)
#+++++txName block +++++++++++++++++
T2tt = dataset.addTxName('T2tt')
T2tt.checked = ''
T2tt.constraint = "[[['t']],[['t']]]"
from utilsPath import utilsPath
utilsPath = databaseRoot + utilsPath
if args.smodelsPath:
sys.path.append(os.path.abspath(args.smodelsPath))
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-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']]]"
if args.smodelsPath:
sys.path.append(os.path.abspath(args.smodelsPath))
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('ATLAS-SUSY-2017-01')
info.url = 'https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PAPERS/SUSY-2017-01/'
info.sqrts = 13
info.lumi = 36.1
info.prettyName = 'MET'
info.private = False
info.arxiv = 'arXiv:1812.09432'
info.contact = 'ATLAS collaboration'
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']}
if args.smodelsPath:
sys.path.append(os.path.abspath(args.smodelsPath))
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
DataSetInput.ntoys = args.ntoys
#+++++++ 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)
utilsPath = databaseRoot + utilsPath
if args.smodelsPath:
sys.path.append(os.path.abspath(args.smodelsPath))
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-EXO-13-006')
info.url = 'http://cms-results.web.cern.ch/cms-results/public-results/publications/EXO-13-006/index.html'
info.sqrts = 8
info.lumi = 18.8
info.prettyName ='hscp search'
info.private = False
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):