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)
T2 = dataset.addTxName('T2')
T2.checked = 'NO'
T2.constraint = "[[['jet']],[['jet']]]"
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
T1tttt.conditionDescription = None
T1tttt.source = "ATLAS"
T1tttt.massConstraint = [['dm>=338.0'], ['dm>=338.0']]
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']]]"
TChiWH.condition = None
TChiWH.conditionDescription = None
TChiWH.source = "CMS"
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']]]"
T1tttt.conditionDescription = None
T1tttt.condition = None
T1tttt.source = "ATLAS"
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']
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'
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 = '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'
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]]
## corr > .7, removing 0 points (18,28): 26 agg regions. agreement 90%. around 2.2s per point on commodore.
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?)
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], [24], [26], [27], [31], [35], [39], [40], [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], [78], [79], [80], [81], [83], [84]]
from smodels_utils.dataPreparation.datasetCreation import DatasetsFromLatex
from smodels_utils.dataPreparation.massPlaneObjects import x, y, z
# DataSetInput.ntoys = 200
#+++++++ global info block ++++++++++++++
info = MetaInfoInput('CMS-PAS-SUS-16-052')
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/....'
info.comment = 'EPS 2017. https://cds.cern.ch/record/2273394, http://inspirehep.net/record/1609006'
info.implementedBy = 'WW'
info.contact = 'CMS'
dataset = DataSetInput('data')
dataset.setInfo(dataType='upperLimit', dataId=None)
#+++++++ next txName block ++++++++++++++
T2bbWWoff = dataset.addTxName('T2bbWWoff')
T2bbWWoff.checked = ''
T2bbWWoff.constraint = "3.47*([[['b','l','nu']],[['b','jet','jet']]])"
#T2bbWWoff.constraint = "3.47*([[['b','mu','nu']],[['b','jet','jet']]]+[[['b','e','nu']],[['b','jet','jet']]])"
T2bbWWoff.conditionDescription = None
# T2bbWWoff.massConstraint = [['dm <= 76.']]*2
T2bbWWoff.condition = None
T2bbWWoff.source = 'CMS' #+++++++ next mass plane block ++++++++++++++
T2bbWWoff_1 = T2bbWWoff.addMassPlane([[x, x - y]] * 2)
#----exclusion source----
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,