26 from modularAnalysis
import inputMdst
27 from modularAnalysis
import reconstructDecay
28 from modularAnalysis
import matchMCTruth
29 from vertex
import kFit
30 from modularAnalysis
import variablesToNtuple
31 from stdPi0s
import stdPi0s
36 my_path = b2.create_path()
39 inputMdst(environmentType=
'default',
40 filename=b2.find_file(
'B02pi0D0_D2kpi_B2Dstarpi_Dstar2Dpi_D2kpi.root',
'examples',
False),
46 stdPi0s(
'eff40_Jan2020Fit', path=my_path)
50 reconstructDecay(
'D0:pi0pi0 -> pi0:eff40_Jan2020Fit pi0:eff40_Jan2020Fit',
'1.7 < M < 2.0', path=my_path)
54 kFit(
'D0:pi0pi0', 0.0,
'mass', path=my_path)
59 reconstructDecay(
'B0:all -> D0:pi0pi0 pi0:eff40_Jan2020Fit',
'5.24 < Mbc < 5.29 and abs(deltaE) < 1.0', path=my_path)
62 matchMCTruth(
'B0:all', path=my_path)
65 B0_vars = vc.inv_mass + vc.mc_truth + \
66 vu.create_aliases_for_selected(
67 vc.inv_mass + vc.mc_truth,
69 vu.create_aliases_for_selected(
70 vc.cluster,
'B0 -> D0 [pi0 -> ^gamma ^gamma]')
72 pi0_vars = vc.mc_truth + vc.kinematics + \
73 [
'extraInfo(BDT)',
'decayAngle(0)'] + \
74 vu.create_aliases_for_selected(
75 vc.cluster + vc.kinematics,
'pi0 -> ^gamma ^gamma')
78 output_file =
'B2A401-KFit-MassFit.root'
80 filename=output_file, treename=
'b0', path=my_path)
82 filename=output_file, treename=
'pi0', path=my_path)