Page 1 of 1

Spectrum convolution in pepper

Posted: Wed Jun 05, 2024 12:44 am
by ovlaguta

Hi all,
It seems there's a convolution problem (or limitation) in pepper (EasySpin version: 6.0.0 (2024-04-01)). Basically, it sets intensity to zero for points below 6x10^-4 relative intensity. The convspec function though works properly. Please see a simple example below.

Sys.S = 1/2;
Sys.g = 2;
Sys.lw = [0 1];

Exp.mwFreq = 9.8;
Exp.Range = [300 400];
Exp.SampleFrame = [0 0 0];
Exp.MolFrame = [0 0 0];
Exp.Harmonic = 0;
Exp.nPoints = 1000;

[m,spc] = pepper(Sys,Exp);
spc = spc/max(spc);
spec = spc*0;
[~, i] = max(spc);
spec(i) = 1;
w = convspec(spec,0.1,1,0,0);
w = w/max(w);
plot(m,spc, m,w)
ylim([0 1e-3])
xlabel('Magnetic field (mT)')
legend('pepper', 'convspec')


Re: Spectrum convolution in pepper

Posted: Wed Jun 05, 2024 4:42 pm
by Stefan Stoll

This is in purpose. The Lorentzian lineshape is truncated in pepper for perfomance reasons. Is this truncation visible in any practically relevant situation?


Re: Spectrum convolution in pepper

Posted: Wed Jun 05, 2024 10:27 pm
by Stefan Stoll

It should be fixed in 6.0.1, which is available now.