12 from modularAnalysis
import *
13 from variables
import variables
14 from ROOT
import gSystem
18 ap = argparse.ArgumentParser()
19 ap.add_argument(
"--mc", help=
"Input file is MC", action=
'store_true', default=
False)
20 ap.add_argument(
"--vfit", help=
"Do vertex fit", action=
'store_true', default=
False)
21 ap.add_argument(
"--tag", help=
"Global tag (data only)", default=
"data_reprocessing_proc7")
22 ap.add_argument(
"--out", help=
"Output file", default=
"Dstar.root")
23 args = ap.parse_args()
29 use_central_database(args.tag, LogLevel.WARNING)
32 gSystem.Load(
'libtop.so')
35 main = basf2.create_path()
38 inputMdstList(
'default',
'Input.root', path=main)
41 fillParticleList(
'K-:all',
'-2.0 < d0 < 2.0 and -4.0 < z0 < 4.0', path=main)
42 fillParticleList(
'pi+:all',
'-2.0 < d0 < 2.0 and -4.0 < z0 < 4.0', path=main)
45 reconstructDecay(
'D0:kpi -> K-:all pi+:all',
'1.8 < M < 1.95', path=main)
47 fitVertex(
'D0:kpi', 0.001,
'',
'rave',
'vertex',
'',
True, path=main)
50 reconstructDecay(
'D*+ -> D0:kpi pi+:all',
'0.0 < Q < 0.020', path=main)
52 fitVertex(
'D*+', 0.001,
'',
'rave',
'vertex',
'ipprofile',
True, path=main)
56 matchMCTruth(
'D*+', path=main)
59 variables.addAlias(
'isBunch',
'isTopRecBunchReconstructed')
60 variables.addAlias(
'bunchOffset',
'topRecBunchCurrentOffset')
61 variables.addAlias(
'M_D0',
'daughter(0, M)')
62 variables.addAlias(
'p_cms',
'useCMSFrame(p)')
64 variables.addAlias(
'p_K',
'daughter(0, daughter(0, p))')
65 variables.addAlias(
'flag_K',
'daughter(0, daughter(0, topFlag))')
66 variables.addAlias(
'kaonLL_K',
'daughter(0, daughter(0, topKaonLogL))')
67 variables.addAlias(
'pionLL_K',
'daughter(0, daughter(0, topPionLogL))')
68 variables.addAlias(
'slotID_K',
'daughter(0, daughter(0, topSlotID))')
70 variables.addAlias(
'p_pi',
'daughter(0, daughter(1, p))')
71 variables.addAlias(
'flag_pi',
'daughter(0, daughter(1, topFlag))')
72 variables.addAlias(
'kaonLL_pi',
'daughter(0, daughter(1, topKaonLogL))')
73 variables.addAlias(
'pionLL_pi',
'daughter(0, daughter(1, topPionLogL))')
74 variables.addAlias(
'slotID_pi',
'daughter(0, daughter(1, topSlotID))')
76 varlist = [
'M_D0',
'Q',
'dQ',
'p_cms',
'isBunch',
'bunchOffset',
'p_K',
'p_pi',
77 'kaonLL_K',
'pionLL_K',
'flag_K',
'slotID_K',
78 'kaonLL_pi',
'pionLL_pi',
'flag_pi',
'slotID_pi']
86 print(basf2.statistics)