# 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
# (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
# Define the SLHA file name
filename="inputFiles/slha/gluino_squarks.slha"
model = Model(BSMparticles = BSMList, SMparticles = SMList)
model.updateParticles(inputFile=filename)
# Perform the decomposition:
listOfTopologies = decomposer.decompose(model, sigmacut = 0.5 * fb,
doCompress=True, doInvisible=True,minmassgap = 5* GeV)
# Print a summary of all the topologies generated:
for top in listOfTopologies:
print (top.describe())
number of vertices: [0, 0], number of vertex particles: [[], []], number of elements: 1 number of vertices: [0, 1], number of vertex particles: [[], [1]], number of elements: 1 number of vertices: [0, 1], number of vertex particles: [[], [2]], number of elements: 2 number of vertices: [1, 1], number of vertex particles: [[1], [1]], number of elements: 8 number of vertices: [1, 1], number of vertex particles: [[1], [2]], number of elements: 10 number of vertices: [1, 1], number of vertex particles: [[2], [2]], number of elements: 3 number of vertices: [0, 2], number of vertex particles: [[], [2, 1]], number of elements: 9 number of vertices: [1, 2], number of vertex particles: [[1], [1, 1]], number of elements: 28 number of vertices: [1, 2], number of vertex particles: [[1], [1, 2]], number of elements: 4 number of vertices: [1, 2], number of vertex particles: [[1], [2, 1]], number of elements: 50 number of vertices: [1, 2], number of vertex particles: [[2], [1, 1]], number of elements: 14 number of vertices: [1, 2], number of vertex particles: [[2], [1, 2]], number of elements: 4 number of vertices: [1, 2], number of vertex particles: [[2], [2, 1]], number of elements: 20 number of vertices: [2, 2], number of vertex particles: [[1, 1], [1, 1]], number of elements: 31 number of vertices: [2, 2], number of vertex particles: [[1, 1], [1, 2]], number of elements: 8 number of vertices: [2, 2], number of vertex particles: [[1, 1], [2, 1]], number of elements: 87 number of vertices: [2, 2], number of vertex particles: [[1, 2], [2, 1]], number of elements: 42 number of vertices: [2, 2], number of vertex particles: [[2, 1], [2, 1]], number of elements: 72 number of vertices: [0, 3], number of vertex particles: [[], [1, 2, 1]], number of elements: 6 number of vertices: [1, 3], number of vertex particles: [[1], [1, 2, 1]], number of elements: 37 number of vertices: [1, 3], number of vertex particles: [[2], [1, 2, 1]], number of elements: 42 number of vertices: [2, 3], number of vertex particles: [[1, 1], [1, 2, 1]], number of elements: 82 number of vertices: [2, 3], number of vertex particles: [[2, 1], [1, 2, 1]], number of elements: 276
# To print specific information about othe i-th topology:
i = 3
top = listOfTopologies[i]
print ("Number of vertices = ",top.vertnumb)
print ("Number of final states = ",top.vertparts)
print ("Number of elements = ",len(top.elementList))
Number of vertices = [1, 1] Number of final states = [[1], [1]] Number of elements = 8
# We can also print information for each element in the topology:
for iel,element in enumerate(top.elementList):
print ('Element',iel,':',element.evenParticles)
print ('masses=',element.mass)
print ('weight=',element.weight,'\n')
Element 0 : [[[higgs]], [[higgs]]] masses= [[2.69E+02 [GeV], 1.29E+02 [GeV]], [2.69E+02 [GeV], 1.29E+02 [GeV]]] weight= ['8.00E+00 [TeV]:5.63E-04 [pb] (1000023, 1000023)', '1.30E+01 [TeV]:1.85E-03 [pb] (1000023, 1000023)'] Element 1 : [[[higgs]], [[W-]]] masses= [[2.69E+02 [GeV], 1.29E+02 [GeV]], [2.69E+02 [GeV], 1.29E+02 [GeV]]] weight= ['8.00E+00 [TeV]:2.74E-02 [pb] (-1000024, 1000023)', '1.30E+01 [TeV]:7.11E-02 [pb] (-1000024, 1000023)'] Element 2 : [[[higgs]], [[W+]]] masses= [[2.69E+02 [GeV], 1.29E+02 [GeV]], [2.69E+02 [GeV], 1.29E+02 [GeV]]] weight= ['8.00E+00 [TeV]:6.84E-02 [pb] (1000023, 1000024)', '1.30E+01 [TeV]:1.36E-01 [pb] (1000023, 1000024)'] Element 3 : [[[Z]], [[W-]]] masses= [[2.69E+02 [GeV], 1.29E+02 [GeV]], [2.69E+02 [GeV], 1.29E+02 [GeV]]] weight= ['8.00E+00 [TeV]:5.29E-03 [pb] (-1000024, 1000023)', '1.30E+01 [TeV]:1.37E-02 [pb] (-1000024, 1000023)'] Element 4 : [[[Z]], [[W+]]] masses= [[2.69E+02 [GeV], 1.29E+02 [GeV]], [2.69E+02 [GeV], 1.29E+02 [GeV]]] weight= ['8.00E+00 [TeV]:1.32E-02 [pb] (1000023, 1000024)', '1.30E+01 [TeV]:2.62E-02 [pb] (1000023, 1000024)'] Element 5 : [[[W-]], [[W+]]] masses= [[2.69E+02 [GeV], 1.29E+02 [GeV]], [2.69E+02 [GeV], 1.29E+02 [GeV]]] weight= ['8.00E+00 [TeV]:5.96E-02 [pb] (-1000024, 1000024)', '1.30E+01 [TeV]:1.45E-01 [pb] (-1000024, 1000024)'] Element 6 : [[[q]], [[q]]] masses= [[9.91E+02 [GeV], 1.29E+02 [GeV]], [9.91E+02 [GeV], 1.29E+02 [GeV]]] weight= ['8.00E+00 [TeV]:1.36E-03 [pb] (None, None)', '1.30E+01 [TeV]:7.31E-03 [pb] (None, None)'] Element 7 : [[[q]], [[q]]] masses= [[9.91E+02 [GeV], 1.29E+02 [GeV]], [9.92E+02 [GeV], 1.29E+02 [GeV]]] weight= ['8.00E+00 [TeV]:4.15E-04 [pb] (2000001, 2000002)', '1.30E+01 [TeV]:2.27E-03 [pb] (2000001, 2000002)']