HVC_DA/default.py at main · pearsonlab/HVC_DA · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
from brian2 import *

# use a large number for inf

eq_dV = '''
dV/dt = (-gL * (V - EL) + psi - w + I_ext_%s(t,i)) / Cm : volt (unless refractory)
'''

eq_dwE = '''
dw/dt = (a * clip(V - EL, 0 * mV, 1000 * mV) - w) / tauw : amp
'''

eq_dwI = '''
dw/dt = (a * clip(V - EL, -1000 * mV, 0 * mV) - w) / tauw : amp
'''

eq_psi = '''
psi = gL * DeltaT * exp((V - VT) / DeltaT) : amp
'''

eqs_reset = 'V = Vr; w += b'

pop_params = '''
gL : siemens
EL : volt
Cm : farad
Vr : volt
DeltaT : volt
VT : volt
a : siemens
b : amp
tauw : second
'''

def assign_default_params_Exc(E_pop):
    E_pop.Cm = 150 * pF
    E_pop.gL = 10 * nS
    E_pop.EL = -80 * mV
    E_pop.VT = -50 * mV
    E_pop.DeltaT = 2 * mV
    E_pop.a  = -6 * nS
    E_pop.tauw = 120 * ms
    E_pop.b  = 0.05 * nA
    E_pop.Vr = -65 * mV

def assign_default_params_Inh(I_pop):
    I_pop.Cm = 160 * pF
    I_pop.gL = 8 * nS
    I_pop.EL = -60 * mV
    I_pop.VT = -50 * mV
    I_pop.DeltaT = 2 * mV
    I_pop.a  = 7 * nS
    I_pop.tauw = 50 * ms
    I_pop.b  = 0 * nA
    I_pop.Vr = -70 * mV

def smooth_AP(v, th=-35):
    ''' To smooth the uneven AP heights due to finite time step
    v needs to be unitless
    '''
    v[v > th] = -20
    return v