cosmicsAxialHoughTrackingValidation.py

#!/usr/bin/env python3
 
 
 
"""
<header>
  <contact>software-tracking@belle2.org</contact>
  <input>CosmicsSimNoBkg.root</input>
  <description>Validates the axial hough finder working on segments in cosmics events.</description>
</header>
"""
 
from tracking.validation.run import TrackingValidationRun
import logging
import trackfindingcdc.cdcdisplay as cdcdisplay
import os
import basf2
VALIDATION_OUTPUT_FILE = 'CosmicsAxialHoughTrackingValidation.root'
N_EVENTS = 1000
ACTIVE = False
 
 
class CosmicsAxialHough(TrackingValidationRun):
    """Validate the CDC axial Hough track finder with cosmic rays"""
    
    n_events = N_EVENTS
    
    generator_module = 'Cosmics'
    
    root_input_file = '../CosmicsSimNoBkg.root'
    
    components = None
 
    def finder_module(self, path):
        """Add the CDC track-finding module to the basf2 path"""
        path.add_module('TFCDC_WireHitPreparer')
        path.add_module("TFCDC_ClusterPreparer")
        path.add_module('TFCDC_SegmentFinderFacetAutomaton',
                        SegmentOrientation="downwards")
        path.add_module('TFCDC_AxialTrackCreatorSegmentHough')
        path.add_module('TFCDC_TrackExporter')
 
        interactive_display = False
        if interactive_display:
            cdc_display_module = cdcdisplay.CDCSVGDisplayModule(os.getcwd(), interactive=True)
            cdc_display_module.draw_recotracks = True
            cdc_display_module.draw_recotrack_seed_trajectories = True
            path.add_module(cdc_display_module)
 
    
    tracking_coverage = {
        'WhichParticles': ['CDC'],  # Include all particles seen in CDC, also secondaries
        'UsePXDHits': False,
        'UseSVDHits': False,
        'UseCDCHits': True,
        'UseOnlyAxialCDCHits': True,
        "UseReassignedHits": True,
    }
 
    
    pulls = True
    
    output_file_name = VALIDATION_OUTPUT_FILE
 
 
def main():
    basf2.set_random_seed(1337)
    validation_run = CosmicsAxialHough()
    validation_run.configure_and_execute_from_commandline()
 
 
if __name__ == '__main__':
    logging.basicConfig(level=logging.INFO)
    if ACTIVE:
        main()
    else:
        print("This validation deactivated and thus basf2 is not executed.\n"
              "If you want to run this validation, please set the 'ACTIVE' flag above to 'True'.\n"
              "Exiting.")