KP_2PVR_obtain_field_distToRef_multiAcq_organizer_singleTrialRe.m

function sOrg = KP_2PVR_obtain_field_distToRef_multiAcq_organizer_singleTrialRe(neuHistInputStruct, in_struct_sOrg, bordersSectors, VREnd, outPath)
transfGroups = fieldnames(neuHistInputStruct.Transformations);
for g = 1:length(transfGroups)
    group = transfGroups{g};
    fprintf('Processing %s\n', group)
    neurons = fieldnames(in_struct_sOrg.sOrg.(group));
    numNeurons = length(neurons);
    refNeuron = neurons{1};
    maxNumFields = size(in_struct_sOrg.sOrg.(group).(refNeuron),1);
    maxNumTrials = size(in_struct_sOrg.sOrg.(group).(refNeuron),3);
    if VREnd
        numRefs = size(bordersSectors.(group),1);
    else
        numRefs = size(bordersSectors.(group),1)-1;
    end
    for n = 1:numNeurons
        neuron = neurons{n};
        sOrg.(group).(neuron) = NaN(maxNumFields,numRefs,maxNumTrials);
        if ~(all(all(all(isnan(in_struct_sOrg.sOrg.(group).(neuron))))))
            fprintf('Processing %s\n', neuron)
            for t = 1:maxNumTrials
                fprintf('trial: %i\n', t)
                for r = 1:numRefs
                    fprintf('ref: %i\n', r)
                    for f = 1:maxNumFields
                        fprintf('field: %i\n', f)
                        if ~isnan(in_struct_sOrg.sOrg.(group).(neuron)(f,r,t))
                            sOrg.(group).(neuron)(f,r,t) = bordersSectors.(group)(r,2) - in_struct_sOrg.sOrg.(group).(neuron)(f,r,t);
                        end
                    end
                end
            end
        end
    end
end
filenameTrialSuperOrg = 'sOrg';
save([outPath, filenameTrialSuperOrg],'sOrg')
end