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.share.models.mssm import BSMList
from smodels.share.models.SMparticles import SMList
from smodels.theory.model import Model
from smodels.theory import decomposer
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="inputFiles/slha/gluino_squarks.slha"
model = Model(BSMparticles = BSMList, SMparticles = SMList)
model.updateParticles(inputFile=filename)
In [5]:
# Perform the decomposition:
topList = decomposer.decompose(model, sigmacut = 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.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.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 = 4.21E-03 [pb] Conditions violation (if any) = None Dataset: SR2jm Txnames = ['T1', 'T2'] Experimental Result: ATLAS-SUSY-2016-07 (upperLimit-type) Theory prediction xsec = 1.03E-01 [pb] Conditions violation (if any) = None Dataset: None Txnames = ['T5WW'] Experimental Result: ATLAS-SUSY-2016-07 (upperLimit-type) Theory prediction xsec = 4.91E-02 [pb] Conditions violation (if any) = None Dataset: None Txnames = ['T6WW'] Experimental Result: ATLAS-SUSY-2016-07 (upperLimit-type) Theory prediction xsec = 1.16E-02 [pb] Conditions violation (if any) = None Dataset: None Txnames = ['T2'] Experimental Result: ATLAS-SUSY-2016-07 (upperLimit-type) Theory prediction xsec = 7.53E-03 [pb] Conditions violation (if any) = None Dataset: None Txnames = ['T1'] Experimental Result: ATLAS-SUSY-2016-07 (efficiencyMap-type) Theory prediction xsec = 9.72E-03 [pb] Conditions violation (if any) = None Dataset: 5j_Meff_1600 Txnames = ['T5GQ', 'T6WW', 'T5WW', 'T1', 'T2', 'T5ZZ', 'TGQ'] Experimental Result: ATLAS-SUSY-2016-14 (upperLimit-type) Theory prediction xsec = 3.50E-05 [pb] Conditions violation (if any) = None Dataset: None Txnames = ['T1tttt'] Experimental Result: ATLAS-SUSY-2016-24 (upperLimit-type) Theory prediction xsec = 3.99E-02 [pb] Conditions violation (if any) = None Dataset: None Txnames = ['TChiWZ'] Experimental Result: ATLAS-SUSY-2016-24 (efficiencyMap-type) Theory prediction xsec = 7.64E-06 [pb] Conditions violation (if any) = None Dataset: WZ-0Jb Txnames = ['TChiWZ'] Experimental Result: ATLAS-SUSY-2016-26 (upperLimit-type) Theory prediction xsec = 2.49E-04 [pb] Conditions violation (if any) = None Dataset: None Txnames = ['T2cc'] Experimental Result: ATLAS-SUSY-2017-03 (upperLimit-type) Theory prediction xsec = 3.99E-02 [pb] Conditions violation (if any) = None Dataset: None Txnames = ['TChiWZ'] Experimental Result: ATLAS-SUSY-2018-06 (upperLimit-type) Theory prediction xsec = 3.99E-02 [pb] Conditions violation (if any) = None Dataset: None Txnames = ['TChiWZ'] Experimental Result: ATLAS-SUSY-2018-06 (efficiencyMap-type) Theory prediction xsec = 1.26E-05 [pb] Conditions violation (if any) = None Dataset: SR_low Txnames = ['TChiWZ'] Experimental Result: ATLAS-SUSY-2018-10 (upperLimit-type) Theory prediction xsec = 2.46E-02 [pb] Conditions violation (if any) = None Dataset: None Txnames = ['T6WW'] Experimental Result: ATLAS-SUSY-2018-10 (efficiencyMap-type) Theory prediction xsec = 2.47E-04 [pb] Conditions violation (if any) = None Dataset: 4j0blowx_3 Txnames = ['T6WW'] Experimental Result: ATLAS-SUSY-2018-22 (upperLimit-type) Theory prediction xsec = 1.16E-02 [pb] Conditions violation (if any) = None Dataset: None Txnames = ['T2'] Experimental Result: ATLAS-SUSY-2018-22 (upperLimit-type) Theory prediction xsec = 2.46E-02 [pb] Conditions violation (if any) = None Dataset: None Txnames = ['T6WW'] Experimental Result: ATLAS-SUSY-2018-22 (efficiencyMap-type) Theory prediction xsec = 2.28E-03 [pb] Conditions violation (if any) = None Dataset: SR2j_1600 Txnames = ['T2'] Experimental Result: ATLAS-SUSY-2018-23 (upperLimit-type) Theory prediction xsec = 2.07E-01 [pb] Conditions violation (if any) = None Dataset: None Txnames = ['TChiWH'] Experimental Result: ATLAS-SUSY-2018-31 (upperLimit-type) Theory prediction xsec = 2.87E-04 [pb] Conditions violation (if any) = None Dataset: None Txnames = ['T6bbHH'] Experimental Result: ATLAS-SUSY-2018-31 (efficiencyMap-type) Theory prediction xsec = 2.05E-06 [pb] Conditions violation (if any) = None Dataset: SRA_M Txnames = ['T6bbHH'] Experimental Result: ATLAS-SUSY-2018-31 (combined-type) Theory prediction xsec = 4.14E-06 [pb] Conditions violation (if any) = None Dataset: (combined) Txnames = ['T6bbHH', 'T6bbHH', 'T6bbHH', 'T6bbHH'] Experimental Result: ATLAS-SUSY-2018-32 (upperLimit-type) Theory prediction xsec = 1.45E-01 [pb] Conditions violation (if any) = None Dataset: None Txnames = ['TChiWW'] Experimental Result: ATLAS-SUSY-2019-08 (upperLimit-type) Theory prediction xsec = 2.07E-01 [pb] Conditions violation (if any) = None Dataset: None Txnames = ['TChiWH'] Experimental Result: ATLAS-SUSY-2019-08 (efficiencyMap-type) Theory prediction xsec = 9.06E-06 [pb] Conditions violation (if any) = None Dataset: SR_LM_High_MCT Txnames = ['TChiWH'] Experimental Result: ATLAS-SUSY-2019-08 (combined-type) Theory prediction xsec = 2.94E-05 [pb] Conditions violation (if any) = None Dataset: (combined) Txnames = ['TChiWH', 'TChiWH', 'TChiWH', 'TChiWH', 'TChiWH', 'TChiWH', 'TChiWH', 'TChiWH'] Experimental Result: CMS-SUS-16-009 (upperLimit-type) Theory prediction xsec = 3.50E-05 [pb] Conditions violation (if any) = 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 Dataset: None Txnames = ['T1'] Experimental Result: CMS-SUS-16-033 (upperLimit-type) Theory prediction xsec = 3.50E-05 [pb] Conditions violation (if any) = None Dataset: None Txnames = ['T1tttt'] Experimental Result: CMS-SUS-16-033 (upperLimit-type) Theory prediction xsec = 1.16E-02 [pb] Conditions violation (if any) = None Dataset: None Txnames = ['T2'] Experimental Result: CMS-SUS-16-033 (efficiencyMap-type) Theory prediction xsec = 6.19E-03 [pb] Conditions violation (if any) = None Dataset: SR3_Njet5_Nb0_HT500_MHT_500 Txnames = ['T1', 'T2', 'T1tttt', 'T5GQ', 'TGQ'] Experimental Result: CMS-SUS-16-034 (upperLimit-type) Theory prediction xsec = 3.99E-02 [pb] Conditions violation (if any) = 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 Dataset: None Txnames = ['T1tttt'] Experimental Result: CMS-SUS-16-036 (upperLimit-type) Theory prediction xsec = 1.16E-02 [pb] Conditions violation (if any) = None Dataset: None Txnames = ['T2'] Experimental Result: CMS-SUS-16-036 (upperLimit-type) Theory prediction xsec = 7.53E-03 [pb] Conditions violation (if any) = 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 Dataset: None Txnames = ['T1tttt'] Experimental Result: CMS-SUS-16-037 (upperLimit-type) Theory prediction xsec = 3.50E-05 [pb] Conditions violation (if any) = 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 Dataset: None Txnames = ['TChiWZ'] Experimental Result: CMS-SUS-16-039 (upperLimit-type) Theory prediction xsec = 2.07E-01 [pb] Conditions violation (if any) = 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 Dataset: None Txnames = ['T1tttt'] Experimental Result: CMS-SUS-16-042 (upperLimit-type) Theory prediction xsec = 3.50E-05 [pb] Conditions violation (if any) = 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 Dataset: None Txnames = ['TChiWH'] Experimental Result: CMS-SUS-16-045 (upperLimit-type) Theory prediction xsec = 2.07E-01 [pb] Conditions violation (if any) = 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 Dataset: None Txnames = ['T1tttt'] Experimental Result: CMS-SUS-17-004 (upperLimit-type) Theory prediction xsec = 2.07E-01 [pb] Conditions violation (if any) = None Dataset: None Txnames = ['TChiWH'] Experimental Result: CMS-SUS-17-004 (upperLimit-type) Theory prediction xsec = 3.99E-02 [pb] Conditions violation (if any) = None Dataset: None Txnames = ['TChiWZ'] Experimental Result: CMS-SUS-19-006 (upperLimit-type) Theory prediction xsec = 3.50E-05 [pb] Conditions violation (if any) = None Dataset: None Txnames = ['T1tttt'] Experimental Result: CMS-SUS-19-006 (upperLimit-type) Theory prediction xsec = 7.53E-03 [pb] Conditions violation (if any) = None Dataset: None Txnames = ['T1'] Experimental Result: CMS-SUS-19-006 (upperLimit-type) Theory prediction xsec = 1.16E-02 [pb] Conditions violation (if any) = None Dataset: None Txnames = ['T2'] Experimental Result: ATLAS-SUSY-2013-02 (upperLimit-type) Theory prediction xsec = 6.74E-03 [pb] Conditions violation (if any) = None Dataset: None Txnames = ['T6WW'] Experimental Result: ATLAS-SUSY-2013-02 (upperLimit-type) Theory prediction xsec = 2.11E-03 [pb] Conditions violation (if any) = None Dataset: None Txnames = ['T2'] Experimental Result: ATLAS-SUSY-2013-02 (upperLimit-type) Theory prediction xsec = 9.07E-03 [pb] Conditions violation (if any) = None Dataset: None Txnames = ['T5WW'] Experimental Result: ATLAS-SUSY-2013-02 (upperLimit-type) Theory prediction xsec = 6.62E-04 [pb] Conditions violation (if any) = None Dataset: None Txnames = ['T1'] Experimental Result: ATLAS-SUSY-2013-02 (efficiencyMap-type) Theory prediction xsec = 7.86E-04 [pb] Conditions violation (if any) = 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 Dataset: None Txnames = ['T1tttt'] Experimental Result: ATLAS-SUSY-2013-04 (efficiencyMap-type) Theory prediction xsec = 1.53E-04 [pb] Conditions violation (if any) = None Dataset: GtGrid_SR_7ej80_0bjet Txnames = ['T5WW', 'T6bbWW', 'T1tttt', 'T5ZZ'] Experimental Result: ATLAS-SUSY-2013-09 (upperLimit-type) Theory prediction xsec = 9.07E-03 [pb] Conditions violation (if any) = None Dataset: None Txnames = ['T5WW'] Experimental Result: ATLAS-SUSY-2013-09 (upperLimit-type) Theory prediction xsec = 3.07E-06 [pb] Conditions violation (if any) = None Dataset: None Txnames = ['T1tttt'] Experimental Result: ATLAS-SUSY-2013-09 (efficiencyMap-type) Theory prediction xsec = 5.00E-05 [pb] Conditions violation (if any) = None Dataset: SR0B Txnames = ['T1tttt', 'T5WW'] Experimental Result: ATLAS-SUSY-2013-11 (upperLimit-type) Theory prediction xsec = 1.85E-02 [pb] Conditions violation (if any) = 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 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 Dataset: None Txnames = ['TChiWZ'] Experimental Result: ATLAS-SUSY-2013-18 (upperLimit-type) Theory prediction xsec = 3.07E-06 [pb] Conditions violation (if any) = None Dataset: None Txnames = ['T1tttt'] Experimental Result: ATLAS-SUSY-2013-18 (efficiencyMap-type) Theory prediction xsec = 1.99E-07 [pb] Conditions violation (if any) = None Dataset: SR-1l-6j-B Txnames = ['T1tttt'] Experimental Result: ATLAS-SUSY-2013-20 (upperLimit-type) Theory prediction xsec = 6.74E-03 [pb] Conditions violation (if any) = None Dataset: None Txnames = ['T6WW'] Experimental Result: CMS-PAS-SUS-13-016 (upperLimit-type) Theory prediction xsec = 3.07E-06 [pb] Conditions violation (if any) = 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 Dataset: sr0 Txnames = ['T1tttt'] Experimental Result: CMS-SUS-12-024 (upperLimit-type) Theory prediction xsec = 3.07E-06 [pb] Conditions violation (if any) = None Dataset: None Txnames = ['T1tttt'] Experimental Result: CMS-SUS-12-024 (efficiencyMap-type) Theory prediction xsec = 3.77E-08 [pb] Conditions violation (if any) = None Dataset: MET4_HT4_nb3 Txnames = ['T1tttt'] Experimental Result: CMS-SUS-12-028 (upperLimit-type) Theory prediction xsec = 6.62E-04 [pb] Conditions violation (if any) = None Dataset: None Txnames = ['T1'] Experimental Result: CMS-SUS-12-028 (upperLimit-type) Theory prediction xsec = 3.07E-06 [pb] Conditions violation (if any) = None Dataset: None Txnames = ['T1tttt'] Experimental Result: CMS-SUS-12-028 (upperLimit-type) Theory prediction xsec = 2.11E-03 [pb] Conditions violation (if any) = 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 Dataset: None Txnames = ['T1tttt'] Experimental Result: CMS-SUS-13-004 (upperLimit-type) Theory prediction xsec = 3.07E-06 [pb] Conditions violation (if any) = None Dataset: None Txnames = ['T1tttt'] Experimental Result: CMS-SUS-13-006 (upperLimit-type) Theory prediction xsec = 1.85E-02 [pb] Conditions violation (if any) = None Dataset: None Txnames = ['TChiWZ'] Experimental Result: CMS-SUS-13-006 (upperLimit-type) Theory prediction xsec = 9.59E-02 [pb] Conditions violation (if any) = None Dataset: None Txnames = ['TChiWH'] Experimental Result: CMS-SUS-13-007 (upperLimit-type) Theory prediction xsec = 3.07E-06 [pb] Conditions violation (if any) = 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 Dataset: mu350Nb3 Txnames = ['T1tttt'] Experimental Result: CMS-SUS-13-012 (upperLimit-type) Theory prediction xsec = 3.07E-06 [pb] Conditions violation (if any) = 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 Dataset: None Txnames = ['T1'] Experimental Result: CMS-SUS-13-012 (upperLimit-type) Theory prediction xsec = 2.11E-03 [pb] Conditions violation (if any) = None Dataset: None Txnames = ['T2'] Experimental Result: CMS-SUS-13-012 (efficiencyMap-type) Theory prediction xsec = 3.45E-04 [pb] Conditions violation (if any) = None Dataset: 6NJet8_1500HTinf_300MHTinf Txnames = ['T6bbWW', 'TChiWW', 'TChiZZ', 'T5WW', 'T1', 'T2', 'T5ZZ', 'T6WW', 'T1tttt'] Experimental Result: CMS-SUS-13-013 (upperLimit-type) Theory prediction xsec = 3.07E-06 [pb] Conditions violation (if any) = None Dataset: None Txnames = ['T1tttt'] Experimental Result: CMS-SUS-13-013 (efficiencyMap-type) Theory prediction xsec = 4.17E-08 [pb] Conditions violation (if any) = None Dataset: SR28_HighPt Txnames = ['T1tttt'] Experimental Result: CMS-SUS-13-019 (upperLimit-type) Theory prediction xsec = 6.62E-04 [pb] Conditions violation (if any) = None Dataset: None Txnames = ['T1'] Experimental Result: CMS-SUS-13-019 (upperLimit-type) Theory prediction xsec = 2.11E-03 [pb] Conditions violation (if any) = None Dataset: None Txnames = ['T2'] Experimental Result: CMS-SUS-13-019 (upperLimit-type) Theory prediction xsec = 3.07E-06 [pb] Conditions violation (if any) = None Dataset: None Txnames = ['T1tttt'] Experimental Result: CMS-SUS-14-010 (upperLimit-type) Theory prediction xsec = 3.07E-06 [pb] Conditions violation (if any) = None Dataset: None Txnames = ['T1tttt']
In [ ]: