""" comparison.py Compare Jonas 1993 cell with Schmidt-Hieber 2007 """ from neuron import h from CA3.library import jonas1993 # global variables tstop = 10 Vrest = -70 jonas = jonas1993() # create VClamps on cells jonas_vclamp = h.SEClamp(0.5, sec = jonas.soma) # fill same VClamp parameters for both jonas_vclamp.rs = 0.1 jonas_vclamp.amp1 = Vrest jonas_vclamp.dur1 = tstop # set one synapse originally at the soma jonas_syn = h.synapse(0.5, sec = jonas.soma) jonas_syn.tonset = 0 jonas_syn.tau0 = 0.2 # in ms (tau onset in Jonas et al., 1993) jonas_syn.tau1 = 2.5 # in ms (tau decay in Jonas et al., 1993) jonas_syn.gmax = 300e-6 # 300 pS in Jonas et al., 1993 # run simulation
#!/usr/bin/env python """ Fig18A.py Reproduces figure 18 A in Jonas et al., 1993 """ import numpy as np import matplotlib.pyplot as plt from CA3.library import jonas1993 from neuron import h cell = jonas1993() # initial parameters # voltage-clamp at the soma VC_patch = h.SEClamp(0.5, sec = cell.soma) VC_patch.rs = 0.1 VC_patch.amp1 = -70 # set one synapse originally at the soma # for parameters check page 624 and 649 syn = h.synapse(0.5, sec = cell.soma) syn.tonset = 0 syn.tau0 = 0.2 # in ms (tau onset in Jonas et al., 1993) syn.tau1 = 2.5 # in ms (tau decay in Jonas et al., 1993) syn.gmax = 300e-6 # 300 pS in Jonas et al., 1993 def simulation(tstop, with_time = False):