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)