30 import modularAnalysis
as ma
33 my_path = b2.create_path()
36 ma.inputMdst(environmentType=
'default',
37 filename=b2.find_file(
'JPsi2ee_e2egamma.root',
'examples',
False),
42 ma.fillParticleList(decayString=
'e+:uncorrected',
43 cut=
'electronID > 0.2 and dr < 2 and abs(dz) < 4',
45 ma.fillParticleList(decayString=
'gamma:bremsinput',
50 ma.matchMCTruth(list_name=
'e+:uncorrected', path=my_path)
51 ma.matchMCTruth(list_name=
'gamma:bremsinput', path=my_path)
55 ma.correctBrems(outputList=
'e+:corrected',
56 inputList=
'e+:uncorrected',
57 gammaList=
'gamma:bremsinput',
59 ma.matchMCTruth(list_name=
'e+:corrected',
63 ma.reconstructDecay(decayString=
'J/psi:uncorrected -> e+:uncorrected e-:uncorrected',
66 ma.reconstructDecay(decayString=
'J/psi:corrected -> e+:corrected e-:corrected',
71 ma.matchMCTruth(list_name=
'J/psi:uncorrected', path=my_path)
72 ma.matchMCTruth(list_name=
'J/psi:corrected', path=my_path)
75 ma.fillParticleListFromMC(decayString=
'J/psi:MC', cut=
"", path=my_path)
100 'extraInfo(bremsCorrected)'
110 'daughter(0, extraInfo(bremsCorrected))',
111 'daughter(1, extraInfo(bremsCorrected))'
114 ma.variablesToNtuple(decayString=
'e+:uncorrected',
116 filename=
'e_uncorrected.root',
118 ma.variablesToNtuple(decayString=
'e+:corrected',
120 filename=
'e_corrected.root',
122 ma.variablesToNtuple(decayString=
'J/psi:uncorrected',
124 filename=
'Jpsi_uncorrected.root',
126 ma.variablesToNtuple(decayString=
'J/psi:corrected',
128 filename=
'Jpsi_corrected.root',
130 ma.variablesToNtuple(decayString=
'J/psi:MC',
132 filename=
'Jpsi_MC.root',