17 from modularAnalysis
import inputMdst
18 from modularAnalysis
import fillParticleList
19 from modularAnalysis
import reconstructDecay
20 from modularAnalysis
import matchMCTruth
21 from kinfit
import UnmeasuredfitKinematic1C
22 from modularAnalysis
import variablesToNtuple
27 mypath = b2.create_path()
29 b2.set_log_level(b2.LogLevel.WARNING)
32 inputMdst(environmentType=
'default',
33 filename=b2.find_file(
'darkphotonmumu_mdst.root',
'examples',
False),
37 fillParticleList(
'mu-:z0',
'chiProb > 0.001 and p > 1.0', path=mypath)
38 fillParticleList(
'mu-:z0fit',
'chiProb > 0.001 and p > 1.0', path=mypath)
41 reconstructDecay(
'Z0:mm_rec -> mu+:z0 mu-:z0',
'9.0 < M < 11.0', path=mypath)
42 reconstructDecay(
'Z0:mm_kinfit -> mu+:z0fit mu-:z0fit',
'9.0 < M < 11.0', path=mypath)
45 matchMCTruth(
'Z0:mm_rec', path=mypath)
46 matchMCTruth(
'Z0:mm_kinfit', path=mypath)
49 UnmeasuredfitKinematic1C(
'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',
60 'OrcaKinFitErrorCode',
61 'OrcaKinFitUnmeasuredTheta',
62 'OrcaKinFitUnmeasuredPhi',
63 'OrcaKinFitUnmeasuredE'],
'extraInfo(variable)',
"")
66 output_file =
'B2A421-Orcakinfit_1CFit.root'
68 filename=output_file, treename=
'Z0_mm_rec', path=mypath)
70 filename=output_file, treename=
'Z0_mm_kinfit', path=mypath)