Pulse Mode Analysis

This repository contains a Python implementation for analyzing pulse mode measurements. The main class, PulseModeAnalysis, provides methods for processing waveform data and extracting various parameters, such as baseline, pulse shape, and charge.

To install the library use pip:

pip install git+https://github.com/martinunland/pulse_mode_analysis.git

Or download the wheel in repositories and pip:

pip install pulse_mode_analysis-0.1.1-py3-none-any.whl

Usage

First, import the PulseModeAnalysis class:

from pulse_mode_analysis import PulseModeAnalysis

Next, create an instance of PulseModeAnalysis with the desired sampling interval, baseline minimum time, and baseline maximum time:

analysis = PulseModeAnalysis(sampling_interval=1e-9, baseline_tmin=1e-9, baseline_tmax=10e-9)
(...)
analysis.update_time_axis(waveform)

Now, you can use the various methods provided by the PulseModeAnalysis class to process and analyze the waveform data, for example:

data =... #data is a block of waveforms you got from the DAQ
baseline, baseline_error = analysis.get_baseline(data)
for waveform in data:
  waveform -= baseline
  pedestal_charge, transit_time, charge, amplitude, FWHM, RT, FT = analysis.process_waveform(waveform)

The get_simple_intensity method calculates a first approximation of the PMT intensity and baseline of the waveform data. It takes the waveform data as input and an optional boolean mask for the signal region. If no mask is provided, the entire waveform is integrated.

(mean, error), (baseline, baseline_error) = analysis.get_simple_intensity(data)