bec/ophyd_devices
1
0
Fork 0
Code Issues 9 Pull Requests 10 Actions Packages Projects Releases 126 Wiki Activity

Testing, fixing and blacking

Browse Source
This commit is contained in:
Mohacsi Istvan 2022-12-01 11:01:01 +01:00
parent 6b2011318d
commit 0bffecc201
1 changed files with 19 additions and 12 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

23
ophyd_devices/epics/devices/specMotors.py
View File

@ -179,10 +179,15 @@ class MonoTheta2(VirtualEpicsSignalRO):
) )
return theta2 return theta2
MONO_THETA2_OFFSETS_FILENAME = "/import/work/sls/spec/local/X12SA/macros/spec_data/mono_th2_offsets.txt"
MONO_THETA2_OFFSETS_FILENAME = (
"/import/work/sls/spec/local/X12SA/macros/spec_data/mono_th2_offsets.txt"
)
class EnergyKev(VirtualEpicsSignalRO): class EnergyKev(VirtualEpicsSignalRO):
"""Converts the pusher motor position to energy in keV""" """Converts the pusher motor position to energy in keV"""
_mono_add_offs = True _mono_add_offs = True
_mono_a3_enc_offs2 = -19.7072 _mono_a3_enc_offs2 = -19.7072
_mono_a2_pusher_offs2 = 5.93905 _mono_a2_pusher_offs2 = 5.93905
@ -199,15 +204,17 @@ class EnergyKev(VirtualEpicsSignalRO):
if self._th2_offsets is None: if self._th2_offsets is None:
return 0.0 return 0.0
max_offs = 0.2 max_offs = np.max(self._th2_offsets[:, 1])
if np.max(self._th2_offsets) > max_offs:
raise ValueError(f"\nThe empirical moth2 corrections are as high as {max_offs} deg\nThis is unreasonable and the corrections will not be used.\n\n***PLEASE INFORM BEAMLINE SCIENTISTS***\n") if max_offs > 0.2:
raise ValueError(
f"\nThe empirical moth2 corrections are as high as {max_offs} deg\nThis is unreasonable and the corrections will not be used.\n\n***PLEASE INFORM BEAMLINE SCIENTISTS***\n"
)
offs = np.interp(energy_keV, self._th2_offsets[:, 0], self._th2_offsets[:, 1]) offs = np.interp(energy_keV, self._th2_offsets[:, 0], self._th2_offsets[:, 1])
print(offs) print(offs)
return offs return offs
def calc(self, val): def calc(self, val):
_mono_sintheta2_to_Ekev = -self._mono_hce / self._mono_2d2 _mono_sintheta2_to_Ekev = -self._mono_hce / self._mono_2d2
asin_arg = (val - self._mono_a2_pusher_offs2) / self._mono_a1_lever_length2 asin_arg = (val - self._mono_a2_pusher_offs2) / self._mono_a1_lever_length2
@ -217,7 +224,7 @@ class EnergyKev(VirtualEpicsSignalRO):
E_keV = _mono_sintheta2_to_Ekev / sin(theta2_deg / 180.0 * 3.141592) E_keV = _mono_sintheta2_to_Ekev / sin(theta2_deg / 180.0 * 3.141592)
if self._mono_add_offs: if self._mono_add_offs:
theta2_deg -= _mono_get_th2_offs(E_keV) theta2_deg -= self._mono_get_th2_offs(E_keV)
E_keV = _mono_sintheta2_to_Ekev / sin(theta2_deg / 180.0 * 3.141592) E_keV = _mono_sintheta2_to_Ekev / sin(theta2_deg / 180.0 * 3.141592)
return E_keV return E_keV
@ -270,7 +277,7 @@ class Bpm4i(Device):
read_attrs=None, read_attrs=None,
configuration_attrs=None, configuration_attrs=None,
parent=None, parent=None,
**kwargs **kwargs,
): ):
super().__init__( super().__init__(
prefix, prefix,
@ -279,7 +286,7 @@ class Bpm4i(Device):
read_attrs=read_attrs, read_attrs=read_attrs,
configuration_attrs=configuration_attrs, configuration_attrs=configuration_attrs,
parent=parent, parent=parent,
**kwargs **kwargs,
) )
self.sum.name = self.name self.sum.name = self.name
# Ensure the scaler counts automatically # Ensure the scaler counts automatically