/*
This file is part of libszdist.

libszdist is free software: you can redistribute it and/or modify
it under the terms of the GNU General Public License as published by
the Free Software Foundation, either version 3 of the License, or
(at your option) any later version.

libszdist is distributed in the hope that it will be useful,
but WITHOUT ANY WARRANTY; without even the implied warranty of
MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
GNU General Public License for more details.

You should have received a copy of the GNU General Public License
along with libszdist.  If not, see <http://www.gnu.org/licenses/>.
*/


// Compares results of inversion algorithms : libszdist vs. IfT with data
// from ACTRIS SMPS workshop, Oct. 2013, Leipzig.
//

fileIft="../files/ift2013/smps131002.in2";
fileTst="../files/out/smps20131002.csv";
info="ACTRIS SMPS workshop, 2013-10-02, eq len = 9.8 m, flow ratio = 1:5";

exec('plotResults.sci');


scf();
//-----------------------------------------------------------------------------
subplot(2,1,1);

[scRef,concRef]=avgScan(fileIft,"ift2013");   // no total conc in these files
plotScan(scRef,color("green"),"Ift");

errL=[scRef(:,1), scRef(:,2)*0.90];
plotScan(errL,color("grey"),"IfT -10%");
errH=[scRef(:,1), scRef(:,2)*1.10];
plotScan(errH,color("grey"),"IfT +10%");

[scTst,concTst]=avgScan(fileTst,"lamp");
lbl=msprintf("LaMP - %.0f #/cm3",concTst);
plotScan(scTst,color("red"),lbl);

xgrid(2);
a=gca();
a.grid=[color("grey"),color("grey"),-1];
a.data_bounds(1,1)=10;
a.data_bounds(2,1)=1000;
xtitle(info,"Dp (nm)","dN/dLog");

//-----------------------------------------------------------------------------
subplot(2,1,2);
plotRatio(scTst,scRef,color("red"),"LaMP / IfT");
xgrid(2);
a=gca();
a.grid=[color("grey"),color("grey"),-1];
a.data_bounds(1,1)=10;        // xmin
a.data_bounds(2,1)=1000;      // xmax
a.data_bounds(1,2)=0.8;       // ymin
a.data_bounds(2,2)=1.2;       // ymax
xtitle("","Dp (nm)","ratio");