cait.simulate

cait.simulate.generate_ps_par(h)[source]

Simulate PS parameters modelled after the TUM40 phonon detector.

Parameters:

h (1D array) – Array of pulse heights.

Returns:

(t0, An, At, tau_n, tau_in, tau_t)

Return type:

list of 6 1D arrays

cait.simulate.simulate_baselines(path_h5, size, rms_thresholds, lamb=0.01, kde=True, sim_poly=True, verb=False)[source]

Creates fake baselines with given nps and drift structure.

Parameters:

  • path_h5 (string) – Path to the file from which nps and bl drifts come.

  • size (int) – Nmbr of baselines to simulate.

  • rms_thresholds (list of two ints) – Threshold for the bl rms fit error above which they get not included in the nps and drifts.

  • lamb (float) – Parameter for the noise simulation method.

  • kde (bool) – If True we sample the coefficients of the bl fit with a kernel density estimation.

  • sim_poly (pool) – If True we simulate the polynomials for the baselines.

  • verb (bool) – If true feedback about progress in code.

Returns:

(ch_nmbr, size bl, rec_len) of the simulated baselines, (ch_nmbr, size bl, rec_len) of the simulated polynomials.

Type:

tuple (3D array, 3D array)

cait.simulate.simulate_events(path_h5, type, name_appendix, size, record_length, nmbr_channels, ph_intervals=[(0, 1), (0, 1)], discrete_ph=None, exceptional_sev_naming=None, channels_exceptional_sev=[0], t0_interval=(-20, 20), fake_noise=True, use_bl_from_idx=0, take_idx=None, rms_thresholds=[1, 1], lamb=0.01, sample_length=0.04, saturation=False, reuse_bl=False, ps_dev=False)[source]

Simulates pulses on a noise baseline.

Options are to take measured noise baselines or fake bl, sev of events or testpulses, pulse heights, onset intervals.

Parameters:

  • path_h5 (string) – Path of the dataset on which the sim is based on.

  • type (string) – Either events, testpulses or noise - what type of event is simulated; if testpulses the other channels than the first will beignored for the PHs.

  • size (integer > 0) – Size of the simulated dataset, if fake_noise is False it must be smaller than the number of noise baselines in the hdf5 file.

  • record_length (integer) – Length of the record window in the hdf5 datset.

  • nmbr_channels (integer) – The number of channels in the hdf5 dataset.

  • ph_intervals (list of c intervals) – Intervals for the simulated pulse heights, with c the nmbr_channels; the heights are samples uniformly from these intervals.

  • discrete_ph (list of intervals) – c Lists of values for the pulse heights for the c channels or None for ph intervals; the heights are sampled uniformly for the list.

  • exceptional_sev_naming (string) – If set, this is full group name in the HDF5 set for the sev used for the simulation of events - by setting this, e.g. carrier events can be simulated.

  • channel_exceptional_sev (list of ints) – The channels for that the exceptional sev is used, e.g. if only for phonon channel, choose [0], if for botch phonon and light, choose [0,1].

  • t0_interval (tuple) – Interval (l,u) with l lower bound for onset and u upper bound. The onsets are sampled uniformly from the interval. In ms.

  • fake_noise (bool) – If true use simulated noise baselines, otherwise measured ones.

  • use_bl_from_idx (int) – The start index of the baselines that are used.

  • take_idx (list) – The event indices which we want to use for the simulation.

  • rms_threshold (float) – Above which value noise baselines are excluded for the distribution of polynomial coefficients; also, a cut value for the baselines if not the fake ones but the ones from the h5 set are taken.

  • lamb (float) – The parameter for the bl simulation method.

  • sample_length (float) – The length in ms of one sample from an event.

  • saturation (bool) – If True the logistic curve is applied to the pulses.

  • reuse_bl (bool) – If True the same baselines are used multiple times to have enough of them (use this with care to not have identical copies of events).

  • ps_dev (bool) – If True the pulse shape parameters are modelled with deviations.

Returns:

The simulated events, the true pulse heights, the onsets of the events.

Return type:

(3D array of size (nmbr channels, size, record_length), the simulated events, 2D array of size (nmbr channels, size), the true pulse heights, 1D array (size), the onsets of the events)