How To: Print out the theory predictions

In [1]:
#Set up the path to SModelS installation folder if running on a different folder
import sys,os
sys.path.append(os.path.join(os.getenv("HOME"),"smodels/"))
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 theoryPredictionFor
from smodels.experiment import smsAnalysisFactory, smsHelpers
In [3]:
## define where the database resides
smsHelpers.base=os.path.join(os.getenv("HOME"),"smodels-database/")
#and load the analyses:
listofanalyses = smsAnalysisFactory.load()
12:22:48.117 INFO     smodels.experiment.smsAnalysisFactory:39  useSuperseded is not set, skipping superseded results

In [4]:
#Define the SLHA input file name
filename="%s/inputFiles/slha/gluino_squarks.slha" % installDirectory()
In [5]:
#Perform the decomposition:
listOfTopologies = slhaDecomposer.decompose (filename, sigcut = 0.03 * fb, doCompress=True, doInvisible=True,minmassgap = 5* GeV)
12:22:51.690 INFO     smodels.theory.slhaDecomposer:132 Ignoring t+ decays
12:22:51.695 INFO     smodels.theory.slhaDecomposer:132 Ignoring higgs decays
12:22:51.696 INFO     smodels.theory.slhaDecomposer:132 Ignoring H0 decays
12:22:51.696 INFO     smodels.theory.slhaDecomposer:132 Ignoring A0 decays
12:22:51.696 INFO     smodels.theory.slhaDecomposer:132 Ignoring H+ decays
12:22:51.730 INFO     smodels.theory.crossSection:512 Ignoring 76 lower order cross-sections

In [6]:
#Compute the theory prediction for each analysis using the results from the decomposition:
analysesPredictions = [theoryPredictionFor(analysis, listOfTopologies) for analysis in listofanalyses]
12:23:05.784 INFO     smodels.experiment.smsInterpolation:174 Masses out of range for ATLAS-SUSY-2013-19/T6bbWW (no extrapolation)
12:23:05.860 INFO     smodels.experiment.smsInterpolation:174 Masses out of range for CMS-SUS-13-011/T6bbWW (no extrapolation)
12:23:06.338 INFO     smodels.experiment.smsInterpolation:174 Masses out of range for CMS-PAS-SUS-13-008/T6ttWW (no extrapolation)

In [7]:
#Print information about each theory prediction (cluster) for each analysis:
#(Since this is a test database, there are very few results)
for anaPrediction in analysesPredictions:
    if not anaPrediction: continue #skip analyses without results
    for theoryPred in anaPrediction:
        print "Analysis name = ",theoryPred.analysis.label
        print "Theory prediction = ",theoryPred.value
        print "Conditions violation (if any) = ",theoryPred.conditions
        print "Mass of cluster = ",theoryPred.mass
Analysis name =  CMS-SUS-12-028:T2
Theory prediction =  ['8.00E+00 [TeV]:1.77E-03 [pb]']
Conditions violation (if any) =  {'None': None}
Mass of cluster =  [[9.91E+02 [GeV], 1.29E+02 [GeV]], [9.91E+02 [GeV], 1.29E+02 [GeV]]]
Analysis name =  CMS-SUS-12-028:T1
Theory prediction =  ['8.00E+00 [TeV]:3.92E-04 [pb]']
Conditions violation (if any) =  {'None': None}
Mass of cluster =  [[8.65E+02 [GeV], 1.29E+02 [GeV]], [8.65E+02 [GeV], 1.29E+02 [GeV]]]
Analysis name =  ATLAS-SUSY-2013-12:TChiWZ
Theory prediction =  ['8.00E+00 [TeV]:1.85E-02 [pb]']
Conditions violation (if any) =  {'None': None}
Mass of cluster =  [[2.69E+02 [GeV], 1.29E+02 [GeV]], [2.69E+02 [GeV], 1.29E+02 [GeV]]]
Analysis name =  ATLAS-SUSY-2013-11:TChiWZ
Theory prediction =  ['8.00E+00 [TeV]:1.85E-02 [pb]']
Conditions violation (if any) =  {'None': None}
Mass of cluster =  [[2.69E+02 [GeV], 1.29E+02 [GeV]], [2.69E+02 [GeV], 1.29E+02 [GeV]]]
Analysis name =  CMS-SUS-13-012:T2
Theory prediction =  ['8.00E+00 [TeV]:1.77E-03 [pb]']
Conditions violation (if any) =  {'None': None}
Mass of cluster =  [[9.91E+02 [GeV], 1.29E+02 [GeV]], [9.91E+02 [GeV], 1.29E+02 [GeV]]]
Analysis name =  CMS-SUS-13-012:T1
Theory prediction =  ['8.00E+00 [TeV]:3.92E-04 [pb]']
Conditions violation (if any) =  {'None': None}
Mass of cluster =  [[8.65E+02 [GeV], 1.29E+02 [GeV]], [8.65E+02 [GeV], 1.29E+02 [GeV]]]
Analysis name =  CMS-PAS-SUS-13-019:T2
Theory prediction =  ['8.00E+00 [TeV]:1.77E-03 [pb]']
Conditions violation (if any) =  {'None': None}
Mass of cluster =  [[9.91E+02 [GeV], 1.29E+02 [GeV]], [9.91E+02 [GeV], 1.29E+02 [GeV]]]
Analysis name =  CMS-PAS-SUS-13-019:T1
Theory prediction =  ['8.00E+00 [TeV]:3.92E-04 [pb]']
Conditions violation (if any) =  {'None': None}
Mass of cluster =  [[8.65E+02 [GeV], 1.29E+02 [GeV]], [8.65E+02 [GeV], 1.29E+02 [GeV]]]
Analysis name =  ATLAS-SUSY-2013-02:T2
Theory prediction =  ['8.00E+00 [TeV]:1.77E-03 [pb]']
Conditions violation (if any) =  {'None': None}
Mass of cluster =  [[9.91E+02 [GeV], 1.29E+02 [GeV]], [9.91E+02 [GeV], 1.29E+02 [GeV]]]
Analysis name =  ATLAS-SUSY-2013-02:T1
Theory prediction =  ['8.00E+00 [TeV]:3.92E-04 [pb]']
Conditions violation (if any) =  {'None': None}
Mass of cluster =  [[8.65E+02 [GeV], 1.29E+02 [GeV]], [8.65E+02 [GeV], 1.29E+02 [GeV]]]
Analysis name =  CMS-SUS-13-006:TChiWZ
Theory prediction =  ['8.00E+00 [TeV]:1.85E-02 [pb]']
Conditions violation (if any) =  {'None': None}
Mass of cluster =  [[2.69E+02 [GeV], 1.29E+02 [GeV]], [2.69E+02 [GeV], 1.29E+02 [GeV]]]

In []: