def create_models(self):
self.generic_cell = IafCell(id="generic_iaf_cell",
C = self.get_bioparameter("iaf_C").value,
thresh = self.get_bioparameter("iaf_thresh").value,
reset = self.get_bioparameter("iaf_reset").value,
leak_conductance = self.get_bioparameter("iaf_conductance").value,
leak_reversal = self.get_bioparameter("iaf_leak_reversal").value)
self.exc_syn = ExpTwoSynapse(id="exc_syn",
gbase = self.get_bioparameter("chem_exc_syn_gbase").value,
erev = self.get_bioparameter("chem_exc_syn_erev").value,
tau_decay = self.get_bioparameter("chem_exc_syn_decay").value,
tau_rise = self.get_bioparameter("chem_exc_syn_rise").value)
self.inh_syn = ExpTwoSynapse(id="inh_syn",
gbase = self.get_bioparameter("chem_inh_syn_gbase").value,
erev = self.get_bioparameter("chem_inh_syn_erev").value,
tau_decay = self.get_bioparameter("chem_inh_syn_decay").value,
tau_rise = self.get_bioparameter("chem_inh_syn_rise").value)
self.elec_syn = ExpTwoSynapse(id="elec_syn",
gbase = self.get_bioparameter("elec_syn_gbase").value,
erev = self.get_bioparameter("elec_syn_erev").value,
tau_decay = self.get_bioparameter("elec_syn_decay").value,
tau_rise = self.get_bioparameter("elec_syn_rise").value)
self.offset_current = PulseGenerator(id="offset_current",
delay= self.get_bioparameter("unphysiological_offset_current_del").value,
duration= self.get_bioparameter("unphysiological_offset_current_dur").value,
amplitude= self.get_bioparameter("unphysiological_offset_current").value)
def create_generic_neuron_cell(self):
self.generic_neuron_cell = IafCell(id="generic_neuron_iaf_cell",
C = self.get_bioparameter("neuron_iaf_C").value,
thresh = self.get_bioparameter("neuron_iaf_thresh").value,
reset = self.get_bioparameter("neuron_iaf_reset").value,
leak_conductance = self.get_bioparameter("neuron_iaf_conductance").value,
leak_reversal = self.get_bioparameter("neuron_iaf_leak_reversal").value)
def s(doc, net, n):
i = IafCell(id=n.name(),
C="1.0 nF",
thresh="-50mV",
reset="-65mV",
leak_conductance="10 nS",
leak_reversal="-65mV")
pop = Population(id="popOf" + n.name(), component=n.name(), size=1)
doc.iaf_cells.append(i)
net.populations.append(pop)
from neuroml import NeuroMLDocument
from neuroml import IafCell
from neuroml import Network
from neuroml import ExpOneSynapse
from neuroml import Population
from neuroml import PulseGenerator
from neuroml import ExplicitInput
from neuroml import SynapticConnection
import neuroml.writers as writers
from random import random
nml_doc = NeuroMLDocument(id="IafNet")
IafCell0 = IafCell(id="iaf0",
C="1.0 nF",
thresh="-50mV",
reset="-65mV",
leak_conductance="10 nS",
leak_reversal="-65mV")
nml_doc.iaf_cells.append(IafCell0)
IafCell1 = IafCell(id="iaf1",
C="1.0 nF",
thresh="-50mV",
reset="-65mV",
leak_conductance="20 nS",
leak_reversal="-65mV")
nml_doc.iaf_cells.append(IafCell1)
syn0 = ExpOneSynapse(id="syn0", gbase="65nS", erev="0mV", tau_decay="3ms")
chem_exc_syn_rise = BioParameter("chem_exc_syn_rise", "3ms", "BlindGuess", "0.1")
chem_exc_syn_decay = BioParameter("chem_exc_syn_decay", "10ms", "BlindGuess", "0.1")
chem_inh_syn_gbase = BioParameter("chem_inh_syn_gbase", "0.2nS", "BlindGuess", "0.1")
chem_inh_syn_erev = BioParameter("chem_inh_syn_erev", "-80mV", "BlindGuess", "0.1")
chem_inh_syn_rise = BioParameter("chem_inh_syn_rise", "3ms", "BlindGuess", "0.1")
chem_inh_syn_decay = BioParameter("chem_inh_syn_decay", "10ms", "BlindGuess", "0.1")
unphysiological_offset_current = BioParameter("unphysiological_offset_current", "0.21nA", "KnownError", "0")
unphysiological_offset_current_dur = BioParameter("unphysiological_offset_current_dur", "200ms", "KnownError", "0")
generic_cell = IafCell(id="generic_iaf_cell",
C = iaf_C.value,
thresh = iaf_thresh.value,
reset = iaf_reset.value,
leak_conductance = iaf_conductance.value,
leak_reversal = iaf_leak_reversal.value)
exc_syn = ExpTwoSynapse(id="exc_syn",
gbase = chem_exc_syn_gbase.value,
erev = chem_exc_syn_erev.value,
tau_decay = chem_exc_syn_decay.value,
tau_rise = chem_exc_syn_rise.value)
inh_syn = ExpTwoSynapse(id="inh_syn",
gbase = chem_inh_syn_gbase.value,
erev = chem_inh_syn_erev.value,
tau_decay = chem_inh_syn_decay.value,