function algorParam = sepia_config()

% General algorithm parameters
algorParam = struct();
algorParam.general.isBET = 0;
algorParam.general.isInvert = 0;
algorParam.general.isRefineBrainMask = 1;

% Total field recovery algorithm parameters
algorParam.unwrap.echoCombMethod = 'ROMEO total field calculation';
algorParam.unwrap.offsetCorrect = 'On';
algorParam.unwrap.mask = 'SEPIA mask';
algorParam.unwrap.qualitymaskThreshold = 0.5;
algorParam.unwrap.useRomeoMask = 0;
algorParam.unwrap.isEddyCorrect = 1;
algorParam.unwrap.isSaveUnwrappedEcho = 0;
algorParam.unwrap.excludeMaskThreshold = 0.3;
algorParam.unwrap.excludeMethod = 'Weighting map';

% Background field removal algorithm parameters
algorParam.bfr.refine_method = 'None';
algorParam.bfr.refine_order = 4;
algorParam.bfr.erode_radius = 0;
algorParam.bfr.erode_before_radius = 0;
algorParam.bfr.method = 'VSHARP';
algorParam.bfr.radius = 12:-1:1;

% QSM algorithm parameters
algorParam.qsm.reference_tissue = 'Brain mask';
algorParam.qsm.method = 'LPCNN';