TestTreeFits Class Reference

Inheritance diagram for TestTreeFits:

Collaboration diagram for TestTreeFits:

Public Member Functions
def	testFit (self)

Detailed Description

The unit test case for TreeFitter

Definition at line 19 of file treeFit_referenced_ignore_charged.py.

Member Function Documentation

testFit()

def testFit

(

self

)

Run the test fit

Definition at line 22 of file treeFit_referenced_ignore_charged.py.

    def testFit(self):
        """Run the test fit"""
 
        testFile = tempfile.NamedTemporaryFile()
 
        main = basf2.create_path()
 
        inputfile = b2test_utils.require_file('analysis/tests/100_noBKG_B0ToPiPiPi0.root', py_case=self)
        ma.inputMdst(inputfile, path=main)
 
        ma.fillParticleList('pi+:a', 'pionID > 0.5', path=main)
 
        ma.fillParticleList('gamma:a', '', path=main)
        ma.reconstructDecay('pi0:a -> gamma:a gamma:a', '0.125 < InvM < 0.145', 0, path=main)
 
        ma.reconstructDecay('B0:rec -> pi-:a pi+:a pi0:a', '', 0, path=main)
        ma.matchMCTruth('B0:rec', path=main)
 
        conf = 0
        main.add_module('TreeFitter',
                        particleList='B0:rec',
                        confidenceLevel=conf,
                        massConstraintList=[],
                        massConstraintListParticlename=[],
                        expertUseReferencing=True,
                        ipConstraint=True,
                        updateAllDaughters=True,
                        ignoreFromVertexFit='B0 -> ^pi- pi+ pi0')
 
        ntupler = basf2.register_module('VariablesToNtuple')
        ntupler.param('fileName', testFile.name)
        ntupler.param('variables', ['chiProb', 'M', 'isSignal', 'Mbc'])
        ntupler.param('particleList', 'B0:rec')
        main.add_module(ntupler)
 
        basf2.process(main)
 
        ntuplefile = TFile(testFile.name)
        ntuple = ntuplefile.Get('ntuple')
 
        self.assertFalse(ntuple.GetEntries() == 0, "Ntuple is empty.")
 
        allBkg = ntuple.GetEntries("isSignal == 0")
        allSig = ntuple.GetEntries("isSignal > 0")
        SigMbcReasonable = ntuple.GetEntries("isSignal > 0 && Mbc > 5.27")
 
        truePositives = ntuple.GetEntries("(chiProb > 0) && (isSignal > 0)")
        falsePositives = ntuple.GetEntries("(chiProb > 0) && (isSignal == 0)")
 
        mustBeZero = ntuple.GetEntries(f"(chiProb < {conf})")
 
        print(f"True fit survivors: {truePositives} out of {allSig} true candidates")
        print(f"False fit survivors: {falsePositives} out of {allBkg} false candidates")
 
        self.assertTrue(falsePositives < 2156, "Background rejection too small.")
 
        self.assertTrue(truePositives > 34, "Signal rejection too high")
        self.assertFalse(mustBeZero, f"We should have dropped all candidates with confidence level less than {conf}.")
 
        self.assertTrue(SigMbcReasonable > 33, "Signal kinematics is wrongly reconstructed too much")
 
        print("Test passed, cleaning up.")
 
 

---

The documentation for this class was generated from the following file:

• analysis/tests/treeFit_referenced_ignore_charged.py