def test_z_1_cms(self):
'''Check Z mass in ee->Z->ddbar (CMS).
Will fail if physics algorithms are modified,
so should probably be removed from test suite,
or better: be made optional.
'''
from heppy.papas.detectors.CMS import cms
fname = '/'.join([os.environ['HEPPY'], 'test/data/ee_Z_ddbar.root'])
config.components[0].files = [fname]
for s in config.sequence:
if hasattr(s, 'detector'):
s.detector = cms
self.looper = Looper(self.outdir, config, nEvents=100, nPrint=0)
self.looper.loop()
self.looper.write()
rootfile = '/'.join([
self.outdir,
'heppy.analyzers.GlobalEventTreeProducer.GlobalEventTreeProducer_1/tree.root'
])
mean, sigma = plot_ee_mass(rootfile)
self.assertAlmostEqual(mean, 83.58, 1)
self.assertAlmostEqual(sigma, 7.06, 1)
def test_z_mumu_clic(self):
'''Check Z mass in ee->Z->mumu (CLIC).
Will fail if physics algorithms are modified,
so should probably be removed from test suite,
or better: be made optional.
'''
from heppy.papas.detectors.CLIC import clic
fname = '/'.join([os.environ['HEPPY'], 'test/data/ee_Z_mumu.root'])
config.components[0].files = [fname]
for s in config.sequence:
if hasattr(s, 'detector'):
s.detector = clic
self.looper = Looper(self.outdir, config, nEvents=500, nPrint=0)
self.looper.loop()
self.looper.write()
rootfile = '/'.join([
self.outdir,
'heppy.analyzers.GlobalEventTreeProducer.GlobalEventTreeProducer_1/tree.root'
])
mean, sigma = plot_ee_mass(rootfile, nbins=400, xmin=70, xmax=110)
self.assertAlmostEqual(mean, 90.84, 1)
self.assertAlmostEqual(sigma, 1.32, 1)
def test_z_2_clic(self):
'''Check Z mass in ee->Z->ddbar (CLIC).
Will fail if physics algorithms are modified,
so should probably be removed from test suite,
or better: be made optional.
'''
from heppy.papas.detectors.CLIC import clic
fname = '/'.join([os.environ['HEPPY'],
'test/data/ee_Z_ddbar.root'])
config.components[0].files = [fname]
for s in config.sequence:
if hasattr( s,'detector'):
s.detector = clic
self.looper = Looper( self.outdir, config,
nEvents=100,
nPrint=0 )
self.looper.loop()
self.looper.write()
rootfile = '/'.join([self.outdir,
'heppy.analyzers.GlobalEventTreeProducer.GlobalEventTreeProducer_1/tree.root'])
mean, sigma = plot_ee_mass(rootfile)
self.assertAlmostEqual(mean, 90.15, 1)
self.assertAlmostEqual(sigma, 3.3, 1)
