Skip to content
Snippets Groups Projects
Commit 37d83286 authored by Loïc Le Guyader's avatar Loïc Le Guyader
Browse files

Merge branch 'checkTimApdWindow-meanTrace' into 'master'

checkTimApdWindow plots the mean of the MCPraw traces if they are loaded instead…

See merge request !32
parents fd8ae729 b3537574
No related branches found
No related tags found
No related merge requests found
...@@ -836,14 +836,14 @@ def checkTimApdWindow(data, mcp=1, use_apd=True, intstart=None, intstop=None): ...@@ -836,14 +836,14 @@ def checkTimApdWindow(data, mcp=1, use_apd=True, intstart=None, intstop=None):
apdChannels={1:3, 2:2, 3:1, 4:0} apdChannels={1:3, 2:2, 3:1, 4:0}
npulses_max = data['npulses_sase3'].max().values npulses_max = data['npulses_sase3'].max().values
tid = data['npulses_sase3'].where(data['npulses_sase3'] == npulses_max, tid = data['npulses_sase3'].where(data['npulses_sase3'] == npulses_max,
drop=True)[0].trainId.values drop=True).trainId.values
if 'MCP{}raw'.format(mcp) not in data: if 'MCP{}raw'.format(mcp) not in data:
tid, data_from_train = data.attrs['run'].train_from_id(tid) ntid, data_from_train = data.attrs['run'].train_from_id(tid[0])
trace = data_from_train['SCS_UTC1_ADQ/ADC/1:network']['digitizers.channel_1_' trace = data_from_train['SCS_UTC1_ADQ/ADC/1:network']['digitizers.channel_1_'
+'{}.raw.samples'.format(mcpToChannel[mcp])] +'{}.raw.samples'.format(mcpToChannel[mcp])]
print('no raw data for MCP{}. Loading trace from MCP{}'.format(mcp, mcp)) print('no raw data for MCP{}. Loading trace from MCP{}'.format(mcp, mcp))
else: else:
trace = data['MCP{}raw'.format(mcp)].sel(trainId=tid).T trace = data['MCP{}raw'.format(mcp)].sel({'trainId':tid}).mean(dim='trainId')
if use_apd: if use_apd:
pulseStart = data.attrs['run'].get_array( pulseStart = data.attrs['run'].get_array(
'SCS_UTC1_ADQ/ADC/1', 'SCS_UTC1_ADQ/ADC/1',
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment