info.publication = 'Phys. Rev. D 100, 012006 (2019)' SR = {'obsN' : [7, 1, 5, 7, 6], 'expN' : [8, 2.5, 4.6, 2.8, 5.7], 'bgErr' : [4, 1.3, 1.2, 1, 2.3], 'SR' : ['SRHad-Low', 'SRHad-High', 'SR1Lbb-High', 'SR1Lbb-Medium', 'SR1Lbb-Low']} DATA = {'expExcl' : ['orig/Expectedlimit0lbb.csv', 'orig/Expectedlimit0lbb.csv', 'orig/Expectedlimit1lbb.csv', 'orig/Expectedlimit1lbb.csv', 'orig/Expectedlimit1lbb.csv'], 'obsExcl' : ['orig/Observedlimit0lbb.csv', 'orig/Observedlimit0lbb.csv', 'orig/Observedlimit1lbb.csv', 'orig/Observedlimit1lbb.csv', 'orig/Observedlimit1lbb.csv'], 'effMap' : ['orig/EffMap_TChiWH_SRHad-Low.txt', 'orig/EffMap_TChiWH_SRHad-High.txt', 'orig/EffMap_TChiWH_SR1Lbb-High.txt', 'orig/EffMap_TChiWH_SR1Lbb-Medium.txt', 'orig/EffMap_TChiWH_SR1Lbb-Low.txt'], 'fig' : ['12a.png', '12a.png', '12b.png', '12b.png', '12b.png']} 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]],
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 = '' T1ttttoff.constraint = "[[['b','W','b','W']],[['b','W','b','W']]]" T1ttttoff.conditionDescription = None T1ttttoff.condition = None
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]) #+++++++ 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 ++++++++++++++
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' #+++++++ next mass plane block ++++++++++++++ plane = HSCPM1.addMassPlane([[x], [x]]) plane.setSources(dataLabels=['upperLimits'],
from smodels_utils.dataPreparation.databaseCreation import databaseCreator from smodels_utils.dataPreparation.massPlaneObjects import x, y, z #+++++++ global info block ++++++++++++++ info = MetaInfoInput('ATLAS-SUSY-2016-27') info.url = "https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PAPERS/SUSY-2016-27/" info.sqrts = 13 info.lumi = 36.1 info.prettyName = " " info.private = False info.arxiv = 'https://arxiv.org/abs/1802.03158' info.contact = 'ATLAS collaboration' info.publication = 'Phys. Rev. D 97, 092006 (2018)' #+++++++ dataset block ++++++++++++++ dataset = DataSetInput('data') dataset.setInfo(dataType='upperLimit', dataId=None) #+++++++ next txName block ++++++++++++++ T5Gamma = dataset.addTxName('T5Gamma') T5Gamma.checked = 'No' T5Gamma.constraint = "[[['jet','jet'],['photon']],[['jet','jet'],['photon']]]" T5Gamma.conditionDescription = None T5Gamma.condition = None 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(
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') T1ttttoff.checked = '' T1ttttoff.constraint = "[[[b,W,b,W]],[[b,W,b,W]]]" T1ttttoff.condition = None
#+++++++ 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']] T1ttttoff=dataset.addTxName('T1ttttoff') T1ttttoff.checked='' T1ttttoff.constraint="[[['b','W','b','W']],[['b','W','b','W']]]" T1ttttoff.condition=None
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+++++++++ T5gg_1 = T5gg.addMassPlane(2 * [[x, y, 1.0]]) T5gg_1.figure = 'Fig. 6-c'
#+++++++ global info block ++++++++++++++ info = MetaInfoInput('ATLAS-SUSY-2016-24') info.url = 'https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PAPERS/SUSY-2016-24/' info.sqrts = 13 info.lumi = 36.1 info.prettyName = '2+ leptons (e,mu) + jets + Etmiss' info.private = False info.arxiv = 'https://arxiv.org/abs/1803.02762' info.contact = 'ATLAS collaboration' info.publication ='Eur. Phys. J. C 78 (2018) 995' #+++++++ dataset block ++++++++++++++ dataset = DataSetInput('data') dataset.setInfo(dataType = 'upperLimit', dataId = None) #+++++++ next txName block ++++++++++++++ TChipChimSlepSlep = dataset.addTxName('TChipChimSlepSlep') TChipChimSlepSlep.checked = 'No' TChipChimSlepSlep.constraint = "[[['nu'],['e+']],[['nu'],['e-']]]+[[['nu'],['e-']],[['nu'],['e+']]]+[[['nu'],['mu+']],[['nu'],['mu-']]]+[[['nu'],['mu-']],[['nu'],['mu+']]]" TChipChimSlepSlep.conditionDescription = None 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'
from smodels_utils.dataPreparation.databaseCreation import databaseCreator from smodels_utils.dataPreparation.massPlaneObjects import x, y, z #+++++++ global info block ++++++++++++++ info = MetaInfoInput('ATLAS-SUSY-2018-31') info.url = 'https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PAPERS/SUSY-2018-31/' info.sqrts = 13 info.lumi = 139. info.prettyName = 'higgs + b-jets + MET' info.private = False info.arxiv = 'arXiv:1908.03122' 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
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 #++++++next mass plane block+++++++++ T2bb_1 = T2bb.addMassPlane(2 * [[x, y]])
#+++++++ global info block ++++++++++++++ info = MetaInfoInput('ATLAS-SUSY-2015-02') info.url = 'http://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PAPERS/SUSY-2015-02/' info.sqrts = 13 info.lumi = 3.2 info.prettyName = 'Top 1l' 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----
from smodels_utils.dataPreparation.databaseCreation import databaseCreator from smodels_utils.dataPreparation.massPlaneObjects import x, y, z #+++++++ global info block ++++++++++++++ info = MetaInfoInput('ATLAS-SUSY-2018-31') info.url = 'https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PAPERS/SUSY-2018-31/' info.sqrts = 13 info.lumi = 139. info.prettyName = 'higgs + b-jets + MET' info.private = False info.arxiv = 'arXiv:1908.03122' 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
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)' obsN = [0, 0] expN = [0.5, 0.48] bgErr = [0.28, 0.275] SR = ['SRyy-SL', 'SRyy-SH'] for i in range(len(obsN)): #+++++++ dataset block ++++++++++++++ dataset = DataSetInput(SR[i] + '-q') dataset.setInfo(dataType = 'efficiencyMap', dataId = SR[i] + '-q', observedN = obsN[i], expectedBG = expN[i], bgError = bgErr[i]) #+++++++ next txName block ++++++++++++++ T6Gamma = dataset.addTxName('T6Gamma') T6Gamma.checked = 'No' T6Gamma.constraint = "[[['jet'],['photon']],[['jet'],['photon']]]" T6Gamma.conditionDescription = None T6Gamma.condition = None T6Gamma.source = "ATLAS" #+++++++ next mass plane block ++++++++++++++ T6Gamma_1 = T6Gamma.addMassPlane(2*[[x, y, 0.]]) T6Gamma_1.figure = 'Fig.9' T6Gamma_1.figureUrl = 'https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PAPERS/SUSY-2016-27/fig_09.png' T6Gamma_1.dataUrl = "https://www.hepdata.net/record/ins1654357?version=1&table=Cross section UL 2" T6Gamma_1.setSources(dataLabels = ['expExclusion', 'obsExclusion', 'efficiencyMap'], dataFiles = ['orig/HEPData-ins1654357-v1-Exclusion_contour_(expected)_3.csv', 'orig/HEPData-ins1654357-v1-Exclusion_contour_(observed)_4.csv', 'orig/EffMap_T6Gamma_' + SR[i] + '.txt'],
#+++++++ global info block ++++++++++++++ info = MetaInfoInput('ATLAS-SUSY-2016-07') info.url = "https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PAPERS/SUSY-2016-07/" info.sqrts = 13 info.lumi = 36.1 info.prettyName = "hadronic jets + MET" info.private = False info.arxiv = 'https://arxiv.org/abs/1712.02332' info.contact = 'ATLAS collaboration' info.publication = 'Phys. Rev. D 97, 112001 (2018)' #+++++++ dataset block ++++++++++++++ dataset = DataSetInput('data') dataset.setInfo(dataType = 'upperLimit', dataId = None) topos = [ "T1", "T2", "T6WWleft", "T5WW", "T5ZZ", "T5WZh", "T6WZh", "T5WWoff", "T6WWoffleft" ] # topos = [ "T5WZh", "T6WZh" ] constraints = { "T2": "[[[q]],[[q]]]", "T1": "[[[q,q]],[[q,q]]]", "T6WWleft": "[[['q'],['W']],[['q'],['W']]]", "T6WWoffleft": "2.23 * [[['jet'],['jet','jet']],[['jet'],['jet','jet']]]", "T5WW": "[[['q','q'],['W']],[['q','q'],['W']]]", "T5ZZ": "[[['q','q'],['Z']],[['q','q'],['Z']]]", "T5WZh": "[[[q,q],[Z]],[[q,q],[W]]]+[[[q,q],[higgs]],[[q,q],[W]]]+[[[q,q],[W]],[[q,q],[W]]]+[[[q,q],[Z]],[[q,q],[Z]]]+[[[q,q],[higgs]],[[q,q],[higgs]]]+[[[q,q],[higgs]],[[q,q],[Z]]]", "T6WZh": "[[[jet],[Z]],[[jet],[W]]]+[[[jet],[higgs]],[[jet],[W]]]+[[[jet],[W]],[[jet],[W]]]+[[[jet],[Z]],[[jet],[Z]]]+[[[jet],[higgs]],[[jet],[higgs]]]+[[[jet],[Z]],[[jet],[higgs]]]", "T5WWoff": "2.23*[[['jet','jet'],['jet','jet']],[['jet','jet'],['jet','jet']]]" } conditions = { "T5WZh": "Csim([[[q,q],[Z]],[[q,q],[W]]]+[[[q,q],[higgs]],[[q,q],[W]]],2.*([[[q,q],[W]],[[q,q],[W]]]))", "T6WZh": "Csim([[[jet],[Z]],[[jet],[W]]]+[[[jet],[higgs]],[[jet],[W]]],2.*([[[jet],[W]],[[jet],[W]]]))" } figure = { "T2": "Fig. 76a", "T1": "Fig. 76b", "T6WWleft": ('Fig. 77a', 'Fig. 77b' ), "T5WW": ( 'Fig. 77c', 'Fig. 77d' ), "T5ZZ": ('Fig.79', None) , "T5WZh": ('Fig. 78b', None ), "T6WZh": ('Fig. 78a', None ), "T5WWoff": ( 'Fig. 77c', 'Fig. 77d' ), "T6WWoffleft": ('Fig. 77a', 'Fig. 77b' ) } baseUrl = 'https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PAPERS/SUSY-2016-07/' figureUrl = { "T2": 'figaux_076a.png', "T1": 'figaux_076b.png', "T6WWleft": ( 'figaux_077a.png', 'figaux_077b.png' ), "T5WW": ( 'figaux_077c.png', 'figaux_077d.png' ),
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+++++++++ T1bbbb_1 = T1bbbb.addMassPlane(2 * [[x, y]]) T1bbbb_1.figure = 'Fig. 6-a'
#+++++++ global info block ++++++++++++++ info = MetaInfoInput('ATLAS-CONF-2012-166') info.comment = 'superseding publication contain more Data' info.sqrts = '8.0' info.private = False info.lumi = '13.0' info.url = 'https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/CONFNOTES/ATLAS-CONF-2012-166/' info.supersededBy = 'ATLAS-SUSY-2013-15' info.prettyName = '1 lepton + 4(1 b-)jets + Etmiss' info.implementedBy = 'MT' #+++++++ dataset block ++++++++++++++ dataset = DataSetInput('data') dataset.setInfo(dataType = 'upperLimit', dataId = None) #+++++++ next txName block ++++++++++++++ T2tt = dataset.addTxName('T2tt') T2tt.checked ="AL" T2tt.constraint ="[[['t']],[['t']]]" T2tt.conditionDescription ="None" T2tt.condition ="None" T2tt.source = "ATLAS" #+++++++ next mass plane block ++++++++++++++ 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'],
'expExcl': 'orig/HEPData-ins1639856-v4-Table_28.csv', 'obsExcl': 'orig/HEPData-ins1639856-v4-Table_29.csv', 'upLimit': 'orig/HEPData-ins1639856-v4-Table_70.csv' }, 'constraint': [ '[[[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
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 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']
'constraint': '[[[b, L, nu]],[[b, jet,jet]]]', 'massConstr': [['dm < 80'], ['dm < 80']], 'massPlanes': [2 * [[x, y]], 2 * [[x, x - y]]] } Txnames = [T2tt, T2ttoff, T2bbll] # Txnames = [T2tt,T2ttoff] exclusionCurves = { 'obs': ['orig/Table17.csv', 'orig/Table20.csv'], 'exp': ['orig/Table16.csv', 'orig/Table19.csv'] } for i, sr in enumerate(SR['SR']): #+++++++ 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'
"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.]]), "TGQ": ([[x, 0.], [y, 0.]], [[x, 695.], [y, 695.]], [[x, 995.], [y, 995.]]) } datasets = {} dsnames = getSignalRegionsEMBaked("orig/%s.embaked" % "T2") for dsname in dsnames: #+++++++ dataset block ++++++++++++++ if not dsname in stats: print("cannot find stats for %s" % dsname) continue dataset = DataSetInput(dsname) dst = stats[dsname] dataset.setInfo(dataType='efficiencyMap', dataId=dsname, observedN=dst["nobs"], expectedBG=dst["nb"], bgError=dst["deltanb"]) datasets[dsname] = dataset rjrs = [1, 2, 3, 5] for topo in topos: for dsname in dsnames: #+++++++ next txName block ++++++++++++++ Tx = datasets[dsname].addTxName(topo)
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]]) 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'
#+++++++ global info block ++++++++++++++ info = MetaInfoInput('ATLAS-SUSY-2015-01') info.url = "https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PAPERS/SUSY-2015-01/" info.sqrts = 13 info.lumi = 3.2 info.prettyName = "2b" info.private = False info.arxiv = 'https://arxiv.org/abs/1606.08772v2' info.contact = 'ATLAS collaboration' info.publication = 'http://link.springer.com/article/10.1140/epjc/s10052-016-4382-4' #info.comment = #info.supersedes = #info.supersededBy = #+++++++ dataset block ++++++++++++++ dataset = DataSetInput("SRB") dataset.setInfo(dataType='efficiencyMap', dataId="SRB", observedN=6, expectedBG=12, bgError=2.5) 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----
from smodels_utils.dataPreparation.databaseCreation import databaseCreator from smodels_utils.dataPreparation.massPlaneObjects import x, y, z #+++++++ global info block ++++++++++++++ info = MetaInfoInput('ATLAS-SUSY-2016-19') info.url = "https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PAPERS/SUSY-2016-19/" info.sqrts = 13 info.lumi = 36.1 info.prettyName = "stops to staus" info.private = False info.arxiv = 'https://arxiv.org/abs/1803.10178' info.contact = 'ATLAS collaboration' info.publication = 'in process of being published' #+++++++ dataset block ++++++++++++++ dataset = DataSetInput('data') dataset.setInfo(dataType='upperLimit', dataId=None) #+++++++ next txName block ++++++++++++++ T4bnutaubnutau = dataset.addTxName('T4bnutaubnutau') T4bnutaubnutau.checked = '' T4bnutaubnutau.constraint = "[[[b,nu],[ta]],[[b,nu],[ta]]]" T4bnutaubnutau.conditionDescription = None T4bnutaubnutau.condition = None 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'],
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 #++++++++++mass plane++++++++++++++++++ TChiWH_1 = TChiWH.addMassPlane(2 * [[x, y]])
'observedN': 4, 'expectedBG': 6.9, 'bgError': 1.5, 'index': 5 }, "SR6jt": { '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',
#+++++++ 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 T2.source = "ATLAS" #+++++++ next mass plane block ++++++++++++++ T2_1 = T2.addMassPlane([[x, y]] * 2) #---- new efficiency map ----- #----figure----
''' obsN = [18, 9, 11, 206, 53, 38, 20, 22, 22, 1, 0, 3] expN = [18.7, 9.3, 8.6, 179., 52.4, 39.3, 20.6, 27.6, 18.9, 7.7, 0.91, 3.64] bgErr = [2.7, 2.2, 2.1, 26., 7.4, 7.6, 6.5, 4.9, 3.4, 1.2, 0.73, 0.79] SR = [ 'SRAT0', 'SRATW', 'SRATT', 'SRBT0', 'SRBTW', 'SRBTT', 'SRC1', 'SRC2', 'SRC3', 'SRC4', 'SRC5', 'SRE' ] fig = [ '09c', '09b', '09a', '09c', '09b', '09a', '10a', '10b', '10c', '10d', '10e', '11b' ] for i in range(len(obsN)): #+++++++ dataset block ++++++++++++++ dataset = DataSetInput(SR[i]) dataset.setInfo(dataType='efficiencyMap', dataId=SR[i], observedN=obsN[i], expectedBG=expN[i], bgError=bgErr[i]) #+++++++ next txName block ++++++++++++++ T2tt = dataset.addTxName('T2tt') T2tt.checked = 'False' T2tt.constraint = '[[[t]],[[t]]]' T2tt.conditionDescription = None T2tt.condition = None T2tt.massConstraint = [['dm >= 169.0'], ['dm >= 169.0']] T2tt.source = 'ATLAS' #+++++++ next txName block ++++++++++++++ T2ttoff = dataset.addTxName('T2ttoff')