add_simulated_spikes_to_Cells.m

function Cells = add_simulated_spikes_to_Cells(Cells,params,spikes)
    % spikes is the output of fit_glm_to_Cells    
    Cells.Trials = remove_trials(Cells.Trials,params.exclude_trials);
    nTrials = numel(Cells.Trials.pokedR);
    Cells.spike_time_s = remove_trials_from_spikes(Cells.spike_time_s,params.exclude_trials);  
    expt=params.dm.dspec.expt;
    time_s = (expt.binSize:expt.binSize:expt.trial(1).duration_s) - expt.binSize/2 - expt.trial(1).(params.ref_event);       
    for i=1:numel(spikes)
        cellno = spikes(i).cellno;
        if params.distribution=="poisson"        
            spikes(i).samples = poissrnd(spikes(i).Yhat);
        else
            error('not implemented yet');
        end
        Cells.spike_time_s.(params.ref_event){cellno} = cell(nTrials,1);
        for t=1:nTrials
            ndx = buildGLM.getSpikeIndicesforTrial(expt,t);
            these_spikes = spikes(i).samples(ndx);
            lim = max(these_spikes);              
            for k=1:lim
                Cells.spike_time_s.(params.ref_event){cellno}{t} = [Cells.spike_time_s.(params.ref_event){cellno}{t} time_s(these_spikes>=k)];
            end
            Cells.spike_time_s.(params.ref_event){cellno}{t}       = sort(Cells.spike_time_s.(params.ref_event){cellno}{t}');      
        end
    end
        



end