allparticlecombiner.py

#!/usr/bin/env python3
# -*- coding: utf-8 -*-
 
 
import unittest
import tempfile
import basf2
import b2test_utils
import modularAnalysis as ma
from variables import variables
from vertex import kFit
from ROOT import TFile
 
 
class TestAllParticleCombiner(unittest.TestCase):
    """The unit test case for the AllParticleCombiner"""
 
    def testCombination(self):
        """Run the test fit"""
 
        testFile = tempfile.NamedTemporaryFile()
        # make logging more reproducible by replacing some strings.
        # Also make sure the testfile name is replaced if necessary
        b2test_utils.configure_logging_for_tests({testFile.name: "${testfile}"})
 
        basf2.set_random_seed("1234")
        main = basf2.create_path()
        inputfile = b2test_utils.require_file(
            'analysis/tests/mdst.root', py_case=self)
        main.add_module(
            'RootInput',
            inputFileNames=[inputfile],
            logLevel=basf2.LogLevel.ERROR)
 
        ma.fillParticleList('pi+:fromPV', 'dr < 2 and abs(dz) < 5', path=main)
 
        ma.variablesToExtraInfo('pi+:fromPV', {'medianValueInList(pi+:fromPV, dz)': 'medianDZ'}, path=main)
        variables.addAlias('medianDzFromPV', 'extraInfo(medianDZ)')
        variables.addAlias('dzFromMedianDz', 'formula(dz - medianDzFromPV)')
        ma.applyCuts('pi+:fromPV', 'abs(dzFromMedianDz) < 0.05', path=main)
 
        ma.combineAllParticles(['pi+:fromPV'], 'vpho:PVFit', path=main)
 
        kFit('vpho:PVFit', conf_level=-1, fit_type='vertex', constraint='iptube', path=main)
 
        variables.addAlias('PVX', 'x')
        variables.addAlias('PVY', 'y')
        variables.addAlias('PVZ', 'z')
        variables.addAlias('nGoodTracksFromPV', 'nParticlesInList(pi+:fromPV)')
 
        main.add_module(
            'VariablesToNtuple',
            particleList='vpho:PVFit',
            fileName=testFile.name,
            variables=[
                'PVX',
                'PVY',
                'PVZ',
                'chiProb',
                'nGoodTracksFromPV',
                'distance'])
 
        ma.summaryOfLists(['vpho:PVFit'], path=main)
 
        basf2.process(main)
 
        ntuplefile = TFile(testFile.name)
        ntuple = ntuplefile.Get('ntuple')
 
        self.assertFalse(ntuple.GetEntries() == 0, "Ntuple is empty.")
 
        converged = ntuple.GetEntries("chiProb > 0")
 
        self.assertFalse(converged == 0, "No fit converged.")
 
        print("Test passed, cleaning up.")
 
 
if __name__ == '__main__':
    unittest.main()