20 import modularAnalysis
as ma
23 import stdCharged
as stdc
26 my_path = b2.create_path()
29 ma.inputMdst(environmentType=
'default',
30 filename=b2.find_file(
'Dst2D0pi.root',
'examples',
False),
37 stdc.stdPi(listtype=
'all', path=my_path)
39 stdc.stdPi(listtype=
'loose', path=my_path)
41 stdc.stdK(listtype=
'loose', path=my_path)
45 ma.reconstructDecay(decayString=
'D0:kpi -> K-:loose pi+:loose', cut=
'1.8 < M < 1.9', path=my_path)
49 ma.reconstructDecay(decayString=
'D*+ -> D0:kpi pi+:all', cut=
'0.0 < Q < 0.2', path=my_path)
52 ma.matchMCTruth(list_name=
'D*+', path=my_path)
55 dstar_vars = vc.inv_mass + vc.mc_truth
68 fs_hadron_vars = vu.create_aliases_for_selected(
69 list_of_variables=vc.pid + vc.track + vc.mc_truth,
70 decay_string=
'D*+ -> [D0 -> ^K- ^pi+] ^pi+')
72 d0_vars = vu.create_aliases_for_selected(
73 list_of_variables=vc.inv_mass + vc.mc_truth,
74 decay_string=
'D*+ -> ^D0 pi+',
79 output_file =
'B2A301-Dstar2D0Pi-Reconstruction.root'
80 ma.variablesToNtuple(
'D*+', dstar_vars + d0_vars + fs_hadron_vars,
81 filename=output_file, treename=
'dsttree', path=my_path)