How To: Print out the theoretical decomposition

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.theory import slhaDecomposer
from smodels.installation import installDirectory
from smodels.tools.physicsUnits import fb, GeV
In [3]:
# Define the SLHA file name
filename="%s/inputFiles/slha/gluino_squarks.slha" % installDirectory()
In [4]:
# Perform the decomposition:
listOfTopologies = slhaDecomposer.decompose(filename, sigcut = 0.5 * fb, 
                                            doCompress=True, doInvisible=True,minmassgap = 5* GeV)
In [5]:
# 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
In [6]:
# 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
In [7]:
# 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]']