|
| 1 | +#!/usr/bin/env python3 |
| 2 | + |
| 3 | +""" |
| 4 | +.. module:: testCrossSectionComp |
| 5 | + :synopsis: Tests xsec comparison. |
| 6 | +
|
| 7 | +.. moduleauthor:: Andre Lessa <lessa.a.p@gmail.com> |
| 8 | +
|
| 9 | +""" |
| 10 | + |
| 11 | +import sys |
| 12 | +sys.path.insert(0,"../") |
| 13 | +import unittest |
| 14 | +from smodels.base.crossSection import XSection,XSectionInfo,XSectionList |
| 15 | +from smodels.share.models.mssm import BSMList |
| 16 | +from smodels.share.models.SMparticles import SMList |
| 17 | +from smodels.base.model import Model |
| 18 | +from smodels.base.physicsUnits import fb, pb |
| 19 | + |
| 20 | + |
| 21 | +class XSecCompTest(unittest.TestCase): |
| 22 | + |
| 23 | + def testCrossSectionComp(self): |
| 24 | + |
| 25 | + slhafile = './testFiles/slha/lightEWinos.slha' |
| 26 | + model = Model(BSMparticles=BSMList, SMparticles=SMList) |
| 27 | + model.updateParticles(inputFile=slhafile, |
| 28 | + ignorePromptQNumbers=['spin','eCharge','colordim']) |
| 29 | + |
| 30 | + xsecA = model.xsections[0] |
| 31 | + xsecB = model.xsections[1] |
| 32 | + self.assertEqual(xsecA.info.label,'8 TeV (LO)') |
| 33 | + self.assertEqual(xsecB.info.label,'13 TeV (LO)') |
| 34 | + |
| 35 | + self.assertEqual(xsecA.pid,(-1000037, 1000023)) |
| 36 | + self.assertEqual(xsecB.pid,(-1000037, 1000023)) |
| 37 | + |
| 38 | + self.assertTrue(xsecA < xsecB) |
| 39 | + self.assertFalse(xsecA > 5*fb) |
| 40 | + self.assertTrue(xsecA > 4e-3) |
| 41 | + |
| 42 | + maxXsec = model.xsections.getMaxXsec() |
| 43 | + self.assertAlmostEqual(maxXsec.asNumber(pb),11.3,1) |
| 44 | + |
| 45 | + self.assertFalse(model.xsections > 15) |
| 46 | + self.assertTrue(model.xsections > 15.*fb) |
| 47 | + self.assertFalse(model.xsections < xsecB) |
| 48 | + |
| 49 | +if __name__ == "__main__": |
| 50 | + unittest.main() |
0 commit comments