21 import modularAnalysis
as ma
24 sys_tests = [
'def',
'scale',
'scaleUp',
'scaleDown',
'Efficiency']
29 if test
not in sys_tests:
30 print(
"Unknown systematic test {}".format(test))
34 my_path = b2.create_path()
37 ma.inputMdst(environmentType=
'default',
38 filename=b2.find_file(
'JPsi2ee_e2egamma.root',
'examples',
False),
43 ma.fillParticleList(decayString=
'e+:sel',
44 cut=
'electronID > 0.2 and dr < 2 and abs(dz) < 4',
50 ma.trackingMomentum(inputListNames=[
'e+:sel'], scale=1.00056, path=my_path)
51 elif test ==
'scaleUp':
52 ma.trackingMomentum(inputListNames=[
'e+:sel'], scale=1.001, path=my_path)
53 elif test ==
'scaleDown':
54 ma.trackingMomentum(inputListNames=[
'e+:sel'], scale=0.999, path=my_path)
55 elif test ==
'Efficiency':
56 ma.trackingEfficiency(inputListNames=[
'e+:sel'], fraction=0.01, path=my_path)
59 ma.reconstructDecay(decayString=
'J/psi:out -> e+:sel e-:sel',
68 ma.variablesToNtuple(decayString=
'J/psi:out',
70 filename=
'Jpsi_out_{}.root'.format(test),