Belle II Software  release-05-01-25
toCDCfromEclCKFTrackingValidation.py
1 #!/usr/bin/env python3
2 # -*- coding: utf-8 -*-
3 
4 """
5 <header>
6  <contact>simon.kurz@desy.de</contact>
7  <input>EvtGenSimNoBkg.root</input>
8  <output>toCDCfromEclCKFTrackingValidation.root</output>
9  <description>This module validates toCDCfromEclCKF module using Y(4S) runs.</description>
10 </header>
11 """
12 
13 VALIDATION_OUTPUT_FILE = 'toCDCfromEclCKFTrackingValidation_expert.root'
14 N_EVENTS = 10000
15 ACTIVE = True
16 
17 import basf2
18 basf2.set_random_seed(1337)
19 
20 import logging
21 import tracking
22 from tracking.validation.run import TrackingValidationRun
23 import reconstruction
24 
25 
27  """Validate the ecl seeded ckf"""
28 
29  n_events = N_EVENTS
30 
31  generator_module = 'generic'
32 
33  root_input_file = '../EvtGenSimNoBkg.root'
34 
35  @staticmethod
36  def finder_module(path):
37  """Only run the following (necessary) modules"""
38  path.add_module('SetupGenfitExtrapolation',
39  energyLossBrems=False, noiseBrems=False)
40 
42 
43  # needed for truth matching
44  path.add_module('TrackFinderMCTruthRecoTracks',
45  RecoTracksStoreArrayName="MCRecoTracks",
46  WhichParticles=[],
47  UsePXDHits=False,
48  UseSVDHits=False,
49  UseCDCHits=True)
50 
51  path.add_module("TFCDC_WireHitPreparer",
52  wirePosition="aligned",
53  useSecondHits=False,
54  flightTimeEstimation="outwards")
55 
56  path.add_module("ToCDCFromEclCKF",
57  inputWireHits="CDCWireHitVector",
58  minimalEnRequirementCluster=0.3,
59  eclSeedRecoTrackStoreArrayName='EclSeedRecoTracks',
60  hitFindingDirection="backward",
61  outputRecoTrackStoreArrayName="CKFCDCRecoTracksFromEcl",
62  outputRelationRecoTrackStoreArrayName="EclSeedRecoTracks",
63  writeOutDirection="forward",
64  stateBasicFilterParameters={"maximalHitDistance": 7.5, "maximalHitDistanceEclSeed": 75.0},
65  stateExtrapolationFilterParameters={"direction": "backward"},
66  pathFilter="arc_length_fromEcl",
67  inputECLshowersStoreArrayName="ECLShowers",
68  trackFindingDirection="backward",
69  seedComponent="ECL"
70  )
71 
72  path.add_module("ToCDCCKF",
73  inputWireHits="CDCWireHitVector",
74  inputRecoTrackStoreArrayName="CKFCDCRecoTracksFromEcl",
75  relatedRecoTrackStoreArrayName="RecoTracks",
76  relationCheckForDirection="backward",
77  outputRecoTrackStoreArrayName="RecoTracks",
78  outputRelationRecoTrackStoreArrayName="CKFCDCRecoTracksFromEcl",
79  writeOutDirection="backward",
80  stateBasicFilterParameters={"maximalHitDistance": 0.75},
81  stateExtrapolationFilterParameters={"direction": "forward"},
82  pathFilter="arc_length",
83  seedComponent="ECL"
84  )
85 
86  path.add_module("DAFRecoFitter", recoTracksStoreArrayName="RecoTracks")
87 
88  path.add_module('TrackCreator', recoTrackColName='RecoTracks')
89 
90  path.add_module("MCRecoTracksMatcher",
91  mcRecoTracksStoreArrayName="MCRecoTracks",
92  prRecoTracksStoreArrayName="RecoTracks",
93  UseCDCHits=True,
94  UseSVDHits=False,
95  UsePXDHits=False)
96 
97 
98  tracking_coverage = {
99  'UsePXDHits': False,
100  'UseSVDHits': False,
101  'UseCDCHits': True,
102  'WhichParticles': [],
103  }
104 
105 
106  fit_tracks = False
107 
108  use_fit_information = True
109 
110  pulls = True
111 
112  resolution = True
113 
114  output_file_name = VALIDATION_OUTPUT_FILE
115 
116  extended = True
117 
118 
119 def main():
120  validation_run = toCDCfromEclCKF()
121  validation_run.configure_and_execute_from_commandline()
122 
123 
124 if __name__ == '__main__':
125  logging.basicConfig(level=logging.INFO)
126 
127  if ACTIVE:
128  main()
tracking.validation.run
Definition: run.py:1
toCDCfromEclCKFTrackingValidation.toCDCfromEclCKF
Definition: toCDCfromEclCKFTrackingValidation.py:26
main
int main(int argc, char **argv)
Run all tests.
Definition: test_main.cc:77
tracking.run.tracked_event_generation.ReadOrGenerateTrackedEventsRun.finder_module
finder_module
Name of the finder module to be used - can be everything that is accepted by tracking....
Definition: tracked_event_generation.py:25
reconstruction.add_ecl_modules
def add_ecl_modules(path, components=None)
Definition: reconstruction.py:618
tracking.validation.run.TrackingValidationRun
Definition: run.py:18