5505-public/nappMATLAB-toolbox
0
0 Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

Updated comment to nappBackground to make more clear that pos 2 is always cutoff in linear fit

Browse Source
This commit is contained in:
bartels-rausch committed 2026-06-01 13:13:43 +02:00
1 parent caec12ce3d
commit 9d6d271461
4 files changed
+259 -122
Download Patch File Download Diff File Expand all files Collapse all files

No files matched your search

beamtime/documentedExampleScripts/TestAndDevelopShirley.m
+45
View File
@@ -78,6 +78,25 @@ FitList(tmpDataSelection,:).spectrum_counts{:} = y;
nappPlot(FitList,tmpDataSelection,'tableVariable','spectrum_countsNorm','energyScale','kinetic');
nappPlot(FitList,tmpDataSelection,'tableVariable','spectrum_counts','energyScale','binding');
%% TEST 0.1: Linear Background
% exclude one central area
nappBackground(FitList,109,116,'tableSelection',tmpDataSelection,'mode', 'linfit','tableVariable','spectrum_countsNorm','energyScale','kinetic','showFigures','on');
% exclude one central area and start
nappBackground(FitList,[109 106],[116 119],'tableSelection',tmpDataSelection,'mode', 'linfit','tableVariable','spectrum_countsNorm','energyScale','kinetic','showFigures','on');
% exclude 2 central area and start
nappBackground(FitList,[109 106 107],[116 119 108],'tableSelection',tmpDataSelection,'mode', 'linfit','tableVariable','spectrum_countsNorm','energyScale','kinetic','showFigures','on');
% exclude 2 central area and start
nappBackground(FitList,[109 106 107],[116 NaN 108],'tableSelection',tmpDataSelection,'mode', 'linfit','tableVariable','spectrum_countsNorm','energyScale','kinetic','showFigures','on');
% exclude one central area and start
nappBackground(FitList,[109 106],116,'tableSelection',tmpDataSelection,'mode', 'linfit','tableVariable','spectrum_countsNorm','energyScale','kinetic','showFigures','on');
% exclude one central area
nappBackground(FitList,35,39,'tableSelection',tmpDataSelection,'mode', 'linfit','tableVariable','spectrum_counts','energyScale','binding','showFigures','on');
% exclude one central area and start
nappBackground(FitList,[35 31],[39 40],'tableSelection',tmpDataSelection,'mode', 'linfit','tableVariable','spectrum_counts','energyScale','binding','showFigures','on');
%% TEST 1.1: Unconstrained kinetic Gauss linear
% Gaussian + Linear, individual area/width/position per peak
@@ -97,6 +116,32 @@ tmpFitParameters.peakPosLowerBound = [112 109];
tmpFitParameters.peakPosUpperBound = [114 111];
tmpFitParameters.backgroundRange = [12 6]; % <- selects fit range around peak1
FitList = nappFitGaussXps(FitList,tmpFitParameters,tmpDataSelection,...
'tableVariable','spectrum_countsNorm','energyScale','kinetic',...
'mode','unconstrained','backgroundMode','linear','showFigures','on');
%% TEST 1.1B: Unconstrained kinetic Gauss linear
% Gaussian + Linear, individual area/width/position per peak
clearvars tmpFitParameters;
clc
tmpFitParameters.peakArea = [5e5 100]; % <- 2 values => 2 peaks
tmpFitParameters.peakWidth = [1.6 0.6];
tmpFitParameters.peakPos = [113 110];
% Bounds for areas (lower and upper)
tmpFitParameters.peakAreaLowerBound = [1e3 1e3];
tmpFitParameters.peakAreaUpperBound = [Inf Inf];
% Bounds for widths
tmpFitParameters.peakWidthLowerBound = [1.0 0.6];
tmpFitParameters.peakWidthUpperBound = [2.0 1.2];
% Bounds for positions
tmpFitParameters.peakPosLowerBound = [112 109];
tmpFitParameters.peakPosUpperBound = [114 111];
tmpFitParameters.backgroundRange = [12 6]; % <- selects fit range around peak1
% Bounds for background
tmpFitParameters.bkgrd = [0 0];
tmpFitParameters.pbkgrdLowerBound = [0 0];
tmpFitParameters.bkgrdUpperBound = [0 0];
FitList = nappFitGaussXps(FitList,tmpFitParameters,tmpDataSelection,...
'tableVariable','spectrum_countsNorm','energyScale','kinetic',...
'mode','unconstrained','backgroundMode','linear','showFigures','on');