Belle II Software light-2406-ragdoll
neutralHadron4MomentumCalculation.py
1
8
9from basf2 import Path, process, B2INFO, statistics, find_file
10from stdCharged import stdK, stdPi, stdPr
11import modularAnalysis as ma
12
13mypath = Path()
14ma.inputMdst(
15 environmentType='default',
16 filename=find_file(
17 'B02Dpnbar_D2kpipi.root',
18 data_type='examples',
19 silent=False),
20 path=mypath)
21
22stdK('higheff', path=mypath)
23stdPi('higheff', path=mypath)
24stdPr('higheff', path=mypath)
25ma.reconstructDecay('D-:sig -> K+:higheff pi-:higheff pi-:higheff', 'abs(dM) < 0.05', path=mypath)
26ma.fillParticleList('anti-n0:good', 'isFromECL > 0', path=mypath)
27ma.reconstructDecayWithNeutralHadron('B0:sig -> D-:sig p+:higheff ^anti-n0:good', '', chargeConjugation=False, path=mypath)
28
29fname = 'neutralHadron4MomentumCalculation.root'
30ma.variablesToNtuple('B0:sig', ['deltaE', 'M', 'daughter(2,p)', 'daughter(2,clusterE)'], treename='B0', filename=fname, path=mypath)
31ma.variablesToNtuple('anti-n0:good', ['p'], treename='nbar', filename=fname, path=mypath)
32
33process(mypath)
34B2INFO(statistics)