QxrdDetectorGeometry Class Reference

#include <qxrddetectorgeometry.h>

Inheritance diagram for QxrdDetectorGeometry:

Collaboration diagram for QxrdDetectorGeometry:

Public Slots
static double	convertEnergyToWavelength (double energy)
static double	convertWavelengthToEnergy (double wavelength)
static double	convertTwoThetaToQ (double twoTheta, double wavelength)
static double	convertQToTwoTheta (double Q, double wavelength)
static double	getRadius (double xCenter, double yCenter, double distance, double xPixel, double yPixel, double pixelLength, double pixelHeight, double cos_beta, double sin_beta, double cos_rotation, double sin_rotation)
static double	getTwoTheta (double xCenter, double yCenter, double distance, double xPixel, double yPixel, double pixelLength, double pixelHeight, double cos_beta, double sin_beta, double cos_rotation, double sin_rotation)
static void	getTwoThetaChi (double xCenter, double yCenter, double distance, double xPixel, double yPixel, double pixelLength, double pixelHeight, double rotation, double cos_beta, double sin_beta, double cos_alpha, double sin_alpha, double cos_rotation, double sin_rotation, double *twoTheta, double *chi)
static void	getQChi (double xCenter, double yCenter, double distance, double energy, double xPixel, double yPixel, double pixelLength, double pixelHeight, double rotation, double cos_beta, double sin_beta, double cos_alpha, double sin_alpha, double cos_rotation, double sin_rotation, double *q, double *chi)
static void	getXY (double xCenter, double yCenter, double distance, double energy, double q, double chi, double pixelLength, double pixelHeight, double rotation, double cos_beta, double sin_beta, double cos_alpha, double sin_alpha, double cos_rotation, double sin_rotation, double *xPixel, double *yPixel)
Public Slots inherited from QcepObject
virtual void	printLine (QString line)
virtual void	printMessage (QString msg, QDateTime dt=QDateTime::currentDateTime()) const
virtual void	criticalMessage (QString msg, QDateTime ts=QDateTime::currentDateTime()) const
virtual void	statusMessage (QString msg, QDateTime ts=QDateTime::currentDateTime()) const
virtual QString	settingsScript ()
QString	scriptValueLiteral (QVariant v)

Public Member Functions
	QxrdDetectorGeometry (QString name, QcepObject *parent)
Public Member Functions inherited from QcepObject
	QcepObject (QString name, QcepObject *parent)
virtual	~QcepObject ()
virtual void	writeSettings (QSettings *set, QString section)
virtual void	readSettings (QSettings *set, QString section)
QString	get_Name () const
void	set_Name (QString name)

Additional Inherited Members
Static Public Member Functions inherited from QcepObject
static int	allocatedObjects ()
static int	deletedObjects ()
static QSet< QcepObject * >	allocatedObjectsSet ()
static QString	addSlashes (QString str)
Properties inherited from QcepObject
QString	name

Detailed Description

Definition at line 8 of file qxrddetectorgeometry.h.

Constructor & Destructor Documentation

QxrdDetectorGeometry::QxrdDetectorGeometry ( QString  name, QcepObject *  parent  )

explicit

Definition at line 5 of file qxrddetectorgeometry.cpp.

                                                                           :
    QcepObject(name, parent)
{
}

Member Function Documentation

double QxrdDetectorGeometry::convertEnergyToWavelength ( double  energy )

staticslot

Definition at line 22 of file qxrddetectorgeometry.cpp.

Referenced by QxrdCenterFinder::getXY().

{
  return 12398.4172/energy;
}

Here is the caller graph for this function:

double QxrdDetectorGeometry::convertQToTwoTheta ( double  Q, double  wavelength  )

staticslot

Definition at line 52 of file qxrddetectorgeometry.cpp.

{
  // convert twoTheta to degrees before returning
  return 2.0*asin(Q*wavelength/(4*M_PI))*180.0/M_PI;
}

double QxrdDetectorGeometry::convertTwoThetaToQ ( double  twoTheta, double  wavelength  )

staticslot

Definition at line 45 of file qxrddetectorgeometry.cpp.

Referenced by getQChi(), and QxrdCenterFinder::getXY().

{
  // convert twoTheta from degrees to radians first
  return 4.0*M_PI*sin( (twoTheta*M_PI/180.0) /2.0)/wavelength;
}

Here is the caller graph for this function:

double QxrdDetectorGeometry::convertWavelengthToEnergy ( double  wavelength )

staticslot

Definition at line 36 of file qxrddetectorgeometry.cpp.

{
  return 12398.4172/wavelength;
}

void QxrdDetectorGeometry::getQChi ( double  xCenter, double  yCenter, double  distance, double  energy, double  xPixel, double  yPixel, double  pixelLength, double  pixelHeight, double  rotation, double  cos_beta, double  sin_beta, double  cos_alpha, double  sin_alpha, double  cos_rotation, double  sin_rotation, double *  q, double *  chi  )

staticslot

Definition at line 226 of file qxrddetectorgeometry.cpp.

References convertTwoThetaToQ(), and getTwoThetaChi().

Referenced by QxrdIntegratorCache::getChi(), QxrdCenterFinder::getChi(), QxrdIntegratorCache::getQ(), and QxrdCenterFinder::getQ().

{

  double wavelength;
  double twoTheta;
  wavelength = 12398.4172/energy;

  getTwoThetaChi(xCenter,yCenter,distance,xPixel,yPixel,
                 pixelLength,pixelHeight,rotation,cos_beta,sin_beta,
                 cos_alpha,sin_alpha,cos_rotation,sin_rotation,
                 &twoTheta,chi);

  *q=convertTwoThetaToQ(twoTheta,wavelength);
}

Here is the call graph for this function:

Here is the caller graph for this function:

double QxrdDetectorGeometry::getRadius ( double  xCenter, double  yCenter, double  distance, double  xPixel, double  yPixel, double  pixelLength, double  pixelHeight, double  cos_beta, double  sin_beta, double  cos_rotation, double  sin_rotation  )

staticslot

cos_alpha

Definition at line 58 of file qxrddetectorgeometry.cpp.

{
  double pixelLength_mm,pixelHeight_mm;
  double xMeasured,yMeasured;
  double bottom;
  double xPhysical,yPhysical;

  // pixellength comes in in units of micron
  // We convert pixelLength & pixelHeight into mm units so
  // that they are comparable with distance (in units of
  //  millimeters)
  pixelLength_mm = pixelLength/1000.0;
  pixelHeight_mm = pixelHeight/1000.0;

  xMeasured = (xPixel-xCenter)*pixelLength_mm;
  yMeasured = (yPixel-yCenter)*pixelHeight_mm;

  bottom = distance+(xMeasured*cos_rotation+
                     yMeasured*sin_rotation)*sin_beta/*+
      (-xMeasured*sin_rotation+
       yMeasured*cos_rotation)*sin_alpha*cos_beta*/;

  // calculate the x y coordinates on the imaginary
  // detector using fancy math
  xPhysical = distance*((xMeasured*cos_rotation+
                         yMeasured*sin_rotation)*cos_beta
                        /*-(-xMeasured*sin_rotation+
                          yMeasured*cos_rotation)*sin_alpha*/)/bottom;

  yPhysical = distance*(-xMeasured*sin_rotation+
                        yMeasured*cos_rotation)/bottom;

  return sqrt(xPhysical*xPhysical + yPhysical*yPhysical);
}

double QxrdDetectorGeometry::getTwoTheta ( double  xCenter, double  yCenter, double  distance, double  xPixel, double  yPixel, double  pixelLength, double  pixelHeight, double  cos_beta, double  sin_beta, double  cos_rotation, double  sin_rotation  )

staticslot

cos_alpha

Definition at line 97 of file qxrddetectorgeometry.cpp.

Referenced by QxrdIntegratorCache::getTTH(), and QxrdCenterFinder::getTTH().

{
  double pixelLength_mm,pixelHeight_mm;
  double xMeasured,yMeasured;
  double bottom;
  double xPhysical,yPhysical;
  double twoTheta;

  // pixellength comes in in units of micron
  // We convert pixelLength & pixelHeight into mm units so
  // that they are comparable with distance (in units of
  //  millimeters)
  pixelLength_mm = pixelLength/1000.0;
  pixelHeight_mm = pixelHeight/1000.0;

  xMeasured = (xPixel-xCenter)*pixelLength_mm;
  yMeasured = (yPixel-yCenter)*pixelHeight_mm;

  bottom = distance+(xMeasured*cos_rotation+
                     yMeasured*sin_rotation)*sin_beta/*+
      (-xMeasured*sin_rotation+
       yMeasured*cos_rotation)*sin_alpha*cos_beta*/;

  // calculate the x y coordinates on the imaginary
  // detector using fancy math
  xPhysical = distance*((xMeasured*cos_rotation+
                         yMeasured*sin_rotation)*cos_beta
                        /*-(-xMeasured*sin_rotation+
                          yMeasured*cos_rotation)*sin_alpha*/)/bottom;

  yPhysical = distance*(-xMeasured*sin_rotation+
                        yMeasured*cos_rotation)/bottom;

  twoTheta = atan2(sqrt(xPhysical*xPhysical+
                         yPhysical*yPhysical),distance);
  // Convert to radians
  twoTheta = twoTheta * 180.0/M_PI;

  return twoTheta;
//
//  // explicitly convert chi to degrees
//  *chi=atan2(yPhysical,xPhysical)*180.0/M_PI;
//
//  // then add rotation to it so that chi always points
//  // to the right. Also, we have to multiply chi by -1
//  // because we have been defining our angles the inverse
//  // of the way they should be. There is a probably a
//  // better way to do this if I really thought through
//  // exactly how chi is defined. For the moment, through,
//  // this does exactly the right thing.
//  *chi = (*chi + rotation)*(-1);
//
//  // make sure that chi is b/n 0 and 360
//  *chi = mod(*chi, 360.0);
}

Here is the caller graph for this function:

void QxrdDetectorGeometry::getTwoThetaChi ( double  xCenter, double  yCenter, double  distance, double  xPixel, double  yPixel, double  pixelLength, double  pixelHeight, double  rotation, double  cos_beta, double  sin_beta, double  cos_alpha, double  sin_alpha, double  cos_rotation, double  sin_rotation, double *  twoTheta, double *  chi  )

staticslot

Definition at line 166 of file qxrddetectorgeometry.cpp.

References mod().

Referenced by getQChi().

{
  double pixelLength_mm,pixelHeight_mm;
  double xMeasured,yMeasured;
  double bottom;
  double xPhysical,yPhysical;

  // pixellength comes in in units of micron
  // We convert pixelLength & pixelHeight into mm units so
  // that they are comparable with distance (in units of
  //  millimeters)
  pixelLength_mm = pixelLength/1000.0;
  pixelHeight_mm = pixelHeight/1000.0;

  xMeasured = (xPixel-xCenter)*pixelLength_mm;
  yMeasured = (yPixel-yCenter)*pixelHeight_mm;

  bottom = distance+(xMeasured*cos_rotation+
                     yMeasured*sin_rotation)*sin_beta+
      (-xMeasured*sin_rotation+
       yMeasured*cos_rotation)*sin_alpha*cos_beta;

  // calculate the x y coordinates on the imaginary
  // detector using fancy math
  xPhysical = distance*((xMeasured*cos_rotation+
                         yMeasured*sin_rotation)*cos_beta
                        -(-xMeasured*sin_rotation+
                          yMeasured*cos_rotation)*sin_alpha)/bottom;

  yPhysical = distance*(-xMeasured*sin_rotation+
                        yMeasured*cos_rotation)*cos_alpha/bottom;

  *twoTheta = atan2(sqrt(xPhysical*xPhysical+
                         yPhysical*yPhysical),distance);
  // Convert to radians
  *twoTheta = *twoTheta * 180.0/M_PI;

  // explicitly convert chi to degrees
  *chi=atan2(yPhysical,xPhysical)*180.0/M_PI;

  // then add rotation to it so that chi always points
  // to the right. Also, we have to multiply chi by -1
  // because we have been defining our angles the inverse
  // of the way they should be. There is a probably a
  // better way to do this if I really thought through
  // exactly how chi is defined. For the moment, through,
  // this does exactly the right thing.
  *chi = (*chi + rotation*180.0/M_PI)*(-1);

  // make sure that chi is b/n 0 and 360
  *chi = mod(*chi, 360.0);
}

Here is the call graph for this function:

Here is the caller graph for this function:

void QxrdDetectorGeometry::getXY ( double  xCenter, double  yCenter, double  distance, double  energy, double  q, double  chi, double  pixelLength, double  pixelHeight, double  rotation, double  cos_beta, double  sin_beta, double  cos_alpha, double  sin_alpha, double  cos_rotation, double  sin_rotation, double *  xPixel, double *  yPixel  )

staticslot

Definition at line 257 of file qxrddetectorgeometry.cpp.

References mod().

Referenced by QxrdCenterFinder::getXY().

{
  double wavelength;
  double twoTheta;
  double xPhysical,yPhysical;
  double bottom;
  double pixelLength_mm,pixelHeight_mm;
  double xMeasured,yMeasured;

  double tan_chi;

  wavelength = 12398.4172/energy;

  // rotate chi back to point whatever way it is supposed
  // to point
  chi = chi*(-1) - rotation;

  chi = mod(chi, 360.0);

  // explicitly convert chi to radians
  chi*=M_PI/180.0;

  twoTheta = 2.0*asin(wavelength*q/(4.0*M_PI));

  tan_chi = tan(chi);
  xPhysical = fabs(distance*tan(twoTheta)/sqrt(1.0+
                                               tan_chi*tan_chi));

  // one must determine explicitly the sign of xPhysical
  // by inspecting the diagram. This is b/c at one point
  // we take a sqrt in our derivation.
  if (chi>M_PI/2.0 && chi<3.0*M_PI/2.0)
    xPhysical = -1.0*xPhysical;

  yPhysical=fabs(xPhysical*tan_chi);

  // set the sign of y explicitly
  if (chi > M_PI && chi < 2.0*M_PI)
    yPhysical = -1.0*fabs(yPhysical);

  // I should worry about cos_alpha being 0
  bottom = distance*cos_beta-xPhysical*sin_beta-
           cos_beta*(xPhysical*cos_beta+distance)/(
               (xPhysical*cos_alpha)/(yPhysical*sin_alpha)+1);

  xMeasured = (distance*xPhysical*cos_rotation-
               distance*xPhysical*cos_beta*sin_rotation/
               (xPhysical*cos_alpha/yPhysical+sin_alpha))/
      bottom;

  yMeasured = (distance*xPhysical*sin_rotation+
               distance*xPhysical*cos_beta*cos_rotation/
               (xPhysical*cos_alpha/yPhysical+sin_alpha))/
      bottom;

  // convert pixelLength & pixelHeight into mm units so
  // that they are comparable with distance (in units of
  // millimeters)
  pixelLength_mm = pixelLength/1000.0;
  pixelHeight_mm = pixelHeight/1000.0;

  *xPixel = xMeasured/pixelLength_mm + xCenter;
  *yPixel = yMeasured/pixelHeight_mm + yCenter;
}

Here is the call graph for this function:

Here is the caller graph for this function:

---

The documentation for this class was generated from the following files:

• qxrddetectorgeometry.h (Commit a65ccc9... : jennings : 2016-03-15 14:00:18 -0500)
• qxrddetectorgeometry.cpp (Commit a65ccc9... : jennings : 2016-03-15 14:00:18 -0500)