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,