// genesis
int stc_number
int msgcount
int i, j, GABAtarget_index
str list, GABAtarget, indexstr
str output, name
int startchar, endchar
// plots for each stc how many of its PF synapses are stimulated by the stimulus
stc_number = 212 // 211 // 400 // 500 // 840 // 84 // 50 // 84 // 50
list = ""
output = "stc-boutonsEz_" @ {stc_number} @ ".dat"
echo {output}
ce /molecular_layer/Stellate[{stc_number}]/soma/spike
msgcount = {getmsg . -out -count}
echo {stc_number} {msgcount}
for (j = 0; {j < msgcount - 1}; j = j + 1)
// echo {j}
GABAtarget = {getmsg . -out -destination {j}}
startchar = {findchar {GABAtarget} [}
endchar = {findchar {GABAtarget} ]}
echo {GABAtarget} {startchar} {endchar}
if ({startchar > 0})
GABAtarget_index = {{substring {GABAtarget} {startchar + 1} {endchar - 1}}}
else GABAtarget_index = 0
end // if
echo GABAtarget_index {GABAtarget_index}
list = (list) @ {j} @ " " @ {GABAtarget_index} @ " " @ \
{getfield {GABAtarget} x} @ " " @ \
{getfield {GABAtarget} y} @ " " @ \
{getfield {GABAtarget} z} @ " " @ {chr 10}
end // for
// the soma
list = (list) @ {j} @ " " @ {stc_number} @ " " @ \
{getfield . x} @ " " @ \
{getfield . y} @ " " @ \
{getfield . z} @ " " @ {chr 10}
// and now the dendrite
foreach name ({el /molecular_layer/Stellate[{stc_number}]/d#[][TYPE=compartment]})
list = (list) @ {j} @ " " @ {stc_number} @ " " @ \
{getfield {name} x} @ " " @ \
{getfield {name} y} @ " " @ \
{getfield {name} z} @ " " @ {chr 10}
end // foreach
echo {list} > {output}