genDoublePulse.py

#!/usr/bin/env python3
# -*- coding: utf-8 -*-
 
from basf2 import *
import os
import math
 
# ----------------------------------------------------------------------------
# Example of generating double calibration pulses with TOPDoublePulseGenerator
# ----------------------------------------------------------------------------
 
# Create path
main = create_path()
 
# Set number of events to generate
eventinfosetter = register_module('EventInfoSetter')
eventinfosetter.param('evtNumList', [100])
main.add_module(eventinfosetter)
 
# Gearbox: access to database (xml files)
gearbox = register_module('Gearbox')
main.add_module(gearbox)
 
# Geometry
geometry = register_module('Geometry')
geometry.param('useDB', False)
geometry.param('components', ['TOP'])
main.add_module(geometry)
 
# pulse generator
calpulse = register_module('TOPDoublePulseGenerator')
calpulse.param('asicChannels', [0])
calpulse.param('moduleIDs', [1])
# intervals = [math.exp(-float(i)/99.0)+1 for i in range(128)]
# intervals = intervals + intervals
# calpulse.param('sampleTimeIntervals', intervals)
calpulse.param('useDatabase', True)
calpulse.param('outputFileName', 'usedSampleTimes.root')
main.add_module(calpulse)
 
# output
main.add_module('RootOutput')
 
# Show progress of processing
progress = register_module('Progress')
main.add_module(progress)
 
# Process events
process(main)
 
# Print call statistics
print(statistics)