efish/in_vivo_1
1
0
Fork 0
Code Issues Pull Requests Releases Wiki Activity

[04/2020] add juxta settings

Browse Source
This commit is contained in:
efish
2020-04-30 14:54:13 +02:00
parent 3acc4fa942
commit a2a6d8312b
16 changed files with 1185 additions and 98 deletions
Download Patch File Download Diff File Expand all files Collapse all files

16
torus/code/multisensory_stimuli.py
View File

@@ -33,12 +33,12 @@ class MultichannelStimulus(object):
def create_contrast_stim(self):
t = np.arange(0, self._contrast_sweep_duration, self._dt)
sweep_up = np.linspace(0., 1.0, len(t))
am = np.sin((self._eod_freq + self._df) * 2 * np.pi * t)
dc = np.sin(2*np.pi*self._dcf*t)
am = np.sin((self._eod_freq + self._df) * 2 * np.pi * t) * 0.5
dc = np.sin(2*np.pi*self._dcf*t) * 0.5
am_sweep = am * sweep_up
dc_sweep = dc * sweep_up
self._contrast_stimulus = np.hstack((am_sweep, dc_sweep, am + dc_sweep, dc + am_sweep))
self._contrast_stimulus = np.hstack((am_sweep, dc_sweep, (am + dc_sweep), (dc + am_sweep)))
header = {" ": "Contrast sweeps for ampullary and tuberous pathways",
"sd": 1.0,
"deltat": "%.5f s" % self._dt,
@@ -55,7 +55,7 @@ class MultichannelStimulus(object):
am = np.sin((self._eod_freq + self._df) * 2 * np.pi * t)
dc = np.sin(2*np.pi*self._df*t + phase_sweep)
self._phase_stimulus = am + dc
self._phase_stimulus = (am + dc)/2.
header = {" ": "Phase sweep stimulating ampullary and tuberous pathways",
"sd": 1.0,
"deltat": "%.5f s" % self._dt,
@@ -68,9 +68,9 @@ class MultichannelStimulus(object):
t = np.arange(0, self._zap_duration, self._dt)
m = self._maxfreq/self._zap_duration
dc = np.sin(2*np.pi*m*t*t)
am = np.sin(2*np.pi*(m*t + self._eod_freq)*t)
self._zap_stimulus = np.hstack((dc, am, dc+am))
dc = np.sin(2*np.pi*m*t*t) * 0.5
am = np.sin(2*np.pi*(m*t + self._eod_freq)*t) * 0.5
self._zap_stimulus = np.hstack((dc, am, (dc + am)))
header = {" ": "Zap stimulus for ampullary and tuberous pathways",
"sd": 1.0,
"deltat": "%.5f s" % self._dt,
@@ -112,7 +112,7 @@ if __name__ == "__main__":
help=args.dt_help)
parser.add_argument('-mf', '--maxfreq', type=float, default=100.,
help=args.maxfreq_help)
parser.add_argument('outfile prefix', help="name of the output file may include the path")
parser.add_argument('outfile', help="prefix name of the output file may include the path")
args = parser.parse_args()