cdcCombinedTrackingValidation.py

#!/usr/bin/env python3
# -*- coding: utf-8 -*-
 
"""
<header>
  <contact>software-tracking@belle2.org</contact>
  <input>EvtGenSimNoBkg.root</input>
  <output>CDCCombinedTrackingValidation.root</output>
  <description>This script validates the combined CDC tracking chain with a legendre step first
  and a cellular automaton step second in Y(4S) runs.</description>
</header>
"""
 
VALIDATION_OUTPUT_FILE = 'CDCCombinedTrackingValidation.root'
N_EVENTS = 1000
ACTIVE = True
 
import basf2
basf2.set_random_seed(1337)
 
import logging
import tracking
 
from tracking.validation.run import TrackingValidationRun
from tracking.adjustments import adjust_module
 
 
class CDCCombined(TrackingValidationRun):
    """Validate the combined CDC track-finding chain include Legendre step"""
    
    n_events = N_EVENTS
    
    generator_module = 'generic'
    
    root_input_file = '../EvtGenSimNoBkg.root'
    
    contact = 'software-tracking@belle2.org - why are we not founding this?'
 
    def finder_module(self, path):
        """Add the CDC track-finding module to the basf2 path"""
        tracking.add_cdc_track_finding(path, with_ca=True)
        # adjust_module(path, "TFCDC_WireHitPreparer",
        #               UseNLoops=1)
 
    
    tracking_coverage = {
        'WhichParticles': ['CDC'],  # Include all particles seen in CDC, also secondaries
        'UsePXDHits': False,
        'UseSVDHits': False,
        'UseCDCHits': True,
        'UseOnlyAxialCDCHits': False,
        "UseReassignedHits": True,
        "UseNLoops": 1,
        "UseOnlyBeforeTOP": True,
        'MinCDCAxialHits': 8,
        'MinCDCStereoHits': 6,
        "AllowFirstCDCSuperLayerOnly": True,
        "MergeDecayInFlight": True,
        'EnergyCut': 0,
    }
    
    pulls = True
    
    output_file_name = VALIDATION_OUTPUT_FILE
 
 
def main():
    validation_run = CDCCombined()
    validation_run.configure_and_execute_from_commandline()
 
 
if __name__ == '__main__':
    logging.basicConfig(level=logging.INFO)
    if ACTIVE:
        main()