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')
T1ttttoff.checked = ''
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
T1tttt.source = "CMS"
T1tttt.massConstraint = None
T1ttttoff = dataset.addTxName('T1ttttoff')
T1ttttoff.constraint = "[[['b','W','b','W']],[['b','W','b','W']]]"
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
T1tttt.conditionDescription = None
T1tttt.source = "ATLAS"
T1tttt.massConstraint = [['dm>=338.0'], ['dm>=338.0']]
T1ttttoff = dataset.addTxName('T1ttttoff')
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 = 'Andre Lessa <*****@*****.**>; Jan Heisig <*****@*****.**>; SModelS'
info.comment = 'Search for long-lived charged particles implemented using a correction factor for the 8 TeV efficiencies.'
#++Define list of datasets++
datasetNames = ['c000', 'c100', 'c200', 'c300']
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
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']]]"
T1tttt.condition=None
T1tttt.conditionDescription = None
T1tttt.source="CMS"
T1tttt.massConstraint=[['dm>=338.0'],['dm>=338.0']]
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']]]"
T2tt.condition = None
T2tt.conditionDescription = None
T2tt.massConstraint = None
T2tt.source = "CMS"
sys.path.append(os.path.abspath(databaseRoot))
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-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')
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']]]"
TChiWZ.condition = None
TChiWZ.conditionDescription = None
TChiWZ.source = "CMS"
#offshell txName block
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']]]"
T2tt.condition = None
T2tt.conditionDescription = None
T2tt.source = "CMS"
T2tt.massConstraint = None
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-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+++++++++
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-039')
info.url = 'http://cms-results.web.cern.ch/cms-results/public-results/publications/SUS-16-039/index.html'
info.sqrts = 13
info.lumi = 35.9
info.prettyName = 'Multilepton EWK searches'
info.private = False
info.arxiv = 'https://arxiv.org/abs/1709.05406'
info.contact = '*****@*****.**'
info.publication = 'Submitted to J. High Energy Phys.'
info.comment = 'Moriond 2017. Negative values in TChiChipmSlepL and TChiChipmSlepStau root files. Added TChiChipmStauStau.'
#+++++++ 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
TChiWZ.conditionDescription = None
TChiWZ.source = "CMS"
#offshell txName block
from smodels_utils.dataPreparation.databaseCreation import databaseCreator
from smodels_utils.dataPreparation.massPlaneObjects import x, y, z
#+++++++ global info block ++++++++++++++
info = MetaInfoInput('CMS-PAS-SUS-15-002')
info.url = 'http://cms-results.web.cern.ch/cms-results/public-results/superseded/SUS-15-002/index.html'
info.sqrts = 13
info.lumi = 2.2
info.prettyName = '>= 4jets + Etmiss, HT, HTmiss'
info.private = False
info.arxiv = ''
info.contact = ''
info.publication = ''
info.comment = 'Only CDS entry https://cds.cern.ch/record/2114817'
info.supersededBy = 'CMS-SUS-15-002'
#+++++++ dataset block ++++++++++++++
dataset = DataSetInput('data')
dataset.setInfo(dataType = 'upperLimit', dataId = None)
#+++++++ next txName block ++++++++++++++
T1bbbb = dataset.addTxName('T1bbbb')
T1bbbb.checked = ''
T1bbbb.constraint = "[[['b','b']],[['b','b']]]"
T1bbbb.conditionDescription = None
T1bbbb.condition = None
T1bbbb.source = "CMS"
#+++++++ next mass plane block ++++++++++++++
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-049')
info.url = 'http://cms-results.web.cern.ch/cms-results/public-results/publications/SUS-16-049/index.html'
info.sqrts = 13
info.lumi = 35.9
info.prettyName = 'All hadronic stop'
info.private = False
info.arxiv = 'http://arxiv.org/abs/1707.03316'
info.contact = '*****@*****.**'
info.publication = 'J. High Energy Phys. 10 (2017) 005, http://dx.doi.org/10.1007/JHEP10(2017)005'
info.comment = 'Moriond 2017. Paper contains 6 SMS interpretations, 5 of which are implemented here. Not implemented: Fig. 9 with stop -> b chi+1 on one branch and stop -> t chi01 on the other branch, because both the BRs and the 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']]]"
T2tt.condition = None
T2tt.conditionDescription = None
T2tt.massConstraint = [[
'dm > 169.'
]] * 2 #Use only on-shell region to avoid interpolating in the excluded band
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-016')
info.url = "http://cms-results.web.cern.ch/cms-results/public-results/preliminary-results/SUS-16-016/index.html"
info.sqrts = 13
info.lumi = 12.9
info.prettyName = '>= 1 jet + Etmiss, alpha_T'
info.private = False
info.comment = "Only CDS entry https://cds.cern.ch/record/2205163. Superseded by CMS-SUS-16-033 and CMS-SUS-16-036."
info.supersededBy = 'CMS-SUS-16-033'
#+++++++ dataset block ++++++++++++++
dataset = DataSetInput('data')
dataset.setInfo(dataType='upperLimit', dataId=None)
#+++++++ next txName block ++++++++++++++
T1bbbb = dataset.addTxName('T1bbbb')
T1bbbb.checked = ''
T1bbbb.constraint = "[[['b','b']],[['b','b']]]"
T1bbbb.conditionDescription = None
T1bbbb.condition = None
T1bbbb.source = "CMS"
#+++++++ next mass plane block ++++++++++++++
T1bbbb_1 = T1bbbb.addMassPlane(2 * [[x, y]])
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-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
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']]]"
T1bbbb.condition = None
T1bbbb.conditionDescription = None
T1bbbb.source = "CMS"
#++++++next mass plane block+++++++++
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
T1.conditionDescription = None
T1.source = "CMS"
#++++++next mass plane block+++++++++
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"
T1tttt.massConstraint = None
T1ttttoff = dataset.addTxName('T1ttttoff')
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]]
## corr > .7: 28 agg regions. agreement 94%. around X.Xs per point on commodore.
## aggregate = [[1], [2, 3, 4, 19, 21, 22], [5], [6, 7, 8], [9], [10, 11, 12], [13], [14], [15], [16], [17], [18], [20], [23], [24, 25, 26], [27], [28], [29, 30], [31], [32, 33, 34], [35], [36, 37, 38], [39], [40], [41], [42], [43], [44]]
else:
databaseRoot = '../../../'
sys.path.append(os.path.abspath(databaseRoot))
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-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+']]]"
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
T1tttt.conditionDescription = None
T1tttt.source = "CMS"
T1tttt.massConstraint = [['dm>=338.0'], ['dm>=338.0']]
T1ttttoff = dataset.addTxName('T1ttttoff')
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
TChiWZ.conditionDescription = None
TChiWZ.source = "CMS"
TChiWZ.massConstraint = None
#++++++ mass plane block+++++++++
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']]]"
T1tttt.conditionDescription = None
T1tttt.condition = None
T1tttt.source = "ATLAS"
T1tttt.massConstraint = None
T1ttttoff = dataset.addTxName('T1ttttoff')
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.
# aggregation = [[1, 6, 10, 11, 22, 38], [2], [3], [4], [5], [7, 12, 13, 16, 17, 18, 21, 28, 29, 30, 32, 33, 34, 36, 37, 58, 68, 69, 70, 82], [8], [9], [14], [15], [19], [20], [23], [24], [25], [26], [27], [31], [35], [39], [40], [41], [42], [43], [44], [45], [46], [47], [48], [49], [50], [51], [52], [53], [54], [55], [56], [57], [59], [60], [61], [62], [63], [64], [65], [66], [67], [71], [72], [73], [74], [75], [76], [77], [78], [79], [80], [81], [83], [84]]
## corr>.4: remove 0-points (23,25,41,77), 56 aggregate regions, agreements 98/97/97. Takes around 10s per point on higgs. overexcludes for one case (T2tt?)
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']
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'
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):
#+++++++ dataset block ++++++++++++++
dataset = DataSetInput(name)
dataset.setInfo(dataType='efficiencyMap',
dataId=name,
observedN=observedNs[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
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
T1tttt.source = "CMS"
T1tttt.massConstraint = 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('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]])
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']]]"
TChiWH.condition = None
TChiWH.conditionDescription = None
TChiWH.source = "CMS"
TChiWH.massConstraint = None
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 = 'Efficiency maps created using Pythia 8.230 and Rivet 2.6.0'
#info.supersedes =
#info.supersededBy =
datasets = {
"SR2jl": {
'observedN': 263,
'expectedBG': 283,
'bgError': 24,
'index': 0
},
"SR2jm": {
'observedN': 191,
'expectedBG': 191,
'bgError': 21,
'index': 1
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)
T2 = dataset.addTxName('T2')
T2.checked = 'NO'
T2.constraint = "[[['jet']],[['jet']]]"
T2.conditionDescription = None
T2.condition = None
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):
#+++++++ 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])