pxdDataReductionAnalysis.py

#!/usr/bin/env python3
 
 
 
import basf2 as b2
from tracking import add_tracking_reconstruction
from simulation import add_simulation
 
b2.set_random_seed(1)
 
# Create paths
main = b2.create_path()
 
# Add modules to paths
main.add_module('EventInfoSetter', expList=[0], runList=[1], evtNumList=[100])
main.add_module('EvtGenInput', logLevel=b2.LogLevel.INFO)
 
add_simulation(main, forceSetPXDDataReduction=True, usePXDDataReduction=False)
 
add_tracking_reconstruction(main, ['SVD', 'CDC'], mcTrackFinding=True)
 
main.add_module('PXDROIFinder',
                recoTrackListName='RecoTracks',
                PXDInterceptListName='PXDIntercepts',
                ROIListName='ROIs',
                logLevel=b2.LogLevel.DEBUG)
 
main.add_module('PXDROIFinderAnalysis',
                recoTrackListName='RecoTracks',
                PXDInterceptListName='PXDIntercepts',
                ROIListName='ROIs',
                writeToRoot=True,
                rootFileName='pxdDataRedAnalysis_SVDCDC_MCTF_test',
                logLevel=b2.LogLevel.RESULT)
main.add_module('Progress')
 
# Process events
b2.process(main)
 
print(b2.statistics)