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') T1ttttoff.checked = '' T1ttttoff.constraint = "[[['b','W','b','W']],[['b','W','b','W']]]" T1ttttoff.condition = None T1ttttoff.source = "ATLAS" T1ttttoff.massConstraint = [['160 < dm < 338.0'], ['160 < dm < 338.0']] T1ttofftt = dataset.addTxName('T1ttofftt')
from smodels_utils.dataPreparation.massPlaneObjects import x, y, z #+++++++ global info block ++++++++++++++ info = MetaInfoInput('CMS-PAS-SUS-16-022') info.url = 'http://cms-results.web.cern.ch/cms-results/public-results/preliminary-results/SUS-16-022/' info.sqrts = 13 info.lumi = 12.9 info.prettyName = '>= 3 leptons' info.comment = 'https://cds.cern.ch/record/2205165' #+++++++ dataset block ++++++++++++++ dataset = DataSetInput('data') dataset.setInfo(dataType='upperLimit', dataId=None) #+++++++ next txName block ++++++++++++++ T1tttt = dataset.addTxName('T1tttt') T1tttt.checked = '' T1tttt.constraint = "[[['t','t']],[['t','t']]]" T1tttt.conditionDescription = None T1tttt.condition = None T1tttt.source = "CMS" T1tttt.massConstraint = None T1ttttoff = dataset.addTxName('T1ttttoff') T1ttttoff.checked = '' T1ttttoff.constraint = "[[['b','W','b','W']],[['b','W','b','W']]]" T1ttttoff.conditionDescription = None T1ttttoff.condition = None T1ttttoff.massConstraint = [['dm <= 338.0'], ['dm <= 338.0']] T1ttttoff.source = "CMS" #+++++++ next mass plane block ++++++++++++++ T1tttt_1 = T1tttt.addMassPlane(2 * [[x, y]])
] #I've artificially increase the expUL a bit, so the SR with the largest Mreco is selected for i, name in enumerate(datasetNames): #+++++++ dataset block ++++++++++++++ dataset = DataSetInput(name) dataset.setInfo(dataType='efficiencyMap', dataId=name, observedN=observedNs[i], expectedBG=expectedBGs[i], bgError=bgErrors[i], upperLimit=obsUpperLimits[i], expectedUpperLimit=expUpperLimits[i]) #+++++++ txnames ++++++++++++++++++++ #+++++++ next txName block ++++++++++++++ HSCPM1 = dataset.addTxName('THSCPM1b') HSCPM1.checked = '' HSCPM1.constraint = "[[],[]]" HSCPM1.condition = None HSCPM1.finalState = ['HSCP', 'HSCP'] HSCPM1.massConstraints = None HSCPM1.dataUrl = None HSCPM1.source = 'SModelS' #+++++++ next mass plane block ++++++++++++++ plane = HSCPM1.addMassPlane([[x], [x]]) plane.setSources( dataLabels=['efficiencyMap'], dataFiles=['orig/effmap_M1_stau_cons_mre' + name + '_trim.txt'], dataFormats=['txt']) # plane.addSource(dataLabels='obsExclusion',dataFiles='orig/CMS-PAS-EXO-16-036_Figure_003-b.dat', dataFormats='txt', unit='pb') #+++++++ next txName block ++++++++++++++
info.private = False info.arxiv = 'https://arxiv.org/abs/1606.08772v2' info.contact = 'ATLAS collaboration' info.publication = 'http://link.springer.com/article/10.1140/epjc/s10052-016-4382-4' #info.comment = #info.supersedes = #info.supersededBy = #+++++++ dataset block ++++++++++++++ dataset = DataSetInput("SRB") dataset.setInfo(dataType='efficiencyMap', dataId="SRB", observedN=6, expectedBG=12, bgError=2.5) T2bb = dataset.addTxName('T2bb') T2bb.checked = 'NO' T2bb.constraint = "[[['b']],[['b']]]" T2bb.conditionDescription = None T2bb.condition = None T2bb.source = "ATLAS" #+++++++ next mass plane block ++++++++++++++ T2bb_1 = T2bb.addMassPlane([[x, y]] * 2) #---- new efficiency map ----- #----figure---- T2bb_1.figureUrl = "https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PAPERS/SUSY-2015-01/figaux_06e.png; https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PAPERS/SUSY-2015-01/figaux_07e.png" T2bb_1.figure = "Fig.6e,Fig.7e" #----exclusion source---- T2bb_1.addSource('obsExclusion', 'orig/T2bb_Obs_Excl.dat', 'txt',
SR = {'obsN' : [7, 1, 5, 7, 6], 'expN' : [8, 2.5, 4.6, 2.8, 5.7], 'bgErr' : [4, 1.3, 1.2, 1, 2.3], 'SR' : ['SRHad-Low', 'SRHad-High', 'SR1Lbb-High', 'SR1Lbb-Medium', 'SR1Lbb-Low']} DATA = {'expExcl' : ['orig/Expectedlimit0lbb.csv', 'orig/Expectedlimit0lbb.csv', 'orig/Expectedlimit1lbb.csv', 'orig/Expectedlimit1lbb.csv', 'orig/Expectedlimit1lbb.csv'], 'obsExcl' : ['orig/Observedlimit0lbb.csv', 'orig/Observedlimit0lbb.csv', 'orig/Observedlimit1lbb.csv', 'orig/Observedlimit1lbb.csv', 'orig/Observedlimit1lbb.csv'], 'effMap' : ['orig/EffMap_TChiWH_SRHad-Low.txt', 'orig/EffMap_TChiWH_SRHad-High.txt', 'orig/EffMap_TChiWH_SR1Lbb-High.txt', 'orig/EffMap_TChiWH_SR1Lbb-Medium.txt', 'orig/EffMap_TChiWH_SR1Lbb-Low.txt'], 'fig' : ['12a.png', '12a.png', '12b.png', '12b.png', '12b.png']} for i in range(len(SR['obsN'])): #+++++++ dataset block ++++++++++++++ dataset = DataSetInput(SR['SR'][i]) dataset.setInfo(dataType = 'efficiencyMap', dataId = SR['SR'][i], observedN = SR['obsN'][i], expectedBG = SR['expN'][i], bgError = SR['bgErr'][i]) #+++++++ next txName block ++++++++++++++ newTx = dataset.addTxName('TChiWH') newTx.checked = 'No' newTx.constraint = "[[['W']],[['higgs']]]" newTx.conditionDescription = None newTx.condition = None newTx.source = 'ATLAS' #+++++++ next mass plane block ++++++++++++++ newPlane = newTx.addMassPlane(2*[[x, y]]) newPlane.figure = 'Fig.12a' newPlane.figureUrl = 'https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PAPERS/SUSY-2017-01/fig_' + DATA['fig'][i] newPlane.dataUrl = 'https://www.hepdata.net/record/ins1711261?version=1&table=Upper limit 0lbb' newPlane.setSources(dataLabels = ['expExclusion', 'obsExclusion', 'efficiencyMap'], dataFiles = [DATA['expExcl'][i], DATA['obsExcl'][i], DATA['effMap'][i]], dataFormats = ['csv', 'csv', 'txt']) databaseCreator.create()
info.comment = 'superseding publication contain more Data' info.sqrts = '8.0' info.private = False info.lumi = '13.0' info.url = 'https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/CONFNOTES/ATLAS-CONF-2012-166/' info.supersededBy = 'ATLAS-SUSY-2013-15' info.prettyName = '1 lepton + 4(1 b-)jets + Etmiss' info.implementedBy = 'MT' #+++++++ dataset block ++++++++++++++ dataset = DataSetInput('data') dataset.setInfo(dataType = 'upperLimit', dataId = None) #+++++++ next txName block ++++++++++++++ T2tt = dataset.addTxName('T2tt') T2tt.checked ="AL" T2tt.constraint ="[[['t']],[['t']]]" T2tt.conditionDescription ="None" T2tt.condition ="None" T2tt.source = "ATLAS" #+++++++ next mass plane block ++++++++++++++ T2tt = T2tt.addMassPlane(2*[[x, y]]) T2tt.figure = 'Fig.(aux) 3' T2tt.figureUrl = 'https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/CONFNOTES/ATLAS-CONF-2012-166//figaux_03.png' T2tt.dataUrl = 'Not defined' T2tt.setSources(dataLabels= ['obsExclusion', 'upperLimits'], dataFiles= ['orig/exclusion_T2tt.txt', 'orig/T2tt.txt'], dataFormats= ['txt', 'txt']) #+++++++ next txName block ++++++++++++++
info.url = 'https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PAPERS/SUSY-2016-24/' info.sqrts = 13 info.lumi = 36.1 info.prettyName = '2+ leptons (e,mu) + jets + Etmiss' info.private = False info.arxiv = 'https://arxiv.org/abs/1803.02762' info.contact = 'ATLAS collaboration' info.publication ='Eur. Phys. J. C 78 (2018) 995' #+++++++ dataset block ++++++++++++++ dataset = DataSetInput('data') dataset.setInfo(dataType = 'upperLimit', dataId = None) #+++++++ next txName block ++++++++++++++ TChipChimSlepSlep = dataset.addTxName('TChipChimSlepSlep') TChipChimSlepSlep.checked = 'No' TChipChimSlepSlep.constraint = "[[['nu'],['e+']],[['nu'],['e-']]]+[[['nu'],['e-']],[['nu'],['e+']]]+[[['nu'],['mu+']],[['nu'],['mu-']]]+[[['nu'],['mu-']],[['nu'],['mu+']]]" TChipChimSlepSlep.conditionDescription = None TChipChimSlepSlep.condition = None TChipChimSlepSlep.source = 'ATLAS' #+++++++ next mass plane block ++++++++++++++ TChipChimSlepSlep_1 = TChipChimSlepSlep.addMassPlane(2*[[x, 0.5*(x+y), y]]) TChipChimSlepSlep_1.figure = 'Fig.8a' TChipChimSlepSlep_1.figureUrl = 'https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PAPERS/SUSY-2016-24/fig_08a.png' TChipChimSlepSlep_1.dataUrl = 'https://www.hepdata.net/record/ins1658902?version=1&table=Table78' TChipChimSlepSlep_1.setSources(dataLabels = ['expExclusion', 'obsExclusion', 'upperLimits'], dataFiles = ['orig/HEPData-ins1658902-v1-Table_13.csv', 'orig/HEPData-ins1658902-v1-Table_14.csv', 'orig/HEPData-ins1658902-v1-Table_78.csv'], units = [ None, None, 'fb' ], coordinates = [ {x: 0, y: 1, 'value': None}, {x: 0, y: 1, 'value': None}, {x : 1, y: 0, 'value' :2} ], dataFormats = ['csv', 'csv', 'csv'])
'observedN': 3, 'expectedBG': 4.2, 'bgError': 1.2, 'index': 6 } } for SR, data in datasets.items(): #+++++++ dataset block ++++++++++++++ dataset = DataSetInput(SR) dataset.setInfo(dataType='efficiencyMap', dataId=SR, observedN=data['observedN'], expectedBG=data['expectedBG'], bgError=data['bgError']) T1 = dataset.addTxName('T1') T1.checked = 'NO' T1.constraint = "[[['jet','jet']],[['jet','jet']]]" T1.conditionDescription = None T1.condition = None T1.source = "SModelS" #+++++++ next mass plane block ++++++++++++++ T1_1 = T1.addMassPlane([[x, y]] * 2) T1_1.addSource('efficiencyMap', 'orig/ATLAS_2016_I1458270_eff.dat', 'txt', coordinateMap={ x: 0, y: 1, 'eff': data['index'] + 2 })
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 T2tt.source = "CMS" T2ttoff = dataset.addTxName('T2ttoff') T2ttoff.checked = '' T2ttoff.constraint = "[[['b','W']],[['b','W']]]" T2ttoff.condition = None T2ttoff.conditionDescription = None T2ttoff.massConstraint = [['80<dm<169'], ['80<dm<169']]
info.sqrts = 13 info.lumi = 35.9 info.prettyName = 'Photon + HT' info.private = False info.arxiv = 'http://arxiv.org/abs/1707.06193' info.contact = '*****@*****.**' info.publication = 'to appear in J. High Energy Phys.' info.comment = 'Moriond 2017. Several UL maps in the context of gauge-mediated SSB. Implemented T5gg and T6gg; others with mixed BRs are not useable. NB results only for high gluino mass.' #+++++++ dataset block ++++++++++++++ dataset = DataSetInput('data') dataset.setInfo(dataType='upperLimit', dataId=None) #+++++txName block +++++++++++++++++ T5gg = dataset.addTxName('T5gg') T5gg.checked = '' T5gg.constraint = "[[['jet','jet'],['photon']],[['jet','jet'],['photon']]]" T5gg.condition = None T5gg.conditionDescription = None T5gg.source = "CMS" #++++++next mass plane block+++++++++ T5gg_1 = T5gg.addMassPlane(2 * [[x, y, 1.0]]) T5gg_1.figure = 'Fig. 6-c' T5gg_1.figureUrl = 'http://cms-results.web.cern.ch/cms-results/public-results/publications/SUS-16-047/CMS-SUS-16-047_Figure_006-c.png' T5gg_1.dataUrl = 'http://cms-results.web.cern.ch/cms-results/public-results/publications/SUS-16-047/CMS-SUS-16-047_Figure_006-c.root' T5gg_1.histoDataUrl = 'http://cms-results.web.cern.ch/cms-results/public-results/publications/SUS-16-047/CMS-SUS-16-047_Figure_006-c.root' T5gg_1.exclusionDataUrl = 'http://cms-results.web.cern.ch/cms-results/public-results/publications/SUS-16-047/CMS-SUS-16-047_Figure_006-c.root' T5gg_1.setSources(
info.contact = '*****@*****.**' info.publication = '' # +++++++ dataset block ++++++++++++++ SRA_incl = DataSetInput('SRA_incl') SRA_incl.setInfo( dataType='efficiencyMap', dataId="SRA_incl", observedN=17, expectedBG=17.1, bgError=2.8, # upperLimit = '2.637E-01*fb', expectedUpperLimit = '2.639E-01*fb', ## FIXME this line is nonsense jsonfile="orig/BkgOnlyA.json") # +++++++ next txName block ++++++++++++++ T6bbHH = SRA_incl.addTxName("T6bbHH") T6bbHH.checked = 'NO' T6bbHH.constraint = "[[['b'],['higgs']],[['b'],['higgs']]]" T6bbHH.conditionDescription = None T6bbHH.condition = None T6bbHH.source = 'ATLAS' # +++++++ next mass plane block ++++++++++++++ T6bbHHM60 = T6bbHH.addMassPlane(2 * [[x, y, 60.]]) T6bbHHM60.figure = "figaux_05b (acceptance) and figaux_9b (efficiency)" T6bbHHM60.figureUrl = "https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PAPERS/SUSY-2018-31/figaux_05b.png" T6bbHHM60.dataUrl = "https://www.hepdata.net/record/ins1748602?version=1&table=Efficiency_SRB and others" T6bbHHM60.addSource("obsExclusion", "orig/HEPData-ins1748602-v1-M60_Obs.csv", unit=None, dataFormat="csv")
info.private = False info.arxiv = 'https://arxiv.org/abs/1606.03903' info.source = "ATLAS" info.publication = 'http://journals.aps.org/prd/abstract/10.1103/PhysRevD.94.052009' info.comment = 'Very weird exclusion due to best CLs selection, and fluctuation in the Obs data' #info.supersedes = #info.supersededBy = #+++++++ dataset block ++++++++++++++ dataset = DataSetInput("SR2") dataset.setInfo(dataType='efficiencyMap', dataId="SR2", observedN=1, expectedBG=1.25, bgError=0.26) T2tt = dataset.addTxName('T2tt') T2tt.checked = 'NO' T2tt.constraint = "[[['t']],[['t']]]" T2tt.conditionDescription = None T2tt.condition = None T2tt.source = "ATLAS" #+++++++ next mass plane block ++++++++++++++ T2tt_1 = T2tt.addMassPlane([[x, y]] * 2) #---- new efficiency map ----- #----figure---- T2tt_1.figureUrl = "http://hepdata.cedar.ac.uk/resource/1469069/figAuxiliaryFigure10b.png; http://hepdata.cedar.ac.uk/resource/1469069/figAuxiliaryFigure13b.png" T2tt_1.figure = "Aux.Fig.10b,Aux.Fig.13b" #----exclusion source---- T2tt_1.addSource('obsExclusion', 'orig/T2tt_ATLAS_SUSY_2015_02_Obs.txt', 'txt',
info.contact = 'ATLAS collaboration' info.publication = '' # +++++++ dataset block ++++++++++++++ SRA_L = DataSetInput('SRA_L') SRA_L.setInfo( dataType='efficiencyMap', dataId="SRA_L", observedN=12, expectedBG=8.4, bgError=1.7, # upperLimit = '2.637E-01*fb', expectedUpperLimit = '2.639E-01*fb', ## FIXME this line is nonsense jsonfile="orig/BkgOnlyA.json") # +++++++ next txName block ++++++++++++++ T6bbHH = SRA_L.addTxName("T6bbHH") T6bbHH.checked = 'NO' T6bbHH.constraint = "[[['b'],['higgs']],[['b'],['higgs']]]" T6bbHH.conditionDescription = None T6bbHH.condition = None T6bbHH.source = 'ATLAS' # +++++++ next mass plane block ++++++++++++++ T6bbHHM60 = T6bbHH.addMassPlane(2 * [[x, y, 60.]]) T6bbHHM60.figure = "figaux_02d(acceptance) and figaux_06d (efficiency)" T6bbHHM60.figureUrl = "https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PAPERS/SUSY-2018-31/figaux_02d.png" T6bbHHM60.dataUrl = "https://www.hepdata.net/record/ins1748602?version=1&table=Efficiency_SRB and others" T6bbHHM60.addSource("obsExclusion", "orig/HEPData-ins1748602-v1-M60_Obs.csv", unit=None, dataFormat="csv")
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---- T2bbWWoff_1.addSource('upperLimits', 'orig/CMS-PAS-SUS-16-052_Figure_004.root', 'canvas', objectName='cCONT_', index=2) T2bbWWoff_1.addSource( 'obsExclusion', 'orig/CMS-PAS-SUS-16-052_Figure_004.root',\
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 TChiWZoff = dataset.addTxName('TChiWZoff') TChiWZoff.checked = '' TChiWZoff.constraint = "71.*([[['mu+','mu-']],[['l','nu']]] + [[['e+','e-']],[['l','nu']]])" TChiWZoff.condition = "cGtr([[['mu+','mu-']],[['l','nu']]],[[['e+','e-']],[['l','nu']]])" TChiWZoff.massConstraint = [['dm < 86.0'], ['dm < 76.0']] TChiWZoff.conditionDescription = None TChiWZoff.source = "CMS"
info = MetaInfoInput('ATLAS-SUSY-2016-19') info.url = "https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PAPERS/SUSY-2016-19/" info.sqrts = 13 info.lumi = 36.1 info.prettyName = "stops to staus" info.private = False info.arxiv = 'https://arxiv.org/abs/1803.10178' info.contact = 'ATLAS collaboration' info.publication = 'in process of being published' #+++++++ dataset block ++++++++++++++ dataset = DataSetInput('data') dataset.setInfo(dataType='upperLimit', dataId=None) #+++++++ next txName block ++++++++++++++ T4bnutaubnutau = dataset.addTxName('T4bnutaubnutau') T4bnutaubnutau.checked = '' T4bnutaubnutau.constraint = "[[[b,nu],[ta]],[[b,nu],[ta]]]" T4bnutaubnutau.conditionDescription = None T4bnutaubnutau.condition = None T4bnutaubnutau.source = "ATLAS" #+++++++ next mass plane block ++++++++++++++ T4bnutaubnutau_1 = T4bnutaubnutau.addMassPlane(2 * [[x, y, 0.]]) T4bnutaubnutau_1.figure = 'Fig.7' T4bnutaubnutau_1.figureUrl = 'https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PAPERS/SUSY-2016-19/fig_07.png' T4bnutaubnutau_1.setSources( dataLabels=['expExclusion', 'obsExclusion', 'upperLimits'], dataFiles=[ 'orig/Exp_Excl.csv', 'orig/Obs_Excl.csv', 'orig/Obs_UL_fixed.csv' ], units=[None, None, 'fb'],
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 TChiWZoff = dataset.addTxName('TChiWZoff') TChiWZoff.checked = '' TChiWZoff.constraint = "71.*([[['mu+','mu-']],[['l','nu']]] + [[['e+','e-']],[['l','nu']]])" TChiWZoff.condition = "cGtr([[['mu+','mu-']],[['l','nu']]],[[['e+','e-']],[['l','nu']]])" TChiWZoff.massConstraint = [['dm < 86.0'], ['dm < 76.0']] TChiWZoff.conditionDescription = None TChiWZoff.source = "CMS"
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+++++++++ T1_1 = T1.addMassPlane(2 * [[x, y]]) T1_1.figure = 'Fig. 12-c' T1_1.figureUrl = 'http://cms-results.web.cern.ch/cms-results/public-results/publications/SUS-16-033/CMS-SUS-16-033_Figure_012-c.png' T1_1.dataUrl = 'http://cms-results.web.cern.ch/cms-results/public-results/publications/SUS-16-033/CMS-SUS-16-033_Figure_012-c.root' T1_1.histoDataUrl = 'http://cms-results.web.cern.ch/cms-results/public-results/publications/SUS-16-033/CMS-SUS-16-033_Figure_012-c.root' T1_1.exclusionDataUrl = 'http://cms-results.web.cern.ch/cms-results/public-results/publications/SUS-16-033/CMS-SUS-16-033_Figure_012-c.root' names = [
info.sqrts = 13 info.lumi = 35.9 info.prettyName = '0L + jets + Etmiss (using MT2)' info.private = False info.arxiv = 'https://arxiv.org/abs/1705.04650' info.contact = '*****@*****.**' info.publication = 'Eur. Phys. J. C 77 (2017) 710, http://dx.doi.org/10.1140/epjc/s10052-017-5267-x' info.comment = 'Moriond 2017. The mixed decay topology for stop pair production with stop->t \chi1_0 (50% BR) and stop->b \chi^{\pm} -> b W^{\pm} \chi1_0 (50%), has not been implemented since the branching ratio is fixed and the mass difference between the chargino and neutralino is fixed to 5 GeV. ' #+++++++ dataset block ++++++++++++++ dataset = DataSetInput('data') dataset.setInfo(dataType='upperLimit', dataId=None) #+++++txName block +++++++++++++++++ T1bbbb = dataset.addTxName('T1bbbb') T1bbbb.checked = '' T1bbbb.constraint = "[[['b','b']],[['b','b']]]" T1bbbb.condition = None T1bbbb.conditionDescription = None T1bbbb.source = "CMS" #++++++next mass plane block+++++++++ T1bbbb_1 = T1bbbb.addMassPlane(2 * [[x, y]]) T1bbbb_1.figure = 'Fig. 6-a' T1bbbb_1.figureUrl = 'http://cms-results.web.cern.ch/cms-results/public-results/publications/SUS-16-036/CMS-SUS-16-036_Figure_006-a.png' T1bbbb_1.dataUrl = 'http://cms-results.web.cern.ch/cms-results/public-results/publications/SUS-16-036/CMS-SUS-16-036_Figure_006-a.root' T1bbbb_1.histoDataUrl = 'http://cms-results.web.cern.ch/cms-results/public-results/publications/SUS-16-036/CMS-SUS-16-036_Figure_006-a.root' T1bbbb_1.exclusionDataUrl = 'http://cms-results.web.cern.ch/cms-results/public-results/publications/SUS-16-036/CMS-SUS-16-036_Figure_006-a.root' labels = [
"T5ZZ": (2*[[x, y, 1.]], None), "T5WW": ( 2*[[x, 0.5*(x+y), y]], 2*[[x, y , 60.]] ), "T5WWoff": ( 2*[[x, 0.5*(x+y), y]], 2*[[x, y, 60.]] ) # "T5WW": ( 2*[[x, 0.5*(x+y), y]], 2*[[x, 60. + y * (x - 60.), 60.]] ), # "T5WWoff": ( 2*[[x, 0.5*(x+y), y]], 2*[[x, 60. + y * (x - 60.), 60.]] ) # "T6WW": ( 2*[[x, 0.5*(x+y), y]], 2*[[x, 60. + y * (x - 60.), 60.]] ), # "T6WWoffleft": ( 2*[[x, 0.5*(x+y), y]], 2*[[x, 60. + y * (x - 60.), 60.]] ), # "T5WZh": (2*[[x, 60. + y * (x - 60.), 60.]],None), # "T6WZh": (2*[[x, 60. + y * (x - 60.), 60.]], None ), } rjrs = [ 1, 2, 3, 5 ] for topo in topos: #+++++++ next txName block ++++++++++++++ Tx = dataset.addTxName( topo ) Tx.checked = '' Tx.constraint = constraints[topo] Tx.conditionDescription = None Tx.condition = None if topo in conditions.keys(): Tx.condition = conditions[topo] Tx.source = "ATLAS" mPlanes = massPlanes[topo] if type (mPlanes) == list: Tx_1 = Tx.addMassPlane(massPlanes[topo]) Tx_1.figure = figure[topo] Tx_1.figureUrl = baseUrl + figureUrl[topo] eIdx = exclusionIndex[topo] meffrjr = "Meff" meffrjrurl = "Meff"
#+++++++ 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') T6ttWWoff.constraint = "3.5*([[['t+'],['l-','nu']],[['t-'],['jet','jet']]] + [[['t-'],['l+','nu']],[['t+'],['jet','jet']]])" T6ttWWoff.conditionDescription = "[[['t'],['mu','nu']],[['t'],['jet','jet']]]>[[['t'],['e','nu']],[['t'],['jet','jet']]]" T6ttWWoff.condition = "Cgtr([[['t'],['mu','nu']],[['t'],['jet','jet']]],[[['t'],['e','nu']],[['t'],['jet','jet']]])" T6ttWWoff.massConstraint = [['dm >= 169.0', 'dm <= 76.0'], ['dm >= 169.0', 'dm <= 76.0']] T6ttWWoff.source = 'ATLAS' #+++++++ next mass plane block ++++++++++++++ T6ttWWLSP060 = T6ttWW.addMassPlane(2 * [[x, y, 60.0]])
info.url = 'http://cms-results.web.cern.ch/cms-results/public-results/publications/SUS-16-035/index.html' info.sqrts = 13 info.lumi = 35.9 info.prettyName = '2 SS leptons' info.private = False info.arxiv = 'https://arxiv.org/abs/1704.07323' info.contact = '*****@*****.**' info.publication = 'Eur.Phys.J. C77 (2017) no.9, 578, http://dx.doi.org/10.1140/epjc/s10052-017-5079-z' info.comment = 'Moriond 2017. Seven SMS interpretations in paper, all implemented.' #+++++++ dataset block ++++++++++++++ dataset = DataSetInput('data') dataset.setInfo(dataType='upperLimit', dataId=None) #+++++++txName block++++++++++++++++++++ T1tttt = dataset.addTxName('T1tttt') T1tttt.checked = '' T1tttt.constraint = "[[['t','t']],[['t','t']]]" T1tttt.condition = None T1tttt.conditionDescription = None T1tttt.source = "CMS" T1tttt.massConstraint = [['dm>=338.0'], ['dm>=338.0']] T1ttttoff = dataset.addTxName('T1ttttoff') T1ttttoff.checked = '' T1ttttoff.constraint = "[[[b,W,b,W]],[[b,W,b,W]]]" T1ttttoff.condition = None T1ttttoff.source = "CMS" T1ttttoff.massConstraint = [['dm < 338.0'], ['dm < 338.0']] #++++++next mass plane block+++++++++
'[[[b],[W]],[[b],[W]]]', '[[[b],[L,nu]],[[b],[L,nu]]]+[[[b],[L,nu]],[[b],[q,q]]]+[[[b],[q,q]],[[b],[q,q]]]' ], 'massPlane': 2 * [[x, x - 10, y]] } DATA = [T2tt] #+++++++ dataset block ++++++++++++++ dataset = DataSetInput('data') dataset.setInfo(dataType='upperLimit', dataId=None) for TX in DATA: #+++++++ next txName block ++++++++++++++ newTx = dataset.addTxName(TX['name'][0]) newTx.checked = 'False' newTx.constraint = TX['constraint'][0] newTx.conditionDescription = None newTx.condition = None newTx.source = 'ATLAS' #+++++++ next txName block ++++++++++++++ newTxOff1 = dataset.addTxName(TX['name'][1]) newTxOff1.checked = 'False' newTxOff1.constraint = TX['constraint'][1] newTxOff1.conditionDescription = None newTxOff1.condition = None newTxOff1.massConstraint = [['80 <= dm < 169.0'], ['80 <= dm < 169.0']] newTxOff1.source = 'ATLAS' #+++++++ next mass plane block ++++++++++++++ newTxOff2 = dataset.addTxName(TX['name'][2])
info = MetaInfoInput('CMS-PAS-EXO-16-036') info.url = 'http://cms-results.web.cern.ch/cms-results/public-results/preliminary-results/EXO-16-036/index.html' info.sqrts = 13 info.lumi = 12.9 info.prettyName = 'hscp search' info.private = False info.contact = '*****@*****.**' info.comment = 'Upper limits digitized from Track+TOF png plots. Used conservative gluino bounds (50% gluino-gluon probability) and curve for the cloud model for the squark R-hadron constraints.' #+++++++ dataset block ++++++++++++++ dataset = DataSetInput('data') dataset.setInfo(dataType='upperLimit', dataId=None) #+++++++ txnames ++++++++++++++++++++ #+++++++ next txName block ++++++++++++++ HSCPM1 = dataset.addTxName('THSCPM1b') HSCPM1.checked = '' HSCPM1.constraint = "[[],[]]" HSCPM1.condition = None HSCPM1.finalState = ['HSCP', 'HSCP'] HSCPM1.massConstraints = None HSCPM1.dataUrl = 'http://cms-results.web.cern.ch/cms-results/public-results/preliminary-results/EXO-16-036/CMS-PAS-EXO-16-036_Figure_003-b.png' HSCPM1.source = 'CMS' #+++++++ next mass plane block ++++++++++++++ plane = HSCPM1.addMassPlane([[x], [x]]) plane.setSources(dataLabels=['upperLimits'], dataFiles=['orig/CMS-PAS-EXO-16-036_Figure_003-b_stauDP.dat'], dataFormats=['txt'], units=['pb']) #+++++++ next txName block ++++++++++++++
#+++++++ dataset block ++++++++++++++ dataset = DataSetInput(sr) dataset.setInfo(dataType='efficiencyMap', dataId=sr, observedN=SR['obsN'][i], expectedBG=SR['expN'][i], bgError=SR['bgErr'][i], upperLimit=str(SR['NsigUL'][i] / 36.1).strip() + '*fb', expectedUpperLimit=str(SR['NsigULexp'][i] / 36.1).strip() + '*fb') #+++++++ next txName block ++++++++++++++ for TX in Txnames: #The SRs 'bffN' and 'tN_diag_low' only apply to T2tt and T2bbll if sr in ['bffN', 'tN_diag_low'] and TX['name'] != 'T2bbll': continue newTx = dataset.addTxName(TX['name']) newTx.checked = 'False' newTx.constraint = TX['constraint'] newTx.conditionDescription = None newTx.condition = None newTx.source = 'ATLAS' newTx.massConstraint = TX['massConstr'] for ip, plane in enumerate(TX['massPlanes']): if not os.path.isfile(TX['sources'][ip].replace('SR', sr)): continue #+++++++ next mass plane block ++++++++++++++ newPlane = newTx.addMassPlane(plane) newPlane.figure = TX['info']['figure'] newPlane.figureUrl = TX['info']['figureUrl'] newPlane.dataUrl = TX['info']['dataUrl'] newPlane.addSource('efficiencyMap',
info = MetaInfoInput('ATLAS-SUSY-2016-27') info.url = "https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PAPERS/SUSY-2016-27/" info.sqrts = 13 info.lumi = 36.1 info.prettyName = " " info.private = False info.arxiv = 'https://arxiv.org/abs/1802.03158' info.contact = 'ATLAS collaboration' info.publication = 'Phys. Rev. D 97, 092006 (2018)' #+++++++ dataset block ++++++++++++++ dataset = DataSetInput('data') dataset.setInfo(dataType='upperLimit', dataId=None) #+++++++ next txName block ++++++++++++++ T5Gamma = dataset.addTxName('T5Gamma') T5Gamma.checked = 'No' T5Gamma.constraint = "[[['jet','jet'],['photon']],[['jet','jet'],['photon']]]" T5Gamma.conditionDescription = None T5Gamma.condition = None T5Gamma.source = "ATLAS" #+++++++ next mass plane block ++++++++++++++ T5Gamma_1 = T5Gamma.addMassPlane(2 * [[x, y, 1.]]) T5Gamma_1.figure = 'Fig.8' T5Gamma_1.figureUrl = 'https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PAPERS/SUSY-2016-27/fig_08.png' T5Gamma_1.dataUrl = "https://www.hepdata.net/record/ins1654357?version=1&table=Cross section UL 1" T5Gamma_1.setSources( dataLabels=['expExclusion', 'obsExclusion', 'upperLimits'], dataFiles=[ 'orig/HEPData-ins1654357-v1-Exclusion_contour_(expected)_1.csv', 'orig/HEPData-ins1654357-v1-Exclusion_contour_(observed)_2.csv',
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+++++++++ TChiWZ_1 = TChiWZ.addMassPlane(2 * [[x, y]]) TChiWZ_1.figure = 'Fig. 8' TChiWZ_1.figureUrl = 'http://cms-results.web.cern.ch/cms-results/public-results/publications/SUS-16-034/CMS-SUS-16-034_Figure_008.png' TChiWZ_1.dataUrl = 'http://cms-results.web.cern.ch/cms-results/public-results/publications/SUS-16-034/CMS-SUS-16-034_Figure_008.root' TChiWZ_1.histoDataUrl = 'http://cms-results.web.cern.ch/cms-results/public-results/publications/SUS-16-034/CMS-SUS-16-034_Figure_008.root' TChiWZ_1.exclusionDataUrl = 'http://cms-results.web.cern.ch/cms-results/public-results/publications/SUS-16-034/CMS-SUS-16-034_Figure_008.root'
info.sqrts = 13 info.lumi = 35.9 info.prettyName = 'EWK WH(bb)' info.private = False info.arxiv = 'http://arxiv.org/abs/1706.09933' info.contact = '*****@*****.**' info.publication = 'JHEP 11 (2017) 029 , http://dx.doi.org/10.1007/JHEP11(2017)029' info.comment = 'Moriond 2017' #+++++++ dataset block ++++++++++++++ dataset = DataSetInput('data') dataset.setInfo(dataType='upperLimit', dataId=None) #+++++txName block +++++++++++++++++ TChiWH = dataset.addTxName('TChiWH') TChiWH.checked = '' TChiWH.constraint = "[[['W']],[['higgs']]]" TChiWH.condition = None TChiWH.conditionDescription = None TChiWH.source = "CMS" TChiWH.massConstraint = None #++++++++++mass plane++++++++++++++++++ TChiWH_1 = TChiWH.addMassPlane(2 * [[x, y]]) TChiWH_1.figure = 'Fig. 6' TChiWH_1.figureUrl = 'http://cms-results.web.cern.ch/cms-results/public-results/publications/SUS-16-043/CMS-SUS-16-043_Figure_006-b.png' TChiWH_1.dataUrl = 'http://cms-results.web.cern.ch/cms-results/public-results/publications/SUS-16-043/CMS-SUS-16-043_Figure_006-b.root' TChiWH_1.histoDataUrl = 'http://cms-results.web.cern.ch/cms-results/public-results/publications/SUS-16-043/CMS-SUS-16-043_Figure_006-b.root' TChiWH_1.exclusionDataUrl = 'http://cms-results.web.cern.ch/cms-results/public-results/publications/SUS-16-043/CMS-SUS-16-043_Figure_006-b.root'
#+++++++ global info block ++++++++++++++ info = MetaInfoInput('CMS-PAS-SUS-16-024') info.url = 'http://cms-results.web.cern.ch/cms-results/public-results/preliminary-results/SUS-16-024/' info.sqrts = 13 info.lumi = 12.9 info.private = False info.comment = 'https://cds.cern.ch/record/2205168. Superseded by CMS-SUS-16-039.' info.supersededBy = 'CMS-SUS-16-039' #+++++++ dataset block ++++++++++++++ dataset = DataSetInput('data') dataset.setInfo(dataType='upperLimit', dataId=None) #+++++++ next txName block ++++++++++++++ TChiChipmStauL = dataset.addTxName('TChiChipmStauL') TChiChipmStauL.checked = '' TChiChipmStauL.constraint = "[[['ta'],['ta']],[['nu'],['L']]]" TChiChipmStauL.conditionDescription = None TChiChipmStauL.condition = None TChiChipmStauL.source = "CMS" #+++++++ next mass plane block ++++++++++++++ TChiChipmStauL_1 = TChiChipmStauL.addMassPlane(2 * [[x, .5 * (x + y), y]]) TChiChipmStauL_1.figure = 'Figure 9' TChiChipmStauL_1.figureUrl = 'http://cms-results.web.cern.ch/cms-results/public-results/preliminary-results/SUS-16-024/CMS-PAS-SUS-16-024_Figure_009.png' TChiChipmStauL_1.dataUrl = 'http://cms-results.web.cern.ch/cms-results/public-results/preliminary-results/SUS-16-024/CMS-PAS-SUS-16-024_Figure_009.root' TChiChipmStauL_1.histoDataUrl = 'http://cms-results.web.cern.ch/cms-results/public-results/preliminary-results/SUS-16-024/CMS-PAS-SUS-16-024_Figure_009.root' TChiChipmStauL_1.dataUrl = 'http://cms-results.web.cern.ch/cms-results/public-results/preliminary-results/SUS-16-024/CMS-PAS-SUS-16-024_Figure_009.root' TChiChipmStauL_1.exclusionDataUrl = 'http://cms-results.web.cern.ch/cms-results/public-results/preliminary-results/SUS-16-024/CMS-PAS-SUS-16-024_Figure_009.root' TChiChipmStauL_1.setSources(dataLabels=[ 'expExclusion', 'expExclusionM1', 'expExclusionP1', 'expectedUpperLimits',
info.private = False info.arxiv = 'https://arxiv.org/abs/1605.03814' info.contact = 'ATLAS collaboration' info.publication = 'http://link.springer.com/article/10.1140/epjc/s10052-016-4184-8' info.comment = 'UL analyses does not have digital data available' #info.supersedes = #info.supersededBy = #+++++++ dataset block ++++++++++++++ dataset = DataSetInput("SR6jm") dataset.setInfo(dataType='efficiencyMap', dataId="SR6jm", observedN=4, expectedBG=6.9, bgError=1.5) T2 = dataset.addTxName('T2') T2.checked = 'NO' T2.constraint = "[[['jet']],[['jet']]]" T2.conditionDescription = None T2.condition = None T2.source = "ATLAS" #+++++++ next mass plane block ++++++++++++++ T2_1 = T2.addMassPlane([[x, y]] * 2) #---- new efficiency map ----- #----figure---- T2_1.figureUrl = "http://hepdata.cedar.ac.uk/resource/1458270/figAuxiliaryFigure24b.png" T2_1.figure = "Fig. Aux. 24b" #----exclusion source---- T2_1.addSource('obsExclusion', 'orig/Obs_Line_T2.dat', 'txt',
datasetNames = ['c000','c100','c200','c300'] observedNs = [42,7,0,0] expectedBGs = [44.,5.6,0.56,0.02] bgErrors = [9.,1.1,0.11,0.004] obsUpperLimits = ['1.15*fb','0.441*fb','0.16*fb','0.159*fb'] expUpperLimits = ['1.23*fb','0.338*fb','0.16*fb','0.159*fb'] for i,name in enumerate(datasetNames): #+++++++ dataset block ++++++++++++++ dataset = DataSetInput(name) dataset.setInfo(dataType = 'efficiencyMap', dataId = name, observedN=observedNs[i], expectedBG=expectedBGs[i], bgError=bgErrors[i], upperLimit = obsUpperLimits[i], expectedUpperLimit = expUpperLimits[i]) #+++++++ txnames ++++++++++++++++++++ #+++++++ next txName block ++++++++++++++ HSCPM1 = dataset.addTxName('THSCPM1') HSCPM1.checked ='' HSCPM1.constraint = "[[],[]]" HSCPM1.condition =None HSCPM1.finalState = ['HSCP','HSCP'] HSCPM1.massConstraints = None HSCPM1.dataUrl = None HSCPM1.source = 'SModelS' #+++++++ next mass plane block ++++++++++++++ plane = HSCPM1.addMassPlane([[x],[x]]) plane.setSources(dataLabels= ['efficiencyMap'],dataFiles=['orig/eff_HSCPM1_'+name+'.txt'], dataFormats=['txt']) #+++++++ next txName block ++++++++++++++ HSCPM3 = dataset.addTxName('THSCPM3') HSCPM3.checked ='' HSCPM3.constraint = "[[['*']],[['*']]]" ##Here '*' represents any (single) even particle HSCPM3.condition =None