18 from modularAnalysis
import inputMdst
19 from modularAnalysis
import fillParticleList
20 from modularAnalysis
import reconstructDecay
21 from modularAnalysis
import matchMCTruth
22 from kinfit
import fitKinematic3C
23 from modularAnalysis
import variablesToNtuple
26 from stdPhotons
import stdPhotons
30 mypath = b2.create_path()
32 b2.set_log_level(b2.LogLevel.WARNING)
35 inputMdst(environmentType=
'default',
36 filename=b2.find_file(
'darkphotonmumu_mdst.root',
'examples',
False),
43 fillParticleList(
'mu-:z0',
'chiProb > 0.001 and p > 1.0', path=mypath)
46 reconstructDecay(
'Z0:mm -> mu+:z0 mu-:z0',
'', path=mypath)
50 reconstructDecay(
"beam:rec -> gamma:loose Z0:mm",
"", path=mypath)
51 reconstructDecay(
"beam:kinfit -> gamma:loose Z0:mm",
"", path=mypath)
55 matchMCTruth(
'beam:rec', path=mypath)
56 matchMCTruth(
'beam:kinfit', path=mypath)
59 fitKinematic3C(
'beam:kinfit', path=mypath)
63 mugvars = vc.inv_mass + vc.kinematics + vc.mc_truth + vc.mc_kinematics + vc.momentum_uncertainty
64 z0vars = vc.inv_mass + vc.kinematics + vc.mc_kinematics + vc.mc_truth + \
65 vu.create_aliases_for_selected(mugvars,
'beam -> ^gamma [^Z0 -> ^mu+ ^mu-]')
68 vu.create_aliases([
'OrcaKinFitProb',
70 'OrcaKinFitErrorCode'],
'extraInfo(variable)',
"")
73 output_file =
'B2A422-Orcakinfit_3CFit.root'
75 filename=output_file, treename=
'Z0_mm_rec', path=mypath)
77 filename=output_file, treename=
'Z0_mm_kinfit', path=mypath)