Belle II Software  release-05-01-25
SVDValidationRecoDigitPerformance.py
1 # !/usr/bin/env python3
2 # -*- coding: utf-8 -*-
3 
4 """
5 <header>
6  <input>SVDValidationTTreeRecoDigit.root</input>
7  <input>SVDValidationTTreeStrip.root</input>
8  <output>SVDRecoDigitPerformance.root</output>
9  <description>
10  Validation plots related to RecoDigit performance.
11  </description>
12  <contact>
13  SVD Software Group, svd-software@belle2.org
14  </contact>
15 </header>
16 """
17 
18 from plotUtils import *
19 
20 inputRD = R.TFile.Open("../SVDValidationTTreeRecoDigit.root")
21 inputS = R.TFile.Open("../SVDValidationTTreeStrip.root")
22 
23 treeRD = inputRD.Get("tree")
24 treeS = inputS.Get("tree")
25 
26 histsRDP = R.TFile.Open("SVDRecoDigitPerformance.root", "recreate")
27 
28 ploter(
29  name='TimeResolution',
30  title='Resolution of the reconstructed time of a single strip',
31  nbins=200,
32  xmin=-100,
33  xmax=100,
34  x_label='RecoDigit time resolution (ns)',
35  y_label='counts',
36  granules=gD,
37  tree=treeRD,
38  expr='recodigit_time - truehit_time',
39  cut='',
40  descr='Definition: (reconstructed time of the RecoDigit) - (time of the Truehit);\
41  for all RecoDigit related to one TrueHit.',
42  check='Distribution peaks around 0.',
43  isShifter=True)
44 
45 ploter(
46  name='StripCharge',
47  title='Charge of the strip',
48  nbins=50,
49  xmin=0,
50  xmax=120000,
51  x_label='Strip charge (# of electrons)',
52  y_label='counts',
53  granules=gD,
54  tree=treeS,
55  expr='strip_charge',
56  cut='',
57  descr='Strip charge as the highest charge among its digits. For all RecoDigit related to one TrueHit.',
58  check='Distribution peak around 15-30 ke.',
59  isShifter=False)