def synapses(self):
syn_ = h.MyExp2Syn(self.dend2(0.5))
syn_.tau1 = 0.5
syn_.tau2 = 3
self.prelist.append(syn_)
syn_ = h.MyExp2Syn(self.dend1(0.5))
syn_.tau1 = 0.5
syn_.tau2 = 3
self.prelist.append(syn_)
syn_ = h.MyExp2Syn(self.soma(0.5))
syn_.tau1 = 1
syn_.tau2 = 8
syn_.e = -75
self.prelist.append(syn_)
syn_ = h.MyExp2Syn(self.soma(0.5))
syn_.tau1 = 35
syn_.tau2 = 100
syn_.e = -75
def addSynapses(self,myTauValue):
# Define synapses in various areas of the cell
# and set the connection parameters (synapse
# rise time, decay time, reversal potential,
# conductance/weight).
# Also create lists of synapses onto the model cell:
# excsyn_list: excitatory synapses onto the cell
# inhdendsyn_list: inhibitory synapses onto the cell dendrites
# inhsomasyn_list: inhibitory synapses onto the cell body objref preInhSoma_list, preInhDend_list, preExcDend_list
self.preInhSoma_list = []
self.preInhDend_list = []
self.preExcDend_list = []
s=0
self.recInhSomaCurrent = []
for sec in self.inhsomasyn_list:
syn_ = h.MyExp2Syn(sec(0.5))
self.preInhSoma_list.append(syn_) # AMPA EC
syn_.tau1 = myTauValue #0.5
syn_.tau2 = 3
syn_.e = -70
self.recInhSomaCurrent.append(h.Vector())
self.recInhSomaCurrent[s].record(self.preInhSoma_list[s]._ref_i)
#sprint(cmdstr,"objref recInhSomaCurrent%d", s)
#{execute(cmdstr)}
#sprint(cmdstr,"recInhSomaCurrent%d = new Vector()", s)
#{execute(cmdstr)}
#sprint(cmdstr,"recInhSomaCurrent%d.record(&preInhSoma_list.object(%d).i)", s, s)
#{execute(cmdstr)}
s = s + 1
totInhSoma = s
s = 0
self.recInhDendCurrent = []
for sec in self.inhdendsyn_list:
syn_ = h.MyExp2Syn(sec(0.5))
self.preInhDend_list.append(syn_) # AMPA EC
syn_.tau1 = 0.5
syn_.tau2 = 3
syn_.e = -70
self.recInhDendCurrent.append(h.Vector())
self.recInhDendCurrent[s].record(self.preInhDend_list[s]._ref_i)
#sprint(cmdstr,"recInhDendCurrent%d = new Vector()", s)
#sprint(cmdstr,"recInhDendCurrent%d.record(&preInhDend_list.object(%d).i)", s, s)
s = s + 1
totInhDend = s
s = 0
self.recExcCurrent = []
for sec in self.excsyn_list:
syn_ = h.MyExp2Syn(sec(0.5))
self.preExcDend_list.append(syn_) # AMPA EC
syn_.tau1 = 1
syn_.tau2 = 5
syn_.e = 0
self.recExcCurrent.append(h.Vector())
self.recExcCurrent[s].record(self.preExcDend_list[s]._ref_i)
#sprint(cmdstr,"objref recExcCurrent%d", s)
#{execute(cmdstr)}
#sprint(cmdstr,"recExcCurrent%d = new Vector()", s)
#{execute(cmdstr)}
#sprint(cmdstr,"recExcCurrent%d.record(&preExcDend_list.object(%d).i)", s, s)
#{execute(cmdstr)}
s = s + 1
self.totExc = s
self.excitatory_syn_weight = 0.005 # the maximum synaptic conductance in microSiemens, aka the synaptic amplitude, of the excitatory connections
self.inhDend_syn_weight = 0.003 # the maximum synaptic conductance of the inhibitory connections onto the dendrites
self.inhSoma_syn_weight = 0.006 # the maximum synaptic conductance of the inhibitory connections onto the soma
def addSynapses(self):
self.pre_list = []
# E0
syn_ = h.MyExp2Syn(self.lmM1(0.5))
self.pre_list.append(syn_) # AMPA EC
syn_.tau1 = 0.5
syn_.tau2 = 3
syn_.e = 0
# E1
syn_ = h.MyExp2Syn(self.lmM2(0.5))
self.pre_list.append(syn_) # AMPA EC (not used)
syn_.tau1 = 0.5
syn_.tau2 = 3
syn_.e = 0
# E2
syn_ = h.MyExp2Syn(self.radM1(0.5))
self.pre_list.append(syn_) # AMPA CA3 Shaffer collateral
syn_.tau1 = 0.5
syn_.tau2 = 3
syn_.e = 0
# E3
syn_ = h.MyExp2Syn(self.radM2(0.5))
self.pre_list.append(syn_) # AMPA CA3 Shaffer collateral
syn_.tau1 = 0.5
syn_.tau2 = 3
syn_.e = 0
# E4
syn_ = h.MyExp2Syn(self.radT1(0.5))
self.pre_list.append(syn_) # AMPA CA3 Shaffer collateral
syn_.tau1 = 0.5
syn_.tau2 = 3
syn_.e = 0
# E5
syn_ = h.MyExp2Syn(self.radT2(0.5))
self.pre_list.append(syn_) # AMPA CA3 Shaffer collateral
syn_.tau1 = 0.5
syn_.tau2 = 3
syn_.e = 0
# E6
syn_ = h.MyExp2Syn(self.oriT1(0.5))
self.pre_list.append(syn_) # AMPA PC
syn_.tau1 = 0.5
syn_.tau2 = 3
syn_.e = 0
# E7
syn_ = h.MyExp2Syn(self.oriT2(0.5))
self.pre_list.append(syn_) # AMPA PC
syn_.tau1 = 0.5
syn_.tau2 = 3
syn_.e = 0
# I8
syn_ = h.MyExp2Syn(self.soma(0.5))
self.pre_list.append(syn_) # GABA-A Neighboring axo-axonic cell
syn_.tau1 = 1
syn_.tau2 = 8
syn_.e = -75
# I9
syn_ = h.MyExp2Syn(self.soma(0.6))
self.pre_list.append(syn_) # GABA-A Bistratified cell
syn_.tau1 = 1
syn_.tau2 = 8
syn_.e = -75
# I10
syn_ = h.MyExp2Syn(self.oriT1(0.6))
self.pre_list.append(syn_) # GABA-A Septum
syn_.tau1 = 1
syn_.tau2 = 8
syn_.e = -75
# I11
syn_ = h.MyExp2Syn(self.oriT2(0.6))
self.pre_list.append(syn_) # GABA-A Septum
syn_.tau1 = 1
syn_.tau2 = 8
syn_.e = -75
# I12
syn_ = h.MyExp2Syn(self.oriT1(0.6))
self.pre_list.append(syn_) # GABA-B Septum
syn_.tau1 = 35
syn_.tau2 = 100
syn_.e = -75
# I13
syn_ = h.MyExp2Syn(self.oriT2(0.6))
self.pre_list.append(syn_) # GABA-B Septum
syn_.tau1 = 35
syn_.tau2 = 100
syn_.e = -75
def addSynapses(self):
self.pre_list = []
# E0
syn_ = h.MyExp2Syn(self.lm_thick1(0.5))
self.pre_list.append(syn_) # AMPA EC
syn_.tau1 = 0.5
syn_.tau2 = 3
syn_.e = 0
# E1
syn_ = h.MyExp2Syn(self.lm_thick2(0.5))
self.pre_list.append(syn_) # AMPA EC
syn_.tau1 = 0.5
syn_.tau2 = 3
syn_.e = 0
# E2
syn_ = h.MyExp2Syn(self.radTmed(0.5))
self.pre_list.append(syn_) # AMPA CA3 Shaffer collateral
syn_.tau1 = 0.5
syn_.tau2 = 3
syn_.e = 0
# E3
syn_ = h.NMDA(self.radTmed(0.5))
self.pre_list.append(syn_) # NMDA CA3 Shaffer collateral
syn_.tcon = 2.3
syn_.tcoff = 100
syn_.gNMDAmax = 1 # use connection weight to determine max cond
# E4
syn_ = h.MyExp2Syn(self.radTprox(0.5))
self.pre_list.append(syn_) # AMPA PC Recurrent collateral
syn_.tau1 = 0.5
syn_.tau2 = 3
syn_.e = 0
# I5
syn_ = h.MyExp2Syn(self.soma(0.5))
self.pre_list.append(syn_) # GABA-A basket cell
syn_.tau1 = 1
syn_.tau2 = 8
syn_.e = -75
# I6
syn_ = h.MyExp2Syn(self.axon(0.1))
self.pre_list.append(syn_) # GABA-A AA cell
syn_.tau1 = 1
syn_.tau2 = 8
syn_.e = -75
# I7
syn_ = h.MyExp2Syn(self.lm_thick1(0.5))
self.pre_list.append(syn_) # GABA-A OLM cell
syn_.tau1 = 1
syn_.tau2 = 8
syn_.e = -75
# I8
syn_ = h.MyExp2Syn(self.lm_thick2(0.5))
self.pre_list.append(syn_) # GABA-A OLM cell
syn_.tau1 = 1
syn_.tau2 = 8
syn_.e = -75
# I9
syn_ = h.MyExp2Syn(self.lm_thick1(0.5))
self.pre_list.append(syn_) # GABA-B OLM cell
syn_.tau1 = 35
syn_.tau2 = 100
syn_.e = -75
# I10
syn_ = h.MyExp2Syn(self.lm_thick2(0.5))
self.pre_list.append(syn_) # GABA-B OLM Cell
syn_.tau1 = 35
syn_.tau2 = 100
syn_.e = -75
# I11
syn_ = h.MyExp2Syn(self.radTmed(0.8))
self.pre_list.append(syn_) # GABA-A Bistratified
syn_.tau1 = 1
syn_.tau2 = 8
syn_.e = -75
# I12
syn_ = h.MyExp2Syn(self.radTmed(0.7))
self.pre_list.append(syn_) # GABA-A Bistratified
syn_.tau1 = 1
syn_.tau2 = 8
syn_.e = -75
# I13
syn_ = h.MyExp2Syn(self.radTmed(0.6))
self.pre_list.append(syn_) # GABA-A Bistratified
syn_.tau1 = 1
syn_.tau2 = 8
syn_.e = -75
# I14
syn_ = h.MyExp2Syn(self.radTmed(0.4))
self.pre_list.append(syn_) # GABA-A Bistratified
syn_.tau1 = 1
syn_.tau2 = 8
syn_.e = -75
# I15
syn_ = h.MyExp2Syn(self.radTmed(0.3))
self.pre_list.append(syn_) # GABA-A Bistratified
syn_.tau1 = 1
syn_.tau2 = 8
syn_.e = -75
# I16
syn_ = h.MyExp2Syn(self.radTmed(0.2))
self.pre_list.append(syn_) # GABA-A Bistratified
syn_.tau1 = 1
syn_.tau2 = 8
syn_.e = -75
# I17
syn_ = h.MyExp2Syn(self.radTmed(0.8))
self.pre_list.append(syn_) # GABA-B Bistratified
syn_.tau1 = 35
syn_.tau2 = 100
syn_.e = -75
# I18
syn_ = h.MyExp2Syn(self.radTmed(0.7))
self.pre_list.append(syn_) # GABA-B Bistratified
syn_.tau1 = 35
syn_.tau2 = 100
syn_.e = -75
# I19
syn_ = h.MyExp2Syn(self.radTmed(0.6))
self.pre_list.append(syn_) # GABA-B Bistratified
syn_.tau1 = 35
syn_.tau2 = 100
syn_.e = -75
# I20
syn_ = h.MyExp2Syn(self.radTmed(0.4))
self.pre_list.append(syn_) # GABA-B Bistratified
syn_.tau1 = 35
syn_.tau2 = 100
syn_.e = -75
# I21
syn_ = h.MyExp2Syn(self.radTmed(0.3))
self.pre_list.append(syn_) # GABA-B Bistratified
syn_.tau1 = 35
syn_.tau2 = 100
syn_.e = -75
# I22
syn_ = h.MyExp2Syn(self.radTmed(0.2))
self.pre_list.append(syn_) # GABA-B Bistratified
syn_.tau1 = 35
syn_.tau2 = 100
syn_.e = -75
# I23
syn_ = h.STDPE2(self.radTmed(0.5))
self.pre_list.append(syn_) # AMPA modifiable CA3 Schaffer collaterals
syn_.tau1 = 0.5
syn_.tau2 = 3
syn_.e = 0
def synapses(self):
syn_ = h.MyExp2Syn(self.radM1(0.5))
syn_.tau1 = 0.5
syn_.tau2 = 3
syn_.e = 0
self.prelist.append(syn_)
syn_ = h.MyExp2Syn(self.radM2(0.5))
syn_.tau1 = 0.5
syn_.tau2 = 3
syn_.e = 0
self.prelist.append(syn_)
syn_ = h.MyExp2Syn(self.radT1(0.5))
syn_.tau1 = 0.5
syn_.tau2 = 3
syn_.e = 0
self.prelist.append(syn_)
syn_ = h.MyExp2Syn(self.radT2(0.5))
syn_.tau1 = 0.5
syn_.tau2 = 3
syn_.e = 0
self.prelist.append(syn_)
syn_ = h.MyExp2Syn(self.oriT1(0.5))
syn_.tau1 = 0.5
syn_.tau2 = 3
syn_.e = 0
self.prelist.append(syn_)
syn_ = h.MyExp2Syn(self.oriT2(0.5))
syn_.tau1 = 0.5
syn_.tau2 = 3
syn_.e = 0
self.prelist.append(syn_)
syn_ = h.MyExp2Syn(self.soma(0.5))
syn_.tau1 = 1
syn_.tau2 = 8
syn_.e = -75
self.prelist.append(syn_)
syn_ = h.MyExp2Syn(self.soma(0.6))
syn_.tau1 = 1
syn_.tau2 = 8
syn_.e = -75
self.prelist.append(syn_)
syn_ = h.MyExp2Syn(self.oriT1(0.6))
syn_.tau1 = 1
syn_.tau2 = 8
syn_.e = -75
self.prelist.append(syn_)
syn_ = h.MyExp2Syn(self.oriT2(0.6))
syn_.tau1 = 1
syn_.tau2 = 8
syn_.e = -75
self.prelist.append(syn_)
syn_ = h.MyExp2Syn(self.oriT1(0.6))
syn_.tau1 = 35
syn_.tau2 = 100
syn_.e = -75
self.prelist.append(syn_)
syn_ = h.MyExp2Syn(self.oriT2(0.6))
syn_.tau1 = 35
syn_.tau2 = 100
syn_.e = -75
self.prelist.append(syn_)
def synapses(self):
syn_ = h.MyExp2Syn(self.lm_thick1(0.5))
syn_.tau1 = 0.5
syn_.tau2 = 3.0
syn_.e = 0.0
self.prelist.append(syn_)
syn_ = h.MyExp2Syn(self.lm_thick2(0.5))
syn_.tau1 = 0.5
syn_.tau2 = 3.0
syn_.e = 0.0
self.prelist.append(syn_)
syn_ = h.MyExp2Syn(self.radTmed(0.5))
syn_.tau1 = 0.5
syn_.tau2 = 3.0
syn_.e = 0.0
self.prelist.append(syn_)
syn_ = h.NMDA(self.radTmed(0.5))
syn_.tcon = 2.3
syn_.tcoff = 100.
syn_.gNMDAmax = 1.0
self.prelist.append(syn_)
syn_ = h.MyExp2Syn(self.radTprox(0.5))
syn_.tau1 = 0.5
syn_.tau2 = 3.0
syn_.e = 0.0
self.prelist.append(syn_)
syn_ = h.MyExp2Syn(self.soma(0.5))
syn_.tau1 = 1.0
syn_.tau2 = 8.0
syn_.e = -75
self.prelist.append(syn_)
syn_ = h.MyExp2Syn(self.axon(0.5))
syn_.tau1 = 1.0
syn_.tau2 = 8.0
syn_.e = -75.0
self.prelist.append(syn_)
syn_ = h.MyExp2Syn(self.lm_thick1(0.5))
syn_.tau1 = 1.0
syn_.tau2 = 8.0
syn_.e = -75
self.prelist.append(syn_)
syn_ = h.MyExp2Syn(self.lm_thick2(0.5))
syn_.tau1 = 1.0
syn_.tau2 = 8.0
syn_.e = -75
self.prelist.append(syn_)
syn_ = h.MyExp2Syn(self.lm_thick1(0.5))
syn_.tau1 = 35
syn_.tau2 = 100
syn_.e = -75
self.prelist.append(syn_)
syn_ = h.MyExp2Syn(self.lm_thick2(0.5))
syn_.tau1 = 35
syn_.tau2 = 100
syn_.e = -75
self.prelist.append(syn_)
for l in np.arange(0.2, 0.9, 0.1):
if l == 0.5:
continue
syn_ = h.MyExp2Syn(self.radTmed(l))
syn_.tau1 = 1.0
syn_.tau2 = 8.0
syn_.e = -75
self.prelist.append(syn_)
for l in np.arange(0.2, 0.9, 0.1):
if l == 0.5:
continue
syn_ = h.MyExp2Syn(self.radTmed(l))
syn_.tau1 = 35
syn_.tau2 = 100
syn_.e = -75
self.prelist.append(syn_)
syn_ = h.STDPE2(self.radTmed(0.5))
syn_.tau1 = 0.5
syn_.tau2 = 3.0
syn_.e = 0.0