import neuron
import plotting_convention
plt.close('all')
###############################################################################
# Main script, set parameters and create cell, synapse and electrode objects
###############################################################################
lambdaf = 200.
# folder = "morphologies/cell_hallermann_myelin"
folder = "morphologies/EyalEtAl2016"
# folder = "morphologies/L23_PC_cADpyr229_5"
# get the template name
f = file(join(folder, "ActiveModels/model_0603_cell08_cm045.hoc"), 'r')
templatename = utils.get_templatename(f)
f.close()
neuron.load_mechanisms(join(folder, "mechanisms"))
add_synapses = False
cell = LFPy.TemplateCell(
morphology=join(folder, "morphs/2013_03_06_cell08_876_H41_05_Cell2.ASC"),
# cell = LFPy.TemplateCell(morphology=join(folder, "morphology/dend-C260897C-P3_axon-C220797A-P3_-_Clone_0.asc"),
templatefile=join(folder, 'ActiveModels/model_0603_cell08_cm045.hoc'),
templatename=templatename,
templateargs=1 if add_synapses else 0,
tstop=50.,
dt=2.**-4,
extracellular=True,
tstart=-50,
lambda_f=lambdaf,
# morph = join('morphologies', 'axon.hoc') # Mainen&Sejnowski, 1996
# morph = join(folder, 'cell_simple.hoc')
# morph = join(folder, 'cell_simple.hoc')
custom_code = [
join(folder, 'Cell parameters.hoc'),
join(folder, 'charge.hoc')
]
# join(folder, 'pruning_full.hoc')]
if RANK == 0:
# folder = "morphologies/cell_hallermann_myelin"
folder = "morphologies/L23_PC_cADpyr229_1"
# folder = "morphologies/L23_PC_cADpyr229_5"
# get the template name
f = file(join(folder, "template.hoc"), 'r')
templatename = utils.get_templatename(f)
f.close()
# get biophys template name
f = file(join(folder, "biophysics.hoc"), 'r')
biophysics = utils.get_templatename(f)
f.close()
# get morphology template name
f = file(join(folder, "morphology.hoc"), 'r')
morphology = utils.get_templatename(f)
f.close()
# f = file(join(folder, "morphology/dend-C170897A-P3_axon-C260897C-P4_-_Clone_4.asc"), 'r')
# morphasc = utils.get_templatename(f)
# f.close()
# get synapses template name