Belle II Software  release-05-01-25
test2_BtoDh_Kshh.py
1 #!/usr/bin/env/python3
2 # -*-coding: utf-8-*-
3 
4 """
5 <header>
6  <input>14140200.udst.root</input>
7  <output>BtoDh_Kshh_Validation.root</output>
8  <contact>niharikarout@physics.iitm.ac.in</contact>
9 </header>
10 """
11 
12 import basf2 as b2
13 import modularAnalysis as ma
14 from stdCharged import stdK, stdPi
15 from stdV0s import stdKshorts
16 
17 kshhPath = b2.Path()
18 
19 ma.inputMdst('default', '14140200.udst.root', path=kshhPath)
20 
21 stdPi('all', path=kshhPath)
22 stdK('all', path=kshhPath)
23 stdKshorts(path=kshhPath)
24 
25 ma.reconstructDecay('D0 -> K_S0:merged pi+:all pi-:all', '1.84 < M < 1.89', path=kshhPath)
26 ma.reconstructDecay('B-:ch3 ->D0 K-:all', '5.24 < Mbc < 5.3 and abs(deltaE) < 0.15', path=kshhPath)
27 
28 # the variables that are printed out are: Mbc, deltaE and the daughter particle invariant masses.
29 
30 ma.variablesToHistogram(
31  filename='BtoDh_Kshh_Validation.root',
32  decayString='B-:ch3',
33  variables=[
34  ('Mbc', 100, 5.2, 5.3),
35  ('deltaE', 100, -1, 1),
36  ('daughter(0, InvM)', 100, 1.8, 1.9)], # D0 invariant mass
37  variables_2d=[
38  ('Mbc', 50, 5.23, 5.31, 'deltaE', 50, -0.7, 0.7)], path=kshhPath)
39 
40 b2.process(kshhPath)
41 print(b2.statistics)