diana_tutorial_averagepulse.cfg

Go to the documentation of this file.

# /** @example diana_tutorial_averagepulse.cfg 
#     @brief Tutorial 2. 
#     @author Marco Vignati */
# 
# Sequence to compute the average signal pulse of each channel.
# It includes the filter  MTutorialAmplitudeFilter and  the module MTutorialAveragePulse.
# Create the directory "avg" if it does not exist then execute as:
#
# shell> diana -C cfg/tutorial/diana_tutorial_averagepulse.cfg
#
# Average pulses are saved as global quantities in the file
# produced by MTutorialAveragePulse. To look at this output,
# open a TBrowser whitin root, open the file, go into
# the "Global" directory, right click on an average pulse
# and then "Draw".
 
 
sequence TutorialAveragePulseSequence
 
# Reads a root file containing events. The file being read
# is produced by diana_tutorial.cfg
reader RootFileReader
verbosity = info
enable = true
InputFile = data/Analyzed.list
EventNumberPrintStep = 10000
endmod
 
# Load QRunData and QDetChannelCollection from the file read by the reader,
# and save it to the output file written by the writer.
module RunDataLoader
verbosity = info
enable = true
Input = CurrentReader
Output = CurrentWriter
endmod
 
# Filter events using raw event information. In this case 
# events marked as "signal" by the DAQ are selected
# The Logic=CASE options tells to ignore the result of previous
# filters. Since there are no filters before this filter,
# the Logic option is useless in this case.
# Since RawDataFilter information is already present in the
# event from the diana_tutorial.cfg sequence, we specify
# an ExtraLabel. The label of the bool variable in the QEvent will be
# RawDataFilter_AveragePulse@Passed instead of RawDataFilter@Passed.
filter RawDataFilter
enable = true
verbosity = info
Logic = CASE
KeepSignal = true
KeepHeater = false
KeepNoise = false
KeepThermometers = false
ExtraLabel = AveragePulse
endmod
 
# select waveforms with only 1 identified peak and
# with small baseline slope (removes pileups).
filter BadPulse
enable = true
verbosity = info
Logic = AND
NumberOfPeaks = 1
MaxBaselineSlope = 0.02
endmod
 
# Select events with amplitude in the range
# 50-1000. The amplitude variable can be chosen by
# the user. In this case we use PulseBasicParameters@MaxBaseline,
# the Maximum-Baseline value in mV computed by module PulseBasicParameters.
# The Logic is set to AND, i.e.  the output of this filter is
# combined with the output of previouse filters.
# As result only events which are marked as Signal (selcted by
# RawDataFilter above) with amplitude in the 50-1000 mV range
# will be available for MTutorialAveragePulse
filter TutorialAmplitudeFilter
enable = true
verbosity = info
Logic = AND
AmplitudeMinimum = 50
AmplitudeMaximum = 3000
AmplitudeLabel = PulseBasicParameters@MaxBaseline
endmod
 
# Average pulses. The average is performed on each
# channel separately.
module TutorialAveragePulse
enable = true
verbosity = info
Output = avg/averagepulses.root
endmod
 
# Write events to file with file prefix "AveragePulses"
writer RootFileWriter
verbosity = info
enable = true
Description =  data with amplitude
WriteFilesList = true
OutputFilesList = AveragePulses.list
OutputFilePrefix = AveragePulses
OutputDir = data
AliasFileName = ${DIANA_INSTALL}/cfg/aliases.txt
endmod
 
endseq