Belle II Software  release-05-01-25
phase2_sim_reco_tracking.py
1 #!/usr/bin/env python3
2 # -*- coding: utf-8 -*-
3 
4 
8 
9 import os
10 from basf2 import *
11 from tracking import *
12 from svd import *
13 from rawdata import *
14 from ROOT import Belle2
15 from reconstruction import add_reconstruction
16 import os.path
17 import sys
18 
19 print('***')
20 print('*** Used steering script:')
21 with open(sys.argv[0], 'r') as fin:
22  print(fin.read(), end="")
23 print('*** end of the script.')
24 print('***')
25 
26 fileIN = sys.argv[1]
27 fileOUT = sys.argv[2]
28 
29 main = create_path()
30 main.add_module('RootInput', inputFileName=fileIN)
31 
32 gearbox = register_module('Gearbox')
33 geomfile = '/geometry/Beast2_phase2.xml'
34 if geomfile != 'None':
35  gearbox.param('fileName', geomfile)
36 
37 main.add_module(gearbox)
38 geometry = register_module('Geometry')
39 geometry.param('useDB', False)
40 geometry.param('components', ['SVD'])
41 main.add_module(geometry)
42 
43 # SVD reconstruction
44 main.add_module('SVDCoGTimeEstimator')
45 main.add_module('SVDSimpleClusterizer', Clusters="SVDClusters")
46 
47 # add tracking recontruction
48 add_tracking_reconstruction(main, components=['SVD'])
49 # remove cut on CoG in SpacePointCreator
50 for module in main.modules():
51  if module.name() == 'SVDSpacePointCreator':
52  module.param("MinClusterTime", -999)
53 
54 main.add_module('VXDDedxPID', enableDebugOutput=True)
55 
56 main.add_module(
57  'RootOutput',
58  outputFileName=fileOUT,
59  branchNames=[
60  "SVDShaperDigits",
61  "SVDClusters",
62  "SVDRecoDigits",
63  "Tracks",
64  "RecoTracks",
65  "TracksToRecoTracks",
66  "TrackFitResults",
67  "MCParticles",
68  "SVDTrueHits",
69  "TracksToVXDDedxLikelihoods",
70  "TracksToVXDDedxTracks",
71  "VXDDedxLikelihoods",
72  "VXDDedxTracks"])
73 print_path(main)
74 main.add_module('ProgressBar')
75 # Process events
76 process(main)
77 print(statistics)