18 from modularAnalysis
import inputMdst
19 from modularAnalysis
import fillParticleList
20 from modularAnalysis
import reconstructDecay
21 from modularAnalysis
import matchMCTruth
22 from kinfit
import fitKinematic4C
23 from modularAnalysis
import variablesToNtuple
28 mypath = b2.create_path()
30 b2.set_log_level(b2.LogLevel.WARNING)
33 inputMdst(environmentType=
'default',
34 filename=b2.find_file(
'B2A424-SimulateMuonPairs.root',
'examples',
False),
38 fillParticleList(
'mu-:z0',
'chiProb > 0.001 and p > 1.0', path=mypath)
39 fillParticleList(
'mu-:z0fit',
'chiProb > 0.001 and p > 1.0', path=mypath)
42 reconstructDecay(
'Z0:mm_rec -> mu+:z0 mu-:z0',
'9.0 < M < 11.0', path=mypath)
43 reconstructDecay(
'Z0:mm_kinfit -> mu+:z0fit mu-:z0fit',
'9.0 < M < 11.0', path=mypath)
46 matchMCTruth(
'Z0:mm_rec', path=mypath)
47 matchMCTruth(
'Z0:mm_kinfit', path=mypath)
50 fitKinematic4C(
'Z0:mm_kinfit', path=mypath)
53 muvars = vc.kinematics + vc.mc_truth + vc.mc_kinematics + vc.momentum_uncertainty
54 z0vars = vc.inv_mass + vc.kinematics + vc.mc_kinematics + vc.mc_truth + \
55 vu.create_aliases_for_selected(muvars,
'Z0 -> ^mu+ ^mu-')
58 vu.create_aliases([
'OrcaKinFitProb',
'OrcaKinFitChi2',
'OrcaKinFitErrorCode'],
'extraInfo(variable)',
"")
62 output_file =
'B2A420-OrcaKinFit.root'
64 filename=output_file, treename=
'Z0_mm_rec', path=mypath)
66 filename=output_file, treename=
'Z0_mm_kinfit', path=mypath)