# Set up the path to SModelS installation folder
import sys; sys.path.append("."); import smodels_paths
from smodels.tools.physicsUnits import GeV
from smodels.tools import databaseBrowser
# First load the browser and check the total number of results loaded
browser = databaseBrowser.Browser("official")
print(len(browser))
# Check the upper limit for the HM200 signal region:
print(browser.getULForSR(expid='CMS-SUS-13-011', datasetID='HM200'))
# Check the upper limit for the LM300 signal region:
print(browser.getULForSR(expid='CMS-SUS-13-011', datasetID='LM300'))
# Check the upper limit for T1tttt simplified model with gluino mass 1 TeV and LSP mass 100 GeV:
print(browser.getULFor(expid='CMS-SUS-16-037', txname='T1tttt', massarray = [[1000.*GeV,100.*GeV],[1000.*GeV,100.*GeV]]))
# Check the upper limit for T1tttt simplified model with gluino mass 1.5 TeV and LSP mass 100 GeV:
print(browser.getULFor(expid='CMS-SUS-16-037', txname='T1tttt', massarray = [[1500.*GeV,100.*GeV],[1500.*GeV,100.*GeV]]))
# No results will be given if masses outside the data grid are used:
print(browser.getULFor(expid='CMS-SUS-16-037', txname='T1tttt', massarray = [[5000.*GeV,100.*GeV],[5000.*GeV,100.*GeV]]))