@@ -278,7 +278,7 @@ def __init__(self, filename, nsx_override=None, nev_override=None,
278
278
exts = ['nsx' , 'nev' , 'sif' , 'ccf' ]
279
279
ext_overrides = [nsx_override , nev_override , sif_override , ccf_override ]
280
280
for ext , ext_override in zip (exts , ext_overrides ):
281
- if ext_override is not None and self ._avail_files [ext ] == False :
281
+ if ext_override is not None and self ._avail_files [ext ] is False :
282
282
raise ValueError ('Specified {} file {} could not be '
283
283
'found.' .format (ext , ext_override ))
284
284
@@ -1108,8 +1108,8 @@ def __get_nsx_rec_times_variant_a(self, nsx_nb):
1108
1108
bytes_in_headers = self .__nsx_params [self .__nsx_spec [nsx_nb ]](
1109
1109
'bytes_in_headers' , nsx_nb )
1110
1110
nb_data_points = int (
1111
- (self .__get_file_size (filename ) - bytes_in_headers ) /
1112
- (2 * self .__nsx_basic_header [nsx_nb ]['channel_count' ]) - 1 )
1111
+ (self .__get_file_size (filename ) - bytes_in_headers )
1112
+ / (2 * self .__nsx_basic_header [nsx_nb ]['channel_count' ]) - 1 )
1113
1113
1114
1114
# add n_start
1115
1115
n_starts = [(0 * t_unit ).rescale (highest_res )]
@@ -1267,8 +1267,8 @@ def __get_nsx_param_variant_a(self, param_name, nsx_nb):
1267
1267
nsx_parameters = {
1268
1268
'labels' : labels ,
1269
1269
'units' : np .array (
1270
- [b'uV' ] *
1271
- self .__nsx_basic_header [nsx_nb ]['channel_count' ]),
1270
+ [b'uV' ]
1271
+ * self .__nsx_basic_header [nsx_nb ]['channel_count' ]),
1272
1272
'min_analog_val' : - 1 * np .array (dig_factor ),
1273
1273
'max_analog_val' : np .array (dig_factor ),
1274
1274
'min_digital_val' : np .array (
@@ -1277,9 +1277,9 @@ def __get_nsx_param_variant_a(self, param_name, nsx_nb):
1277
1277
[1000 ] * self .__nsx_basic_header [nsx_nb ]['channel_count' ]),
1278
1278
'timestamp_resolution' : 30000 ,
1279
1279
'bytes_in_headers' :
1280
- self .__nsx_basic_header [nsx_nb ].dtype .itemsize +
1281
- self .__nsx_ext_header [nsx_nb ].dtype .itemsize *
1282
- self .__nsx_basic_header [nsx_nb ]['channel_count' ],
1280
+ self .__nsx_basic_header [nsx_nb ].dtype .itemsize
1281
+ + self .__nsx_ext_header [nsx_nb ].dtype .itemsize
1282
+ * self .__nsx_basic_header [nsx_nb ]['channel_count' ],
1283
1283
'sampling_rate' :
1284
1284
30000 / self .__nsx_basic_header [nsx_nb ]['period' ] * pq .Hz ,
1285
1285
'time_unit' : pq .CompoundUnit ("1.0/{0}*s" .format (
@@ -1310,11 +1310,11 @@ def __get_nsx_param_variant_b(self, param_name, nsx_nb):
1310
1310
'bytes_in_headers' :
1311
1311
self .__nsx_basic_header [nsx_nb ]['bytes_in_headers' ],
1312
1312
'sampling_rate' :
1313
- self .__nsx_basic_header [nsx_nb ]['timestamp_resolution' ] /
1314
- self .__nsx_basic_header [nsx_nb ]['period' ] * pq .Hz ,
1313
+ self .__nsx_basic_header [nsx_nb ]['timestamp_resolution' ]
1314
+ / self .__nsx_basic_header [nsx_nb ]['period' ] * pq .Hz ,
1315
1315
'time_unit' : pq .CompoundUnit ("1.0/{0}*s" .format (
1316
- self .__nsx_basic_header [nsx_nb ]['timestamp_resolution' ] /
1317
- self .__nsx_basic_header [nsx_nb ]['period' ]))}
1316
+ self .__nsx_basic_header [nsx_nb ]['timestamp_resolution' ]
1317
+ / self .__nsx_basic_header [nsx_nb ]['period' ]))}
1318
1318
1319
1319
return nsx_parameters [param_name ]
1320
1320
@@ -1336,8 +1336,8 @@ def __get_nsx_databl_param_variant_a(
1336
1336
# extract parameters from nsx basic extended and data header
1337
1337
data_parameters = {
1338
1338
'nb_data_points' : int (
1339
- (self .__get_file_size (filename ) - bytes_in_headers ) /
1340
- (2 * self .__nsx_basic_header [nsx_nb ]['channel_count' ]) - 1 ),
1339
+ (self .__get_file_size (filename ) - bytes_in_headers )
1340
+ / (2 * self .__nsx_basic_header [nsx_nb ]['channel_count' ]) - 1 ),
1341
1341
'databl_idx' : 1 ,
1342
1342
'databl_t_start' : t_starts [0 ],
1343
1343
'databl_t_stop' : t_stops [0 ]}
@@ -1410,8 +1410,8 @@ def __get_nonneural_evtypes_variant_a(self, data):
1410
1410
'name' : 'serial_input_port' ,
1411
1411
'field' : 'digital_input' ,
1412
1412
'mask' :
1413
- self .__is_set (data ['packet_insertion_reason' ], 0 ) &
1414
- self .__is_set (data ['packet_insertion_reason' ], 7 ),
1413
+ self .__is_set (data ['packet_insertion_reason' ], 0 )
1414
+ & self .__is_set (data ['packet_insertion_reason' ], 7 ),
1415
1415
'desc' : "Events of the serial input port" }}
1416
1416
1417
1417
# analog input events via threshold crossings
@@ -1451,8 +1451,8 @@ def __get_nonneural_evtypes_variant_b(self, data):
1451
1451
'name' : 'serial_input_port' ,
1452
1452
'field' : 'digital_input' ,
1453
1453
'mask' :
1454
- self .__is_set (data ['packet_insertion_reason' ], 0 ) &
1455
- self .__is_set (data ['packet_insertion_reason' ], 7 ),
1454
+ self .__is_set (data ['packet_insertion_reason' ], 0 )
1455
+ & self .__is_set (data ['packet_insertion_reason' ], 7 ),
1456
1456
'desc' : "Events of the serial input port" }}
1457
1457
1458
1458
return event_types
@@ -1658,7 +1658,7 @@ def __merge_time_ranges(
1658
1658
"" .format (user_n_stops [start_stop_id ]))
1659
1659
1660
1660
if (user_n_stops [start_stop_id ] < min_time
1661
- or user_n_starts [start_stop_id ] > max_time ):
1661
+ or user_n_starts [start_stop_id ] > max_time ):
1662
1662
user_n_stops .pop (start_stop_id )
1663
1663
user_n_starts .pop (start_stop_id )
1664
1664
self ._print_verbose (
@@ -1721,7 +1721,7 @@ def __merge_time_ranges(
1721
1721
mask = (n_stops_files > start ) & (n_stops_files < stop )
1722
1722
red_n_stops_files = np .array (n_stops_files )[mask ]
1723
1723
merged_n_stops = [
1724
- tp * highest_res for tp in red_n_stops_files ] + [stop ]
1724
+ tp * highest_res for tp in red_n_stops_files ] + [stop ]
1725
1725
else :
1726
1726
merged_n_stops = [stop ]
1727
1727
# Define combined user and file n_starts and n_stops
@@ -1747,7 +1747,7 @@ def __merge_time_ranges(
1747
1747
n_stops .extend (merged_n_stops )
1748
1748
1749
1749
if len (n_starts ) > len (user_n_starts ) and \
1750
- len (n_stops ) > len (user_n_stops ):
1750
+ len (n_stops ) > len (user_n_stops ):
1751
1751
self ._print_verbose (
1752
1752
"Additional recording pauses were detected. There will be "
1753
1753
"more segments than the user expects." )
@@ -1826,9 +1826,9 @@ def __read_spiketrain(
1826
1826
1827
1827
if scaling == 'voltage' :
1828
1828
st .waveforms = (
1829
- waveforms [mask ] * self .__nev_params ('waveform_unit' ) *
1830
- self .__nev_params ('digitization_factor' )[channel_id ] /
1831
- 1000. )
1829
+ waveforms [mask ] * self .__nev_params ('waveform_unit' )
1830
+ * self .__nev_params ('digitization_factor' )[channel_id ]
1831
+ / 1000. )
1832
1832
elif scaling == 'raw' :
1833
1833
st .waveforms = waveforms [mask ] * pq .dimensionless
1834
1834
else :
@@ -2008,11 +2008,11 @@ def __read_channelindex(
2008
2008
nev_dig_factor = self .__nev_ext_header [
2009
2009
b'NEUEVWAV' ]['digitization_factor' ][get_idx ],
2010
2010
nev_energy_threshold = self .__nev_ext_header [
2011
- b'NEUEVWAV' ]['energy_threshold' ][get_idx ] * pq .uV ,
2011
+ b'NEUEVWAV' ]['energy_threshold' ][get_idx ] * pq .uV ,
2012
2012
nev_hi_threshold = self .__nev_ext_header [
2013
- b'NEUEVWAV' ]['hi_threshold' ][get_idx ] * pq .uV ,
2013
+ b'NEUEVWAV' ]['hi_threshold' ][get_idx ] * pq .uV ,
2014
2014
nev_lo_threshold = self .__nev_ext_header [
2015
- b'NEUEVWAV' ]['lo_threshold' ][get_idx ] * pq .uV ,
2015
+ b'NEUEVWAV' ]['lo_threshold' ][get_idx ] * pq .uV ,
2016
2016
nb_sorted_units = self .__nev_ext_header [
2017
2017
b'NEUEVWAV' ]['nb_sorted_units' ][get_idx ],
2018
2018
waveform_size = self .__waveform_size [self .__nev_spec ](
@@ -2027,15 +2027,15 @@ def __read_channelindex(
2027
2027
# filter type codes (extracted from blackrock manual)
2028
2028
chidx .annotate (
2029
2029
nev_hi_freq_corner = self .__nev_ext_header [b'NEUEVFLT' ][
2030
- 'hi_freq_corner' ][get_idx ] /
2031
- 1000. * pq .Hz ,
2030
+ 'hi_freq_corner' ][get_idx ]
2031
+ / 1000. * pq .Hz ,
2032
2032
nev_hi_freq_order = self .__nev_ext_header [b'NEUEVFLT' ][
2033
2033
'hi_freq_order' ][get_idx ],
2034
2034
nev_hi_freq_type = flt_type [self .__nev_ext_header [
2035
2035
b'NEUEVFLT' ]['hi_freq_type' ][get_idx ]],
2036
2036
nev_lo_freq_corner = self .__nev_ext_header [
2037
- b'NEUEVFLT' ]['lo_freq_corner' ][get_idx ] /
2038
- 1000. * pq .Hz ,
2037
+ b'NEUEVFLT' ]['lo_freq_corner' ][get_idx ]
2038
+ / 1000. * pq .Hz ,
2039
2039
nev_lo_freq_order = self .__nev_ext_header [
2040
2040
b'NEUEVFLT' ]['lo_freq_order' ][get_idx ],
2041
2041
nev_lo_freq_type = flt_type [self .__nev_ext_header [
@@ -2051,9 +2051,9 @@ def __read_channelindex(
2051
2051
channel_id )
2052
2052
chidx .annotate (
2053
2053
nsx_hi_freq_corner = self .__nsx_ext_header [k ][
2054
- 'hi_freq_corner' ][get_idx ] / 1000. * pq .Hz ,
2054
+ 'hi_freq_corner' ][get_idx ] / 1000. * pq .Hz ,
2055
2055
nsx_lo_freq_corner = self .__nsx_ext_header [k ][
2056
- 'lo_freq_corner' ][get_idx ] / 1000. * pq .Hz ,
2056
+ 'lo_freq_corner' ][get_idx ] / 1000. * pq .Hz ,
2057
2057
nsx_hi_freq_order = self .__nsx_ext_header [k ][
2058
2058
'hi_freq_order' ][get_idx ],
2059
2059
nsx_lo_freq_order = self .__nsx_ext_header [k ][
@@ -2558,8 +2558,8 @@ def __str__(self):
2558
2558
'sampling_rate' , nsx_nb )
2559
2559
avail_el = [
2560
2560
el for el in self .__nsx_ext_header [nsx_nb ]['electrode_id' ]]
2561
- output += "\n Analog Parameters (NS" + str ( nsx_nb ) + ")" \
2562
- " \n ===================================="
2561
+ output += "\n Analog Parameters (NS" \
2562
+ + str ( nsx_nb ) + ") \n ===================================="
2563
2563
output += "\n Resolution (Hz): %i" % analog_res
2564
2564
output += "\n Available channel IDs: " + \
2565
2565
", " .join (["%i" % a for a in avail_el ]) + "\n "
0 commit comments