subjVec = {'ss040804','mho040625','bg040719'};
threshVec = [5000 2000 1000 500];
for ss = 1:length(subjVec)
    subjDir = ['/teal/scr1/dti/sisr/' subjVec{ss}];
    fgDir = [subjDir '/conTrack/resamp_LDOCC'];
    cd fgDir;
    mtrComputeManyFiberDensities(subjDir,'paths*100k*',);
end