function [brweights, branch_syns, nrnweights, nrn_syns, brstrengths, brsynratio] = getsynstate(fn)
        defaults
        brweights   = zeros(ninputs, npyrs*nbranches);
        nrnweights  = zeros(ninputs, npyrs);
        branch_syns = zeros(ninputs, npyrs*nbranches);
        nrn_syns    = zeros(ninputs, npyrs);
        
        %ff = sprintf('./data/%s_%d_%d/synstate.dat', CONDITION, ncase, run-1)
        ss = load(fn);
        for i=1:size(ss,1)
            bid=ss(i,2);
            nid=ss(i,3);
            srcid=ss(i,5);
            bstrength = ss(i,6);
            w=ss(i,7);
            if (srcid >=0 && bid <= npyrs*nbranches)
                brweights(srcid+1, bid+1) = brweights(srcid+1,  bid+1) + w;
                brstrengths(srcid+1, bid+1)=bstrength;
                nrnweights(srcid+1,   nid+1) = nrnweights(srcid+1,  nid+1) + w;
            end
            if (srcid >=0 && bid <= npyrs*nbranches &&  w > 0.7)
                branch_syns(srcid+1,  bid+1) = branch_syns(srcid+1, bid+1)+1;
                nrn_syns(srcid+1,  nid+1) = nrn_syns(srcid+1, nid+1)+1;
            end

        end
        
        %for n=1:ninputs
        %    mg = find(ss(:,5)==n && ss(:,2) <= npyrs*nbranches)
        %    brweights(n, )
        %end

end