toCDCCKFTrackingValidationBkg.py

#!/usr/bin/env python3
 

 
"""
<header>
  <contact>software-tracking@belle2.org</contact>
  <input>EvtGenSim.root</input>
  <output>toCDCCKFTrackingValidationBkg.root</output>
  <description>This module validates ToCDCCKF module using Y(4S) runs with background.</description>
</header>
"""
 
from tracking.path_utils import add_vxd_track_finding_vxdtf2
from tracking.validation.run import TrackingValidationRun
import logging
import basf2
import svd
 
VALIDATION_OUTPUT_FILE = 'toCDCCKFTrackingValidationBkg.root'
N_EVENTS = 1000
ACTIVE = True
 
 
class toCDCCKFValidationBkg(TrackingValidationRun):
    """Validate the ToCDCCKF Kalman finder/filter algorithm with Y(4S) events"""
    
    n_events = N_EVENTS
    
    generator_module = 'generic'
    
    root_input_file = '../EvtGenSim.root'
 
    @staticmethod
    def finder_module(path):
        """Add the ToCDCCKF module and related modules to the basf2 path"""
 
        path.add_module('SetupGenfitExtrapolation',
                        energyLossBrems=False, noiseBrems=False)
 
        svd.add_svd_reconstruction(path)
        add_vxd_track_finding_vxdtf2(path, reco_tracks="RecoTracksSVD", components=["SVD"])
        path.add_module("DAFRecoFitter", recoTracksStoreArrayName="RecoTracksSVD")
 
        path.add_module("TFCDC_WireHitPreparer",
                        wirePosition="aligned",
                        useSecondHits=False,
                        flightTimeEstimation="outwards")
        path.add_module("ToCDCCKF",
                        inputWireHits="CDCWireHitVector",
                        inputRecoTrackStoreArrayName="RecoTracksSVD",
                        relatedRecoTrackStoreArrayName="CKFCDCRecoTracks",
                        relationCheckForDirection="backward",
                        outputRecoTrackStoreArrayName="CKFCDCRecoTracks",
                        outputRelationRecoTrackStoreArrayName="RecoTracksSVD",
                        writeOutDirection="backward",
                        stateBasicFilterParameters={"maximalHitDistance": 0.75},
                        stateExtrapolationFilterParameters={"direction": "forward"},
                        pathFilter="arc_length")
 
        path.add_module("RelatedTracksCombiner",
                        CDCRecoTracksStoreArrayName="CKFCDCRecoTracks",
                        VXDRecoTracksStoreArrayName="RecoTracksSVD",
                        recoTracksStoreArrayName="RecoTracks")
 
        path.add_module("DAFRecoFitter", recoTracksStoreArrayName="RecoTracks")
 
        path.add_module('TrackCreator', recoTrackColName='RecoTracks')
 
    
    tracking_coverage = {
        'UsePXDHits': False,
        'UseSVDHits': True,
        'UseCDCHits': True,
        'WhichParticles': [],
    }
 
    
    fit_tracks = False
    
    use_fit_information = True
    
    pulls = True
    
    resolution = True
    
    output_file_name = VALIDATION_OUTPUT_FILE
 
 
def main():
    basf2.set_random_seed(1337)
    validation_run = toCDCCKFValidationBkg()
    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.")