def run_simulation_CL1(tauE_, wE_, tauI_, wI_, Group, Synapses, Target, time,
dt_, delta, rate_correction):
restore('initialised')
# set the parameters
Group.set_states({'tauE': tauE_ * ms, 'tauI': tauI_ * ms})
print(f'taueE: {tauE_} - tauI {tauI_}')
Synapses.set_states({'wE': wE_ * nS, 'wI': wI_ * nS})
print(f'wE: {wE_} - wI {wI_}')
_, model_spike_monitor = nc.add_monitors(Group)
target_spike_monitor = SpikeMonitor(Target,
name='CL1_target_spike_monitor')
run(time)
gf = metric.compute_gamma_factor(model_spike_monitor,
target_spike_monitor,
time,
dt_=dt_,
delta=delta,
rate_correction=rate_correction)
print(f'Gamma factor: {gf}')
return gf
# params_TL2['tauI'] = 1 * ms
# synapses_TL2['wI'] = 300 * nS
print(params_TL2)
print(synapses_TL2)
# Neuron group
G_TL2 = nc.generate_neuron_groups(N_TL2,
eqs,
threshold_eqs,
reset_eqs,
params_TL2,
name='TL2_test')
# Add monitors
STM_TL2, SPM_TL2 = nc.add_monitors(G_TL2, name='TL2_test')
# Connect heading to TL2
S_P_HEADING_TL2 = nc.connect_synapses(P_HEADING,
G_TL2,
W_HEADING_TL2,
params=synapses_TL2,
model=synapses_model,
on_pre=synapses_eqs_ex)
# Run simulation
run(T_outbound * time_step * ms)
cx_spiking.plotting.plot_rate_cx_log_spikes(
cx_log.tl2,
TL2_spike_rates,
# params_TN2['tauI'] = 1 * ms
# synapses_TN2['wI'] = 300 * nS
print(params_TN2)
print(synapses_TN2)
# Neuron group
G_TN2 = nc.generate_neuron_groups(N_TN2,
eqs,
threshold_eqs,
reset_eqs,
params_TN2,
name='TN2_test')
# Add monitors
STM_TN2, SPM_TN2 = nc.add_monitors(G_TN2, name='TN2_test')
# Connect heading to TN2
S_P_FLOW_TN2 = nc.connect_synapses(P_FLOW,
G_TN2,
W_FLOW_TN2,
params=synapses_TN2,
model=synapses_model,
on_pre=synapses_eqs_ex)
# Run simulation
run(T_outbound * time_step * ms)
cx_spiking.plotting.plot_rate_cx_log_spikes(
cx_log.tn2,
TN2_spike_rates_max,
# Neuron group
G_TL2 = nc.generate_neuron_groups(N_TL2,
eqs,
threshold_eqs,
reset_eqs,
TL2_neuron_params,
name='TL2_test')
G_CL1 = nc.generate_neuron_groups(N_CL1,
eqs,
threshold_eqs,
reset_eqs,
neuron_params,
name='CL1_test')
# Add monitors
STM_TL2, SPM_TL2 = nc.add_monitors(G_TL2, name='TL2_test')
STM_CL1, SPM_CL1 = nc.add_monitors(G_CL1, name='CL1_test')
# Connect heading to TL2
S_P_HEADING_TL2 = nc.connect_synapses(P_HEADING,
G_TL2,
W_HEADING_TL2,
model=synapses_model,
params=TL2_synapses_params,
on_pre=synapses_eqs_ex)
S_TL2_CL1 = nc.connect_synapses(G_TL2,
G_CL1,
W_TL2_CL1,
model=synapses_model,
params=synapses_params,
on_pre=synapses_eqs_ex)
