Each segment of every section that has the stochastic mechanism must have its own instance of the point process. That's easy to do. If every section in your model is supposed to have the stochastic mechanism, then
- Code: Select all
objref smlist // will be a List that contains all instances of the stochastic mechanism
smlist = new List()
// change SMechName to the name of your stochastic mechanism
forall for (x,0) smlist.append(new SMechName(x))
You must also make sure that each instance gets the "correct" number of channels. The correct number will be the product
segment_area * channel_density
If channel_density is in (number of channels/square micron) then you can just do
- Code: Select all
ii = 0
forall for (x,0) {
smlist.o(ii).Nsingle = area(x)*channel_density
ii += 1
}
If only a subset of sections is supposed to have the stochastic mechanism, then the first step is to create a SectionList and append all the appropriate sections to it. Then instead of forall, which iterates over all sections in your model, substitute
forsec slname
where slname is the name of your SectionList.
Be sure to read the Programmer's Reference entries about forall, for (x,0), area(), List, and SectionList.