BornAgain  1.19.79
Open-source research software to simulate and fit neutron and x-ray reflectometry and grazing-incidence small-angle scattering
Interference1DLattice Class Reference

Description

Interference function of a 1D lattice.

Definition at line 25 of file Interference1DLattice.h.

Inheritance diagram for Interference1DLattice:
[legend]
Collaboration diagram for Interference1DLattice:
[legend]

Public Member Functions

 Interference1DLattice (double length, double xi)
 Constructor of interference function of one-dimensional lattice. More...
 
 ~Interference1DLattice () override
 
void checkNodeArgs () const
 Raises exception if a parameter value is invalid. More...
 
std::string className () const final
 Returns the class name, to be hard-coded in each leaf class that inherits from INode. More...
 
Interference1DLatticeclone () const override
 
double DWfactor (R3 q) const
 structureFactors the Debye-Waller factor for a given wavevector transfer More...
 
double length () const
 
std::vector< const INode * > nodeChildren () const override
 Returns all children. More...
 
std::vector< const INode * > nodeOffspring () const
 Returns all descendants. More...
 
std::vector< ParaMetaparDefs () const final
 Returns the parameter definitions, to be hard-coded in each leaf class. More...
 
virtual double particleDensity () const
 If defined by this interference function's parameters, returns the particle density (per area). Otherwise, returns zero or a user-defined value. More...
 
double positionVariance () const
 Returns the position variance. More...
 
void setDecayFunction (const IProfile1D &decay)
 Sets one-dimensional decay function. More...
 
void setPositionVariance (double var)
 Sets the variance of the position for the calculation of the DW factor It is defined as the variance in each relevant dimension. More...
 
virtual double structureFactor (R3 q, double outer_iff=1.0) const
 The interference function for a given wavevector transfer. More...
 
virtual bool supportsMultilayer () const
 Indicates if this interference function can be used with a sample (DWBA mode) More...
 
virtual void transferToCPP ()
 Used for Python overriding of clone (see swig/tweaks.py) More...
 
double xi () const
 

Protected Member Functions

double iff_no_inner (R3 q, double outer_iff) const
 Calculates the structure factor in the absence of extra inner structure. More...
 

Protected Attributes

std::vector< double > m_P
 
double m_position_var
 

Private Member Functions

double iff_without_dw (R3 q) const override
 Calculates the structure factor without Debye-Waller factor. More...
 

Private Attributes

std::unique_ptr< IProfile1Dm_decay
 
double m_length
 
int m_na
 determines the number of reciprocal lattice points to use More...
 
double m_xi
 

Constructor & Destructor Documentation

◆ Interference1DLattice()

Interference1DLattice::Interference1DLattice ( double  length,
double  xi 
)

Constructor of interference function of one-dimensional lattice.

Parameters
lengthlattice constant in nanometers
xirotation of lattice with respect to x-axis in radians

Definition at line 34 of file Interference1DLattice.cpp.

35  : IInterference(0)
36  , m_length(length)
37  , m_xi(xi)
38  , m_na{0}
39 {
41 }
IInterference(const std::vector< double > &PValues)
int m_na
determines the number of reciprocal lattice points to use
void check(const std::string &name, double value) const
Throws if value is outside limits. Parameter 'name' is for exception message.
Definition: RealLimits.cpp:170
static RealLimits nonnegative()
Creates an object which can have only positive values with 0. included.
Definition: RealLimits.cpp:124

References RealLimits::check(), m_length, and RealLimits::nonnegative().

Referenced by clone().

Here is the call graph for this function:

◆ ~Interference1DLattice()

Interference1DLattice::~Interference1DLattice ( )
overridedefault

Member Function Documentation

◆ checkNodeArgs()

void INode::checkNodeArgs ( ) const
inherited

Raises exception if a parameter value is invalid.

Definition at line 27 of file INode.cpp.

28 {
29  size_t nP = m_P.size();
30  if (parDefs().size() != nP) {
31  std::cerr << "BUG in class " << className() << std::endl;
32  std::cerr << "#m_P = " << nP << std::endl;
33  std::cerr << "#PDf = " << parDefs().size() << std::endl;
34  for (const ParaMeta& pm : parDefs())
35  std::cerr << " PDf: " << pm.name << std::endl;
36  ASSERT(0);
37  }
38  ASSERT(parDefs().size() == nP);
39  for (size_t i = 0; i < nP; ++i) {
40  const ParaMeta pm = parDefs()[i];
41 
43  if (pm.vMin == -INF) {
44  ASSERT(pm.vMax == +INF);
45  // nothing to do
46  } else if (pm.vMax == +INF) {
47  ASSERT(pm.vMin == 0);
48  limits = RealLimits::nonnegative();
49  } else {
50  limits = RealLimits::limited(pm.vMin, pm.vMax);
51  }
52  limits.check(pm.name, m_P[i]);
53  }
54 }
#define ASSERT(condition)
Definition: Assert.h:45
const double INF
Definition: INode.h:26
virtual std::vector< ParaMeta > parDefs() const
Returns the parameter definitions, to be hard-coded in each leaf class.
Definition: INode.h:51
std::vector< double > m_P
Definition: INode.h:63
virtual std::string className() const =0
Returns the class name, to be hard-coded in each leaf class that inherits from INode.
Limits for a real fit parameter.
Definition: RealLimits.h:24
static RealLimits limitless()
Creates an object without bounds (default)
Definition: RealLimits.cpp:139
static RealLimits limited(double left_bound_value, double right_bound_value)
Creates an object bounded from the left and right.
Definition: RealLimits.cpp:134
Metadata of one model parameter.
Definition: INode.h:29
double vMin
Definition: INode.h:33
double vMax
Definition: INode.h:34
std::string name
Definition: INode.h:30

References ASSERT, RealLimits::check(), INode::className(), INF, RealLimits::limited(), RealLimits::limitless(), INode::m_P, ParaMeta::name, RealLimits::nonnegative(), INode::parDefs(), ParaMeta::vMax, and ParaMeta::vMin.

Referenced by BarGauss::BarGauss(), BarLorentz::BarLorentz(), Bipyramid4::Bipyramid4(), Box::Box(), CantellatedCube::CantellatedCube(), Cone::Cone(), ConstantBackground::ConstantBackground(), CosineRippleBox::CosineRippleBox(), CosineRippleGauss::CosineRippleGauss(), CosineRippleLorentz::CosineRippleLorentz(), Cylinder::Cylinder(), DistributionCosine::DistributionCosine(), DistributionGate::DistributionGate(), DistributionGaussian::DistributionGaussian(), DistributionLogNormal::DistributionLogNormal(), DistributionLorentz::DistributionLorentz(), DistributionTrapezoid::DistributionTrapezoid(), Dodecahedron::Dodecahedron(), EllipsoidalCylinder::EllipsoidalCylinder(), FootprintGauss::FootprintGauss(), FootprintSquare::FootprintSquare(), FuzzySphere::FuzzySphere(), GaussSphere::GaussSphere(), HemiEllipsoid::HemiEllipsoid(), HollowSphere::HollowSphere(), HorizontalCylinder::HorizontalCylinder(), Icosahedron::Icosahedron(), LongBoxGauss::LongBoxGauss(), LongBoxLorentz::LongBoxLorentz(), PlatonicOctahedron::PlatonicOctahedron(), PlatonicTetrahedron::PlatonicTetrahedron(), Prism3::Prism3(), Prism6::Prism6(), Profile1DCauchy::Profile1DCauchy(), Profile1DCosine::Profile1DCosine(), Profile1DGate::Profile1DGate(), Profile1DGauss::Profile1DGauss(), Profile1DTriangle::Profile1DTriangle(), Profile1DVoigt::Profile1DVoigt(), Profile2DCauchy::Profile2DCauchy(), Profile2DCone::Profile2DCone(), Profile2DGate::Profile2DGate(), Profile2DGauss::Profile2DGauss(), Profile2DVoigt::Profile2DVoigt(), Pyramid2::Pyramid2(), Pyramid3::Pyramid3(), Pyramid4::Pyramid4(), Pyramid6::Pyramid6(), RotationEuler::RotationEuler(), RotationX::RotationX(), RotationY::RotationY(), RotationZ::RotationZ(), SawtoothRippleBox::SawtoothRippleBox(), SawtoothRippleGauss::SawtoothRippleGauss(), SawtoothRippleLorentz::SawtoothRippleLorentz(), Sphere::Sphere(), Spheroid::Spheroid(), TruncatedCube::TruncatedCube(), TruncatedSphere::TruncatedSphere(), and TruncatedSpheroid::TruncatedSpheroid().

Here is the call graph for this function:

◆ className()

std::string Interference1DLattice::className ( ) const
inlinefinalvirtual

Returns the class name, to be hard-coded in each leaf class that inherits from INode.

Implements INode.

Definition at line 31 of file Interference1DLattice.h.

31 { return "Interference1DLattice"; }

◆ clone()

Interference1DLattice * Interference1DLattice::clone ( ) const
overridevirtual

Implements IInterference.

Definition at line 45 of file Interference1DLattice.cpp.

46 {
47  auto* result = new Interference1DLattice(m_length, m_xi);
48  result->setPositionVariance(m_position_var);
49  result->m_na = m_na;
50  if (m_decay)
51  result->setDecayFunction(*m_decay);
52  return result;
53 }
double m_position_var
Definition: IInterference.h:53
std::unique_ptr< IProfile1D > m_decay
Interference1DLattice(double length, double xi)
Constructor of interference function of one-dimensional lattice.

References Interference1DLattice(), m_decay, m_length, m_na, IInterference::m_position_var, and m_xi.

Here is the call graph for this function:

◆ DWfactor()

double IInterference::DWfactor ( R3  q) const
inherited

structureFactors the Debye-Waller factor for a given wavevector transfer

Definition at line 48 of file IInterference.cpp.

49 {
50  // remove z component for two-dimensional interference functions:
51  if (supportsMultilayer())
52  q.setZ(0.0);
53  return std::exp(-q.mag2() * m_position_var);
54 }
virtual bool supportsMultilayer() const
Indicates if this interference function can be used with a sample (DWBA mode)
Definition: IInterference.h:47

References IInterference::m_position_var, and IInterference::supportsMultilayer().

Referenced by IInterference::iff_no_inner(), and Interference2DSuperLattice::interferenceForXi().

Here is the call graph for this function:

◆ iff_no_inner()

double IInterference::iff_no_inner ( R3  q,
double  outer_iff 
) const
protectedinherited

Calculates the structure factor in the absence of extra inner structure.

Definition at line 56 of file IInterference.cpp.

57 {
58  return DWfactor(q) * (iff_without_dw(q) * outer_iff - 1.0) + 1.0;
59 }
double DWfactor(R3 q) const
structureFactors the Debye-Waller factor for a given wavevector transfer
virtual double iff_without_dw(R3 q) const =0
Calculates the structure factor without Debye-Waller factor.

References IInterference::DWfactor(), and IInterference::iff_without_dw().

Referenced by IInterference::structureFactor().

Here is the call graph for this function:

◆ iff_without_dw()

double Interference1DLattice::iff_without_dw ( R3  q) const
overrideprivatevirtual

Calculates the structure factor without Debye-Waller factor.

Implements IInterference.

Definition at line 71 of file Interference1DLattice.cpp.

72 {
73  ASSERT(m_decay);
74  double result = 0.0;
75  double qxr = q.x();
76  double qyr = q.y();
77  double qx_frac;
78  double xi = m_xi;
79  double a = m_length;
80  double a_rec = M_TWOPI / a;
81 
82  // rotate the q vector to xi angle
83  // so that qx_prime is along the a axis of lattice
84  double qx_prime = qxr * std::cos(xi) + qyr * std::sin(xi);
85 
86  // calculate reciprocal vector fraction
87  int qa_int = static_cast<int>(qx_prime / a_rec);
88  qx_frac = qx_prime - qa_int * a_rec;
89 
90  for (int i = -m_na; i < m_na + 1; ++i) {
91  double qx = qx_frac + i * a_rec;
92  result += m_decay->decayFT(qx);
93  }
94  return result / a;
95 }
#define M_TWOPI
Definition: Constants.h:54

References ASSERT, m_decay, m_length, m_na, M_TWOPI, m_xi, and xi().

Here is the call graph for this function:

◆ length()

double Interference1DLattice::length ( ) const
inline

Definition at line 40 of file Interference1DLattice.h.

40 { return m_length; }

References m_length.

◆ nodeChildren()

std::vector< const INode * > Interference1DLattice::nodeChildren ( ) const
overridevirtual

Returns all children.

Reimplemented from INode.

Definition at line 66 of file Interference1DLattice.cpp.

67 {
68  return std::vector<const INode*>() << m_decay;
69 }

References m_decay.

◆ nodeOffspring()

std::vector< const INode * > INode::nodeOffspring ( ) const
inherited

Returns all descendants.

Definition at line 61 of file INode.cpp.

62 {
63  std::vector<const INode*> result;
64  result.push_back(this);
65  for (const auto* child : nodeChildren()) {
66  for (const auto* p : child->nodeOffspring())
67  result.push_back(p);
68  }
69  return result;
70 }
virtual std::vector< const INode * > nodeChildren() const
Returns all children.
Definition: INode.cpp:56

References INode::nodeChildren().

Here is the call graph for this function:

◆ parDefs()

std::vector<ParaMeta> Interference1DLattice::parDefs ( ) const
inlinefinalvirtual

Returns the parameter definitions, to be hard-coded in each leaf class.

Reimplemented from INode.

Definition at line 32 of file Interference1DLattice.h.

33  {
34  return {{"Length", "nm", "lattice constant", 0, +INF, 0},
35  {"Xi", "rad", "angle x axis and lattice vector??", -INF, +INF, 0}};
36  }

References INF.

◆ particleDensity()

virtual double IInterference::particleDensity ( ) const
inlinevirtualinherited

If defined by this interference function's parameters, returns the particle density (per area). Otherwise, returns zero or a user-defined value.

Reimplemented in InterferenceHardDisk, InterferenceFinite2DLattice, Interference2DParaCrystal, and Interference2DLattice.

Definition at line 44 of file IInterference.h.

44 { return 0.0; }

◆ positionVariance()

double IInterference::positionVariance ( ) const
inlineinherited

Returns the position variance.

Definition at line 40 of file IInterference.h.

40 { return m_position_var; }

References IInterference::m_position_var.

◆ setDecayFunction()

void Interference1DLattice::setDecayFunction ( const IProfile1D decay)

Sets one-dimensional decay function.

Parameters
decayone-dimensional decay function in reciprocal space

Definition at line 57 of file Interference1DLattice.cpp.

58 {
59  m_decay.reset(decay.clone());
60  double decay_length = m_decay->decayLength();
61  double qa_max = m_length * nmax / decay_length / M_TWOPI;
62  m_na = static_cast<int>(std::lround(std::abs(qa_max) + 0.5));
63  m_na = std::max(m_na, min_points);
64 }
IProfile1D * clone() const override=0

References IProfile1D::clone(), m_decay, m_length, m_na, and M_TWOPI.

Referenced by ExemplarySamples::createLattice1D().

Here is the call graph for this function:

◆ setPositionVariance()

void IInterference::setPositionVariance ( double  var)
inherited

Sets the variance of the position for the calculation of the DW factor It is defined as the variance in each relevant dimension.

Definition at line 40 of file IInterference.cpp.

41 {
42  if (var < 0.0)
43  throw std::runtime_error("IInterference::setPositionVariance: "
44  "variance should be positive.");
45  m_position_var = var;
46 }

References IInterference::m_position_var.

Referenced by ExemplarySamples::createFiniteSquareLattice2D(), and ExemplarySamples::createSuperLattice().

◆ structureFactor()

double IInterference::structureFactor ( R3  q,
double  outer_iff = 1.0 
) const
virtualinherited

The interference function for a given wavevector transfer.

Reimplemented in Interference2DSuperLattice.

Definition at line 35 of file IInterference.cpp.

36 {
37  return iff_no_inner(q, outer_iff);
38 }
double iff_no_inner(R3 q, double outer_iff) const
Calculates the structure factor in the absence of extra inner structure.

References IInterference::iff_no_inner().

Here is the call graph for this function:

◆ supportsMultilayer()

virtual bool IInterference::supportsMultilayer ( ) const
inlinevirtualinherited

Indicates if this interference function can be used with a sample (DWBA mode)

Reimplemented in InterferenceFinite3DLattice, and Interference3DLattice.

Definition at line 47 of file IInterference.h.

47 { return true; }

Referenced by IInterference::DWfactor().

◆ transferToCPP()

virtual void ICloneable::transferToCPP ( )
inlinevirtualinherited

Used for Python overriding of clone (see swig/tweaks.py)

Definition at line 32 of file ICloneable.h.

◆ xi()

double Interference1DLattice::xi ( ) const
inline

Definition at line 41 of file Interference1DLattice.h.

41 { return m_xi; }

References m_xi.

Referenced by iff_without_dw().

Member Data Documentation

◆ m_decay

std::unique_ptr<IProfile1D> Interference1DLattice::m_decay
private

Definition at line 50 of file Interference1DLattice.h.

Referenced by clone(), iff_without_dw(), nodeChildren(), and setDecayFunction().

◆ m_length

double Interference1DLattice::m_length
private

◆ m_na

int Interference1DLattice::m_na
private

determines the number of reciprocal lattice points to use

Definition at line 51 of file Interference1DLattice.h.

Referenced by clone(), iff_without_dw(), and setDecayFunction().

◆ m_P

◆ m_position_var

◆ m_xi

double Interference1DLattice::m_xi
private

Definition at line 49 of file Interference1DLattice.h.

Referenced by clone(), iff_without_dw(), and xi().


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