function optionStruct = baseDefinition()
%% test identification of continuously rounded particles
% Base: octahedron
% Roundness: RelAddSph varies from 0 to 1
% Orientation: Set of random orientations

%% orientations to track
rng(1)
nSample = 20;
x = randn(3, nSample);
dirMatrix =  bsxfun(@times, 1./sqrt(sum(x.*x, 1)), x);
optionStruct.dirMatrix = dirMatrix;
% plot(dirMatrix(1,:), dirMatrix(2,:), 'b.')
% plot3(dirMatrix(1,:), dirMatrix(2,:), dirMatrix(3,:), 'b.')

%% relative roughness to track
% rrVector = 0:0.2:1;
% rrVector = [0.1 0.5 0.9];
rrVector = [0 0.05 0.1 0.2 0.4 0.6 0.8 0.9 0.95 1];
optionStruct.rrVector = rrVector;

%% array of random background
optionStruct.bgInfo{1} = rand(length(rrVector), size(dirMatrix, 2))*0;
optionStruct.bgInfo{2} = rand(length(rrVector), size(dirMatrix, 2))*0;

%% base options
  cdef = obtainCrystal('generic_potash_alum', 3);
  hrep.A = cdef.A;
%   hrep.h = cdef.groupMappingMatrix * [1 1.45 1.15]';
  hrep.h = cdef.groupMappingMatrix * [1 1.5 1.2]';
%   hrep.h = cdef.groupMappingMatrix * [1 3 3]';
optionStruct.hrepBase = hrep;
optionStruct.writeBaseImage = false;
optionStruct.writeBlurrImage = false;
optionStruct.writeEdgeImage = false;

%% variable optoins
optionStruct.lambdaLocal = 1;
optionStruct.lambdaError = 'sqRough';
optionStruct.particleSize = 200;
optionStruct.blurrSize = 50;
optionStruct.blurrStd = 12;
optionStruct.noise = 0.0002;
optionStruct.rrDiff = 0.1*0;
optionStruct.aggregate = false;

optionStruct.caseName = 'base';