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
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]
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)