Example #1
0
def get_max_dpp(cell, dpAmp, stimAmp, x_dist_val,y):
	[t,i]=pdl.make_dpp(mydelay,dur1,dur2,mydt,mysimtime,dpAmp, stimAmp)
	sim = pdl.Simulation(cell,mydt,sim_time = mysimtime)
	sim.set_exstim([t,i],x_dist=x_dist_val,y_dist = y,rho = myrho)
	sim.go()
	_,som_rec,ax_rec = sim.get_recording()
	return som_rec,ax_rec
Example #2
0
	save_obj(outaxon,'outaxon-tt23.p')

#looking at waveforms directly dpp (unbalanced)
if TEST_TYPE == 24:
	dur1 = .2
	dur2 = .05
#	dur3 = .2
	amp = -100
	dppamp = -80
	#amp = 0
	#dppamp = 0
	y = 10

	stim_params = [mydelay, dur1, dur2, mysimtime, mydt, myrho]

	[t,i] = pdl.make_dpp(mydelay, dur1, dur2, dppamp, amp, mysimtime, mydt)
	sim = pdl.Simulation(cell,mydt,sim_time = mysimtime)
	sim.set_exstim([t,i],x_dist = x_soma, y_dist = y, rho = myrho)
	sim.go()
	sim.show(showAx=0)
	plt.figure()
	plt.plot(t,i)
	plt.show()

#	thr = ptl.get_thresh(cell, 'dppbal', init_amp, init_step, x_soma, y, stim_params, cnt_lim = cntlim, dpp_amp = dppamp)
#	print thr

#make 3d dpp plot (distance vs. dpp perc) UNBALANCED
if TEST_TYPE == 25:
	cntlim = 30
#	durlist = [.2, .3, .6, 1, 2, 3]
Example #3
0
mydt = .01
#mysimtime = 500
mysimtime = 8
myrho = 7900
mydelay = 1.5
#dur1 = .5
dur1 = 4
dur2 = .1
dpAmp = .1e3
stimAmp = .6e3

x_dist_val = 55


for x_dist_val in [0,15,55,250]:
	[t,i]=pdl.make_dpp(mydelay,dur1,dur2,mydt,mysimtime,dpAmp,stimAmp)
	sim = pdl.Simulation(cell,mydt,sim_time = mysimtime)
	sim.set_exstim([t,i],x_dist=x_dist_val,y_dist=3,rho = myrho) 
	sim.go()
	plt.figure()
	sim.show(showAx=False)

	plt.plot(t,i/max(i)*20)

	[t,i]=pdl.make_square(mydelay+dur1,stimAmp,dur2,mysimtime,mydt)
	sim = pdl.Simulation(cell,mydt,sim_time = mysimtime)
	sim.set_exstim([t,i],x_dist=x_dist_val,y_dist=3,rho = myrho) 
	sim.go()
	sim.show(showAx=False)
	plt.plot(t,i/max(i)*20)