unitTest.m

function success = unitTest()

   success = 0;

   SaveDir = smi_helpers.mkSMITETmpDir('unitTest', 'SimSMLM');

   obj = smi_sim.SimSMLM();
   obj.SZ = 128; %256
   obj.Rho = 50; %30
   obj.NDatasets = 10; %2
   obj.NFrames = 1000; % 10
   obj.ZoomFactor = 2;
   obj.K_OnToOff = 1;
   obj.K_OffToOn = 0.005;
   obj.K_OnToBleach = 0.2;
   obj.EmissionRate = 1000;
   obj.Bg = 15;
   obj.PSFSigma = 1.3;
   obj.StartState = 'Equib';
   obj.LabelingEfficiency = 1;

   NWings = 16;
   obj.simStar(NWings);
   SMD_Data = obj.genNoisySMD(obj.SMD_Model);
   
   [Model, Data] = obj.genImageStack();

   % Display the blobs without noise:
   figure; imagesc(sum(Model, 3)); colormap(gca, gray(256));
   saveas(gcf, fullfile(SaveDir, 'Model.png'));
   %dipshow(Model)
   
    % Display the blobs having Poisson noise:
   figure; imagesc(sum(Data, 3)); colormap(gca, gray(256));
   saveas(gcf, fullfile(SaveDir, 'Data.png'));
   %dipshow(Data)

   success = 1;
    
end