30 from modularAnalysis
import inputMdst
31 from modularAnalysis
import reconstructDecay
32 from modularAnalysis
import matchMCTruth
33 from vertex
import kFit
34 from stdCharged
import stdPi, stdK
35 from modularAnalysis
import variablesToNtuple
40 my_path = b2.create_path()
43 inputMdst(environmentType=
'default',
44 filename=b2.find_file(
'B02pi0D0_D2kpi_B2Dstarpi_Dstar2Dpi_D2kpi.root',
'examples',
False),
51 stdPi(
'all', path=my_path)
53 stdPi(
'loose', path=my_path)
55 stdK(
'loose', path=my_path)
59 reconstructDecay(
'D0:kpi -> K-:loose pi+:loose',
'1.8 < M < 1.9', path=my_path)
63 kFit(
'D0:kpi', 0.0, path=my_path)
68 reconstructDecay(
'D*+ -> D0:kpi pi+:all',
'0.0 <= Q < 0.02 and 2.5 < useCMSFrame(p) < 5.5', path=my_path)
71 matchMCTruth(
'D*+', path=my_path)
75 kFit(
'D*+', 0.0, path=my_path)
79 dstar_vars = vc.inv_mass + vc.mc_truth
81 fs_hadron_vars = vu.create_aliases_for_selected(
82 vc.pid + vc.track + vc.mc_truth,
83 'D*+ -> [D0 -> ^K- ^pi+] ^pi+')
85 d0_vars = vu.create_aliases_for_selected(
86 vc.inv_mass + vc.mc_truth,
87 'D*+ -> ^D0 pi+',
'D0')
91 output_file =
'B2A403-KFit-VertexFit.root'
93 filename=output_file, treename=
'dsttree', path=my_path)