bmtool 0.6.6.2__py3-none-any.whl → 0.6.6.4__py3-none-any.whl

bmtool/synapses.py

@@ -84,6 +84,20 @@ class SynapseTuner:
         h.dt = general_settings['dt']  # Time step (resolution) of the simulation in ms
         h.steps_per_ms = 1 / h.dt
         h.celsius = general_settings['celsius']
+        
+        # get some stuff set up we need for both SingleEvent and Interactive Tuner
+        self._set_up_cell()
+        self._set_up_synapse()
+        
+        self.nstim = h.NetStim()
+        self.nstim2 = h.NetStim()
+        
+        self.vcl = h.VClamp(self.cell.soma[0](0.5))
+        
+        self.nc = h.NetCon(self.nstim, self.syn, self.general_settings['threshold'], self.general_settings['delay'], self.general_settings['weight'])
+        self.nc2 = h.NetCon(self.nstim2, self.syn, self.general_settings['threshold'], self.general_settings['delay'], self.general_settings['weight'])
+        
+        self._set_up_recorders()
 
     def _update_spec_syn_param(self, json_folder_path):
         """
@@ -168,7 +182,7 @@ class SynapseTuner:
         self.ivcl.record(self.vcl._ref_i)
 
 
-    def SingleEvent(self):
+    def SingleEvent(self,plot_and_print=True):
         """
         Simulate a single synaptic event by delivering an input stimulus to the synapse.
         
@@ -176,33 +190,30 @@ class SynapseTuner:
         and then runs the NEURON simulation for a single event. The single synaptic event will occur at general_settings['tstart']
         Will display graphs and synaptic properies works best with a jupyter notebook
         """
-        self._set_up_cell()
-        self._set_up_synapse()
+        self.ispk = None
         
         # user slider values if the sliders are set up
         if hasattr(self, 'dynamic_sliders'):
             syn_props = {var: slider.value for var, slider in self.dynamic_sliders.items()}
             self._set_syn_prop(**syn_props)
+                  
+        # sets values based off optimizer 
+        if hasattr(self,'using_optimizer'):
+            for name, value in zip(self.param_names, self.params):
+                setattr(self.syn, name, value)
 
         # Set up the stimulus
-        self.nstim = h.NetStim()
         self.nstim.start = self.general_settings['tstart']
         self.nstim.noise = 0
-        self.nstim2 = h.NetStim()
         self.nstim2.start = h.tstop
         self.nstim2.noise = 0
-        self.nc = h.NetCon(self.nstim, self.syn, self.general_settings['threshold'], self.general_settings['delay'], self.general_settings['weight'])
-        self.nc2 = h.NetCon(self.nstim2, self.syn, self.general_settings['threshold'], self.general_settings['delay'], self.general_settings['weight'])
         
         # Set up voltage clamp
-        self.vcl = h.VClamp(self.cell.soma[0](0.5))
         vcldur = [[0, 0, 0], [self.general_settings['tstart'], h.tstop, 1e9]]
         for i in range(3):
             self.vcl.amp[i] = self.conn['spec_settings']['vclamp_amp']
             self.vcl.dur[i] = vcldur[1][i]
 
-        self._set_up_recorders()
-
         # Run simulation
         h.tstop = self.general_settings['tstart'] + self.general_settings['tdur']
         self.nstim.interval = self.general_settings['tdur']
@@ -211,15 +222,18 @@ class SynapseTuner:
         h.run()
         
         current = np.array(self.rec_vectors[self.current_name])
-        syn_prop = self._get_syn_prop(short=True)
-        current = (current - syn_prop['baseline']) * 1000  # Convert to pA
+        syn_props = self._get_syn_prop(rise_interval=self.general_settings['rise_interval'],dt=h.dt) 
+        current = (current - syn_props['baseline']) * 1000  # Convert to pA
         current_integral = np.trapz(current, dx=h.dt)  # pA·ms
         
-        self._plot_model([self.general_settings['tstart'] - 5, self.general_settings['tstart'] + self.general_settings['tdur']])
-        syn_props = self._get_syn_prop(rise_interval=self.general_settings['rise_interval'])     
-        for prop in syn_props.items():
-            print(prop)
-        print(f'Current Integral in pA: {current_integral:.2f}')
+        if plot_and_print:
+            self._plot_model([self.general_settings['tstart'] - 5, self.general_settings['tstart'] + self.general_settings['tdur']])    
+            for prop in syn_props.items():
+                print(prop)
+            print(f'Current Integral in pA*ms: {current_integral:.2f}')
+        
+        self.rise_time = syn_props['rise_time']
+        self.decay_time = syn_props['decay_time']
 
 
     def _find_first(self, x):
@@ -263,7 +277,7 @@ class SynapseTuner:
         if self.vclamp:
             isyn = self.ivcl
         else:
-            isyn = self.rec_vectors['i']
+            isyn = self.rec_vectors[self.current_name]
         isyn = np.asarray(isyn)
         tspk = np.asarray(self.tspk)
         if tspk.size:
@@ -336,7 +350,7 @@ class SynapseTuner:
         
         # Plot synaptic current (always included)
         current = self.rec_vectors[self.current_name]
-        syn_prop = self._get_syn_prop(short=True)
+        syn_prop = self._get_syn_prop(short=True,dt=h.dt)
         current = (current - syn_prop['baseline']) 
         current = current * 1000
         
@@ -433,13 +447,12 @@ class SynapseTuner:
         """
         isyn = np.array(self.rec_vectors[self.current_name].to_python())
         tspk = np.append(np.asarray(self.tspk), h.tstop)
-        syn_prop = self._get_syn_prop(short=True)
+        syn_prop = self._get_syn_prop(short=True,dt=h.dt)
         # print("syn_prp[sign] = " + str(syn_prop['sign']))
         isyn = (isyn - syn_prop['baseline']) 
         isyn *= syn_prop['sign']
         ispk = np.floor((tspk + self.general_settings['delay']) / h.dt).astype(int)
 
-        
         try:        
             amp = [isyn[ispk[i]:ispk[i + 1]].max() for i in range(ispk.size - 1)]
             # indexs of where the max of the synaptic current is at. This is then plotted     
@@ -687,6 +700,7 @@ class SynapseTuner:
         display(ui)
         update_ui()
         
+        
     def stp_frequency_response(self, freqs=[1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 15, 20, 35, 50, 100, 200], 
                                 delay=250, plot=True,log_plot=True):
         """
@@ -733,6 +747,7 @@ class SynapseTuner:
         
         return results
 
+
     def _plot_frequency_analysis(self, results,log_plot):
         """
         Plot the frequency-dependent synaptic properties.
@@ -928,15 +943,33 @@ class SynapseOptimizer:
         
     def _calculate_metrics(self) -> Dict[str, float]:
         """Calculate standard metrics from the current simulation using specified frequency"""
-        self.tuner._simulate_model(self.train_frequency, self.train_delay)
-        amp = self.tuner._response_amplitude()
-        ppr, induction, recovery = self.tuner._calc_ppr_induction_recovery(amp, print_math=False)
-        amp = self.tuner._find_max_amp(amp)
+        
+        # Set these to 0 for when we return the dict 
+        induction = 0
+        ppr = 0
+        recovery = 0
+        amp = 0
+        rise_time = 0
+        decay_time = 0
+        
+        if self.run_single_event:
+            self.tuner.SingleEvent(plot_and_print=False)
+            rise_time = self.tuner.rise_time
+            decay_time = self.tuner.decay_time
+        
+        if self.run_train_input:
+            self.tuner._simulate_model(self.train_frequency, self.train_delay)
+            amp = self.tuner._response_amplitude()
+            ppr, induction, recovery = self.tuner._calc_ppr_induction_recovery(amp, print_math=False)
+            amp = self.tuner._find_max_amp(amp)
+        
         return {
             'induction': float(induction),
             'ppr': float(ppr),
             'recovery': float(recovery),
-            'max_amplitude': float(amp)
+            'max_amplitude': float(amp),
+            'rise_time': float(rise_time),
+            'decay_time': float(decay_time)
         }
         
     def _default_cost_function(self, metrics: Dict[str, float], target_metrics: Dict[str, float]) -> float:
@@ -957,7 +990,13 @@ class SynapseOptimizer:
         # Set parameters
         for name, value in zip(param_names, params):
             setattr(self.tuner.syn, name, value)
-        
+            
+        # just do this and have the SingleEvent handle it     
+        if self.run_single_event:
+            self.tuner.using_optimizer = True
+            self.tuner.param_names = param_names
+            self.tuner.params = params
+                    
         # Calculate metrics and error
         metrics = self._calculate_metrics()
         error = float(cost_function(metrics, target_metrics))  # Ensure error is scalar
@@ -974,6 +1013,7 @@ class SynapseOptimizer:
     
     def optimize_parameters(self, target_metrics: Dict[str, float],
                             param_bounds: Dict[str, Tuple[float, float]],
+                            run_single_event:bool = False, run_train_input:bool = True,
                             train_frequency: float = 50,train_delay: float = 250,
                             cost_function: Optional[Callable] = None,
                             method: str = 'SLSQP',init_guess='random') -> SynapseOptimizationResult:
@@ -1005,6 +1045,8 @@ class SynapseOptimizer:
         self.optimization_history = []
         self.train_frequency = train_frequency
         self.train_delay = train_delay
+        self.run_single_event = run_single_event
+        self.run_train_input = run_train_input
         
         param_names = list(param_bounds.keys())
         bounds = [param_bounds[name] for name in param_names]
@@ -1107,10 +1149,14 @@ class SynapseOptimizer:
             print(f"{param}: {float(value):.3f}")
         
         # Plot final model response
-        self.tuner._plot_model([self.tuner.general_settings['tstart'] - self.tuner.nstim.interval / 3, self.tuner.tstop])
-        amp = self.tuner._response_amplitude()
-        self.tuner._calc_ppr_induction_recovery(amp)        
-        
+        if self.run_train_input:
+            self.tuner._plot_model([self.tuner.general_settings['tstart'] - self.tuner.nstim.interval / 3, self.tuner.tstop])
+            amp = self.tuner._response_amplitude()
+            self.tuner._calc_ppr_induction_recovery(amp)
+        if self.run_single_event:
+            self.tuner.ispk=None
+            self.tuner.SingleEvent(plot_and_print=True)
+                    
         
 # dataclass means just init the typehints as self.typehint. looks a bit cleaner
 @dataclass

{bmtool-0.6.6.2.dist-info → bmtool-0.6.6.4.dist-info}/METADATA

@@ -1,6 +1,6 @@
 Metadata-Version: 2.2
 Name: bmtool
-Version: 0.6.6.2
+Version: 0.6.6.4
 Summary: BMTool
 Home-page: https://github.com/cyneuro/bmtool
 Download-URL: 

{bmtool-0.6.6.2.dist-info → bmtool-0.6.6.4.dist-info}/RECORD

@@ -7,7 +7,7 @@ bmtool/graphs.py,sha256=K8BiughRUeXFVvAgo8UzrwpSClIVg7UfmIcvtEsEsk0,6020
 bmtool/manage.py,sha256=_lCU0qBQZ4jSxjzAJUd09JEetb--cud7KZgxQFbLGSY,657
 bmtool/plot_commands.py,sha256=Tqujyf0c0u8olhiHOMwgUSJXIIE1hgjv6otb25G9cA0,12298
 bmtool/singlecell.py,sha256=MQiLucsI6OBIjtcJra3Z9PTFQOE-Zn5ST-R9SmFvrbQ,27049
-bmtool/synapses.py,sha256=KyvDtoL7_mbXLr7TXaaeFTjqUM_ZXoVzSABQ69eTAHI,52876
+bmtool/synapses.py,sha256=gIkfLhKDG2dHHCVJJoKuQrFn_Qut843bfk_-s97wu6c,54553
 bmtool/debug/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
 bmtool/debug/commands.py,sha256=AwtcR7BUUheM0NxvU1Nu234zCdpobhJv5noX8x5K2vY,583
 bmtool/debug/debug.py,sha256=xqnkzLiH3s-tS26Y5lZZL62qR2evJdi46Gud-HzxEN4,207
@@ -16,9 +16,9 @@ bmtool/util/commands.py,sha256=zJF-fiLk0b8LyzHDfvewUyS7iumOxVnj33IkJDzux4M,64396
 bmtool/util/util.py,sha256=00vOAwTVIifCqouBoFoT0lBashl4fCalrk8fhg_Uq4c,56654
 bmtool/util/neuron/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
 bmtool/util/neuron/celltuner.py,sha256=xSRpRN6DhPFz4q5buq_W8UmsD7BbUrkzYBEbKVloYss,87194
-bmtool-0.6.6.2.dist-info/LICENSE,sha256=qrXg2jj6kz5d0EnN11hllcQt2fcWVNumx0xNbV05nyM,1068
-bmtool-0.6.6.2.dist-info/METADATA,sha256=ufAGAV4laD-Jo5B7L77lHrfMgb4_V7iv_6u1Lf3O1Qc,20226
-bmtool-0.6.6.2.dist-info/WHEEL,sha256=In9FTNxeP60KnTkGw7wk6mJPYd_dQSjEZmXdBdMCI-8,91
-bmtool-0.6.6.2.dist-info/entry_points.txt,sha256=0-BHZ6nUnh0twWw9SXNTiRmKjDnb1VO2DfG_-oprhAc,45
-bmtool-0.6.6.2.dist-info/top_level.txt,sha256=gpd2Sj-L9tWbuJEd5E8C8S8XkNm5yUE76klUYcM-eWM,7
-bmtool-0.6.6.2.dist-info/RECORD,,
+bmtool-0.6.6.4.dist-info/LICENSE,sha256=qrXg2jj6kz5d0EnN11hllcQt2fcWVNumx0xNbV05nyM,1068
+bmtool-0.6.6.4.dist-info/METADATA,sha256=loQf1_yqp2RYGGrghu71zSt4TnBtlmBpOSGXlm_21JY,20226
+bmtool-0.6.6.4.dist-info/WHEEL,sha256=In9FTNxeP60KnTkGw7wk6mJPYd_dQSjEZmXdBdMCI-8,91
+bmtool-0.6.6.4.dist-info/entry_points.txt,sha256=0-BHZ6nUnh0twWw9SXNTiRmKjDnb1VO2DfG_-oprhAc,45
+bmtool-0.6.6.4.dist-info/top_level.txt,sha256=gpd2Sj-L9tWbuJEd5E8C8S8XkNm5yUE76klUYcM-eWM,7
+bmtool-0.6.6.4.dist-info/RECORD,,