25 from modularAnalysis
import inputMdst
26 from modularAnalysis
import fillParticleList
27 from modularAnalysis
import reconstructDecay
28 from modularAnalysis
import matchMCTruth
29 from kinfit
import fitKinematic4C
30 from modularAnalysis
import variablesToNtuple
33 from stdPhotons
import stdPhotons
36 mypath = b2.create_path()
39 inputMdst(environmentType=
'default',
40 filename=b2.find_file(
'Y4SEventToetaY1S-evtgen_100.root',
'examples',
False),
45 fillParticleList(
'mu+:pid',
'chiProb > 0.001 and p > 1.0', path=mypath)
48 reconstructDecay(
"eta:gg -> gamma:loose gamma:loose",
"", path=mypath)
50 reconstructDecay(
"Upsilon:uu -> mu+:pid mu-:pid",
"M>2.", path=mypath)
53 reconstructDecay(
"Upsilon(4S) -> eta:gg Upsilon:uu",
"", path=mypath)
54 reconstructDecay(
"Upsilon(4S):4c -> eta:gg Upsilon:uu",
"", path=mypath)
57 matchMCTruth(
'Upsilon(4S)', path=mypath)
58 matchMCTruth(
'Upsilon(4S):4c', path=mypath)
61 fitKinematic4C(
"Upsilon(4S):4c", path=mypath)
64 muvars = vc.mc_truth + vc.pid + vc.kinematics
65 gvars = vc.kinematics + vc.mc_truth + vc.inv_mass
66 etaanduvars = vc.inv_mass + vc.kinematics + vc.mc_truth
67 u4svars = vc.inv_mass + vc.kinematics + vc.mc_truth +\
68 vu.create_aliases([
'FourCFitProb',
'FourCFitChi2'],
'extraInfo(variable)',
"") + \
69 vu.create_aliases_for_selected(etaanduvars,
'Upsilon(4S) -> ^eta ^Upsilon') + \
70 vu.create_aliases_for_selected(muvars,
'Upsilon(4S) -> eta [Upsilon -> ^mu+ ^mu-]') + \
71 vu.create_aliases_for_selected(gvars,
'Upsilon(4S) -> [eta -> ^gamma ^gamma] Upsilon')
73 u4svars_4c = u4svars + vu.create_aliases([
'OrcaKinFitProb',
75 'OrcaKinFitErrorCode'],
'extraInfo(variable)',
"")
77 u4svars_def = u4svars + vu.create_aliases([
'chiProb'],
'extraInfo(variable)',
"")
81 output_file =
'B2A423-Orcakinfit_4CFit.root'
83 filename=output_file, treename=
'Upsilon4s', path=mypath)
85 filename=output_file, treename=
'Upsilon4s_4c', path=mypath)