Response function with RadiationDetectorDSP

[rajxabc]

Hi Felix ,
i am trying to use RadiationDetectorDSP wit SSD , I am getting this error when applying Circuit filter

wfx = wf[1].waveform
flt = IntegratorFilter(gain = 2.0)
output = flt(wfx)
ERROR: DimensionError: e and 2.0 are not dimensionally compatible.
[1] convert(::Type{Quantity{Float64, 𝐈 𝐓, Unitful.FreeUnits{(e,), 𝐈 𝐓, nothing}}} , x::Float64)

Do I need to convert something in this ??
Thanks

[oschulz]

Hi @rajxabc,

RadiationDetectorDSP currently doesn't really support unitful sample values. It's not just a technical issue, it's also a semantic one: With an integrator filter, for example, should the units of the samples changes (multiplied by time)? It would depends on the meaning of filter - whether it's just an "circuit" (so the output of an integrator would have the same unit) or a mathematical operation (integration of current should result in charge, etc.). And biquad filter can do both integration and differentiation, depending on the parameter values. So it's not clear how RadiationDetectorDSP can do the right thing with sample value units (units on the time axis are supported, though).

Just strip the units from the signal before you apply the filter.

[rajxabc]

Gotit Thanks