Output Description¶
A detailed description of the possible output formats generated by running SModelS and their content is given below. For simplicity we will assume that all printer options in the parameters file are set to True, so the output information is maximal [*].
Screen (Stdout) Output¶
The stdout (or log output) is intended to provide extensive information about the database, the decomposition, the theory predictions and the missing topologies. It is most convenient if the input is a single file and not a folder, since the output is quite extensive. If all the options in stdout-printer are set to True (see parameters file), the screen output contains the following information:
- information about the basic input parameters and the status of the run:
Input status: 1
Decomposition output status: 1 #decomposition was successful
# Input File: inputFiles/slha/gluino_squarks.slha
# maxcond = 0.2
# minmassgap = 5.
# ncpus = 1
# sigmacut = 0.01
# Database version: 1.2.0
- a list of all the experimental results considered (if printDatabase = True). Note that this list corresponds to all the results selected in the database options (see parameters file). If addAnaInfo = True, for each experimental result entry a list of all the simplified models (or elements) constrained by the analysis is also shown using the bracket notation including the Z2-odd final state class:
=======================================================
|| ||
|| Selected Experimental Results ||
|| ||
=======================================================
========================================================
Experimental Result ID: ATLAS-SUSY-2015-01
Tx Labels: ['T2bb']
Sqrts: 1.30E+01 [TeV]
-----------------------------
Elements tested by analysis:
[[[b]],[[b]]] (MET, MET)
========================================================
[[[t,t]],[[t,t]]] (MET, MET)
========================================================
Experimental Result ID: ATLAS-SUSY-2013-18
Tx Labels: ['T1bbbb', 'T1tttt']
Sqrts: 8.00E+00 [TeV]
-----------------------------
Elements tested by analysis:
[[[b,b]],[[b,b]]] (MET, MET)
[[[t,t]],[[t,t]]] (MET, MET)
========================================================
- a full list of the topologies generated by the decomposition (if printDecomp = True). Each topology entry contains basic information about the topology as well as the number of elements with this topology and the sum over all the elements weights. If addElementInfo = True, the elements belonging to each topology are also explicitly shown, as well as the element’s mass, final states, weight, the PIDs of the intermediate particles contributing to the element and the element ID:
=======================================================
|| ||
|| Topologies Table ||
|| ||
=======================================================
=====================================================
Topology:
Number of vertices: [0, 0]
Number of vertex parts: [[], []]
Total Global topology weight :
Sqrts: 8.00E+00 [TeV], Weight:4.81E-04 [pb]
Sqrts: 1.30E+01 [TeV], Weight:1.58E-03 [pb]
Total Number of Elements: 1
.........................................................................
Element:
Element ID: 1
Particles in element: [[], []]
Final states in element: ['MET', 'MET']
The element masses are
Branch 0: [1.29E+02 [GeV]]
Branch 1: [1.29E+02 [GeV]]
The element PIDs are
PIDs: [[1000022], [1000022]]
The element weights are:
Sqrts: 8.00E+00 [TeV], Weight:4.81E-04 [pb]
Sqrts: 1.30E+01 [TeV], Weight:1.58E-03 [pb]
=====================================================
Topology:
Number of vertices: [0, 1]
Number of vertex parts: [[], [1]]
Total Global topology weight :
Sqrts: 8.00E+00 [TeV], Weight:1.16E-03 [pb]
Sqrts: 1.30E+01 [TeV], Weight:3.71E-03 [pb]
Total Number of Elements: 7
.........................................................................
Element:
Element ID: 2
Particles in element: [[], [['W+']]]
Final states in element: ['MET', 'MET']
The element masses are
Branch 0: [1.29E+02 [GeV]]
Branch 1: [2.69E+02 [GeV], 1.29E+02 [GeV]]
The element PIDs are
PIDs: [[1000022], [1000024, 1000022]]
The element weights are:
Sqrts: 8.00E+00 [TeV], Weight:2.40E-04 [pb]
Sqrts: 1.30E+01 [TeV], Weight:2.63E-04 [pb]
.........................................................................
Element:
Element ID: 3
Particles in element: [[], [['W-']]]
Final states in element: ['MET', 'MET']
The element masses are
Branch 0: [1.29E+02 [GeV]]
Branch 1: [2.69E+02 [GeV], 1.29E+02 [GeV]]
The element PIDs are
PIDs: [[1000022], [-1000024, -1000022]]
The element weights are:
Sqrts: 8.00E+00 [TeV], Weight:4.01E-05 [pb]
- a list of all the theory predictions obtained and the corresponding experimental result upper limit. For each experimental result, the corresponding id, signal region (data set) and sqrts as well as the constrained simplified models (txnames) are printed. After this basic information, the signal cross section (theory prediction), the list of condition values (if applicable) and the corresponding observed upper limit are shown. Also, if available, the expected upper limit is included. If computeStatistics = True, the \(\chi^2\) and likelihood values are printed (see likelihood calculation). Finally, if printExtendedResults = True, basic information about the elements being constrained, such as their masses, IDs and PIDs, is also shown.
=======================================================
|| ||
|| Theory Predictions and ||
|| Experimental Constraints ||
|| ||
=======================================================
---------------Analysis Label = CMS-SUS-16-036
-------------------Dataset Label = (UL)
-------------------Txname Labels = ['T2']
Analysis sqrts: 1.30E+01 [TeV]
Theory prediction: 1.17E-02 [pb]
Theory conditions:[None]
Observed experimental limit: 7.68E+00 [fb]
Observed r-Value: 1.5212733123965445
Masses in branch 0: [9.91E+02 [GeV], 1.29E+02 [GeV]]
Masses in branch 1: [9.91E+02 [GeV], 1.29E+02 [GeV]]
Contributing elements: [28, 29, 30, 34, 35, 36, 37, 38, 39, 40]
PIDs:[[-2000004, -1000022], [2000004, 1000022]]
PIDs:[[-2000004, -1000022], [2000002, 1000022]]
---------------Analysis Label = CMS-SUS-16-033
-------------------Dataset Label = (UL)
-------------------Txname Labels = ['T2']
Analysis sqrts: 1.30E+01 [TeV]
Theory prediction: 1.17E-02 [pb]
Theory conditions:[None]
Observed experimental limit: 8.64E+00 [fb]
Observed r-Value: 1.35282642984693
Masses in branch 0: [9.91E+02 [GeV], 1.29E+02 [GeV]]
Masses in branch 1: [9.91E+02 [GeV], 1.29E+02 [GeV]]
Contributing elements: [28, 29, 30, 34, 35, 36, 37, 38, 39, 40]
PIDs:[[-2000004, -1000022], [2000004, 1000022]]
- summary information about the missing topologies, if testCoverage = True. The first information corresponds to the total input cross-section considered, which corresponds to the cross-section summed over all elements after decomposition. Note that this value might differ from the total input model cross-section, since it includes the lifetime reweighting and the effect of skipping elements with a weight below the minimum decomposition weight. The total missing topology cross section shown corresponds to the sum of cross sections of all elements which are not tested by any experimental result. If, however, the element is constrained by one or more experimental results, but its mass is outside the efficiency or upper limit grids (see EM-type results and UL-type results), its cross section is included in the total cross section outside the grid. Finally, the elements which contribute to the total missing topology cross section are subdivided into elements with long decays or with asymmetric branches (see coverage tool for more details)
Total cross section considered (fb): 3.635E+03
Total missing topology cross section (fb): 3.037E+03
Total cross section where we are outside the mass grid (fb): 3.376E+00
Total cross section in long cascade decays (fb): 1.317E+03
Total cross section in decays with asymmetric branches (fb): 1.686E+03
- detailed information about the missing topologies with highest cross sections. The element cross section (weight) as well as its description in bracket notation and BSM final state classification are included. If addCoverageID = True, all the elements IDs contributing to the missing topology are shown. These IDs can be traced back to the corresponding elements using the decomposition information obtained with printDecomp = True and addElementInfo = True.
================================================================================
Missing topologies with the highest cross sections (up to 10):
Sqrts (TeV) Weight (fb) Element description
13.0 1.482E+02 # [[[jet],[W]],[[jet,jet],[W]]](MET,MET)
Contributing elements [987, 988, 989, 990, 1001, 1002, 1003, 1004, 1033, 1034, 1035, 1036, 1047, 1048, 1049, 1050, 1061, 1062, 1063, 1064, 1075, 1076, 1077, 1078]
- detailed information about the topologies which are outside the experimental results grid. If addCoverageID = True, all the elements IDs contributing to the missing topology are shown.
================================================================================
Contributions outside the mass grid (up to 10):
Sqrts (TeV) Weight (fb) Element description
13.0 1.438E+00 # [[[t],[W]],[[t],[W]]](MET,MET)
Contributing elements [813, 814, 815]
- information about the missing topologies with long cascade decays. The long cascade decays are classified by the initially produced mother particles. If more than one pair of mothers are contributing to the same class of elements, the full list is given in the comment. For definiteness all lists are sorted. If addCoverageID = True, all the elements IDs contributing to the missing topology are shown.
Mother1 Mother2 Weight (fb) # allMothers
1000021 2000002 3.696E+02 # [[1000021, 2000002]]
Contributing elements [1560, 1561, 1562, 1563, 1564, 1916, 1917, 1921, 1925, 1929, 1930, 1934, 1938, 1942, 1946, 1950, 1951, 1955, 1965, 1969, 1973, 1977, 1981, 1985, 1989, 1993, 1997, 2001, 2005, 2009, 2013, 2017, 2028, 2032, 2036, 2040, 2044, 2048, 2052, 2056, 2060, 2064, 2068, 2072, 2076, 2080, 2084, 2087, 2091, 2095, 2099, 2102, 2106, 2110, 2114, 2118, 2122, 2125, 2129, 2130, 3302, 3305, 3309, 3313, 3317, 3320, 3324, 3328, 3332, 3336, 3340, 3343, 3347, 3348, 3361, 3365, 3369, 3373, 3377, 3381, 3385, 3389, 3393, 3397, 3401, 3405, 3409, 3413, 3433, 3437, 3441, 3445, 3449, 3453, 3457, 3461, 3465, 3469, 3473, 3477, 3481, 3485, 3505, 3509, 3513, 3517, 3521, 3525, 3529, 3533, 3537, 3541, 3545, 3549, 3553, 3557, 3561, 3564, 3568, 3572, 3576, 3579, 3583, 3587, 3591, 3595, 3599, 3602, 3606, 3607, 3618, 3622, 3626, 3630, 3634, 3638, 3642, 3646, 3650, 3654, 3658, 3662, 3666, 3670, 3690, 3694, 3698, 3702, 3706, 3710, 3714, 3718, 3722, 3726, 3730, 3734, 3738, 3742, 3762, 3766, 3770, 3774, 3778, 3782, 3786, 3790, 3794, 3798, 3802, 3806, 3810, 3814, 3834, 3838, 3842, 3846, 3850, 3854, 3858, 3862, 3866, 3870, 3874, 3878, 3882, 3886, 3906, 3910, 3914, 3918, 3922, 3926, 3930, 3934, 3938, 3942, 3946, 3950, 3954, 3958, 3962, 3965, 3969, 3973, 3977, 3980, 3984, 3988, 3992, 3996, 4000, 4003, 4007, 4008, 4019, 4023, 4027, 4031, 4035, 4039, 4043, 4047, 4051, 4055, 4059, 4063, 4067, 4071, 4075, 4076, 4080, 4084, 4088, 4089, 4093, 4097, 4101, 4105, 4109, 4110, 4114, 4123, 4127, 4131, 4135, 4139, 4143, 4147, 4151, 4155, 4159, 4163, 4167, 4171, 4174]
1000021 2000001 2.031E+02 # [[1000021, 2000001], [1000021, 2000003]]
Contributing elements [1971, 1975, 1979, 1987, 1991, 1995, 1999, 2003, 2011, 2019, 2034, 2038, 2042, 2050, 2054, 2058, 2062, 2066, 2074, 2082, 3311, 3315, 3326, 3330, 3334, 3338, 3367, 3371, 3375, 3383, 3387, 3391, 3395, 3399, 3407, 3415, 3435, 3439, 3443, 3447, 3455, 3459, 3463, 3467, 3471, 3479, 3487, 3507, 3511, 3515, 3519, 3527, 3531, 3535, 3539, 3543, 3551, 3559, 3624, 3628, 3632, 3640, 3644, 3648, 3652, 3656, 3664, 3672, 3692, 3696, 3700, 3704, 3712, 3716, 3720, 3724, 3728, 3736, 3744, 3764, 3768, 3772, 3776, 3784, 3788, 3792, 3796, 3800, 3808, 3816, 3836, 3840, 3844, 3848, 3856, 3860, 3864, 3868, 3872, 3880, 3888, 3908, 3912, 3916, 3920, 3928, 3932, 3936, 3940, 3944, 3952, 3960, 4025, 4029, 4033, 4041, 4045, 4049, 4053, 4057, 4065, 4073, 4129, 4133, 4137, 4145, 4149, 4153, 4157, 4161, 4169, 1919, 1923, 1927, 1932, 1936, 1940, 1944, 1948, 1953, 1957, 1967, 1983, 2007, 2015, 2030, 2046, 2070, 2078, 2086, 2089, 2093, 2097, 2101, 2104, 2108, 2112, 2116, 2120, 2124, 2127, 2132, 3304, 3307, 3319, 3322, 3342, 3345, 3350, 3363, 3379, 3403, 3411, 3451, 3475, 3483, 3523, 3547, 3555, 3563, 3566, 3570, 3574, 3578, 3581, 3585, 3589, 3593, 3597, 3601, 3604, 3609, 3620, 3636, 3660, 3668, 3708, 3732, 3740, 3780, 3804, 3812, 3852, 3876, 3884, 3924, 3948, 3956, 3964, 3967, 3971, 3975, 3979, 3982, 3986, 3990, 3994, 3998, 4002, 4005, 4010, 4021, 4037, 4061, 4069, 4078, 4082, 4086, 4091, 4095, 4099, 4103, 4107, 4112, 4116, 4125, 4141, 4165, 4173, 4176]
- information about the missing topologies with asymmetric decays. The asymmetric branch decays are classified by the initially produced mother particles. If more than one pair of mothers are contributing to the same class of elements, the full list is given in the comment. For definiteness all lists are sorted. If addCoverageID = True, all the elements IDs contributing to the missing topology are shown.
Mother1 Mother2 Weight (fb) # allMothers
1000021 1000021 5.297E+02 # [[1000021, 1000021]]
Contributing elements [565, 566, 567, 568, 569, 570, 571, 572, 573, 574, 575, 576, 577, 578, 579, 580, 581, 582, 583, 584, 585, 586, 587, 588, 589, 590, 591, 592, 593, 594, 595, 596, 597, 598, 599, 600, 601, 602, 603, 604, 605, 606, 607, 608, 609, 610, 611, 612, 613, 614, 615, 616, 617, 618, 619, 620, 1414, 1415, 1416, 1417, 1418, 1419, 1420, 1421, 1422, 1423, 1424, 1425, 1426, 1428, 1429, 1430, 1431, 1432, 1433, 1434, 1435, 1436, 1437, 1438, 1439, 1441, 1442, 1443, 1444, 1445, 1446, 1447, 1448, 1449, 1450, 1451, 1453, 1454, 1455, 1456, 1457, 1458, 1459, 1460, 1461, 1462, 1464, 1465, 1466, 1467, 1468, 1470, 1471, 1472, 1474, 1475, 1476, 1477, 1478, 1479, 1480, 1481, 1486, 1487, 1488, 1489, 1493, 1494, 1495, 1496, 1499, 1500, 1501, 1502, 1504, 1505, 1506, 1507, 1509, 1510, 1511, 1513, 1514, 1516]
1000002 1000021 4.149E+02 # [[1000002, 1000021], [1000004, 1000021]]
Contributing elements [466, 478, 480, 492, 494, 504, 506, 1027, 1028, 1029, 1030, 1031, 1032, 1033, 1034, 1035, 1036, 1037, 1038, 1039, 1040, 1055, 1056, 1057, 1058, 1059, 1060, 1061, 1062, 1063, 1064, 1065, 1066, 1068, 422, 423, 424, 425, 426, 427, 428, 429, 430, 431, 432, 433, 434, 435, 468, 470, 472, 482, 484, 496, 498, 508, 510, 514, 518, 523, 527, 531, 535, 539, 543, 547, 551, 556, 560, 1067, 1083, 1084, 1085, 1086, 1087, 1088, 1089, 1090, 1091, 1092, 1093, 1094, 1095, 1096, 1111, 1112, 1113, 1114, 1115, 1116, 1117, 1118, 1119, 1120, 1121, 1122, 1123, 1124, 1222, 1223, 1224, 1225, 1226, 1227, 1228, 1229, 1230, 1231, 1265, 1266, 1267, 1268, 1269, 1270, 1271, 1272, 1273, 1274, 1275, 1276, 1277, 1278, 1321, 1322, 1323, 1324, 1325, 1326, 1327, 1328, 1329, 1330, 1331, 1332, 1333, 1334, 1377, 1378, 1379, 1380, 1381, 1382, 1383, 1384, 1385, 1386, 1387, 1388, 1389, 477, 491, 1009, 1010, 1011, 1012, 1013, 1014, 1015, 1016, 1017, 1018, 1019, 1020, 1021, 1022, 1023, 1024, 1025, 1026]
Log Output¶
The log-type output is identical to the screen output, except that it is redirected to a .log file. The filename is set as the <input file>.log and stored in the output folder (see the runSModelS options).
Summary File Output¶
The summary-type output is similar to the screen output, but restricted to the list of theory predictions and model coverage. The output is printed to the file <input file>.smodels and stored in the output folder (see the runSModelS options).
Below we describe in detail the blocks contained in the summary output:
- information about the basic input parameters and the status of the run:
Input status: 1
Decomposition output status: 1 #decomposition was successful
# Input File: inputFiles/slha/gluino_squarks.slha
# maxcond = 0.2
# minmassgap = 5.
# ncpus = 1
# sigmacut = 0.01
# Database version: 1.2.0
- a list of all the theory predictions obtained and the corresponding experimental result upper limit. If expandedSummary = False only the most constraining experimental result is printed. For each theory prediction entry, the corresponding experimental result id, the signal region (data set) used (only for EM-type results) and the experimental result sqrts is printed. Furthermore, the txnames contributing to the signal cross section, the theory cross section (Theory_Value), the observed upper limit (Exp_limit), the (theory cross section)/(observed upper limit) ratio (r) and, when available, the (theory cross section)/(expected upper limit) ratio (r_expect) are also printed. For UL-type results the condition violation (see upper limit conditions) is also included. Finally, if computeStatistics = True, the \(\chi^2\) and likelihood values (for EM-type results) are printed:
#Analysis Sqrts Cond_Violation Theory_Value(fb) Exp_limit(fb) r r_expected
CMS-SUS-16-036 1.30E+01 0.0 1.169E+01 7.683E+00 1.521E+00 N/A
Signal Region: (UL)
Txnames: T2
--------------------------------------------------------------------------------
CMS-SUS-16-033 1.30E+01 0.0 1.169E+01 8.639E+00 1.353E+00 N/A
Signal Region: (UL)
Txnames: T2
--------------------------------------------------------------------------------
ATLAS-SUSY-2013-02 8.00E+00 0.0 6.270E-01 1.818E+00 3.449E-01 4.146E-01
Signal Region: SR2jt
Txnames: T1, T2
Chi2, Likelihood = 1.302E-01 1.276E-03
- the maximum value for the (theory cross section)/(observed upper limit) ratio. If this value is higher than 1 the input model is likely excluded by one of the experimental results (see confronting predictions)
================================================================================
The highest r value is = 1.5212733123965445
- summary information about the missing topologies, if testCoverage = True. The first information corresponds to the total input cross-section considered, which corresponds to the cross-section summed over all elements after decomposition. Note that this value might differ from the total input model cross-section, since it includes the lifetime reweighting and the effect of skipping elements with a weight below the minimum decomposition weight. The total missing topology cross section corresponds to the sum of all elements cross sections which are not tested by any experimental result. If, however, the element is constrained by one or more experimental results, but its mass is outside the efficiency or upper limit grids (see EM-type results and UL-type results), its cross section is included in the total cross section outside the grid. Finally, the elements which contribute to the total missing topology cross section are subdivided into elements with long decays or with asymmetric branches (see coverage tool for more details)
Total cross section considered (fb): 3.635E+03
Total missing topology cross section (fb): 3.037E+03
Total cross section where we are outside the mass grid (fb): 3.376E+00
Total cross section in long cascade decays (fb): 1.317E+03
Total cross section in decays with asymmetric branches (fb): 1.686E+03
- detailed information about the missing topologies with highest cross sections. The element cross section (weight) as well as its description in bracket notation and final states is included.
================================================================================
Missing topologies with the highest cross sections (up to 10):
Sqrts (TeV) Weight (fb) Element description
13.0 1.482E+02 # [[[jet],[W]],[[jet,jet],[W]]](MET,MET)
13.0 1.445E+02 # [[[jet,jet],[W]],[[jet],[jet,jet],[W]]](MET,MET)
- detailed information about the topologies which are outside the experimental results grid:
================================================================================
Contributions outside the mass grid (up to 10):
Sqrts (TeV) Weight (fb) Element description
13.0 1.438E+00 # [[[t],[W]],[[t],[W]]](MET,MET)
13.0 2.871E-01 # [[[b],[higgs]],[[b],[higgs]]](MET,MET)
- information about the missing topologies with long cascade decays:
Missing topos: long cascade decays (up to 10 entries), sqrts = 13 TeV:
Mother1 Mother2 Weight (fb) # allMothers
1000021 2000002 3.696E+02 # [[1000021, 2000002]]
1000021 2000001 2.031E+02 # [[1000021, 2000001], [1000021, 2000003]]
1000002 1000021 1.612E+02 # [[1000002, 1000021]]
- information about the missing topologies with asymmetric decays:
Missing topos: asymmetric branches (w/o long cascades, up to 10), sqrts = 13 TeV
Mother1 Mother2 Weight (fb) # allMothers
1000021 1000021 5.297E+02 # [[1000021, 1000021]]
1000002 1000021 4.149E+02 # [[1000002, 1000021], [1000004, 1000021]]
1000021 2000002 2.134E+02 # [[1000021, 2000002]]
Python Output¶
The Python-type output is similar to the screen output, however converted to a Python dictionary. If all options are set to True, it includes information about the decomposition, the list of theory predictions and model coverage. The output is printed to the file <input file>.py and stored in the output folder (see the runSModelS options).
Below we describe in detail the dictionary keys and values contained in the Python dictionary output:
- information about the basic input parameters and the status of the run stored under the OutputStatus key:
smodelsOutput = {
'OutputStatus': {'sigmacut': 0.01, 'minmassgap': 5.0, 'maxcond': 0.2, 'ncpus': 1, 'file status': 1, 'decomposition status': 1, 'warnings': 'Input fil
- a full list of the elements generated by the decomposition (if addElementList = True) stored under the Element key. Each list entry contains basic information about the elements. The list can be considerably long, so it is recommended to set addElementList to False, unless the decomposition information is required by the user.
'Element': [{'ID': 1, 'Particles': '[[], []]', 'Masses (GeV)': [[129.0], [129.0]], 'PIDs': [[[1000022], [1000022]]], 'Weights (fb)': {'xsec 8.0 TeV': 0.480, 'xsec 13.0 TeV': 1.580}, 'final states': ['MET', 'MET']},
{'ID': 2, 'Particles': "[[], [['W+']]]", 'Masses (GeV)': [[129.0], [269.0,129.0]], 'PIDs': [[[1000022], [1000024, 1000022]]], 'Weights (fb)': {'xsec 8.0 TeV': 0.2404, 'xsec 13.0 TeV': 0.263}, 'final states': ['MET', 'MET']},..]
- a list of all the theory predictions obtained for the experimental results, stored under the ExptRes key. For each list entry, the corresponding result id, the experimental result type (if UL-type result or EM-type result), the signal region (data set ID), the sqrts and luminosity, the constrained simplified models (txnames), the signal cross section (theory prediction), the corresponding observed upper limit and the maximum condition violation (see upper limit conditions) are shown. The masses of the elements contributing to the signal cross section (if unique) and the \(\chi^2\) likelihood values (if computeStatistics = True) are also included. Finally, if addTxWeights = True, the weight contribution of each txname is also included:
'ExptRes': [{'maxcond': 0.0, 'theory prediction (fb)': 11.68725, 'upper limit (fb)': 7.68,
'expected upper limit (fb)': None, 'TxNames': ['T2'], 'Mass (GeV)': [[991.3084, 129.0], [991.4, 129.0]],
'AnalysisID': 'CMS-SUS-16-036', 'DataSetID': None, 'AnalysisSqrts (TeV)': 13.0,
'lumi (fb-1)': 35.9, 'dataType': 'upperLimit',
'r': 1.52127, 'r_expected': None, 'TxNames weights (fb)': {'T2': 11.687}},...]
- a list of missing topologies (if testCoverage = True), stored under the Missed Topologies key. For each list entry, the element cross section (weight), the element IDs contributing to the topology and the element description in bracket notation and final states is included.
'Missed Topologies': [{'sqrts (TeV)': 13.0, 'weight (fb)': 148.16,
'element': "[[[jet],[W]],[[jet,jet],[W]]](MET,MET)", 'element IDs': [987, 988,...]},...]
- a list of topologies which are outside the experimental results grid (if testCoverage = True), stored under the Outside Grid key. For each list entry, the element cross section (weight) and the element description in bracket notation and final states is included.
'Outside Grid': [{'sqrts (TeV)': 13.0, 'weight (fb)': 1.58, 'element': "[[[jet]],[[t,t]]](MET,MET)"},...]
- a list of topologies with long cascade decays (if testCoverage = True), stored under the Long Cascades key. For each list entry, the element cross section (weight) and the PIDs of the mothers are included. The mother PIDs are given in a nested list, as more than one pair might contribute to the same class of elements.
'Long Cascades': [{'sqrts (TeV)': 13.0, 'weight (fb)': 369.58, 'mother PIDs': [[1000021, 2000002]]},...]
- a list of topologies with asymmetric branch decays (if testCoverage = True), stored under the Asymmetric Branches key. For each list entry, the element cross section (weight) and the PIDs of the mothers are included. The mother PIDs are given in a nested list, as more than one pair might contribute to the same class of elements.
'Asymmetric Branches': [{'sqrts (TeV)': 13.0, 'weight (fb)': 529.73, 'mother PIDs': [[1000021, 1000021]]},...]
XML Output¶
The xml-type output is identical to the python output, however converted to a xml format. The output is printed to the file <input file>.xml and stored in the output folder (see the runSModelS options).
Since the output information and options are the same as described for python output, we simply show below an excerpt of the xml file to illustrate the output format:
<?xml version="1.0" ?>
<smodelsOutput>
<Element_List>
<Element>
<ID>1</ID>
<Masses_GeV_List>
<Masses_GeV_List>
<Masses_GeV>129.0</Masses_GeV>
</Masses_GeV_List>
<Masses_GeV_List>
<Masses_GeV>129.0</Masses_GeV>
</Masses_GeV_List>
</Masses_GeV_List>
<PIDs_List>
<PIDs_List>
<PIDs_List>
<PIDs>1000022</PIDs>
</PIDs_List>
<PIDs_List>
<PIDs>1000022</PIDs>
</PIDs_List>
</PIDs_List>
</PIDs_List>
<Particles>[[], []]</Particles>
<Weights_fb>
<xsec_13.0_TeV>1.5801581999999998</xsec_13.0_TeV>
<xsec_8.0_TeV>0.48082062600000003</xsec_8.0_TeV>
</Weights_fb>
<final_states_List>
<final_states>MET</final_states>
<final_states>MET</final_states>
</final_states_List>
</Element>
<ExptRes_List>
<ExptRes>
<AnalysisID>CMS-SUS-16-036</AnalysisID>
<AnalysisSqrts_TeV>13.0</AnalysisSqrts_TeV>
<DataSetID>None</DataSetID>
<Mass_GeV_List>
<Mass_GeV_List>
<Mass_GeV>991.3084550191742</Mass_GeV>
<Mass_GeV>129.0</Mass_GeV>
</Mass_GeV_List>
<Mass_GeV_List>
<Mass_GeV>991.4592631024188</Mass_GeV>
<Mass_GeV>129.0</Mass_GeV>
</Mass_GeV_List>
</Mass_GeV_List>
<TxNames_List>
<TxNames>T2</TxNames>
</TxNames_List>
<TxNames_weights_fb>
<T2>11.687259769407873</T2>
</TxNames_weights_fb>
<dataType>upperLimit</dataType>
<expected_upper_limit_fb>None</expected_upper_limit_fb>
<lumi_fb-1>35.9</lumi_fb-1>
<maxcond>0.0</maxcond>
<r>1.5212733123965445</r>
<r_expected>None</r_expected>
<theory_prediction_fb>11.687259769407875</theory_prediction_fb>
<upper_limit_fb>7.682550975009415</upper_limit_fb>
</ExptRes>
SLHA Output¶
An SLHA-type output format is also available containing a summary of the theory predictions and missing topologies. The file contains the SLHA-type blocks: SModelS_Settings, SModelS_Exclusion, SModelS_Missing_Topos, SModelS_Outside_Grid, SModelS_Long_Cascade and SModelS_Asymmetric_Branches. Below we give a description of each block together with a sample output.
- information about the main input parameters:
BLOCK SModelS_Settings
0 v1.2.0 #SModelS version
1 1.2.0 #database version
2 0.2 #maximum condition violation
3 1 #compression (0 off, 1 on)
4 5. #minimum mass gap for mass compression [GeV]
5 0.01 #sigmacut [fb]
- information about the status of the input model: excluded (1), not excluded (0) or not tested (-1):
BLOCK SModelS_Exclusion
0 0 1 #output status (-1 not tested, 0 not excluded, 1 excluded)
- followed by the list of experimental results. If the model is excluded, all results with \(r\)-value greater than one are shown. If the point is not excluded, only the result with the highest \(r\)-value is displayed. For each experimental result, the Txname, the \(r\)-value, the condition violation and the experimental result ID are shown. If computeStatistics = True, the \(\chi^2\) and likelihood values for EM-type results are also printed:
1 0 T2 #txname
1 1 1.521E+00 #r value
1 2 1.521E+00 #expected r value
1 3 0.00 #condition violation
1 4 CMS-SUS-16-036 #analysis
1 5 (UL) #signal region
1 6 N/A #Chi2
1 7 N/A #Likelihood
- a list of missing topologies (up to 10) and their weights (if testCoverage = True):
BLOCK SModelS_Missing_Topos #sqrts[TeV] weight[fb] description
0 13 1.482E+02 [[[jet],[W]],[[jet,jet],[W]]] ('MET', 'MET')
1 13 1.445E+02 [[[jet,jet],[W]],[[jet],[jet,jet],[W]]] ('MET', 'MET')
2 13 1.047E+02 [[[b,t],[W]],[[jet,jet],[W]]] ('MET', 'MET')
- a list of topologies which are outside the experimental results grid (if testCoverage = True):
BLOCK SModelS_Outside_Grid #sqrts[TeV] weight[fb] description
0 13 1.590E+00 [[[jet]],[[t,t]]] ('MET', 'MET')
1 13 1.438E+00 [[[t],[W]],[[t],[W]]] ('MET', 'MET')
- a list of topologies with long cascade decays (if testCoverage = True):
BLOCK SModelS_Asymmetric_Branches #Mother1 Mother2 Weight[fb] allMothers
0 1000021 1000021 5.297E+02 [[1000021,1000021]]
1 1000002 1000021 4.149E+02 [[1000002,1000021],[1000004,1000021]]
- a list of topologies with asymmetric branch decays (if testCoverage = True):
9 13 6.214E+01 [[[jet],[W]],[[jet,jet],[higgs]]] ('MET', 'MET')
BLOCK SModelS_Outside_Grid #sqrts[TeV] weight[fb] description
0 13 1.590E+00 [[[jet]],[[t,t]]] ('MET', 'MET')
[*] | Some of the output may change depending on the database version used. |