Generate XDS.INP: Difference between revisions

# revision 0.96 . KD 03/10/2020 make script echo detector serial number if inverted ROTATION_AXIS is detected.

# revision 0.97 . KD 19/10/2020 add UNTRUSTED_RECTANGLEs for Eiger2 (which has a few pixels less than Eiger (thanks to Andreas Förster)

#                                                                                                             

  grep -q "^SOURCE_WAVELENGTH= *1" tmp2 && DET=dtrek

  grep -q BEAM_CENTER_X tmp2      && DET=adsc                                 

  grep -q BRUKER tmp2 && grep -q CBF tmp2  && DET=Bruker-cbf   

  grep -q CMOS1 tmp2 && DET=adsc-CMOS1   

         fi

       fi

elif [ "$DET" == "eiger" ]; then

   if [ "$DLS" == 1 ]; then

     echo Eiger HDF5 from Diamond

     OVERLOAD=`h5dump -d "/entry/instrument/detector/saturation_value" $FIRSTFRAME | awk '/\(0\):/{print $2}'`

   # v0.95: fix the next 2 lines by taking care of negative values with the \- , and stop after first "(0)"

   # If rotation vector set (NeXus)

     rotation_axis=`h5dump -a "/entry/sample/transformations/omega/vector" $FIRSTFRAME 2>/dev/null | grep "(0):" | sed -e "s/^.*://; s/,//g"`

     SN=`h5dump -d "/entry/instrument/detector/detector_number" $FIRSTFRAME | awk '/\(0\): /{print $2}' | sed s/\"//g`

       rotation_axis="-1 0 0"

     fi  

   fi

   #let SKIP=768

   #NX=$(od -t x -j $SKIP -N 4 $FIRSTFRAME |awk 'NR==1{print toupper($2)}'|perl -nle '@array= $_ =~/.{2}/g; print "ibase=16;obase=A;".join("",reverse @array)'|bc)

elif [ "$DET" == "dtrek" ]; then

  DETECTOR="MAR345  MINIMUM_VALID_PIXEL_VALUE=0  OVERLOAD=130000"

  NX=`awk '/FORMAT/{print $2}' tmp2`

  QX=`awk '/PIXEL/{print $3/1000.}' tmp2`

  QY=`awk '/PIXEL/{print $5/1000.}' tmp2`

!                                                  and does not seem to have any downsides.  

eof

# Eiger2 1M ; v0.97 numbers from Andreas Förster

     cat >> XDS.INP << eof

UNTRUSTED_RECTANGLE=    0 1029    512  551

eof

     if ! grep -q Flat_field tmp2 ; then

     cat >> XDS.INP << eof

eof

     fi

# Eiger2 4M ; v0.97 numbers from Andreas Förster

     cat >> XDS.INP << eof

  UNTRUSTED_RECTANGLE=    0 2069  1612 1651

eof

# Pilatus 6M

# FIXME: here we could test if a Flat_field correction was applied like we do for 2M

UNTRUSTED_RECTANGLE=  0 2464  2315 2333

eof

# Eiger 9M

     cat >> XDS.INP << eof

UNTRUSTED_RECTANGLE= 0 3110 2717 2757

eof

# Eiger2 9M ; v0.97 numbers from Andreas Förster

     cat >> XDS.INP << eof

UNTRUSTED_RECTANGLE=    0 3109  2712 2751

eof

# Eiger 16M

     cat >> XDS.INP << eof

  UNTRUSTED_RECTANGLE= 3109 3122      0 4371

eof

# Eiger2 16M ; v0.97 numbers from Andreas Förster

     cat >> XDS.INP << eof

* The script tries to interpret the header of the frames, so is currently limited to Dectris (Pilatus, Eiger), ADSC (Quantum), Rigaku (several types), MAR (CCD and image plate) detectors, and one Bruker detector. Other detectors need some values to be manually filled into XDS.INP - the relevant places are marked with XXX. These are detector properties (type, pixel size and number, min and max counts in a pixel), and experimental parameters like ROTATION_AXIS, OSCILLATION_RANGE, X-RAY_WAVELENGTH, DETECTOR_DISTANCE, and XORG, YORG. For fine-tuning of detector parameters, see the [http://xds.mpimf-heidelberg.mpg.de/html_doc/xds_prepare.html detector-specific templates].

* The authors have made a "best effort" to provide a XDS.INP that results in the correct sign of the anomalous signal. In the case of one detector type (internally called Rigaku SMV) this requires reversal of one detector axis, and a negative DETECTOR_DISTANCE, as is found in some of the [http://xds.mpimf-heidelberg.mpg.de/html_doc/xds_prepare.html detector-specific templates]. '''For an unusual or unknown detector setup, the correct sign of the anomalous signal needs to be established and verified e.g. with a good dataset from a test crystal that has an anomalous signal.''' The authors do not take any responsibility for problems arising from incorrect sign of the anomalous signal, nor - obviously! - for any other mischief arising in or from data processing.

* Sometimes, the x- and y- values of the primary beam position recorded in the header should be used for ORGY and ORGX (i.e reversed) instead of as ORGX and ORGY. For ADSC, this has been implemented in the script for a number of beamlines.

== See also ==