29 from modularAnalysis
import inputMdst
30 from modularAnalysis
import reconstructDecay
31 from modularAnalysis
import matchMCTruth
32 from vertex
import raveFit
33 from stdCharged
import stdPi, stdK
34 from modularAnalysis
import variablesToNtuple
39 my_path = b2.create_path()
42 inputMdst(environmentType=
'default',
43 filename=b2.find_file(
'B02pi0D0_D2kpi_B2Dstarpi_Dstar2Dpi_D2kpi.root',
'examples',
False),
50 stdPi(
'all', path=my_path)
52 stdPi(
'loose', path=my_path)
54 stdK(
'loose', path=my_path)
58 reconstructDecay(
'D0:kpi -> K-:loose pi+:loose',
'1.8 < M < 1.9', path=my_path)
59 reconstructDecay(
'D0:st -> K-:loose pi+:loose',
'1.8 < M < 1.9', path=my_path)
60 reconstructDecay(
'D0:du -> K-:loose pi+:loose',
'1.8 < M < 1.9', path=my_path)
64 raveFit(
'D0:kpi', 0.0, fit_type=
'massvertex', path=my_path)
68 raveFit(
'D0:st', 0.0, path=my_path)
70 raveFit(
'D0:st', 0.0, decay_string=
'^D0 -> K- pi+', constraint=
'ipprofile', path=my_path)
74 raveFit(
'D0:du', 0.0, daughtersUpdate=
True, path=my_path)
79 reconstructDecay(
'D*+:1 -> D0:kpi pi+:all',
80 '0.0 <= Q < 0.02 and 2.5 < useCMSFrame(p) < 5.5', path=my_path)
81 reconstructDecay(
'D*+:2 -> D0:kpi pi+:all',
82 '0.0 <= Q < 0.02 and 2.5 < useCMSFrame(p) < 5.5', path=my_path)
83 reconstructDecay(
'D*+:3 -> D0:kpi pi+:all',
84 '0.0 <= Q < 0.02 and 2.5 < useCMSFrame(p) < 5.5', path=my_path)
87 matchMCTruth(
'D*+:1', path=my_path)
88 matchMCTruth(
'D*+:2', path=my_path)
89 matchMCTruth(
'D*+:3', path=my_path)
93 raveFit(
'D*+:1', 0.0, path=my_path)
97 raveFit(
'D*+:2', 0.0, constraint=
'ipprofile', path=my_path)
101 raveFit(
'D*+:3', 0.0, decay_string=
'D*+ -> D0 ^pi+', constraint=
'ipprofile', path=my_path)
104 dstar_vars = vc.inv_mass + vc.mc_truth + \
105 vc.mc_flight_info + vc.flight_info + vc.vertex
107 fs_hadron_vars = vu.create_aliases_for_selected(
108 vc.pid + vc.track + vc.mc_truth,
109 'D*+ -> [D0 -> ^K- ^pi+] ^pi+')
111 d0_vars = vu.create_aliases_for_selected(
112 vc.inv_mass + vc.mc_truth + vc.vertex,
113 'D*+ -> ^D0 pi+',
'D0')
115 dstt = vc.kinematics + vc.vertex + vc.mc_vertex + vc.flight_info + \
116 vu.create_aliases_for_selected(
120 dstu = vc.kinematics + vu.create_aliases_for_selected(
125 output_file =
'B2A406-Rave-DecayStringVertexFit.root'
127 filename=output_file, treename=
'dsttree1', path=my_path)
129 filename=output_file, treename=
'dsttree2', path=my_path)
131 filename=output_file, treename=
'dsttree3', path=my_path)
133 filename=output_file, treename=
'd0tree1', path=my_path)
135 filename=output_file, treename=
'd0tree2', path=my_path)