diffcalc/startup/i07EH1v.py at master · DiamondLightSource/diffcalc · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
from startup._common_imports import *  # @UnusedWildImport
from diffcalc.settings import NUNAME

if GDA:
    from __main__ import diff1delta, diff1theta, diff1gamma, diff1omega, dcm1energy # @UnresolvedImport

from diffcalc.hkl.you.geometry import YouGeometry, YouPosition


class FourCircleI07EH1v(YouGeometry):
    """For a diffractometer with angles:
          delta, gam, mu, phi
    """
    def __init__(self, beamline_axes_transform=None):
        YouGeometry.__init__(self, 'fourc', {'eta': 0, 'chi': 0}, beamline_axes_transform)

    def physical_angles_to_internal_position(self, physical_angle_tuple):
        delta_phys, gam_phys, mu_phys, phi_phys = physical_angle_tuple
        return YouPosition(mu_phys, delta_phys, gam_phys, 0, 0, phi_phys, 'DEG')

    def internal_position_to_physical_angles(self, internal_position):
        clone_position = internal_position.clone()
        clone_position.changeToDegrees()
        mu_phys, delta_phys, gam_phys, _, _, phi_phys = clone_position.totuple()
        return delta_phys, gam_phys, mu_phys, phi_phys


### Create dummy scannables ###
if GDA:
    ###map GDA scannable to diffcalc axis name###
    _fourc = ScannableGroup('_fourc', (diff1delta, diff1gamma, diff1theta, diff1omega))
    delta=_fourc.diff1delta
    gamma = _fourc.diff1gamma
    omega = _fourc.diff1omega
    chi = diff1chi
    theta = diff1theta
    en=dcm1energy
else:
    #Create dummy axes to run outside GDA in IPython#
    delta = Dummy('delta')
    gamma = Dummy('gamma')
    theta = Dummy('theta')
    omega = Dummy('omega')
    _fourc = ScannableGroup('_fourc', (delta, gamma, theta, omega))
    en = Dummy('en')

en.level = 3

### Configure and import diffcalc objects ###
ESMTGKeV = 1
settings.hardware = ScannableHardwareAdapter(_fourc, en, ESMTGKeV)
settings.geometry = FourCircleI07EH1v()  # @UndefinedVariable
settings.energy_scannable = en
settings.axes_scannable_group= _fourc
settings.energy_scannable_multiplier_to_get_KeV = ESMTGKeV
settings.include_reference = False

from diffcalc.gdasupport.you import *  # @UnusedWildImport

if GDA:
    omega = _fourc.diff1omega
    print "Running in GDA --- aliasing commands"
    alias_commands(globals())

# Load the last ub calculation used
lastub()
# Set reference vector direction returning betain and betaout angles as theta and beta
if ubcalc.name:
    surfnphi('0; 0; 1')

### Set i07 specific limits
def setLimitsAndCuts():
    ''' set motor limits for diffcalc, these are within the actual motor limits
    '''
    setmin(delta, -11.275)
    setmax(delta, 100.0)
    setmin(gamma, -11.0)
    setmax(gamma, 49.934)
    setmin(theta, -6.561)
    setmax(theta, 29.163)
    setcut(omega, -180.0)
    print "Current hardware limits set to:"
    hardware()

if not GDA:
    setLimitsAndCuts()

# TODO: make demo code for (2+2) diffractometer geometry
#if not GDA:
#    demo = startup._demo.Demo(globals(), 'fourc')