# Set up the path to SModelS installation folder
import sys; sys.path.append("."); import smodels_paths
# Import those parts of smodels that are needed for this exercise
from smodels.tools.physicsUnits import TeV, GeV, fb
from smodels.installation import installDirectory
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.theory.theoryPrediction import theoryPredictionsFor
from smodels.experiment.databaseObj import Database
from smodels.tools import coverage
# Define the SLHA input file name
filename="inputFiles/slha/gluino_squarks.slha"
model = Model(BSMparticles = BSMList, SMparticles = SMList)
model.updateParticles(inputFile=filename)
# Load the database, do the decomposition and compute theory predictions:
# (Look at the theory predictions HowTo to learn how to compute theory predictions)
database = Database("official")
expResults = database.getExpResults()
topList = decomposer.decompose(model, sigmacut = 0.03 * fb, doCompress=True, doInvisible=True,minmassgap = 5* GeV)
#(Computing theory predictions will tag the elements which have been tested)
allThPredictions = [theoryPredictionsFor(exp, topList) for exp in expResults]
#Find out missing topologies:
uncovered = coverage.Uncovered(topList)
#First sort coverage groups by label
groups = sorted(uncovered.groups[:], key = lambda g: g.label)
#Print uncovered cross-sections:
for group in groups:
print("\nTotal cross-section for %s (fb): %10.3E" %(group.description,group.getTotalXSec()))
Total cross-section for missing topologies (fb): 2.723E+03 Total cross-section for missing topologies with displaced decays (fb): 0.000E+00 Total cross-section for missing topologies with prompt decays (fb): 2.739E+03 Total cross-section for topologies outside the grid (fb): 2.098E+02
# Get list of topologies which are not tested by any result:
missingTopos = uncovered.getGroup('missing (all)')
# print information about the first few missing topologies and
# the elements contributing to the topology:
for genEl in missingTopos.generalElements[:5]:
print ('\nmissing topology:',str(genEl._outputDescription))
print ('Contributing elements:')
for el in genEl._contributingElements:
print('\t',el)
missing topology: [[[jet],[W]],[[jet,jet],[W]]] (MET,MET) Contributing elements: [[[q],[W+]],[[q,q],[W-]]] [[[q],[W+]],[[q,q],[W+]]] [[[q],[W+]],[[c,q],[W-]]] [[[q],[W+]],[[c,q],[W+]]] [[[q],[W-]],[[q,q],[W-]]] [[[q],[W-]],[[q,q],[W+]]] [[[q],[W-]],[[c,q],[W-]]] [[[q],[W-]],[[c,q],[W+]]] [[[q],[W-]],[[q,q],[W-]]] [[[q],[W-]],[[q,q],[W+]]] [[[q],[W-]],[[c,q],[W-]]] [[[q],[W-]],[[c,q],[W+]]] [[[q],[W+]],[[q,q],[W-]]] [[[q],[W+]],[[q,q],[W+]]] [[[q],[W+]],[[c,q],[W-]]] [[[q],[W+]],[[c,q],[W+]]] [[[c],[W-]],[[q,q],[W-]]] [[[c],[W-]],[[q,q],[W+]]] [[[c],[W-]],[[c,q],[W-]]] [[[c],[W-]],[[c,q],[W+]]] [[[c],[W+]],[[q,q],[W-]]] [[[c],[W+]],[[q,q],[W+]]] [[[c],[W+]],[[c,q],[W-]]] [[[c],[W+]],[[c,q],[W+]]] missing topology: [[[jet,jet],[W]],[[jet],[jet,jet],[W]]] (MET,MET) Contributing elements: [[[q,q],[W-]],[[q],[q,q],[W-]]] [[[q,q],[W-]],[[q],[q,q],[W+]]] [[[q,q],[W-]],[[q],[c,q],[W-]]] [[[q,q],[W-]],[[q],[c,q],[W+]]] [[[q,q],[W+]],[[q],[q,q],[W-]]] [[[q,q],[W+]],[[q],[q,q],[W+]]] [[[q,q],[W+]],[[q],[c,q],[W-]]] [[[q,q],[W+]],[[q],[c,q],[W+]]] [[[c,q],[W-]],[[q],[q,q],[W-]]] [[[c,q],[W-]],[[q],[q,q],[W+]]] [[[c,q],[W-]],[[q],[c,q],[W-]]] [[[c,q],[W-]],[[q],[c,q],[W+]]] [[[c,q],[W+]],[[q],[q,q],[W-]]] [[[c,q],[W+]],[[q],[q,q],[W+]]] [[[c,q],[W+]],[[q],[c,q],[W-]]] [[[c,q],[W+]],[[q],[c,q],[W+]]] [[[q,q],[W-]],[[q],[q,q],[W-]]] [[[q,q],[W-]],[[q],[q,q],[W+]]] [[[q,q],[W-]],[[q],[c,q],[W-]]] [[[q,q],[W-]],[[q],[c,q],[W+]]] [[[q,q],[W+]],[[q],[q,q],[W-]]] [[[q,q],[W+]],[[q],[q,q],[W+]]] [[[q,q],[W+]],[[q],[c,q],[W-]]] [[[q,q],[W+]],[[q],[c,q],[W+]]] [[[c,q],[W-]],[[q],[q,q],[W-]]] [[[c,q],[W-]],[[q],[q,q],[W+]]] [[[c,q],[W-]],[[q],[c,q],[W-]]] [[[c,q],[W-]],[[q],[c,q],[W+]]] [[[c,q],[W+]],[[q],[q,q],[W-]]] [[[c,q],[W+]],[[q],[q,q],[W+]]] [[[c,q],[W+]],[[q],[c,q],[W-]]] [[[c,q],[W+]],[[q],[c,q],[W+]]] [[[q,q],[W-]],[[q],[q,q],[W-]]] [[[q,q],[W-]],[[q],[q,q],[W+]]] [[[q,q],[W-]],[[q],[c,q],[W-]]] [[[q,q],[W-]],[[q],[c,q],[W+]]] [[[q,q],[W+]],[[q],[q,q],[W-]]] [[[q,q],[W+]],[[q],[q,q],[W+]]] [[[q,q],[W+]],[[q],[c,q],[W-]]] [[[q,q],[W+]],[[q],[c,q],[W+]]] [[[c,q],[W-]],[[q],[q,q],[W-]]] [[[c,q],[W-]],[[q],[q,q],[W+]]] [[[c,q],[W-]],[[q],[c,q],[W-]]] [[[c,q],[W-]],[[q],[c,q],[W+]]] [[[c,q],[W+]],[[q],[q,q],[W-]]] [[[c,q],[W+]],[[q],[q,q],[W+]]] [[[c,q],[W+]],[[q],[c,q],[W-]]] [[[c,q],[W+]],[[q],[c,q],[W+]]] [[[q,q],[W-]],[[q],[q,q],[W-]]] [[[q,q],[W-]],[[q],[q,q],[W+]]] [[[q,q],[W-]],[[q],[c,q],[W-]]] [[[q,q],[W-]],[[q],[c,q],[W+]]] [[[q,q],[W+]],[[q],[q,q],[W-]]] [[[q,q],[W+]],[[q],[q,q],[W+]]] [[[q,q],[W+]],[[q],[c,q],[W-]]] [[[q,q],[W+]],[[q],[c,q],[W+]]] [[[c,q],[W-]],[[q],[q,q],[W-]]] [[[c,q],[W-]],[[q],[q,q],[W+]]] [[[c,q],[W-]],[[q],[c,q],[W-]]] [[[c,q],[W-]],[[q],[c,q],[W+]]] [[[c,q],[W+]],[[q],[q,q],[W-]]] [[[c,q],[W+]],[[q],[q,q],[W+]]] [[[c,q],[W+]],[[q],[c,q],[W-]]] [[[c,q],[W+]],[[q],[c,q],[W+]]] [[[q,q],[W-]],[[c],[q,q],[W-]]] [[[q,q],[W-]],[[c],[q,q],[W+]]] [[[q,q],[W-]],[[c],[c,q],[W-]]] [[[q,q],[W-]],[[c],[c,q],[W+]]] [[[q,q],[W+]],[[c],[q,q],[W-]]] [[[q,q],[W+]],[[c],[q,q],[W+]]] [[[q,q],[W+]],[[c],[c,q],[W-]]] [[[q,q],[W+]],[[c],[c,q],[W+]]] [[[c,q],[W-]],[[c],[q,q],[W-]]] [[[c,q],[W-]],[[c],[q,q],[W+]]] [[[c,q],[W-]],[[c],[c,q],[W-]]] [[[c,q],[W-]],[[c],[c,q],[W+]]] [[[c,q],[W+]],[[c],[q,q],[W-]]] [[[c,q],[W+]],[[c],[q,q],[W+]]] [[[c,q],[W+]],[[c],[c,q],[W-]]] [[[c,q],[W+]],[[c],[c,q],[W+]]] missing topology: [[[jet,jet],[W]],[[b,t],[W]]] (MET,MET) Contributing elements: [[[q,q],[W-]],[[b,t-],[W+]]] [[[q,q],[W-]],[[b,t+],[W-]]] [[[q,q],[W+]],[[b,t-],[W+]]] [[[q,q],[W+]],[[b,t+],[W-]]] [[[c,q],[W-]],[[b,t-],[W+]]] [[[c,q],[W-]],[[b,t+],[W-]]] [[[c,q],[W+]],[[b,t-],[W+]]] [[[c,q],[W+]],[[b,t+],[W-]]] missing topology: [[[jet]],[[jet,jet],[W]]] (MET,MET) Contributing elements: [[[q]],[[q,q],[W-]]] [[[q]],[[q,q],[W+]]] [[[q]],[[c,q],[W-]]] [[[q]],[[c,q],[W+]]] [[[q]],[[q,q],[W-]]] [[[q]],[[q,q],[W+]]] [[[q]],[[c,q],[W-]]] [[[q]],[[c,q],[W+]]] [[[q]],[[q,q],[W-]]] [[[q]],[[q,q],[W+]]] [[[q]],[[c,q],[W-]]] [[[q]],[[c,q],[W+]]] [[[c]],[[q,q],[W-]]] [[[c]],[[q,q],[W+]]] [[[c]],[[c,q],[W-]]] [[[c]],[[c,q],[W+]]] [[[q]],[[q,q],[W-]]] [[[q]],[[q,q],[W+]]] [[[q]],[[c,q],[W-]]] [[[q]],[[c,q],[W+]]] missing topology: [[[jet],[W]],[[b,t],[W]]] (MET,MET) Contributing elements: [[[q],[W+]],[[b,t-],[W+]]] [[[q],[W+]],[[b,t+],[W-]]] [[[q],[W-]],[[b,t-],[W+]]] [[[q],[W-]],[[b,t+],[W-]]] [[[q],[W-]],[[b,t-],[W+]]] [[[q],[W-]],[[b,t+],[W-]]] [[[q],[W+]],[[b,t-],[W+]]] [[[q],[W+]],[[b,t+],[W-]]] [[[c],[W-]],[[b,t-],[W+]]] [[[c],[W-]],[[b,t+],[W-]]] [[[c],[W+]],[[b,t-],[W+]]] [[[c],[W+]],[[b,t+],[W-]]]
# Get list of topologies which are not tested by any result:
outsideGrid = uncovered.getGroup('outsideGrid (all)')
# print information about the first few missing topologies and
# the elements contributing to the topology:
for genEl in outsideGrid.generalElements[:5]:
print ('\nmissing topology:',str(genEl._outputDescription))
print ('Contributing elements:')
for el in genEl._contributingElements:
print('\t',el)
print('\t element mass:',el.mass)
missing topology: [[[jet,jet],[higgs]],[[jet,jet],[W]]] (MET,MET) Contributing elements: [[[q,q],[higgs]],[[q,q],[W-]]] element mass: [[8.65E+02 [GeV], 2.69E+02 [GeV], 1.29E+02 [GeV]], [8.65E+02 [GeV], 2.69E+02 [GeV], 1.29E+02 [GeV]]] [[[q,q],[higgs]],[[q,q],[W+]]] element mass: [[8.65E+02 [GeV], 2.69E+02 [GeV], 1.29E+02 [GeV]], [8.65E+02 [GeV], 2.69E+02 [GeV], 1.29E+02 [GeV]]] [[[q,q],[higgs]],[[c,q],[W-]]] element mass: [[8.65E+02 [GeV], 2.69E+02 [GeV], 1.29E+02 [GeV]], [8.65E+02 [GeV], 2.69E+02 [GeV], 1.29E+02 [GeV]]] [[[q,q],[higgs]],[[c,q],[W+]]] element mass: [[8.65E+02 [GeV], 2.69E+02 [GeV], 1.29E+02 [GeV]], [8.65E+02 [GeV], 2.69E+02 [GeV], 1.29E+02 [GeV]]] [[[c,c],[higgs]],[[q,q],[W-]]] element mass: [[8.65E+02 [GeV], 2.69E+02 [GeV], 1.29E+02 [GeV]], [8.65E+02 [GeV], 2.69E+02 [GeV], 1.29E+02 [GeV]]] [[[c,c],[higgs]],[[q,q],[W+]]] element mass: [[8.65E+02 [GeV], 2.69E+02 [GeV], 1.29E+02 [GeV]], [8.65E+02 [GeV], 2.69E+02 [GeV], 1.29E+02 [GeV]]] [[[c,c],[higgs]],[[c,q],[W-]]] element mass: [[8.65E+02 [GeV], 2.69E+02 [GeV], 1.29E+02 [GeV]], [8.65E+02 [GeV], 2.69E+02 [GeV], 1.29E+02 [GeV]]] [[[c,c],[higgs]],[[c,q],[W+]]] element mass: [[8.65E+02 [GeV], 2.69E+02 [GeV], 1.29E+02 [GeV]], [8.65E+02 [GeV], 2.69E+02 [GeV], 1.29E+02 [GeV]]] missing topology: [[[jet],[higgs]],[[jet],[W]]] (MET,MET) Contributing elements: [[[q],[higgs]],[[q],[W+]]] element mass: [[9.92E+02 [GeV], 2.69E+02 [GeV], 1.29E+02 [GeV]], [9.92E+02 [GeV], 2.69E+02 [GeV], 1.29E+02 [GeV]]] [[[q],[W+]],[[q],[higgs]]] element mass: [[9.92E+02 [GeV], 2.69E+02 [GeV], 1.29E+02 [GeV]], [9.94E+02 [GeV], 2.69E+02 [GeV], 1.29E+02 [GeV]]] [[[q],[higgs]],[[q],[W-]]] element mass: [[9.92E+02 [GeV], 2.69E+02 [GeV], 1.29E+02 [GeV]], [9.94E+02 [GeV], 2.69E+02 [GeV], 1.29E+02 [GeV]]] [[[q],[higgs]],[[q],[W-]]] element mass: [[9.94E+02 [GeV], 2.69E+02 [GeV], 1.29E+02 [GeV]], [9.94E+02 [GeV], 2.69E+02 [GeV], 1.29E+02 [GeV]]] [[[c],[higgs]],[[q],[W+]]] element mass: [[9.92E+02 [GeV], 2.69E+02 [GeV], 1.29E+02 [GeV]], [9.92E+02 [GeV], 2.69E+02 [GeV], 1.29E+02 [GeV]]] [[[q],[higgs]],[[q],[W-]]] element mass: [[9.92E+02 [GeV], 2.69E+02 [GeV], 1.29E+02 [GeV]], [9.92E+02 [GeV], 2.69E+02 [GeV], 1.29E+02 [GeV]]] [[[q],[higgs]],[[c],[W-]]] element mass: [[9.92E+02 [GeV], 2.69E+02 [GeV], 1.29E+02 [GeV]], [9.94E+02 [GeV], 2.69E+02 [GeV], 1.29E+02 [GeV]]] [[[q],[higgs]],[[c],[W-]]] element mass: [[9.94E+02 [GeV], 2.69E+02 [GeV], 1.29E+02 [GeV]], [9.94E+02 [GeV], 2.69E+02 [GeV], 1.29E+02 [GeV]]] [[[q],[higgs]],[[q],[W+]]] element mass: [[9.92E+02 [GeV], 2.69E+02 [GeV], 1.29E+02 [GeV]], [9.94E+02 [GeV], 2.69E+02 [GeV], 1.29E+02 [GeV]]] [[[c],[higgs]],[[q],[W-]]] element mass: [[9.92E+02 [GeV], 2.69E+02 [GeV], 1.29E+02 [GeV]], [9.92E+02 [GeV], 2.69E+02 [GeV], 1.29E+02 [GeV]]] [[[q],[higgs]],[[q],[W+]]] element mass: [[9.94E+02 [GeV], 2.69E+02 [GeV], 1.29E+02 [GeV]], [9.94E+02 [GeV], 2.69E+02 [GeV], 1.29E+02 [GeV]]] [[[q],[higgs]],[[c],[W+]]] element mass: [[9.94E+02 [GeV], 2.69E+02 [GeV], 1.29E+02 [GeV]], [9.94E+02 [GeV], 2.69E+02 [GeV], 1.29E+02 [GeV]]] [[[q],[higgs]],[[c],[W+]]] element mass: [[9.92E+02 [GeV], 2.69E+02 [GeV], 1.29E+02 [GeV]], [9.94E+02 [GeV], 2.69E+02 [GeV], 1.29E+02 [GeV]]] [[[q],[W-]],[[q],[higgs]]] element mass: [[9.92E+02 [GeV], 2.69E+02 [GeV], 1.29E+02 [GeV]], [9.94E+02 [GeV], 2.69E+02 [GeV], 1.29E+02 [GeV]]] [[[c],[higgs]],[[q],[W-]]] element mass: [[9.92E+02 [GeV], 2.69E+02 [GeV], 1.29E+02 [GeV]], [9.94E+02 [GeV], 2.69E+02 [GeV], 1.29E+02 [GeV]]] missing topology: [[[jet,jet],[higgs]],[[jet,jet],[higgs]]] (MET,MET) Contributing elements: [[[q,q],[higgs]],[[q,q],[higgs]]] element mass: [[8.65E+02 [GeV], 2.69E+02 [GeV], 1.29E+02 [GeV]], [8.65E+02 [GeV], 2.69E+02 [GeV], 1.29E+02 [GeV]]] [[[q,q],[higgs]],[[c,c],[higgs]]] element mass: [[8.65E+02 [GeV], 2.69E+02 [GeV], 1.29E+02 [GeV]], [8.65E+02 [GeV], 2.69E+02 [GeV], 1.29E+02 [GeV]]] [[[c,c],[higgs]],[[c,c],[higgs]]] element mass: [[8.65E+02 [GeV], 2.69E+02 [GeV], 1.29E+02 [GeV]], [8.65E+02 [GeV], 2.69E+02 [GeV], 1.29E+02 [GeV]]] missing topology: [[[jet,jet],[W]],[[jet,jet],[Z]]] (MET,MET) Contributing elements: [[[q,q],[Z]],[[q,q],[W-]]] element mass: [[8.65E+02 [GeV], 2.69E+02 [GeV], 1.29E+02 [GeV]], [8.65E+02 [GeV], 2.69E+02 [GeV], 1.29E+02 [GeV]]] [[[q,q],[Z]],[[q,q],[W+]]] element mass: [[8.65E+02 [GeV], 2.69E+02 [GeV], 1.29E+02 [GeV]], [8.65E+02 [GeV], 2.69E+02 [GeV], 1.29E+02 [GeV]]] [[[q,q],[Z]],[[c,q],[W-]]] element mass: [[8.65E+02 [GeV], 2.69E+02 [GeV], 1.29E+02 [GeV]], [8.65E+02 [GeV], 2.69E+02 [GeV], 1.29E+02 [GeV]]] [[[q,q],[Z]],[[c,q],[W+]]] element mass: [[8.65E+02 [GeV], 2.69E+02 [GeV], 1.29E+02 [GeV]], [8.65E+02 [GeV], 2.69E+02 [GeV], 1.29E+02 [GeV]]] [[[c,c],[Z]],[[q,q],[W-]]] element mass: [[8.65E+02 [GeV], 2.69E+02 [GeV], 1.29E+02 [GeV]], [8.65E+02 [GeV], 2.69E+02 [GeV], 1.29E+02 [GeV]]] [[[c,c],[Z]],[[q,q],[W+]]] element mass: [[8.65E+02 [GeV], 2.69E+02 [GeV], 1.29E+02 [GeV]], [8.65E+02 [GeV], 2.69E+02 [GeV], 1.29E+02 [GeV]]] [[[c,c],[Z]],[[c,q],[W-]]] element mass: [[8.65E+02 [GeV], 2.69E+02 [GeV], 1.29E+02 [GeV]], [8.65E+02 [GeV], 2.69E+02 [GeV], 1.29E+02 [GeV]]] [[[c,c],[Z]],[[c,q],[W+]]] element mass: [[8.65E+02 [GeV], 2.69E+02 [GeV], 1.29E+02 [GeV]], [8.65E+02 [GeV], 2.69E+02 [GeV], 1.29E+02 [GeV]]] missing topology: [[[jet],[higgs]],[[jet],[higgs]]] (MET,MET) Contributing elements: [[[q],[higgs]],[[q],[higgs]]] element mass: [[9.92E+02 [GeV], 2.69E+02 [GeV], 1.29E+02 [GeV]], [9.94E+02 [GeV], 2.69E+02 [GeV], 1.29E+02 [GeV]]] [[[q],[higgs]],[[q],[higgs]]] element mass: [[9.92E+02 [GeV], 2.69E+02 [GeV], 1.29E+02 [GeV]], [9.92E+02 [GeV], 2.69E+02 [GeV], 1.29E+02 [GeV]]] [[[q],[higgs]],[[q],[higgs]]] element mass: [[9.94E+02 [GeV], 2.69E+02 [GeV], 1.29E+02 [GeV]], [9.94E+02 [GeV], 2.69E+02 [GeV], 1.29E+02 [GeV]]] [[[q],[higgs]],[[c],[higgs]]] element mass: [[9.92E+02 [GeV], 2.69E+02 [GeV], 1.29E+02 [GeV]], [9.92E+02 [GeV], 2.69E+02 [GeV], 1.29E+02 [GeV]]] [[[c],[higgs]],[[c],[higgs]]] element mass: [[9.92E+02 [GeV], 2.69E+02 [GeV], 1.29E+02 [GeV]], [9.92E+02 [GeV], 2.69E+02 [GeV], 1.29E+02 [GeV]]] [[[c],[higgs]],[[q],[higgs]]] element mass: [[9.92E+02 [GeV], 2.69E+02 [GeV], 1.29E+02 [GeV]], [9.94E+02 [GeV], 2.69E+02 [GeV], 1.29E+02 [GeV]]]