Load the browser

In [1]:
#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)
44

Check the upper limit of a efficiency map result (95% C.L. limit for a given signal region)

In [2]:
#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')
1.30E+00 [fb]
1.35E+00 [fb]

Check the upper limit for a upper limit result (95% C.L. limit for a given simplified model/txname and mass)

In [3]:
#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]])
1.70E-01 [pb]
In [4]:
#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]])
8.00E-02 [pb]
In [5]:
#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]])
17:56:56.501 INFO     smodels.experiment.txnameObject:322 projected value is None. Projected point not in convex hull? original point=[[4.00E+03 [GeV], 1.50E+02 [GeV]], [4.00E+03 [GeV], 1.50E+02 [GeV]]]
None
In [ ]: