unpackToTOPDigitsWithTBC.py

#!/usr/bin/env python
 
 
 
# ---------------------------------------------------------------------------------------
# Unpack raw data to TOPDigits using time base calibration
#
# Usage: basf2 unpackToTOPDigitsWithTBC.py -i <input_file.root> -o <output_file.root>
# ---------------------------------------------------------------------------------------
 
import basf2 as b2
 
# Database
b2.conditions.override_globaltags()
b2.conditions.append_globaltag('online')
# b2.conditions.append_testing_payloads('localDB-FEMaps/localDB.txt')  # SCROD mapping from local database
# b2.conditions.append_testing_payloads('localDB-timebase/localDB.txt')  # timebase calibration from local database
 
# Create path
main = b2.create_path()
 
# input
# roinput = b2.register_module('SeqRootInput')  # sroot files
roinput = b2.register_module('RootInput')  # root files
main.add_module(roinput)
 
# conversion from RawCOPPER or RawDataBlock to RawDetector objects (uncomment for pocketDAQ)
# main.add_module('Convert2RawDet')
 
# Initialize TOP geometry parameters (creation of Geant geometry is not needed)
main.add_module('TOPGeometryParInitializer')
 
# Unpacking (format auto detection works now)
unpack = b2.register_module('TOPUnpacker')
main.add_module(unpack)
 
# Convert to TOPDigits
converter = b2.register_module('TOPRawDigitConverter')
converter.param('useSampleTimeCalibration', True)  # enable using calibration
converter.param('useChannelT0Calibration', False)
converter.param('useModuleT0Calibration', False)
converter.param('useCommonT0Calibration', False)
converter.param('calibrationChannel', 0)  # if set, cal pulses will be flagged
converter.param('calpulseHeightMin', 450)  # in [ADC counts]
converter.param('calpulseHeightMax', 900)  # in [ADC counts]
converter.param('calpulseWidthMin', 2.0)  # in [ns]
converter.param('calpulseWidthMax', 6.0)  # in [ns]
main.add_module(converter)
 
# output
output = b2.register_module('RootOutput')
output.param('branchNames', ['TOPDigits', 'TOPRawDigits', 'TOPProductionEventDebugs',
                             'TOPProductionHitDebugs', 'TOPRawDigitsToTOPProductionHitDebugs',
                             # 'TOPRawWaveforms', 'TOPRawDigitsToTOPRawWaveforms',
                             ])
main.add_module(output)
 
# Print progress
progress = b2.register_module('Progress')
main.add_module(progress)
 
# Process events
b2.process(main)
 
# Print statistics
print(b2.statistics)