#Make sure smodels can be found
import sys,os
sys.path.append(os.path.join(os.getenv("HOME"),'smodels'))
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(os.path.join(os.getenv("HOME"),'smodels-database'))
print len(browser)
#Check the upper limit for the ANA10-CUT1 signal region:
print browser.getULForSR(expid='ATLAS-CONF-2013-053', datasetID='ANA10-CUT1')
#Check the upper limit for the ANA6-CUT5 signal region:
print browser.getULForSR(expid='ATLAS-CONF-2013-053', datasetID='ANA6-CUT5')
#Check the upper limit for T2bb simplified model with sbottom mass 300 GeV and LSP mass 100 GeV:
print browser.getULFor(expid='ATLAS-CONF-2013-053', txname='T2bb', massarray = [[300.*GeV,100.*GeV],[300.*GeV,100.*GeV]])
#Check the upper limit for T2bb simplified model with sbottom mass 400 GeV and LSP mass 150 GeV:
print browser.getULFor(expid='ATLAS-CONF-2013-053', txname='T2bb', massarray = [[400.*GeV,150.*GeV],[400.*GeV,150.*GeV]])
#No results will be given if masses outside the data grid are uses:
print browser.getULFor(expid='ATLAS-CONF-2013-053', txname='T2bb', massarray = [[4000.*GeV,150.*GeV],[4000.*GeV,150.*GeV]])