# 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.theory import slhaDecomposer
from smodels.installation import installDirectory
from smodels.tools.physicsUnits import fb, GeV
# Define the SLHA file name
filename="inputFiles/slha/gluino_squarks.slha"
# Perform the decomposition:
listOfTopologies = slhaDecomposer.decompose(filename, sigcut = 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.getParticles())
print ('masses=',element.getMasses())
print ('weight=',element.weight,'\n')
Element 0 : [[['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]', '1.30E+01 [TeV]:1.45E-01 [pb]'] Element 1 : [[['W+']], [['Z']]] 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]', '1.30E+01 [TeV]:2.62E-02 [pb]'] Element 2 : [[['W+']], [['higgs']]] 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]', '1.30E+01 [TeV]:1.36E-01 [pb]'] Element 3 : [[['W-']], [['Z']]] 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]', '1.30E+01 [TeV]:1.37E-02 [pb]'] Element 4 : [[['W-']], [['higgs']]] 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]', '1.30E+01 [TeV]:7.11E-02 [pb]'] Element 5 : [[['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]', '1.30E+01 [TeV]:1.85E-03 [pb]'] 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]', '1.30E+01 [TeV]:7.31E-03 [pb]'] 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]', '1.30E+01 [TeV]:2.27E-03 [pb]']