def make_fF_Curve(self,
projFile,
simulator,
simConfig,
nTrains,
simDuration,
analyseStartTime, analyseStopTime,
analyseThreshold,
maxNumSimultaneousSims = 4):
# Load neuroConstruct project
print "Loading project from file: " + projFile.getAbsolutePath()+", exists: "+ str(projFile.exists())
pm = ProjectManager()
self.myProject = pm.loadProject(projFile)
self.simulator = simulator
simConfig = self.myProject.simConfigInfo.getSimConfig(simConfig)
simConfig.setSimDuration(simDuration)
pm.doGenerate(simConfig.getName(), self.neuroConstructSeed)
while pm.isGenerating():
print "Waiting for the project to be generated with Simulation Configuration: "+str(simConfig)
time.sleep(2)
numGenerated = self.myProject.generatedCellPositions.getNumberInAllCellGroups()
print "Number of cells generated: " + str(numGenerated)
simReferences = {}
if numGenerated > 0:
print "Generating scripts for simulator: %s..."%simulator
if simulator == 'NEURON':
self.myProject.neuronFileManager.setQuitAfterRun(1) # Remove this line to leave the NEURON sim windows open after finishing
self.myProject.neuronSettings.setCopySimFiles(1) # 1 copies hoc/mod files to PySim_0 etc. and will allow multiple sims to run at once
self.myProject.neuronSettings.setGraphicsMode(0) # 0 hides graphs during execution
if simulator == 'GENESIS':
self.myProject.genesisFileManager.setQuitAfterRun(1) # Remove this line to leave the NEURON sim windows open after finishing
self.myProject.genesisSettings.setCopySimFiles(1) # 1 copies hoc/mod files to PySim_0 etc. and will allow multiple sims to run at once
self.myProject.genesisSettings.setGraphicsMode(0) # 0 hides graphs during execution
currentTrain = 0
ADD_TO_START_FINITIALIZE = self.myProject.neuronSettings.getNativeBlock(NativeCodeLocation.START_FINITIALIZE)
ADD_TO_RECORD_I = self.myProject.neuronSettings.getNativeBlock(NativeCodeLocation.AFTER_SIMULATION)
while currentTrain <= nTrains:
while (len(self.simsRunning)>=maxNumSimultaneousSims):
print "Sims currently running: "+str(self.simsRunning)
print "Waiting..."
time.sleep(3) # wait a while...
self.updateSimsRunning()
simRef = "PySim_"+str(currentTrain)
print "Going to run simulation: "+simRef
######## Adjusting the amplitude of the Voltage clamp ###############
TEXT_BEFORE_CREATION = """objref vector, gTrainFile\n""" + "objectvar clamp\n"
### TEXT_BEFORE_INIT = "GranuleCell_mod_tonic[0].Soma {\n" + "clampobj = new VClamp(0.5)\n" + "clampobj.dur[0] = " + str(stimDur) + "\n" + "clampobj.amp[0] = " + str(stimAmp) + "\n" + "}\n" # This should do the trick
TEXT_START_FINITIALIZE = "\n" + "gTrainFile = new File() \n" + "vector = new Vector(100000) \n" + "access GranuleCell_mod_tonic[0].Soma \n" + "clamp = new SEClamp(0.5) \n" + "clamp.amp1 = 0 \n" + "clamp.dur1 = 1e9 \n" + """gTrainFile.ropen("E:/neuroConstruct/models/Dan_GranCell/gTrains_nS/gAMPA_""" + str(currentTrain) + """.txt") \n""" + "vector.scanf(gTrainFile) \n" + "gTrainFile.close \n" + "for i = 0, vector.size() -1 { \n" + " if (vector.x[i] <= 0) { \n" + " vector.x[i] = 1e-20 \n" + " } \n" + " } \n" + "for i = 0, vector.size() -1 { \n" + " vector.x[i] = ( 1 / vector.x[i] ) * 1000 \n" + " } \n" + "vector.play(&clamp.rs, 0.03) \n" # This
TEXT_START_FINITIALIZE = TEXT_START_FINITIALIZE + ADD_TO_START_FINITIALIZE
TEXT_TO_RECORD_I = """// currently not used \n"""
TEXT_TO_RECORD_I = TEXT_TO_RECORD_I + ADD_TO_RECORD_I
self.myProject.neuronSettings.setNativeBlock(NativeCodeLocation.BEFORE_CELL_CREATION, TEXT_BEFORE_CREATION)
self.myProject.neuronSettings.setNativeBlock(NativeCodeLocation.START_FINITIALIZE, TEXT_START_FINITIALIZE)
self.myProject.neuronSettings.setNativeBlock(NativeCodeLocation.AFTER_SIMULATION, TEXT_TO_RECORD_I)
print "Next Train: "+ str(currentTrain)
self.myProject.simulationParameters.setReference(simRef)
if simulator == "NEURON":
self.myProject.neuronFileManager.generateTheNeuronFiles(simConfig,
None,
NeuronFileManager.RUN_HOC,
self.simulatorSeed)
print "Generated NEURON files for: "+simRef
compileProcess = ProcessManager(self.myProject.neuronFileManager.getMainHocFile())
compileSuccess = compileProcess.compileFileWithNeuron(0,0)
print "Compiled NEURON files for: "+simRef
if compileSuccess:
pm.doRunNeuron(simConfig)
print "Set running simulation: "+simRef
self.simsRunning.append(simRef)
if simulator == "GENESIS":
compartmentalisation = GenesisCompartmentalisation()
self.myProject.genesisFileManager.generateTheGenesisFiles(simConfig,
None,
compartmentalisation,
self.simulatorSeed)
print "Generated GENESIS files for: "+simRef
pm.doRunGenesis(simConfig)
print "Set running simulation: "+simRef
self.simsRunning.append(simRef)
time.sleep(1) # Wait for sim to be kicked off
simReferences[simRef] = currentTrain
currentTrain = currentTrain + 1
print
print "Finished running "+str(len(simReferences))+" simulations for project "+ projFile.getAbsolutePath()
print "These can be loaded and replayed in the previous simulation browser in the GUI"
print
while (len(self.simsRunning)>0):
print "Sims currently running: "+str(self.simsRunning)
print "Waiting..."
time.sleep(4) # wait a while...
self.updateSimsRunning()
#simReferences = {'PySim_0.3':0.3,'PySim_0.4':0.4,'PySim_0.5':0.5}
plotFrameFf = PlotManager.getPlotterFrame("F-f curve from project: "+str(self.myProject.getProjectFile())+" on "+simulator , 1, 1)
plotFrameVolts = PlotManager.getPlotterFrame("VoltageTraces from project: "+str(self.myProject.getProjectFile())+" on "+simulator , 1, 1)
plotFrameFf.setViewMode(PlotCanvas.INCLUDE_ORIGIN_VIEW)
info = "F-f curve for Simulation Configuration: "+str(simConfig)
dataSet = DataSet(info, info, "Hz", "Hz", "Input_Freq", "Output_Freq")
dataSet.setGraphFormat(PlotCanvas.USE_CIRCLES_FOR_PLOT)
simList = simReferences.keys()
simList.sort()
traininfos = open("E:/neuroConstruct/models/Dan_GranCell/gAMPA_traininfo.txt")
train_info_list = traininfos.readlines()
currentTrain = 0
for sim in simList:
simDir = File(projFile.getParentFile(), "E:/neuroConstruct/models/Dan_GranCell/simulations/"+sim)
print
print "--- Reloading data from simulation in directory: %s"%simDir.getCanonicalPath()
try:
simData = SimulationData(simDir)
simData.initialise()
print "Data loaded: "
print simData.getAllLoadedDataStores()
times = simData.getAllTimes()
cellSegmentRef = simConfig.getCellGroups().get(0)+"_0"
volts = simData.getVoltageAtAllTimes(cellSegmentRef)
time.sleep(2)
traceInfo = "Voltage at: %s in simulation: %s"%(cellSegmentRef, sim)
dataSetV = DataSet(traceInfo, traceInfo, "ms", "mV", "Time", "Membrane Potential")
Spike_List = []
Spike_List.append(0.00)
for i in range(len(times)):
dataSetV.addPoint(times[i], volts[i])
for i in range(len(times)):
if (volts[i] > 0.0):
if (times[i] - Spike_List[len(Spike_List)-1] > 0.50):
Spike_List.append(times[i])
Spike_List.remove(0.00)
plotFrameVolts.addDataSet(dataSetV)
spikeTimes = SpikeAnalyser.getSpikeTimes(volts, times, analyseThreshold, analyseStartTime, analyseStopTime)
currentTrain_Freq = float(train_info_list[currentTrain-1])
print "Number of spikes in sim %s: %i"%(sim, len(spikeTimes))
avgFreq = 0
if len(Spike_List)>0:
avgFreq = 1000.0 / ((Spike_List[len(Spike_List)-1] - Spike_List[0]) / (len(Spike_List) -1))
dataSet.addPoint(currentTrain_Freq,avgFreq)
currentTrain = currentTrain + 1
except:
print "Error analysing simulation data from: %s"%simDir.getCanonicalPath()
print sys.exc_info()[0]
plotFrameFf.addDataSet(dataSet)
times_original = []
voltages_original = []
for line in open ('models/LarkumEtAl2009/ModelDB/IClamp_somaV.txt', 'rt'):
t, v = [float (x) for x in line.split()]
times_original.append (t)
voltages_original.append (v)
print "============================"
print "plot and compare spike times"
print "============================"
analyseStartTime = 0 # ms
analyseStopTime = 600
analyseThreshold = -20 # mV
spikeTimes_nC = SpikeAnalyser.getSpikeTimes(voltages_nC, times_nC, analyseThreshold, analyseStartTime, analyseStopTime)
print "neuroConcsturct spike times:"
print spikeTimes_nC
spikeTimes_original = SpikeAnalyser.getSpikeTimes(voltages_original, times_original, analyseThreshold, analyseStartTime, analyseStopTime)
print "modelDB spike times:"
print spikeTimes_original
plotFrame = PlotManager.getPlotterFrame("test model: "+str(myProject.getProjectFile()) , 1, 1)
plotFrame.setViewMode(PlotCanvas.INCLUDE_ORIGIN_VIEW)
info = "modelDB vs nC: "+str(simConfig)
dataSet_nC = DataSet(info, info, "ms", "nC", "time", "nC")
dataSet_nC.setGraphFormat(PlotCanvas.USE_CIRCLES_FOR_PLOT)
for t1 in spikeTimes_nC:
dataSet_nC.addPoint(t1,1)
cellSegmentRef = simData.getAllLoadedDataStores().get(
0).getCellGroupName() + "_0"
print "Looking for voltages of " + cellSegmentRef
volts = simData.getVoltageAtAllTimes(cellSegmentRef)
if verbose:
print "Got " + str(
len(volts)
) + " data points on cell seg ref: " + cellSegmentRef
analyseStartTime = 0
analyseStopTime = 2000
analyseThreshold = -20 # mV
spikeTimes = SpikeAnalyser.getSpikeTimes(
volts, times, analyseThreshold, analyseStartTime,
analyseStopTime)
#print "Spike times in %s for sim %s: %s"%(cellSegmentRef, simRef, str(spikeTimes))
if totNumSpikes < 0: totNumSpikes = len(spikeTimes)
for spikeIndex in range(max(totNumSpikes - numSpikes, 0),
totNumSpikes):
spikeTimeTrace = "Times_" + sim + "_" + simConfigName + "_" + str(
spikeIndex)
if not allSpikeTimeDataSets.has_key(spikeTimeTrace):
ds = DataSet(spikeTimeTrace, spikeTimeTrace, "",
"ms", "Number of compartments",
simData = SimulationData(simDir)
simData.initialise()
times = simData.getAllTimes()
cellSegmentRef = simData.getAllLoadedDataStores().get(0).getCellGroupName()+"_0"
print "Looking for voltages of "+cellSegmentRef
volts = simData.getVoltageAtAllTimes(cellSegmentRef)
if verbose: print "Got "+str(len(volts))+" data points on cell seg ref: "+cellSegmentRef
analyseStartTime = 0
analyseStopTime = 2000
analyseThreshold = -20 # mV
spikeTimes = SpikeAnalyser.getSpikeTimes(volts, times, analyseThreshold, analyseStartTime, analyseStopTime)
#print "Spike times in %s for sim %s: %s"%(cellSegmentRef, simRef, str(spikeTimes))
if totNumSpikes <0: totNumSpikes = len(spikeTimes)
for spikeIndex in range(max(totNumSpikes-numSpikes, 0), totNumSpikes):
spikeTimeTrace = "Times_"+sim+"_"+simConfigName+"_"+str(spikeIndex)
if not allSpikeTimeDataSets.has_key(spikeTimeTrace):
ds = DataSet(spikeTimeTrace, spikeTimeTrace, "", "ms", "Number of compartments", "Spike time")
allSpikeTimeDataSets[spikeTimeTrace] = ds
spikeTrace = allSpikeTimeDataSets[spikeTimeTrace]
def generateF_ICurve(self,
projFile,
simulator,
simConfig,
preStimAmp, preStimDel, preStimDur,
stimAmpLow, stimAmpInc, stimAmpHigh,
stimDel, stimDur,
simDuration,
analyseStartTime, analyseStopTime,
analyseThreshold,
maxNumSimultaneousSims = 4):
# Load neuroConstruct project
print "Loading project from file: " + projFile.getAbsolutePath()+", exists: "+ str(projFile.exists())
pm = ProjectManager()
self.myProject = pm.loadProject(projFile)
self.simulator = simulator
simConfig = self.myProject.simConfigInfo.getSimConfig(simConfig)
simConfig.setSimDuration(simDuration)
pm.doGenerate(simConfig.getName(), self.neuroConstructSeed)
while pm.isGenerating():
print "Waiting for the project to be generated with Simulation Configuration: "+str(simConfig)
time.sleep(2)
numGenerated = self.myProject.generatedCellPositions.getNumberInAllCellGroups()
print "Number of cells generated: " + str(numGenerated)
simReferences = {}
if numGenerated > 0:
print "Generating scripts for simulator: %s..."%simulator
if simulator == 'NEURON':
self.myProject.neuronFileManager.setQuitAfterRun(1) # Remove this line to leave the NEURON sim windows open after finishing
self.myProject.neuronSettings.setCopySimFiles(1) # 1 copies hoc/mod files to PySim_0 etc. and will allow multiple sims to run at once
self.myProject.neuronSettings.setGraphicsMode(0) # 0 hides graphs during execution
if simulator == 'GENESIS':
self.myProject.genesisFileManager.setQuitAfterRun(1) # Remove this line to leave the NEURON sim windows open after finishing
self.myProject.genesisSettings.setCopySimFiles(1) # 1 copies hoc/mod files to PySim_0 etc. and will allow multiple sims to run at once
self.myProject.genesisSettings.setGraphicsMode(0) # 0 hides graphs during execution
stimAmp = stimAmpLow
while stimAmp <= stimAmpHigh:
while (len(self.simsRunning)>=maxNumSimultaneousSims):
print "Sims currently running: "+str(self.simsRunning)
print "Waiting..."
time.sleep(3) # wait a while...
self.updateSimsRunning()
simRef = "PySim_"+str(float(stimAmp))
print "Going to run simulation: "+simRef
######## Adjusting the amplitude of the current clamp ###############
# preStim = self.myProject.elecInputInfo.getStim(simConfig.getInputs().get(0))
preStim = self.myProject.elecInputInfo.getStim("Input_3")
preStim.setAmp(NumberGenerator(preStimAmp))
preStim.setDel(NumberGenerator(preStimDel))
preStim.setDur(NumberGenerator(preStimDur))
self.myProject.elecInputInfo.updateStim(preStim)
# stim = self.myProject.elecInputInfo.getStim(simConfig.getInputs().get(1))
stim = self.myProject.elecInputInfo.getStim("Input_4")
stim.setAmp(NumberGenerator(stimAmp))
stim.setDel(NumberGenerator(stimDel))
stim.setDur(NumberGenerator(stimDur))
self.myProject.elecInputInfo.updateStim(stim)
print "Next stim: "+ str(stim)
self.myProject.simulationParameters.setReference(simRef)
if simulator == "NEURON":
self.myProject.neuronFileManager.generateTheNeuronFiles(simConfig,
None,
NeuronFileManager.RUN_HOC,
self.simulatorSeed)
print "Generated NEURON files for: "+simRef
compileProcess = ProcessManager(self.myProject.neuronFileManager.getMainHocFile())
compileSuccess = compileProcess.compileFileWithNeuron(0,0)
print "Compiled NEURON files for: "+simRef
if compileSuccess:
pm.doRunNeuron(simConfig)
print "Set running simulation: "+simRef
self.simsRunning.append(simRef)
if simulator == "GENESIS":
compartmentalisation = GenesisCompartmentalisation()
self.myProject.genesisFileManager.generateTheGenesisFiles(simConfig,
None,
compartmentalisation,
self.simulatorSeed)
print "Generated GENESIS files for: "+simRef
pm.doRunGenesis(simConfig)
print "Set running simulation: "+simRef
self.simsRunning.append(simRef)
time.sleep(1) # Wait for sim to be kicked off
simReferences[simRef] = stimAmp
stimAmp = stimAmp +stimAmpInc
print
print "Finished running "+str(len(simReferences))+" simulations for project "+ projFile.getAbsolutePath()
print "These can be loaded and replayed in the previous simulation browser in the GUI"
print
while (len(self.simsRunning)>0):
print "Sims currently running: "+str(self.simsRunning)
print "Waiting..."
time.sleep(4) # wait a while...
self.updateSimsRunning()
#simReferences = {'PySim_0.3':0.3,'PySim_0.4':0.4,'PySim_0.5':0.5}
plotFrameFI = PlotManager.getPlotterFrame("F-I curve from project: "+str(self.myProject.getProjectFile())+" on "+simulator , 1, 1)
plotFrameVolts = PlotManager.getPlotterFrame("VoltageTraces from project: "+str(self.myProject.getProjectFile())+" on "+simulator , 1, 1)
plotFrameFI.setViewMode(PlotCanvas.INCLUDE_ORIGIN_VIEW)
info = "F-I curve for Simulation Configuration: "+str(simConfig)
dataSet = DataSet(info, info, "nA", "Hz", "Current injected", "Firing frequency")
dataSet.setGraphFormat(PlotCanvas.USE_CIRCLES_FOR_PLOT)
simList = simReferences.keys()
simList.sort()
for sim in simList:
simDir = File(projFile.getParentFile(), "/neuroConstruct/models/Dan_GranCell/simulations/"+sim)
print
print "--- Reloading data from simulation in directory: %s"%simDir.getCanonicalPath()
try:
simData = SimulationData(simDir)
simData.initialise()
print "Data loaded: "
print simData.getAllLoadedDataStores()
times = simData.getAllTimes()
cellSegmentRef = simConfig.getCellGroups().get(0)+"_0"
volts = simData.getVoltageAtAllTimes(cellSegmentRef)
traceInfo = "Voltage at: %s in simulation: %s"%(cellSegmentRef, sim)
dataSetV = DataSet(traceInfo, traceInfo, "mV", "ms", "Membrane potential", "Time")
for i in range(len(times)):
dataSetV.addPoint(times[i], volts[i])
plotFrameVolts.addDataSet(dataSetV)
spikeTimes = SpikeAnalyser.getSpikeTimes(volts, times, analyseThreshold, analyseStartTime, analyseStopTime)
stimAmp = simReferences[sim]
print "Number of spikes at %f nA in sim %s: %i"%(stimAmp, sim, len(spikeTimes))
avgFreq = 0
if len(spikeTimes)>1:
avgFreq = len(spikeTimes)/ ((analyseStopTime - analyseStartTime)/1000.0)
dataSet.addPoint(stimAmp,avgFreq)
except:
print "Error analysing simulation data from: %s"%simDir.getCanonicalPath()
print sys.exc_info()[0]
plotFrameFI.addDataSet(dataSet)
cellsToUse.append(pos.cellNumber)
GeneralUtils.timeCheck("After load sim", True)
spikeLists = ArrayList()
times = simData.getAllTimes()
for cellNum in cellsToUse:
cellSegRef = SimulationData.getCellRef(cellGroupName, cellNum)
volts = simData.getVoltageAtAllTimes(cellSegRef)
spikes = SpikeAnalyser.getSpikeTimes(volts,
times,
threshold,
0,
times[-1])
spikeLists.add(spikes)
ds = SpikeAnalyser.getSlidingSpikeSynchrony(spikeLists,
times,
slideSize,
0,
times[-1])
ds.setDescription(str(cellsToUse))
peak_vals = []
peak_times = []
