29 import modularAnalysis
as ma
34 my_path = b2.create_path()
37 ma.inputMdst(environmentType=
'default',
38 filename=b2.find_file(
'B2pi0D_D2hh_D2hhh_B2munu.root',
'examples',
False),
44 photons = (
'gamma',
'')
45 ma.fillParticleListsFromMC(decayStringsWithCuts=[kaons, pions, photons],
50 ma.reconstructDecay(decayString=
'pi0 -> gamma gamma',
51 cut=
'0.1 < InvM < 0.15',
55 ma.reconstructDecay(decayString=
'D0:ch1 -> K- pi+',
59 ma.reconstructDecay(decayString=
'D0:ch2 -> K- pi+ pi0',
63 ma.reconstructDecay(decayString=
'D0:ch3 -> K- pi+ pi+ pi-',
67 ma.reconstructDecay(decayString=
'D0:ch4 -> K- K+',
71 ma.reconstructDecay(decayString=
'D0:ch5 -> pi+ pi-',
77 ma.copyLists(outputListName=
'D0:all',
78 inputListNames=[
'D0:ch1',
'D0:ch2',
'D0:ch3',
'D0:ch4',
'D0:ch5'],
82 ma.reconstructDecay(decayString=
'B+:D0pi -> anti-D0:all pi+',
83 cut=
'5.24 < Mbc < 5.29 and abs(deltaE) < 1.0',
88 ma.matchMCTruth(list_name=
'B+:D0pi',
92 d_vars = vc.inv_mass + vc.kinematics
93 pi_vars = vc.kinematics
94 b_vars = vc.deltae_mbc + \
96 vu.create_aliases_for_selected(list_of_variables=d_vars,
97 decay_string=
'B+ -> ^anti-D0 pi+',
99 vu.create_aliases_for_selected(list_of_variables=pi_vars,
100 decay_string=
'B+ -> anti-D0 ^pi+',
102 vu.create_aliases(list_of_variables=[
'decayModeID'],
103 wrapper=
'daughter(0,extraInfo({variable}))',
107 output_file =
'B2A303-MultipleDecays-Reconstruction.root'
108 ma.variablesToNtuple(decayString=
'B+:D0pi',
110 filename=output_file,