.. DO NOT EDIT.
.. THIS FILE WAS AUTOMATICALLY GENERATED BY SPHINX-GALLERY.
.. TO MAKE CHANGES, EDIT THE SOURCE PYTHON FILE:
.. "auto_examples/brette_et_al_2007/cuba.py"
.. LINE NUMBERS ARE GIVEN BELOW.

.. only:: html

    .. note::
        :class: sphx-glr-download-link-note

        :ref:`Go to the end <sphx_glr_download_auto_examples_brette_et_al_2007_cuba.py>`
        to download the full example code.

.. rst-class:: sphx-glr-example-title

.. _sphx_glr_auto_examples_brette_et_al_2007_cuba.py:


Benchmark 2 of the simulator review (CUBA)
-------------------------------------------

.. only:: html

----

 Run this example as a Jupyter notebook:

  .. card::
    :width: 25%
    :margin: 2
    :text-align: center
    :link: https://lab.ebrains.eu/hub/user-redirect/git-pull?repo=https%3A%2F%2Fgithub.com%2Fnest%2Fnest-simulator-examples&urlpath=lab%2Ftree%2Fnest-simulator-examples%2Fnotebooks%2Fnotebooks%2Fbrette_et_al_2007%2Fcuba.ipynb&branch=main
    :link-alt: JupyterHub service

    .. image:: https://nest-simulator.org/TryItOnEBRAINS.png


.. grid:: 1 1 1 1
   :padding: 0 0 2 0

   .. grid-item::
     :class: sd-text-muted
     :margin: 0 0 3 0
     :padding: 0 0 3 0
     :columns: 4

     See :ref:`our guide <run_jupyter>` for more information and troubleshooting.

----


This script creates a sparsely coupled network of excitatory and inhibitory
neurons which exhibits self-sustained activity after an initial stimulus.
Connections within and across both populations are created at random. Both
neuron populations receive Poissonian background input. The spike output of
500 neurons from each population are recorded. Neurons are modeled as leaky
integrate-and-fire neurons with current-based synapses (exponential functions).
The model is based on the Vogels & Abbott network model [1]_.

This is Benchmark 2 of the FACETS simulator review (Brette et al., 2007) [2]_:
- Neuron model: integrate-and-fire (``iaf_psc_exp``)
- Synapse model: current-based (CUBA)
- Synapse time course: exponential
- Spike times: grid-constrained

Note that the simulation time is set 10 times longer than for benchmark 1
(coba.py) or benchmark 3 (hh_coba.py). This is necessary as the computational
load here is much lower, so a longer simulation time is necessary to make
reasonable measurements.

References
~~~~~~~~~~

.. [1] Vogels TP, Abbott LF. 2005. Signal propagation and logic gating in
       networks of integrate-and-fire neurons. Journal of Neuroscience.
       25(46):10786-10795.
       https://doi.org/10.1523/JNEUROSCI.3508-05.2005

.. [2] Brette R, Rudolph M, Carnevale T, Hines M, Beeman D, Bower JM, et al.
       2007. Simulation of networks of spiking neurons: a review of tools and
       strategies. Journal of Computational Neuroscience. 23(3):349-398.
       https://doi.org/10.1007/s10827-007-0038-6

.. GENERATED FROM PYTHON SOURCE LINES 59-63

.. code-block:: Python


    import nest
    from brette_et_al_2007_benchmark import run_simulation


.. GENERATED FROM PYTHON SOURCE LINES 64-65

Set benchmark parameters

.. GENERATED FROM PYTHON SOURCE LINES 65-109

.. code-block:: Python


    params = {
        "model": "iaf_psc_exp",  # neuron model
        "model_params": {
            "E_L": -49.0,  # resting membrane potential [mV]
            # see Brette et al, J Comput Neurosci 23:349 (2007), p 393
            "V_m": -49.0,  # initial membrane potential [mV]
            "V_th": -50.0,  # threshold [mV]
            "V_reset": -60.0,  # reset potential [mV]
            "C_m": 200.0,  # capacity of the membrane [pF]
            "tau_m": 20.0,  # membrane time constant [ms]
            "tau_syn_ex": 5.0,  # time const. postsynaptic excitatory currents [ms]
            "tau_syn_in": 10.0,  # time const. postsynaptic inhibitory currents [ms]
            "t_ref": 5.0,  # duration of refractory period [ms]
        },
        "delay": 0.1,  # synaptic delay [ms]
        "E_synapse_params": {
            "weight": 16.2,  # excitatory PSC amplitude [pA]
        },
        "I_synapse_params": {
            "weight": -139.5,  # inhibitory PSC amplitude [pA]
        },
        "stimulus": "poisson_generator",
        "stimulus_params": {
            "rate": 300.0,  # rate of initial Poisson stimulus [spikes/s]
            "start": 1.0,  # start of Poisson generator [ms]
            "stop": 51.0,  # stop of Poisson generator [ms]
            "origin": 0.0,  # origin of time [ms]
        },
        "recorder": "spike_recorder",
        "recorder_params": {
            "record_to": "ascii",
            "label": "cuba",
        },
        "Nrec": 500,  # number of neurons per population to record from
        "Nstim": 50,  # number of neurons to stimulate
        "simtime": 10000.0,  # simulated time [ms] (10x longer than coba/hh_coba)
        "dt": 0.1,  # simulation step [ms]
        "NE": 3200,  # number of excitatory neurons
        "NI": 800,  # number of inhibitory neurons
        "epsilon": 0.02,  # connection probability
        "virtual_processes": 1,  # number of virtual processes to use
    }


.. GENERATED FROM PYTHON SOURCE LINES 110-111

Run the simulation

.. GENERATED FROM PYTHON SOURCE LINES 111-115

.. code-block:: Python


    if __name__ == "__main__":
        nest.set_verbosity("M_WARNING")
        results = run_simulation(params)


.. _sphx_glr_download_auto_examples_brette_et_al_2007_cuba.py:

.. only:: html

  .. container:: sphx-glr-footer sphx-glr-footer-example

    .. container:: sphx-glr-download sphx-glr-download-jupyter

      :download:`Download Jupyter notebook: cuba.ipynb <cuba.ipynb>`

    .. container:: sphx-glr-download sphx-glr-download-python

      :download:`Download Python source code: cuba.py <cuba.py>`

    .. container:: sphx-glr-download sphx-glr-download-zip

      :download:`Download zipped: cuba.zip <cuba.zip>`


.. only:: html

 .. rst-class:: sphx-glr-signature

    `Gallery generated by Sphinx-Gallery <https://sphinx-gallery.github.io>`_