.. currentmodule:: brian

.. index::
   pair: example usage; NeuronGroup
   pair: example usage; run
   pair: example usage; raster_plot
   pair: example usage; Connection
   pair: example usage; SpikeMonitor
   pair: example usage; Refractoriness

.. _example-misc_named_threshold:

Example: named_threshold (misc)
===============================

Example with named threshold and reset variables

::

    from brian import *
    eqs = '''
    dge/dt = -ge/(5*ms) : volt
    dgi/dt = -gi/(10*ms) : volt
    dx/dt = (ge+gi-(x+49*mV))/(20*ms) : volt
    '''
    P = NeuronGroup(4000, model=eqs, threshold='x>-50*mV', \
                  reset=Refractoriness(-60 * mV, 5 * ms, state='x'))
    #P=NeuronGroup(4000,model=eqs,threshold=Threshold(-50*mV,state='x'),\
    #              reset=Reset(-60*mV,state='x')) # without refractoriness
    P.x = -60 * mV
    Pe = P.subgroup(3200)
    Pi = P.subgroup(800)
    Ce = Connection(Pe, P, 'ge', weight=1.62 * mV, sparseness=0.02)
    Ci = Connection(Pi, P, 'gi', weight= -9 * mV, sparseness=0.02)
    M = SpikeMonitor(P)
    run(1 * second)
    raster_plot(M)
    show()