How To: Print out the theory predictions

In [1]:
# Set up the path to SModelS installation folder 
import sys; sys.path.append("."); import smodels_paths
In [2]:
# Import those parts of smodels that are needed for this exercise
# (We will assume the input is a SLHA file. For LHE files, use the lheDecomposer instead)
from smodels.theory import slhaDecomposer
from smodels.installation import installDirectory
from smodels.tools.physicsUnits import fb, GeV
from smodels.theory.theoryPrediction import theoryPredictionsFor
from smodels.experiment.databaseObj import Database
In [3]:
# load the official database:
database = Database( "official")
# Get list of desired experimental results (by default all results):
expResults = database.getExpResults()
In [4]:
# Define the SLHA input file name
filename="%s/inputFiles/slha/gluino_squarks.slha" % installDirectory()
In [5]:
# Perform the decomposition:
topList = slhaDecomposer.decompose(filename, sigcut = 0.03 * fb, doCompress=True, doInvisible=True,minmassgap = 5* GeV)
In [6]:
# Compute the theory prediction for each selected experimental result using the output from the decomposition:
allThPredictions = [theoryPredictionsFor(exp, topList) for exp in expResults]
In [7]:
# Print information about each theory prediction for each result:
for thPreds in allThPredictions:
    if not thPreds: continue #skip results with no predictions
    for theoryPred in thPreds:
        print ("\nExperimental Result: ",theoryPred.expResult.globalInfo.id,"(%s-type)" %theoryPred.dataset.dataInfo.dataType) #Result ID
        print ("Theory prediction xsec = ",theoryPred.xsection.value) #Signal xsection*efficiency*BR
        print ("Conditions violation (if any) = ",theoryPred.conditions) #Condition values (for UL-type results)
        print ("Dataset:",theoryPred.dataset.dataInfo.dataId)  #Corresponding signal region (for EM-type results)
        print ("Txnames = ",[str(tx) for tx in theoryPred.txnames]) #List of simplified models (txnames) contributing to the signal xsec        

Experimental Result:  ATLAS-SUSY-2015-06 (efficiencyMap-type)
Theory prediction xsec =  2.79E-04 [pb]
Conditions violation (if any) =  {'None': None}
Dataset: SR5j
Txnames =  ['T1']

Experimental Result:  CMS-PAS-SUS-16-022 (upperLimit-type)
Theory prediction xsec =  3.50E-05 [pb]
Conditions violation (if any) =  {'None': None}
Dataset: None
Txnames =  ['T1tttt']

Experimental Result:  CMS-PAS-SUS-17-004 (upperLimit-type)
Theory prediction xsec =  2.07E-01 [pb]
Conditions violation (if any) =  {'None': None}
Dataset: None
Txnames =  ['TChiWH']

Experimental Result:  CMS-PAS-SUS-17-004 (upperLimit-type)
Theory prediction xsec =  3.99E-02 [pb]
Conditions violation (if any) =  {'None': None}
Dataset: None
Txnames =  ['TChiWZ']

Experimental Result:  CMS-SUS-16-033 (upperLimit-type)
Theory prediction xsec =  1.15E-02 [pb]
Conditions violation (if any) =  {'None': None}
Dataset: None
Txnames =  ['T2']

Experimental Result:  CMS-SUS-16-033 (upperLimit-type)
Theory prediction xsec =  3.50E-05 [pb]
Conditions violation (if any) =  {'None': None}
Dataset: None
Txnames =  ['T1tttt']

Experimental Result:  CMS-SUS-16-033 (upperLimit-type)
Theory prediction xsec =  7.53E-03 [pb]
Conditions violation (if any) =  {'None': None}
Dataset: None
Txnames =  ['T1']

Experimental Result:  CMS-SUS-16-034 (upperLimit-type)
Theory prediction xsec =  3.99E-02 [pb]
Conditions violation (if any) =  {'None': None}
Dataset: None
Txnames =  ['TChiWZ']

Experimental Result:  CMS-SUS-16-035 (upperLimit-type)
Theory prediction xsec =  3.50E-05 [pb]
Conditions violation (if any) =  {'None': None}
Dataset: None
Txnames =  ['T1tttt']

Experimental Result:  CMS-SUS-16-036 (upperLimit-type)
Theory prediction xsec =  7.53E-03 [pb]
Conditions violation (if any) =  {'None': None}
Dataset: None
Txnames =  ['T1']

Experimental Result:  CMS-SUS-16-036 (upperLimit-type)
Theory prediction xsec =  3.50E-05 [pb]
Conditions violation (if any) =  {'None': None}
Dataset: None
Txnames =  ['T1tttt']

Experimental Result:  CMS-SUS-16-036 (upperLimit-type)
Theory prediction xsec =  1.15E-02 [pb]
Conditions violation (if any) =  {'None': None}
Dataset: None
Txnames =  ['T2']

Experimental Result:  CMS-SUS-16-037 (upperLimit-type)
Theory prediction xsec =  3.50E-05 [pb]
Conditions violation (if any) =  {'None': None}
Dataset: None
Txnames =  ['T1tttt']

Experimental Result:  CMS-SUS-16-039 (upperLimit-type)
Theory prediction xsec =  3.99E-02 [pb]
Conditions violation (if any) =  {'None': None}
Dataset: None
Txnames =  ['TChiWZ']

Experimental Result:  CMS-SUS-16-039 (upperLimit-type)
Theory prediction xsec =  2.07E-01 [pb]
Conditions violation (if any) =  {'None': None}
Dataset: None
Txnames =  ['TChiWH']

Experimental Result:  CMS-SUS-16-041 (upperLimit-type)
Theory prediction xsec =  3.50E-05 [pb]
Conditions violation (if any) =  {'None': None}
Dataset: None
Txnames =  ['T1tttt']

Experimental Result:  CMS-SUS-16-042 (upperLimit-type)
Theory prediction xsec =  3.50E-05 [pb]
Conditions violation (if any) =  {'None': None}
Dataset: None
Txnames =  ['T1tttt']

Experimental Result:  CMS-SUS-16-043 (upperLimit-type)
Theory prediction xsec =  2.07E-01 [pb]
Conditions violation (if any) =  {'None': None}
Dataset: None
Txnames =  ['TChiWH']

Experimental Result:  CMS-SUS-16-045 (upperLimit-type)
Theory prediction xsec =  2.07E-01 [pb]
Conditions violation (if any) =  {'None': None}
Dataset: None
Txnames =  ['TChiWH']

Experimental Result:  CMS-SUS-16-050 (upperLimit-type)
Theory prediction xsec =  3.50E-05 [pb]
Conditions violation (if any) =  {'None': None}
Dataset: None
Txnames =  ['T1tttt']

Experimental Result:  ATLAS-CONF-2013-007 (upperLimit-type)
Theory prediction xsec =  3.07E-06 [pb]
Conditions violation (if any) =  {'None': None}
Dataset: None
Txnames =  ['T1tttt']

Experimental Result:  ATLAS-CONF-2013-024 (efficiencyMap-type)
Theory prediction xsec =  1.57E-06 [pb]
Conditions violation (if any) =  {'None': None}
Dataset: SR2: MET > 300
Txnames =  ['T1bbqq', 'T1qqtt', 'T1tttt']

Experimental Result:  ATLAS-CONF-2013-037 (efficiencyMap-type)
Theory prediction xsec =  7.32E-07 [pb]
Conditions violation (if any) =  {'None': None}
Dataset: SRtN3
Txnames =  ['T1qqtt', 'T1tttt']

Experimental Result:  ATLAS-CONF-2013-047 (efficiencyMap-type)
Theory prediction xsec =  2.67E-04 [pb]
Conditions violation (if any) =  {'None': None}
Dataset: B Medium
Txnames =  ['T1', 'T1bbqq', 'T1qqtt', 'T1tttt', 'T2']

Experimental Result:  ATLAS-CONF-2013-048 (efficiencyMap-type)
Theory prediction xsec =  1.34E-07 [pb]
Conditions violation (if any) =  {'None': None}
Dataset: SR M100
Txnames =  ['T1qqtt', 'T1tttt']

Experimental Result:  ATLAS-CONF-2013-053 (efficiencyMap-type)
Theory prediction xsec =  9.67E-08 [pb]
Conditions violation (if any) =  {'None': None}
Dataset: SRA mCT350
Txnames =  ['TGQbbq']

Experimental Result:  ATLAS-CONF-2013-054 (efficiencyMap-type)
Theory prediction xsec =  9.49E-06 [pb]
Conditions violation (if any) =  {'None': None}
Dataset: 7j80 flavor 0 b-jets
Txnames =  ['T1', 'T1bbqq', 'T1qqtt', 'T1tttt', 'T2']

Experimental Result:  ATLAS-CONF-2013-061 (upperLimit-type)
Theory prediction xsec =  3.07E-06 [pb]
Conditions violation (if any) =  {'None': None}
Dataset: None
Txnames =  ['T1tttt']

Experimental Result:  ATLAS-CONF-2013-061 (efficiencyMap-type)
Theory prediction xsec =  1.22E-06 [pb]
Conditions violation (if any) =  {'None': None}
Dataset: SR-0l-4j-A
Txnames =  ['T1bbqq', 'T1qqtt', 'T1tttt', 'TGQbbq']

Experimental Result:  ATLAS-CONF-2013-062 (efficiencyMap-type)
Theory prediction xsec =  8.87E-07 [pb]
Conditions violation (if any) =  {'None': None}
Dataset: incHL6j_m
Txnames =  ['T1qqtt', 'T1tttt']

Experimental Result:  ATLAS-CONF-2013-089 (upperLimit-type)
Theory prediction xsec =  6.66E-03 [pb]
Conditions violation (if any) =  {'None': None}
Dataset: None
Txnames =  ['T6WW']

Experimental Result:  ATLAS-SUSY-2013-02 (upperLimit-type)
Theory prediction xsec =  9.07E-03 [pb]
Conditions violation (if any) =  {'None': None}
Dataset: None
Txnames =  ['T5WW']

Experimental Result:  ATLAS-SUSY-2013-02 (upperLimit-type)
Theory prediction xsec =  6.66E-03 [pb]
Conditions violation (if any) =  {'None': None}
Dataset: None
Txnames =  ['T6WW']

Experimental Result:  ATLAS-SUSY-2013-02 (upperLimit-type)
Theory prediction xsec =  6.62E-04 [pb]
Conditions violation (if any) =  {'None': None}
Dataset: None
Txnames =  ['T1']

Experimental Result:  ATLAS-SUSY-2013-02 (upperLimit-type)
Theory prediction xsec =  1.96E-03 [pb]
Conditions violation (if any) =  {'None': None}
Dataset: None
Txnames =  ['T2']

Experimental Result:  ATLAS-SUSY-2013-02 (efficiencyMap-type)
Theory prediction xsec =  6.23E-04 [pb]
Conditions violation (if any) =  {'None': None}
Dataset: SR2jt
Txnames =  ['T1', 'T2']

Experimental Result:  ATLAS-SUSY-2013-04 (upperLimit-type)
Theory prediction xsec =  3.07E-06 [pb]
Conditions violation (if any) =  {'None': None}
Dataset: None
Txnames =  ['T1tttt']

Experimental Result:  ATLAS-SUSY-2013-04 (efficiencyMap-type)
Theory prediction xsec =  1.53E-04 [pb]
Conditions violation (if any) =  {'None': None}
Dataset: GtGrid_SR_7ej80_0bjet
Txnames =  ['T1tttt', 'T5WW', 'T5ZZ', 'T6bbWW']

Experimental Result:  ATLAS-SUSY-2013-09 (upperLimit-type)
Theory prediction xsec =  3.07E-06 [pb]
Conditions violation (if any) =  {'None': None}
Dataset: None
Txnames =  ['T1tttt']

Experimental Result:  ATLAS-SUSY-2013-11 (upperLimit-type)
Theory prediction xsec =  1.85E-02 [pb]
Conditions violation (if any) =  {'None': None}
Dataset: None
Txnames =  ['TChiWZ']

Experimental Result:  ATLAS-SUSY-2013-11 (efficiencyMap-type)
Theory prediction xsec =  4.39E-05 [pb]
Conditions violation (if any) =  {'None': None}
Dataset: WWc-DF
Txnames =  ['TChiWW']

Experimental Result:  ATLAS-SUSY-2013-12 (upperLimit-type)
Theory prediction xsec =  1.85E-02 [pb]
Conditions violation (if any) =  {'None': None}
Dataset: None
Txnames =  ['TChiWZ']

Experimental Result:  ATLAS-SUSY-2013-18 (upperLimit-type)
Theory prediction xsec =  3.07E-06 [pb]
Conditions violation (if any) =  {'None': None}
Dataset: None
Txnames =  ['T1tttt']

Experimental Result:  ATLAS-SUSY-2013-18 (efficiencyMap-type)
Theory prediction xsec =  1.32E-07 [pb]
Conditions violation (if any) =  {'None': None}
Dataset: SR-1l-6j-C
Txnames =  ['T1tttt']

Experimental Result:  CMS-PAS-SUS-13-016 (upperLimit-type)
Theory prediction xsec =  3.07E-06 [pb]
Conditions violation (if any) =  {'None': None}
Dataset: None
Txnames =  ['T1tttt']

Experimental Result:  CMS-PAS-SUS-13-016 (efficiencyMap-type)
Theory prediction xsec =  2.88E-08 [pb]
Conditions violation (if any) =  {'None': None}
Dataset: sr0
Txnames =  ['T1tttt']

Experimental Result:  CMS-SUS-12-024 (upperLimit-type)
Theory prediction xsec =  3.07E-06 [pb]
Conditions violation (if any) =  {'None': None}
Dataset: None
Txnames =  ['T1tttt']

Experimental Result:  CMS-SUS-12-024 (efficiencyMap-type)
Theory prediction xsec =  3.79E-08 [pb]
Conditions violation (if any) =  {'None': None}
Dataset: MET4_HT4_nb3
Txnames =  ['T1tttt']

Experimental Result:  CMS-SUS-12-028 (upperLimit-type)
Theory prediction xsec =  3.07E-06 [pb]
Conditions violation (if any) =  {'None': None}
Dataset: None
Txnames =  ['T1tttt']

Experimental Result:  CMS-SUS-12-028 (upperLimit-type)
Theory prediction xsec =  6.62E-04 [pb]
Conditions violation (if any) =  {'None': None}
Dataset: None
Txnames =  ['T1']

Experimental Result:  CMS-SUS-12-028 (upperLimit-type)
Theory prediction xsec =  1.96E-03 [pb]
Conditions violation (if any) =  {'None': None}
Dataset: None
Txnames =  ['T2']

Experimental Result:  CMS-SUS-13-002 (upperLimit-type)
Theory prediction xsec =  3.07E-06 [pb]
Conditions violation (if any) =  {'None': None}
Dataset: None
Txnames =  ['T1tttt']

Experimental Result:  CMS-SUS-13-004 (upperLimit-type)
Theory prediction xsec =  3.07E-06 [pb]
Conditions violation (if any) =  {'None': None}
Dataset: None
Txnames =  ['T1tttt']

Experimental Result:  CMS-SUS-13-006 (upperLimit-type)
Theory prediction xsec =  9.59E-02 [pb]
Conditions violation (if any) =  {'None': None}
Dataset: None
Txnames =  ['TChiWH']

Experimental Result:  CMS-SUS-13-006 (upperLimit-type)
Theory prediction xsec =  1.85E-02 [pb]
Conditions violation (if any) =  {'None': None}
Dataset: None
Txnames =  ['TChiWZ']

Experimental Result:  CMS-SUS-13-007 (upperLimit-type)
Theory prediction xsec =  3.07E-06 [pb]
Conditions violation (if any) =  {'None': None}
Dataset: None
Txnames =  ['T1tttt']

Experimental Result:  CMS-SUS-13-007 (efficiencyMap-type)
Theory prediction xsec =  2.07E-08 [pb]
Conditions violation (if any) =  {'None': None}
Dataset: mu350Nb3
Txnames =  ['T1tttt']

Experimental Result:  CMS-SUS-13-012 (upperLimit-type)
Theory prediction xsec =  1.96E-03 [pb]
Conditions violation (if any) =  {'None': None}
Dataset: None
Txnames =  ['T2']

Experimental Result:  CMS-SUS-13-012 (upperLimit-type)
Theory prediction xsec =  3.07E-06 [pb]
Conditions violation (if any) =  {'None': None}
Dataset: None
Txnames =  ['T1tttt']

Experimental Result:  CMS-SUS-13-012 (upperLimit-type)
Theory prediction xsec =  6.62E-04 [pb]
Conditions violation (if any) =  {'None': None}
Dataset: None
Txnames =  ['T1']

Experimental Result:  CMS-SUS-13-012 (efficiencyMap-type)
Theory prediction xsec =  8.71E-05 [pb]
Conditions violation (if any) =  {'None': None}
Dataset: 6NJet8_1250HT1500_450MHTinf
Txnames =  ['T1', 'T1tttt', 'T2', 'T5WW', 'T5ZZ', 'T6bbWW', 'TChiWW', 'TChiWZ', 'TChiZZ']

Experimental Result:  CMS-SUS-13-013 (upperLimit-type)
Theory prediction xsec =  3.07E-06 [pb]
Conditions violation (if any) =  {'None': None}
Dataset: None
Txnames =  ['T1tttt']

Experimental Result:  CMS-SUS-13-013 (efficiencyMap-type)
Theory prediction xsec =  4.16E-08 [pb]
Conditions violation (if any) =  {'None': None}
Dataset: SR28_HighPt
Txnames =  ['T1tttt']

Experimental Result:  CMS-SUS-13-019 (upperLimit-type)
Theory prediction xsec =  3.07E-06 [pb]
Conditions violation (if any) =  {'None': None}
Dataset: None
Txnames =  ['T1tttt']

Experimental Result:  CMS-SUS-13-019 (upperLimit-type)
Theory prediction xsec =  6.62E-04 [pb]
Conditions violation (if any) =  {'None': None}
Dataset: None
Txnames =  ['T1']

Experimental Result:  CMS-SUS-13-019 (upperLimit-type)
Theory prediction xsec =  1.96E-03 [pb]
Conditions violation (if any) =  {'None': None}
Dataset: None
Txnames =  ['T2']

Experimental Result:  CMS-SUS-14-010 (upperLimit-type)
Theory prediction xsec =  3.07E-06 [pb]
Conditions violation (if any) =  {'None': None}
Dataset: None
Txnames =  ['T1tttt']