def returnTranslation(self, logLines):
_vecTmp = []
for line in logLines:
_vecTmp.append(float(line.split()[-1]))
return XSDataVectorDouble(_vecTmp[0], _vecTmp[1], _vecTmp[2])
def generateXSDataInputXDS(_xsDataCollection):
xsDataCollection = _xsDataCollection
xsDataExperimentalCondition = _xsDataCollection.getSubWedge()[0].getExperimentalCondition()
xsDataSubWedgeList = xsDataCollection.getSubWedge()
xsDataInputXDS = XSDataInputXDSGenerateBackgroundImage()
xsDataBeam = xsDataExperimentalCondition.getBeam()
xsDataDetector = xsDataExperimentalCondition.getDetector()
xsDataGoniostat = xsDataExperimentalCondition.getGoniostat()
dWavelength = xsDataBeam.getWavelength().getValue()
dDistance = xsDataDetector.getDistance().getValue()
dBeamPositionX = xsDataDetector.getBeamPositionX().getValue()
dBeamPositionY = xsDataDetector.getBeamPositionY().getValue()
# Start with the detector
xsDataXDSDetector = EDHandlerXSDataXDSv1_0.getXSDataXDSDetector(xsDataDetector)
xsDataInputXDS.setDetector(xsDataXDSDetector)
# Then the beam
xsDataXDSBeam = XSDataXDSBeam()
xsDataVectorDoubleIncidentBeam = XSDataVectorDouble()
xsDataVectorDoubleIncidentBeam.setV1(0.0)
xsDataVectorDoubleIncidentBeam.setV2(0.0)
xsDataVectorDoubleIncidentBeam.setV3(1.0)
xsDataXDSBeam.setIncident_beam_direction(xsDataVectorDoubleIncidentBeam)
xsDataVectorDoublePolarizationPlaneNormal = XSDataVectorDouble()
xsDataVectorDoublePolarizationPlaneNormal.setV1(0.0)
xsDataVectorDoublePolarizationPlaneNormal.setV2(1.0)
xsDataVectorDoublePolarizationPlaneNormal.setV3(0.0)
xsDataXDSBeam.setPolarization_plane_normal(xsDataVectorDoublePolarizationPlaneNormal)
xsDataXDSBeam.setX_ray_wavelength(XSDataWavelength(dWavelength))
xsDataInputXDS.setBeam(xsDataXDSBeam)
# Then the goniostat
xsDataXDSGoniostat = XSDataXDSGoniostat()
xsDataVectorDoubleRotationAxis = XSDataVectorDouble()
xsDataVectorDoubleRotationAxis.setV1(1.0)
xsDataVectorDoubleRotationAxis.setV2(0.0)
xsDataVectorDoubleRotationAxis.setV3(0.0)
xsDataXDSGoniostat.setRotation_axis(xsDataVectorDoubleRotationAxis)
xsDataXDSGoniostat.setOscillation_range(xsDataGoniostat.getOscillationWidth())
xsDataXDSGoniostat.setStarting_angle(xsDataGoniostat.getRotationAxisStart())
xsDataInputXDS.setGoniostat(xsDataXDSGoniostat)
# # Then the Crystal
#
# xsDataXDSCrystal = XSDataXDSCrystal()
#
# xsDataXDSCrystal.setFriedels_law(XSDataString("FALSE"))
#
# # if ( xsDataCrystal is not None ):
# # xsDataSpaceGroup = xsDataCrystal.getSpaceGroup()
# # if ( xsDataSpaceGroup is not None ):
# # xsDataStringName = xsDataSpaceGroup.getName()
# # if ( xsDataStringName is not None ):
# # xsDataInputXDS.setSymmetry( XSDataString( xsDataStringName.getValue() ) )
# xsDataXDSCrystal.setSpace_group_number(XSDataInteger(0))
#
# xsDataXDSCrystal.setStrong_pixel(XSDataInteger(8))
#
# xsDataCell = XSDataCell()
# xsDataCell.setLength_a(XSDataLength(0.0))
# xsDataCell.setLength_b(XSDataLength(0.0))
# xsDataCell.setLength_c(XSDataLength(0.0))
# xsDataCell.setAngle_alpha(XSDataAngle(0.0))
# xsDataCell.setAngle_beta(XSDataAngle(0.0))
# xsDataCell.setAngle_gamma(XSDataAngle(0.0))
# xsDataXDSCrystal.setUnit_cell_constants(xsDataCell)
#
# xsDataInputXDS.setCrystal(xsDataXDSCrystal)
# Finaly the images
xsDataXDSImage = XSDataXDSImage()
xsDataSubWedgeFirst = xsDataSubWedgeList[0]
xsDataImageFirst = xsDataSubWedgeFirst.getImage()[0]
pyStrPath = xsDataImageFirst.getPath().getValue()
pyStrFileName = EDUtilsFile.getBaseName(pyStrPath)
pyStrDirectory = EDUtilsPath.getFolderName(pyStrPath)
pyStrPrefix = EDUtilsImage.getPrefix(pyStrFileName)
pyStrSuffix = EDUtilsImage.getSuffix(pyStrFileName)
pyStrXDSTemplate = "%s_xdslink_?????.%s" % (pyStrPrefix, pyStrSuffix)
xsDataXDSImage.setName_template_of_data_frames(XSDataString(pyStrXDSTemplate))
iXDSLowestImageNumberGlobal = 1
xsDataXDSImage.setStarting_frame(XSDataInteger(iXDSLowestImageNumberGlobal))
# First we have to find the smallest goniostat rotation axis start:
fGonioStatOscillationStartMin = None
for xsDataSubWedge in xsDataSubWedgeList:
xsDataGoniostat = xsDataSubWedge.getExperimentalCondition().getGoniostat()
fGonioStatOscillationStart = xsDataGoniostat.getRotationAxisStart().getValue()
if (fGonioStatOscillationStartMin is None):
fGonioStatOscillationStartMin = fGonioStatOscillationStart
elif (fGonioStatOscillationStartMin > fGonioStatOscillationStart):
fGonioStatOscillationStartMin = fGonioStatOscillationStart
# Loop through the list of sub wedges
for xsDataSubWedge in xsDataSubWedgeList:
xsDataImageList = xsDataSubWedge.getImage()
xsDataGoniostat = xsDataSubWedge.getExperimentalCondition().getGoniostat()
fGonioStatOscillationStart = xsDataGoniostat.getRotationAxisStart().getValue()
fGonioStatOscillationRange = xsDataGoniostat.getOscillationWidth().getValue()
# First find the lowest and highest image numbers
iLowestImageNumber = None
for xsDataImage in xsDataImageList:
iImageNumber = xsDataImage.getNumber().getValue()
if (iLowestImageNumber is None):
iLowestImageNumber = iImageNumber
elif (iImageNumber < iLowestImageNumber):
iLowestImageNumber = iImageNumber
# Loop through the list of images
iLowestXDSImageNumber = None
iHighestXDSImageNumber = None
for xsDataImage in xsDataImageList:
iImageNumber = xsDataImage.getNumber().getValue()
fImageOscillationStart = fGonioStatOscillationStart + (iImageNumber - iLowestImageNumber) * fGonioStatOscillationRange
iXDSImageNumber = iXDSLowestImageNumberGlobal + int((fImageOscillationStart - fGonioStatOscillationStartMin) / fGonioStatOscillationRange)
# print iXDSImageNumber, fImageOscillationStart, fGonioStatOscillationStartMin, fGonioStatOscillationRange
pyStrSourcePath = xsDataImage.getPath()
pyStrTarget = "%s_xdslink_%05d.%s" % (pyStrPrefix, iXDSImageNumber, pyStrSuffix)
xsDataXDSImageLink = XSDataXDSImageLink()
xsDataFileSource = XSDataFile()
xsDataFileSource.setPath(pyStrSourcePath)
xsDataXDSImageLink.setSource(xsDataFileSource)
xsDataXDSImageLink.setTarget(XSDataString(pyStrTarget))
xsDataInputXDS.addImage_link(xsDataXDSImageLink)
if (iLowestXDSImageNumber is None):
iLowestXDSImageNumber = iXDSImageNumber
elif (iLowestXDSImageNumber > iXDSImageNumber):
iLowestXDSImageNumber = iXDSImageNumber
if (iHighestXDSImageNumber is None):
iHighestXDSImageNumber = iXDSImageNumber
elif (iHighestXDSImageNumber < iXDSImageNumber):
iHighestXDSImageNumber = iXDSImageNumber
xsDataXDSIntegerRange = XSDataXDSIntegerRange()
xsDataXDSIntegerRange.setLower(XSDataInteger(iLowestXDSImageNumber))
xsDataXDSIntegerRange.setUpper(XSDataInteger(iHighestXDSImageNumber))
xsDataXDSImage.addBackground_range(xsDataXDSIntegerRange)
xsDataXDSImage.addData_range(xsDataXDSIntegerRange)
xsDataXDSImage.addSpot_range(xsDataXDSIntegerRange)
xsDataInputXDS.setImage(xsDataXDSImage)
return xsDataInputXDS
def getXSDataXDSDetector(_xsDataDetector):
EDFactoryPluginStatic.loadModule("XSDataXDSv1_0")
from XSDataXDSv1_0 import XSDataXDSDetector
from XSDataXDSv1_0 import XSDataXDSIntegerRange
xsDataXDSDetector = XSDataXDSDetector()
strDetectorType = _xsDataDetector.getType().getValue()
if ((strDetectorType == "q4") or \
(strDetectorType == "q4-2x") or \
(strDetectorType == "q210") or \
(strDetectorType == "q210-2x") or \
(strDetectorType == "q315") or \
(strDetectorType == "q315-2x")):
xsDataXDSDetector.setDetector_name(XSDataString("ADSC"))
xsDataXDSDetector.setMinimum_valid_pixel_value(XSDataInteger(1))
xsDataXDSDetector.setOverload(XSDataInteger(65535))
xsDataXDSIntegerRangeTrustedPixel = XSDataXDSIntegerRange()
xsDataXDSIntegerRangeTrustedPixel.setLower(XSDataInteger(6000))
xsDataXDSIntegerRangeTrustedPixel.setUpper(XSDataInteger(30000))
xsDataXDSDetector.setValue_range_for_trusted_detector_pixels(xsDataXDSIntegerRangeTrustedPixel)
# elif ( strDetectorType == "mar165") or \
# strDetectorType == "mar225") ):
# xsDataXDSDetector.setType( XSDataString( "MARCCD" ) )
elif strDetectorType == "pilatus6m":
xsDataXDSDetector.setDetector_name(XSDataString("PILATUS"))
listUntrustedRectangle = \
[[ 487, 495, 0, 2528], \
[ 981, 989, 0, 2528], \
[1475, 1483, 0, 2528], \
[1969, 1977, 0, 2528], \
[ 0, 2464, 195, 213], \
[ 0, 2464, 407, 425], \
[ 0, 2464, 619, 637], \
[ 0, 2464, 831, 849], \
[ 0, 2464, 1043, 1061], \
[ 0, 2464, 1255, 1273], \
[ 0, 2464, 1467, 1485], \
[ 0, 2464, 1679, 1697], \
[ 0, 2464, 1891, 1909], \
[ 0, 2464, 2103, 2121], \
[ 0, 2464, 2315, 2333]]
for listRectangle in listUntrustedRectangle:
xsDataXDSRectangle = XSDataXDSRectangle()
xsDataXDSRectangle.setX1(XSDataInteger(listRectangle[0]))
xsDataXDSRectangle.setX2(XSDataInteger(listRectangle[1]))
xsDataXDSRectangle.setY1(XSDataInteger(listRectangle[2]))
xsDataXDSRectangle.setY2(XSDataInteger(listRectangle[3]))
xsDataXDSDetector.addUntrusted_rectangle(xsDataXDSRectangle)
xsDataXDSDetector.setMinimum_valid_pixel_value(XSDataInteger(0))
xsDataXDSDetector.setOverload(XSDataInteger(1048500))
xsDataXDSIntegerRangeTrustedPixel = XSDataXDSIntegerRange()
xsDataXDSIntegerRangeTrustedPixel.setLower(XSDataInteger(7000))
xsDataXDSIntegerRangeTrustedPixel.setUpper(XSDataInteger(30000))
xsDataXDSDetector.setValue_range_for_trusted_detector_pixels(xsDataXDSIntegerRangeTrustedPixel)
xsDataXDSDoubleRangeTrustedRegion = XSDataXDSDoubleRange()
xsDataXDSDoubleRangeTrustedRegion.setLower(XSDataDouble(0.0))
xsDataXDSDoubleRangeTrustedRegion.setUpper(XSDataDouble(1.41))
xsDataXDSDetector.setTrusted_region(xsDataXDSDoubleRangeTrustedRegion)
xsDataXDSDetector.setSensor_thickness(XSDataDouble(0.32))
elif strDetectorType == "pilatus2m":
xsDataXDSDetector.setDetector_name(XSDataString("PILATUS"))
listUntrustedRectangle = \
[[ 487, 495, 0, 1680], \
[ 981, 989, 0, 1680], \
[ 0, 1476, 195, 213], \
[ 0, 1476, 407, 425], \
[ 0, 1476, 619, 637], \
[ 0, 1476, 831, 849], \
[ 0, 1476, 1043, 1061], \
[ 0, 1476, 1255, 1273], \
[ 0, 1476, 1467, 1485]]
for listRectangle in listUntrustedRectangle:
xsDataXDSRectangle = XSDataXDSRectangle()
xsDataXDSRectangle.setX1(XSDataInteger(listRectangle[0]))
xsDataXDSRectangle.setX2(XSDataInteger(listRectangle[1]))
xsDataXDSRectangle.setY1(XSDataInteger(listRectangle[2]))
xsDataXDSRectangle.setY2(XSDataInteger(listRectangle[3]))
xsDataXDSDetector.addUntrusted_rectangle(xsDataXDSRectangle)
xsDataXDSDetector.setMinimum_valid_pixel_value(XSDataInteger(0))
xsDataXDSDetector.setOverload(XSDataInteger(1048500))
xsDataXDSIntegerRangeTrustedPixel = XSDataXDSIntegerRange()
xsDataXDSIntegerRangeTrustedPixel.setLower(XSDataInteger(7000))
xsDataXDSIntegerRangeTrustedPixel.setUpper(XSDataInteger(30000))
xsDataXDSDetector.setValue_range_for_trusted_detector_pixels(xsDataXDSIntegerRangeTrustedPixel)
xsDataXDSDoubleRangeTrustedRegion = XSDataXDSDoubleRange()
xsDataXDSDoubleRangeTrustedRegion.setLower(XSDataDouble(0.0))
xsDataXDSDoubleRangeTrustedRegion.setUpper(XSDataDouble(1.41))
xsDataXDSDetector.setTrusted_region(xsDataXDSDoubleRangeTrustedRegion)
xsDataXDSDetector.setSensor_thickness(XSDataDouble(0.32))
elif strDetectorType == "eiger4m":
xsDataXDSDetector.setDetector_name(XSDataString("PILATUS"))
listUntrustedRectangle = \
[[ 1029, 1040, 0, 2167], \
[ 0, 2070, 512, 550], \
[ 0, 2070, 1063, 1103], \
[ 0, 2070, 1614, 1654],
]
for listRectangle in listUntrustedRectangle:
xsDataXDSRectangle = XSDataXDSRectangle()
xsDataXDSRectangle.setX1(XSDataInteger(listRectangle[0]))
xsDataXDSRectangle.setX2(XSDataInteger(listRectangle[1]))
xsDataXDSRectangle.setY1(XSDataInteger(listRectangle[2]))
xsDataXDSRectangle.setY2(XSDataInteger(listRectangle[3]))
xsDataXDSDetector.addUntrusted_rectangle(xsDataXDSRectangle)
xsDataXDSDetector.setMinimum_valid_pixel_value(XSDataInteger(0))
xsDataXDSDetector.setOverload(XSDataInteger(1048500))
xsDataXDSIntegerRangeTrustedPixel = XSDataXDSIntegerRange()
xsDataXDSIntegerRangeTrustedPixel.setLower(XSDataInteger(7000))
xsDataXDSIntegerRangeTrustedPixel.setUpper(XSDataInteger(30000))
xsDataXDSDetector.setValue_range_for_trusted_detector_pixels(xsDataXDSIntegerRangeTrustedPixel)
xsDataXDSDoubleRangeTrustedRegion = XSDataXDSDoubleRange()
xsDataXDSDoubleRangeTrustedRegion.setLower(XSDataDouble(0.0))
xsDataXDSDoubleRangeTrustedRegion.setUpper(XSDataDouble(1.41))
xsDataXDSDetector.setTrusted_region(xsDataXDSDoubleRangeTrustedRegion)
xsDataXDSDetector.setSensor_thickness(XSDataDouble(0.32))
elif strDetectorType == "mar225":
xsDataXDSDetector.setDetector_name(XSDataString("CCDCHESS"))
xsDataXDSDetector.setMinimum_valid_pixel_value(XSDataInteger(0))
xsDataXDSDetector.setOverload(XSDataInteger(65000))
xsDataXDSIntegerRangeTrustedPixel = XSDataXDSIntegerRange()
xsDataXDSIntegerRangeTrustedPixel.setLower(XSDataInteger(7000))
xsDataXDSIntegerRangeTrustedPixel.setUpper(XSDataInteger(30000))
xsDataXDSDetector.setValue_range_for_trusted_detector_pixels(xsDataXDSIntegerRangeTrustedPixel)
xsDataXDSDoubleRangeTrustedRegion = XSDataXDSDoubleRange()
xsDataXDSDoubleRangeTrustedRegion.setLower(XSDataDouble(0.0))
xsDataXDSDoubleRangeTrustedRegion.setUpper(XSDataDouble(1.4))
xsDataXDSDetector.setTrusted_region(xsDataXDSDoubleRangeTrustedRegion)
else:
# This is a temporary solution for the exception problem pointed out in bug #43.
# Instead of raising an exception with a known type we send the error message as a string.
pyStrErrorMessage = "EDHandlerXSDataXDSv1_0.getXSDataXDSDetector: Unknown detector type : " + strDetectorType
raise Exception, pyStrErrorMessage
xsDataXDSDetector.setNx(_xsDataDetector.getNumberPixelX())
xsDataXDSDetector.setNy(_xsDataDetector.getNumberPixelY())
xsDataXDSDetector.setQx(_xsDataDetector.getPixelSizeX())
xsDataXDSDetector.setQy(_xsDataDetector.getPixelSizeY())
xsDataXDSDetector.setDetector_distance(_xsDataDetector.getDistance())
fOrgx = _xsDataDetector.getBeamPositionY().getValue() / _xsDataDetector.getPixelSizeY().getValue()
fOrgy = _xsDataDetector.getBeamPositionX().getValue() / _xsDataDetector.getPixelSizeX().getValue()
xsDataXDSDetector.setOrgx(XSDataDouble(fOrgx))
xsDataXDSDetector.setOrgy(XSDataDouble(fOrgy))
xsDataVectorDoubleXAxis = XSDataVectorDouble()
xsDataVectorDoubleXAxis.setV1(1.0)
xsDataVectorDoubleXAxis.setV2(0.0)
xsDataVectorDoubleXAxis.setV3(0.0)
xsDataXDSDetector.setDirection_of_detector_x_axis(xsDataVectorDoubleXAxis)
xsDataVectorDoubleYAxis = XSDataVectorDouble()
xsDataVectorDoubleYAxis.setV1(0.0)
xsDataVectorDoubleYAxis.setV2(1.0)
xsDataVectorDoubleYAxis.setV3(0.0)
xsDataXDSDetector.setDirection_of_detector_y_axis(xsDataVectorDoubleYAxis)
return xsDataXDSDetector
def _extract_infos(lines, output):
# extract values into the output obj modified in place
# stuff we want:
tmp_params = [
'CRYSTAL MOSAICITY (DEGREES)',
'DIRECT BEAM COORDINATES (REC. ANGSTROEM)',
'DETECTOR COORDINATES (PIXELS) OF DIRECT BEAM',
'DETECTOR ORIGIN (PIXELS) AT',
'CRYSTAL TO DETECTOR DISTANCE (mm)',
'COORDINATES OF UNIT CELL A-AXIS',
'COORDINATES OF UNIT CELL B-AXIS',
'COORDINATES OF UNIT CELL C-AXIS',
'UNIT_CELL_CONSTANTS=',
]
params = [(item, len(item.split())) for item in tmp_params]
parsed = dict()
for line in lines:
for param, param_len in params:
if line.find(param) != -1:
values = [float(x) for x in line.split()[param_len:]]
parsed[param] = values
# now fill in the data model with the values we got (and stop if
# any of them is missing)
try:
output.crystal_mosaicity = XSDataFloat(
parsed['CRYSTAL MOSAICITY (DEGREES)'][0])
output.direct_beam_coordinates = XSDataVectorDouble()
output.direct_beam_coordinates.v1 = parsed[
'DIRECT BEAM COORDINATES (REC. ANGSTROEM)'][0]
output.direct_beam_coordinates.v2 = parsed[
'DIRECT BEAM COORDINATES (REC. ANGSTROEM)'][1]
output.direct_beam_coordinates.v3 = parsed[
'DIRECT BEAM COORDINATES (REC. ANGSTROEM)'][2]
output.direct_beam_detector_coordinates = XSData2DCoordinates()
output.direct_beam_detector_coordinates.x = XSDataFloat(
parsed['DETECTOR COORDINATES (PIXELS) OF DIRECT BEAM'][0])
output.direct_beam_detector_coordinates.y = XSDataFloat(
parsed['DETECTOR COORDINATES (PIXELS) OF DIRECT BEAM'][1])
output.detector_origin = XSData2DCoordinates()
output.detector_origin.x = XSDataFloat(
parsed['DETECTOR ORIGIN (PIXELS) AT'][0])
output.detector_origin.y = XSDataFloat(
parsed['DETECTOR ORIGIN (PIXELS) AT'][1])
output.crystal_to_detector_distance = XSDataFloat(
parsed['CRYSTAL TO DETECTOR DISTANCE (mm)'][0])
output.coordinates_of_unit_cell_a_axis = XSDataVectorDouble()
output.coordinates_of_unit_cell_a_axis.v1 = parsed[
'COORDINATES OF UNIT CELL A-AXIS'][0]
output.coordinates_of_unit_cell_a_axis.v2 = parsed[
'COORDINATES OF UNIT CELL A-AXIS'][1]
output.coordinates_of_unit_cell_a_axis.v3 = parsed[
'COORDINATES OF UNIT CELL A-AXIS'][2]
output.coordinates_of_unit_cell_b_axis = XSDataVectorDouble()
output.coordinates_of_unit_cell_b_axis.v1 = parsed[
'COORDINATES OF UNIT CELL B-AXIS'][0]
output.coordinates_of_unit_cell_b_axis.v2 = parsed[
'COORDINATES OF UNIT CELL B-AXIS'][1]
output.coordinates_of_unit_cell_b_axis.v3 = parsed[
'COORDINATES OF UNIT CELL B-AXIS'][2]
output.coordinates_of_unit_cell_c_axis = XSDataVectorDouble()
output.coordinates_of_unit_cell_c_axis.v1 = parsed[
'COORDINATES OF UNIT CELL C-AXIS'][0]
output.coordinates_of_unit_cell_c_axis.v2 = parsed[
'COORDINATES OF UNIT CELL C-AXIS'][1]
output.coordinates_of_unit_cell_c_axis.v3 = parsed[
'COORDINATES OF UNIT CELL C-AXIS'][2]
parsed['UNIT_CELL_CONSTANTS='] = [
XSDataFloat(x) for x in parsed['UNIT_CELL_CONSTANTS=']
]
# there may be trailing information after the 6 floats (the
# string "as used by INTEGRATE")
unit_cells = parsed['UNIT_CELL_CONSTANTS='][:6]
output.cell_a, output.cell_b, output.cell_c, \
output.cell_alpha, output.cell_beta, output.cell_gamma = unit_cells
except KeyError, ke:
EDVerbose.ERROR('Some parameters could not be found!')
EDVerbose.DEBUG('Those found were: %s' % parsed)
raise ValueError('Some parameters missing')
def generateXSDataInputXDS(_xsDataCollection):
xsDataCollection = _xsDataCollection
xsDataExperimentalCondition = _xsDataCollection.getSubWedge(
)[0].getExperimentalCondition()
xsDataSubWedgeList = xsDataCollection.getSubWedge()
xsDataInputXDS = XSDataInputXDSGenerateBackgroundImage()
xsDataBeam = xsDataExperimentalCondition.getBeam()
xsDataDetector = xsDataExperimentalCondition.getDetector()
xsDataGoniostat = xsDataExperimentalCondition.getGoniostat()
dWavelength = xsDataBeam.getWavelength().getValue()
dDistance = xsDataDetector.getDistance().getValue()
dBeamPositionX = xsDataDetector.getBeamPositionX().getValue()
dBeamPositionY = xsDataDetector.getBeamPositionY().getValue()
# Start with the detector
xsDataXDSDetector = EDHandlerXSDataXDSv1_0.getXSDataXDSDetector(
xsDataDetector)
xsDataInputXDS.setDetector(xsDataXDSDetector)
# Then the beam
xsDataXDSBeam = XSDataXDSBeam()
xsDataVectorDoubleIncidentBeam = XSDataVectorDouble()
xsDataVectorDoubleIncidentBeam.setV1(0.0)
xsDataVectorDoubleIncidentBeam.setV2(0.0)
xsDataVectorDoubleIncidentBeam.setV3(1.0)
xsDataXDSBeam.setIncident_beam_direction(
xsDataVectorDoubleIncidentBeam)
xsDataVectorDoublePolarizationPlaneNormal = XSDataVectorDouble()
xsDataVectorDoublePolarizationPlaneNormal.setV1(0.0)
xsDataVectorDoublePolarizationPlaneNormal.setV2(1.0)
xsDataVectorDoublePolarizationPlaneNormal.setV3(0.0)
xsDataXDSBeam.setPolarization_plane_normal(
xsDataVectorDoublePolarizationPlaneNormal)
xsDataXDSBeam.setX_ray_wavelength(XSDataWavelength(dWavelength))
xsDataInputXDS.setBeam(xsDataXDSBeam)
# Then the goniostat
xsDataXDSGoniostat = XSDataXDSGoniostat()
xsDataVectorDoubleRotationAxis = XSDataVectorDouble()
xsDataVectorDoubleRotationAxis.setV1(1.0)
xsDataVectorDoubleRotationAxis.setV2(0.0)
xsDataVectorDoubleRotationAxis.setV3(0.0)
xsDataXDSGoniostat.setRotation_axis(xsDataVectorDoubleRotationAxis)
xsDataXDSGoniostat.setOscillation_range(
xsDataGoniostat.getOscillationWidth())
xsDataXDSGoniostat.setStarting_angle(
xsDataGoniostat.getRotationAxisStart())
xsDataInputXDS.setGoniostat(xsDataXDSGoniostat)
# # Then the Crystal
#
# xsDataXDSCrystal = XSDataXDSCrystal()
#
# xsDataXDSCrystal.setFriedels_law(XSDataString("FALSE"))
#
## if ( xsDataCrystal is not None ):
## xsDataSpaceGroup = xsDataCrystal.getSpaceGroup()
## if ( xsDataSpaceGroup is not None ):
## xsDataStringName = xsDataSpaceGroup.getName()
## if ( xsDataStringName is not None ):
## xsDataInputXDS.setSymmetry( XSDataString( xsDataStringName.getValue() ) )
# xsDataXDSCrystal.setSpace_group_number(XSDataInteger(0))
#
# xsDataXDSCrystal.setStrong_pixel(XSDataInteger(8))
#
# xsDataCell = XSDataCell()
# xsDataCell.setLength_a(XSDataLength(0.0))
# xsDataCell.setLength_b(XSDataLength(0.0))
# xsDataCell.setLength_c(XSDataLength(0.0))
# xsDataCell.setAngle_alpha(XSDataAngle(0.0))
# xsDataCell.setAngle_beta(XSDataAngle(0.0))
# xsDataCell.setAngle_gamma(XSDataAngle(0.0))
# xsDataXDSCrystal.setUnit_cell_constants(xsDataCell)
#
# xsDataInputXDS.setCrystal(xsDataXDSCrystal)
# Finaly the images
xsDataXDSImage = XSDataXDSImage()
xsDataSubWedgeFirst = xsDataSubWedgeList[0]
xsDataImageFirst = xsDataSubWedgeFirst.getImage()[0]
pyStrPath = xsDataImageFirst.getPath().getValue()
pyStrFileName = EDUtilsFile.getBaseName(pyStrPath)
pyStrDirectory = EDUtilsPath.getFolderName(pyStrPath)
pyStrPrefix = EDUtilsImage.getPrefix(pyStrFileName)
pyStrSuffix = EDUtilsImage.getSuffix(pyStrFileName)
pyStrXDSTemplate = "%s_xdslink_?????.%s" % (pyStrPrefix, pyStrSuffix)
xsDataXDSImage.setName_template_of_data_frames(
XSDataString(pyStrXDSTemplate))
iXDSLowestImageNumberGlobal = 1
xsDataXDSImage.setStarting_frame(
XSDataInteger(iXDSLowestImageNumberGlobal))
# First we have to find the smallest goniostat rotation axis start:
fGonioStatOscillationStartMin = None
for xsDataSubWedge in xsDataSubWedgeList:
xsDataGoniostat = xsDataSubWedge.getExperimentalCondition(
).getGoniostat()
fGonioStatOscillationStart = xsDataGoniostat.getRotationAxisStart(
).getValue()
if (fGonioStatOscillationStartMin is None):
fGonioStatOscillationStartMin = fGonioStatOscillationStart
elif (fGonioStatOscillationStartMin > fGonioStatOscillationStart):
fGonioStatOscillationStartMin = fGonioStatOscillationStart
# Loop through the list of sub wedges
for xsDataSubWedge in xsDataSubWedgeList:
xsDataImageList = xsDataSubWedge.getImage()
xsDataGoniostat = xsDataSubWedge.getExperimentalCondition(
).getGoniostat()
fGonioStatOscillationStart = xsDataGoniostat.getRotationAxisStart(
).getValue()
fGonioStatOscillationRange = xsDataGoniostat.getOscillationWidth(
).getValue()
# First find the lowest and highest image numbers
iLowestImageNumber = None
for xsDataImage in xsDataImageList:
iImageNumber = xsDataImage.getNumber().getValue()
if (iLowestImageNumber is None):
iLowestImageNumber = iImageNumber
elif (iImageNumber < iLowestImageNumber):
iLowestImageNumber = iImageNumber
# Loop through the list of images
iLowestXDSImageNumber = None
iHighestXDSImageNumber = None
for xsDataImage in xsDataImageList:
iImageNumber = xsDataImage.getNumber().getValue()
fImageOscillationStart = fGonioStatOscillationStart + (
iImageNumber -
iLowestImageNumber) * fGonioStatOscillationRange
iXDSImageNumber = iXDSLowestImageNumberGlobal + int(
(fImageOscillationStart - fGonioStatOscillationStartMin) /
fGonioStatOscillationRange)
#print iXDSImageNumber, fImageOscillationStart, fGonioStatOscillationStartMin, fGonioStatOscillationRange
pyStrSourcePath = xsDataImage.getPath()
pyStrTarget = "%s_xdslink_%05d.%s" % (
pyStrPrefix, iXDSImageNumber, pyStrSuffix)
xsDataXDSImageLink = XSDataXDSImageLink()
xsDataFileSource = XSDataFile()
xsDataFileSource.setPath(pyStrSourcePath)
xsDataXDSImageLink.setSource(xsDataFileSource)
xsDataXDSImageLink.setTarget(XSDataString(pyStrTarget))
xsDataInputXDS.addImage_link(xsDataXDSImageLink)
if (iLowestXDSImageNumber is None):
iLowestXDSImageNumber = iXDSImageNumber
elif (iLowestXDSImageNumber > iXDSImageNumber):
iLowestXDSImageNumber = iXDSImageNumber
if (iHighestXDSImageNumber is None):
iHighestXDSImageNumber = iXDSImageNumber
elif (iHighestXDSImageNumber < iXDSImageNumber):
iHighestXDSImageNumber = iXDSImageNumber
xsDataXDSIntegerRange = XSDataXDSIntegerRange()
xsDataXDSIntegerRange.setLower(
XSDataInteger(iLowestXDSImageNumber))
xsDataXDSIntegerRange.setUpper(
XSDataInteger(iHighestXDSImageNumber))
xsDataXDSImage.addBackground_range(xsDataXDSIntegerRange)
xsDataXDSImage.addData_range(xsDataXDSIntegerRange)
xsDataXDSImage.addSpot_range(xsDataXDSIntegerRange)
xsDataInputXDS.setImage(xsDataXDSImage)
return xsDataInputXDS
def getXSDataXDSDetector(_xsDataDetector):
EDFactoryPluginStatic.loadModule("XSDataXDSv1_0")
from XSDataXDSv1_0 import XSDataXDSDetector
from XSDataXDSv1_0 import XSDataXDSIntegerRange
xsDataXDSDetector = XSDataXDSDetector()
strDetectorType = _xsDataDetector.getType().getValue()
if ((strDetectorType == "q4") or \
(strDetectorType == "q4-2x") or \
(strDetectorType == "q210") or \
(strDetectorType == "q210-2x") or \
(strDetectorType == "q315") or \
(strDetectorType == "q315-2x")):
xsDataXDSDetector.setDetector_name(XSDataString("ADSC"))
xsDataXDSDetector.setMinimum_valid_pixel_value(XSDataInteger(1))
xsDataXDSDetector.setOverload(XSDataInteger(65535))
xsDataXDSIntegerRangeTrustedPixel = XSDataXDSIntegerRange()
xsDataXDSIntegerRangeTrustedPixel.setLower(XSDataInteger(6000))
xsDataXDSIntegerRangeTrustedPixel.setUpper(XSDataInteger(30000))
xsDataXDSDetector.setValue_range_for_trusted_detector_pixels(
xsDataXDSIntegerRangeTrustedPixel)
# elif ( strDetectorType == "mar165") or \
# strDetectorType == "mar225") ):
# xsDataXDSDetector.setType( XSDataString( "MARCCD" ) )
elif strDetectorType == "pilatus6m":
xsDataXDSDetector.setDetector_name(XSDataString("PILATUS"))
listUntrustedRectangle = \
[[ 487, 495, 0, 2528], \
[ 981, 989, 0, 2528], \
[1475, 1483, 0, 2528], \
[1969, 1977, 0, 2528], \
[ 0, 2464, 195, 213], \
[ 0, 2464, 407, 425], \
[ 0, 2464, 619, 637], \
[ 0, 2464, 831, 849], \
[ 0, 2464, 1043, 1061], \
[ 0, 2464, 1255, 1273], \
[ 0, 2464, 1467, 1485], \
[ 0, 2464, 1679, 1697], \
[ 0, 2464, 1891, 1909], \
[ 0, 2464, 2103, 2121], \
[ 0, 2464, 2315, 2333]]
for listRectangle in listUntrustedRectangle:
xsDataXDSRectangle = XSDataXDSRectangle()
xsDataXDSRectangle.setX1(XSDataInteger(listRectangle[0]))
xsDataXDSRectangle.setX2(XSDataInteger(listRectangle[1]))
xsDataXDSRectangle.setY1(XSDataInteger(listRectangle[2]))
xsDataXDSRectangle.setY2(XSDataInteger(listRectangle[3]))
xsDataXDSDetector.addUntrusted_rectangle(xsDataXDSRectangle)
xsDataXDSDetector.setMinimum_valid_pixel_value(XSDataInteger(0))
xsDataXDSDetector.setOverload(XSDataInteger(1048500))
xsDataXDSIntegerRangeTrustedPixel = XSDataXDSIntegerRange()
xsDataXDSIntegerRangeTrustedPixel.setLower(XSDataInteger(7000))
xsDataXDSIntegerRangeTrustedPixel.setUpper(XSDataInteger(30000))
xsDataXDSDetector.setValue_range_for_trusted_detector_pixels(
xsDataXDSIntegerRangeTrustedPixel)
xsDataXDSDoubleRangeTrustedRegion = XSDataXDSDoubleRange()
xsDataXDSDoubleRangeTrustedRegion.setLower(XSDataDouble(0.0))
xsDataXDSDoubleRangeTrustedRegion.setUpper(XSDataDouble(1.41))
xsDataXDSDetector.setTrusted_region(
xsDataXDSDoubleRangeTrustedRegion)
xsDataXDSDetector.setSensor_thickness(XSDataDouble(0.32))
else:
# This is a temporary solution for the exception problem pointed out in bug #43.
# Instead of raising an exception with a known type we send the error message as a string.
pyStrErrorMessage = "EDHandlerXSDataXDSv1_0.getXSDataXDSDetector: Unknown detector type : " + strDetectorType
raise Exception, pyStrErrorMessage
xsDataXDSDetector.setNx(_xsDataDetector.getNumberPixelX())
xsDataXDSDetector.setNy(_xsDataDetector.getNumberPixelY())
xsDataXDSDetector.setQx(_xsDataDetector.getPixelSizeX())
xsDataXDSDetector.setQy(_xsDataDetector.getPixelSizeY())
xsDataXDSDetector.setDetector_distance(_xsDataDetector.getDistance())
fOrgx = _xsDataDetector.getBeamPositionY().getValue(
) / _xsDataDetector.getPixelSizeY().getValue()
fOrgy = _xsDataDetector.getBeamPositionX().getValue(
) / _xsDataDetector.getPixelSizeX().getValue()
xsDataXDSDetector.setOrgx(XSDataDouble(fOrgx))
xsDataXDSDetector.setOrgy(XSDataDouble(fOrgy))
xsDataVectorDoubleXAxis = XSDataVectorDouble()
xsDataVectorDoubleXAxis.setV1(1.0)
xsDataVectorDoubleXAxis.setV2(0.0)
xsDataVectorDoubleXAxis.setV3(0.0)
xsDataXDSDetector.setDirection_of_detector_x_axis(
xsDataVectorDoubleXAxis)
xsDataVectorDoubleYAxis = XSDataVectorDouble()
xsDataVectorDoubleYAxis.setV1(0.0)
xsDataVectorDoubleYAxis.setV2(1.0)
xsDataVectorDoubleYAxis.setV3(0.0)
xsDataXDSDetector.setDirection_of_detector_y_axis(
xsDataVectorDoubleYAxis)
return xsDataXDSDetector
def getXSDataXDSDetector(_xsDataDetector):
EDFactoryPluginStatic.loadModule("XSDataXDSv1_0")
from XSDataXDSv1_0 import XSDataXDSDetector
from XSDataXDSv1_0 import XSDataXDSIntegerRange
xsDataXDSDetector = XSDataXDSDetector()
strDetectorType = _xsDataDetector.getType().getValue()
if ((strDetectorType == "q4") or \
(strDetectorType == "q4-2x") or \
(strDetectorType == "q210") or \
(strDetectorType == "q210-2x") or \
(strDetectorType == "q315") or \
(strDetectorType == "q315-2x")):
xsDataXDSDetector.setDetector_name(XSDataString("ADSC"))
xsDataXDSDetector.setMinimum_valid_pixel_value(XSDataInteger(1))
xsDataXDSDetector.setOverload(XSDataInteger(65535))
xsDataXDSIntegerRangeTrustedPixel = XSDataXDSIntegerRange()
xsDataXDSIntegerRangeTrustedPixel.setLower(XSDataInteger(6000))
xsDataXDSIntegerRangeTrustedPixel.setUpper(XSDataInteger(30000))
xsDataXDSDetector.setValue_range_for_trusted_detector_pixels(xsDataXDSIntegerRangeTrustedPixel)
# elif ( strDetectorType == "mar165") or \
# strDetectorType == "mar225") ):
# xsDataXDSDetector.setType( XSDataString( "MARCCD" ) )
elif strDetectorType == "pilatus6m":
xsDataXDSDetector.setDetector_name(XSDataString("PILATUS"))
listUntrustedRectangle = \
[[ 487, 495, 0, 2528], \
[ 981, 989, 0, 2528], \
[1475, 1483, 0, 2528], \
[1969, 1977, 0, 2528], \
[ 0, 2464, 195, 213], \
[ 0, 2464, 407, 425], \
[ 0, 2464, 619, 637], \
[ 0, 2464, 831, 849], \
[ 0, 2464, 1043, 1061], \
[ 0, 2464, 1255, 1273], \
[ 0, 2464, 1467, 1485], \
[ 0, 2464, 1679, 1697], \
[ 0, 2464, 1891, 1909], \
[ 0, 2464, 2103, 2121], \
[ 0, 2464, 2315, 2333]]
for listRectangle in listUntrustedRectangle:
xsDataXDSRectangle = XSDataXDSRectangle()
xsDataXDSRectangle.setX1(XSDataInteger(listRectangle[0]))
xsDataXDSRectangle.setX2(XSDataInteger(listRectangle[1]))
xsDataXDSRectangle.setY1(XSDataInteger(listRectangle[2]))
xsDataXDSRectangle.setY2(XSDataInteger(listRectangle[3]))
xsDataXDSDetector.addUntrusted_rectangle(xsDataXDSRectangle)
xsDataXDSDetector.setMinimum_valid_pixel_value(XSDataInteger(0))
xsDataXDSDetector.setOverload(XSDataInteger(1048500))
xsDataXDSIntegerRangeTrustedPixel = XSDataXDSIntegerRange()
xsDataXDSIntegerRangeTrustedPixel.setLower(XSDataInteger(7000))
xsDataXDSIntegerRangeTrustedPixel.setUpper(XSDataInteger(30000))
xsDataXDSDetector.setValue_range_for_trusted_detector_pixels(xsDataXDSIntegerRangeTrustedPixel)
xsDataXDSDoubleRangeTrustedRegion = XSDataXDSDoubleRange()
xsDataXDSDoubleRangeTrustedRegion.setLower(XSDataDouble(0.0))
xsDataXDSDoubleRangeTrustedRegion.setUpper(XSDataDouble(1.41))
xsDataXDSDetector.setTrusted_region(xsDataXDSDoubleRangeTrustedRegion)
xsDataXDSDetector.setSensor_thickness(XSDataDouble(0.32))
elif strDetectorType == "pilatus2m":
xsDataXDSDetector.setDetector_name(XSDataString("PILATUS"))
listUntrustedRectangle = \
[[ 487, 495, 0, 1680], \
[ 981, 989, 0, 1680], \
[ 0, 1476, 195, 213], \
[ 0, 1476, 407, 425], \
[ 0, 1476, 619, 637], \
[ 0, 1476, 831, 849], \
[ 0, 1476, 1043, 1061], \
[ 0, 1476, 1255, 1273], \
[ 0, 1476, 1467, 1485]]
for listRectangle in listUntrustedRectangle:
xsDataXDSRectangle = XSDataXDSRectangle()
xsDataXDSRectangle.setX1(XSDataInteger(listRectangle[0]))
xsDataXDSRectangle.setX2(XSDataInteger(listRectangle[1]))
xsDataXDSRectangle.setY1(XSDataInteger(listRectangle[2]))
xsDataXDSRectangle.setY2(XSDataInteger(listRectangle[3]))
xsDataXDSDetector.addUntrusted_rectangle(xsDataXDSRectangle)
xsDataXDSDetector.setMinimum_valid_pixel_value(XSDataInteger(0))
xsDataXDSDetector.setOverload(XSDataInteger(1048500))
xsDataXDSIntegerRangeTrustedPixel = XSDataXDSIntegerRange()
xsDataXDSIntegerRangeTrustedPixel.setLower(XSDataInteger(7000))
xsDataXDSIntegerRangeTrustedPixel.setUpper(XSDataInteger(30000))
xsDataXDSDetector.setValue_range_for_trusted_detector_pixels(xsDataXDSIntegerRangeTrustedPixel)
xsDataXDSDoubleRangeTrustedRegion = XSDataXDSDoubleRange()
xsDataXDSDoubleRangeTrustedRegion.setLower(XSDataDouble(0.0))
xsDataXDSDoubleRangeTrustedRegion.setUpper(XSDataDouble(1.41))
xsDataXDSDetector.setTrusted_region(xsDataXDSDoubleRangeTrustedRegion)
xsDataXDSDetector.setSensor_thickness(XSDataDouble(0.32))
elif strDetectorType == "eiger4m":
xsDataXDSDetector.setDetector_name(XSDataString("PILATUS"))
listUntrustedRectangle = \
[[ 1029, 1040, 0, 2167], \
[ 0, 2070, 512, 550], \
[ 0, 2070, 1063, 1103], \
[ 0, 2070, 1614, 1654],
]
for listRectangle in listUntrustedRectangle:
xsDataXDSRectangle = XSDataXDSRectangle()
xsDataXDSRectangle.setX1(XSDataInteger(listRectangle[0]))
xsDataXDSRectangle.setX2(XSDataInteger(listRectangle[1]))
xsDataXDSRectangle.setY1(XSDataInteger(listRectangle[2]))
xsDataXDSRectangle.setY2(XSDataInteger(listRectangle[3]))
xsDataXDSDetector.addUntrusted_rectangle(xsDataXDSRectangle)
xsDataXDSDetector.setMinimum_valid_pixel_value(XSDataInteger(0))
xsDataXDSDetector.setOverload(XSDataInteger(1048500))
xsDataXDSIntegerRangeTrustedPixel = XSDataXDSIntegerRange()
xsDataXDSIntegerRangeTrustedPixel.setLower(XSDataInteger(7000))
xsDataXDSIntegerRangeTrustedPixel.setUpper(XSDataInteger(30000))
xsDataXDSDetector.setValue_range_for_trusted_detector_pixels(xsDataXDSIntegerRangeTrustedPixel)
xsDataXDSDoubleRangeTrustedRegion = XSDataXDSDoubleRange()
xsDataXDSDoubleRangeTrustedRegion.setLower(XSDataDouble(0.0))
xsDataXDSDoubleRangeTrustedRegion.setUpper(XSDataDouble(1.41))
xsDataXDSDetector.setTrusted_region(xsDataXDSDoubleRangeTrustedRegion)
xsDataXDSDetector.setSensor_thickness(XSDataDouble(0.32))
elif strDetectorType == "eiger9m":
xsDataXDSDetector.setDetector_name(XSDataString("PILATUS"))
listUntrustedRectangle = \
[[ 1029, 1040, 0, 3269], \
[ 2069, 2082, 0, 3269], \
[ 0, 3110, 513, 553], \
[ 0, 3110, 1064, 1104], \
[ 0, 3110, 1615, 1655], \
[ 0, 3110, 2166, 2206], \
[ 0, 3110, 2717, 2757], \
]
for listRectangle in listUntrustedRectangle:
xsDataXDSRectangle = XSDataXDSRectangle()
xsDataXDSRectangle.setX1(XSDataInteger(listRectangle[0]))
xsDataXDSRectangle.setX2(XSDataInteger(listRectangle[1]))
xsDataXDSRectangle.setY1(XSDataInteger(listRectangle[2]))
xsDataXDSRectangle.setY2(XSDataInteger(listRectangle[3]))
xsDataXDSDetector.addUntrusted_rectangle(xsDataXDSRectangle)
xsDataXDSDetector.setMinimum_valid_pixel_value(XSDataInteger(0))
xsDataXDSDetector.setOverload(XSDataInteger(1048500))
xsDataXDSIntegerRangeTrustedPixel = XSDataXDSIntegerRange()
xsDataXDSIntegerRangeTrustedPixel.setLower(XSDataInteger(7000))
xsDataXDSIntegerRangeTrustedPixel.setUpper(XSDataInteger(30000))
xsDataXDSDetector.setValue_range_for_trusted_detector_pixels(xsDataXDSIntegerRangeTrustedPixel)
xsDataXDSDoubleRangeTrustedRegion = XSDataXDSDoubleRange()
xsDataXDSDoubleRangeTrustedRegion.setLower(XSDataDouble(0.0))
xsDataXDSDoubleRangeTrustedRegion.setUpper(XSDataDouble(1.41))
xsDataXDSDetector.setTrusted_region(xsDataXDSDoubleRangeTrustedRegion)
xsDataXDSDetector.setSensor_thickness(XSDataDouble(0.45))
elif strDetectorType == "eiger16m":
xsDataXDSDetector.setDetector_name(XSDataString("PILATUS"))
# From https://www.psi.ch/sls/pxi/ComputingEN/XDS.INP.E16M.X06SA_2016.txt
# !EXCLUSION OF HORIZONTAL DEAD AREAS OF THE EIGER 16M DETECTOR + ONE PIXEL ON EACH SIDE
# UNTRUSTED_RECTANGLE= 0 4151 513 553
# UNTRUSTED_RECTANGLE= 0 4151 1064 1104
# UNTRUSTED_RECTANGLE= 0 4151 1615 1655
# UNTRUSTED_RECTANGLE= 0 4151 2166 2206
# UNTRUSTED_RECTANGLE= 0 4151 2717 2757
# UNTRUSTED_RECTANGLE= 0 4151 3268 3308
# UNTRUSTED_RECTANGLE= 0 4151 3819 3859
# !EXCLUSION OF VERTICAL DEAD AREAS OF THE EIGER 16M DETECTOR + ONE PIXEL ON EACH SIDE
# UNTRUSTED_RECTANGLE= 1029 1042 0 4372
# UNTRUSTED_RECTANGLE= 2069 2082 0 4372
# UNTRUSTED_RECTANGLE= 3109 3122 0 4372
listUntrustedRectangle = \
[
[ 0, 4151, 513, 553],
[ 0, 4151, 1064, 1104],
[ 0, 4151, 1615, 1655],
[ 0, 4151, 2166, 2206],
[ 0, 4151, 2717, 2757],
[ 0, 4151, 3268, 3308],
[ 0, 4151, 3819, 3859],
[ 1029, 1042, 0, 4372],
[ 2069, 2082, 0, 4372],
[ 3109, 3122, 0, 4372],
]
for listRectangle in listUntrustedRectangle:
xsDataXDSRectangle = XSDataXDSRectangle()
xsDataXDSRectangle.setX1(XSDataInteger(listRectangle[0]))
xsDataXDSRectangle.setX2(XSDataInteger(listRectangle[1]))
xsDataXDSRectangle.setY1(XSDataInteger(listRectangle[2]))
xsDataXDSRectangle.setY2(XSDataInteger(listRectangle[3]))
xsDataXDSDetector.addUntrusted_rectangle(xsDataXDSRectangle)
xsDataXDSDetector.setMinimum_valid_pixel_value(XSDataInteger(0))
xsDataXDSDetector.setOverload(XSDataInteger(1048500))
xsDataXDSIntegerRangeTrustedPixel = XSDataXDSIntegerRange()
xsDataXDSIntegerRangeTrustedPixel.setLower(XSDataInteger(4000))
xsDataXDSIntegerRangeTrustedPixel.setUpper(XSDataInteger(30000))
xsDataXDSDetector.setValue_range_for_trusted_detector_pixels(xsDataXDSIntegerRangeTrustedPixel)
xsDataXDSDoubleRangeTrustedRegion = XSDataXDSDoubleRange()
xsDataXDSDoubleRangeTrustedRegion.setLower(XSDataDouble(0.0))
xsDataXDSDoubleRangeTrustedRegion.setUpper(XSDataDouble(1.21))
xsDataXDSDetector.setTrusted_region(xsDataXDSDoubleRangeTrustedRegion)
xsDataXDSDetector.setSensor_thickness(XSDataDouble(0.32))
elif strDetectorType == "mar225":
xsDataXDSDetector.setDetector_name(XSDataString("CCDCHESS"))
xsDataXDSDetector.setMinimum_valid_pixel_value(XSDataInteger(0))
xsDataXDSDetector.setOverload(XSDataInteger(65000))
xsDataXDSIntegerRangeTrustedPixel = XSDataXDSIntegerRange()
xsDataXDSIntegerRangeTrustedPixel.setLower(XSDataInteger(7000))
xsDataXDSIntegerRangeTrustedPixel.setUpper(XSDataInteger(30000))
xsDataXDSDetector.setValue_range_for_trusted_detector_pixels(xsDataXDSIntegerRangeTrustedPixel)
xsDataXDSDoubleRangeTrustedRegion = XSDataXDSDoubleRange()
xsDataXDSDoubleRangeTrustedRegion.setLower(XSDataDouble(0.0))
xsDataXDSDoubleRangeTrustedRegion.setUpper(XSDataDouble(1.4))
xsDataXDSDetector.setTrusted_region(xsDataXDSDoubleRangeTrustedRegion)
else:
# This is a temporary solution for the exception problem pointed out in bug #43.
# Instead of raising an exception with a known type we send the error message as a string.
pyStrErrorMessage = "EDHandlerXSDataXDSv1_0.getXSDataXDSDetector: Unknown detector type : " + strDetectorType
raise BaseException(pyStrErrorMessage)
xsDataXDSDetector.setNx(_xsDataDetector.getNumberPixelX())
xsDataXDSDetector.setNy(_xsDataDetector.getNumberPixelY())
xsDataXDSDetector.setQx(_xsDataDetector.getPixelSizeX())
xsDataXDSDetector.setQy(_xsDataDetector.getPixelSizeY())
xsDataXDSDetector.setDetector_distance(_xsDataDetector.getDistance())
fOrgx = _xsDataDetector.getBeamPositionY().getValue() / _xsDataDetector.getPixelSizeY().getValue()
fOrgy = _xsDataDetector.getBeamPositionX().getValue() / _xsDataDetector.getPixelSizeX().getValue()
xsDataXDSDetector.setOrgx(XSDataDouble(fOrgx))
xsDataXDSDetector.setOrgy(XSDataDouble(fOrgy))
xsDataVectorDoubleXAxis = XSDataVectorDouble()
xsDataVectorDoubleXAxis.setV1(1.0)
xsDataVectorDoubleXAxis.setV2(0.0)
xsDataVectorDoubleXAxis.setV3(0.0)
xsDataXDSDetector.setDirection_of_detector_x_axis(xsDataVectorDoubleXAxis)
xsDataVectorDoubleYAxis = XSDataVectorDouble()
xsDataVectorDoubleYAxis.setV1(0.0)
xsDataVectorDoubleYAxis.setV2(1.0)
xsDataVectorDoubleYAxis.setV3(0.0)
xsDataXDSDetector.setDirection_of_detector_y_axis(xsDataVectorDoubleYAxis)
return xsDataXDSDetector