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

Description

A peak shape that is Gaussian in the radial direction and uses the Mises-Fisher distribution in the angular direction.

Definition at line 91 of file IPeakShape.h.

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

Public Member Functions

 GaussFisherPeakShape (double max_intensity, double radial_size, double kappa)
 
 ~GaussFisherPeakShape () override
 
bool angularDisorder () const override
 Indicates if the peak shape encodes angular disorder, in which case all peaks in a spherical shell are needed. More...
 
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...
 
GaussFisherPeakShapeclone () const override
 
virtual std::vector< const INode * > nodeChildren () const
 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...
 
double peakDistribution (R3 q, R3 q_lattice_point) const override
 Peak shape at q from a reciprocal lattice point at q_lattice_point. More...
 
virtual void transferToCPP ()
 Used for Python overriding of clone (see swig/tweaks.py) More...
 

Protected Attributes

std::vector< double > m_P
 

Private Attributes

double m_kappa
 
double m_max_intensity
 
double m_radial_size
 

Constructor & Destructor Documentation

◆ GaussFisherPeakShape()

GaussFisherPeakShape::GaussFisherPeakShape ( double  max_intensity,
double  radial_size,
double  kappa 
)

Definition at line 140 of file IPeakShape.cpp.

141  : m_max_intensity(max_intensity)
142  , m_radial_size(radial_size)
143  , m_kappa(kappa)
144 {
145 }

Referenced by clone().

◆ ~GaussFisherPeakShape()

GaussFisherPeakShape::~GaussFisherPeakShape ( )
overridedefault

Member Function Documentation

◆ angularDisorder()

bool GaussFisherPeakShape::angularDisorder ( ) const
inlineoverridevirtual

Indicates if the peak shape encodes angular disorder, in which case all peaks in a spherical shell are needed.

Reimplemented from IPeakShape.

Definition at line 107 of file IPeakShape.h.

107 { return true; }

◆ 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
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
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 GaussFisherPeakShape::className ( ) const
inlinefinalvirtual

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

Implements INode.

Definition at line 97 of file IPeakShape.h.

97 { return "GaussFisherPeakShape"; }

◆ clone()

GaussFisherPeakShape * GaussFisherPeakShape::clone ( ) const
overridevirtual

Implements IPeakShape.

Definition at line 149 of file IPeakShape.cpp.

150 {
152 }
GaussFisherPeakShape(double max_intensity, double radial_size, double kappa)
Definition: IPeakShape.cpp:140

References GaussFisherPeakShape(), m_kappa, m_max_intensity, and m_radial_size.

Here is the call graph for this function:

◆ nodeChildren()

◆ 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> GaussFisherPeakShape::parDefs ( ) const
inlinefinalvirtual

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

Reimplemented from INode.

Definition at line 98 of file IPeakShape.h.

99  {
100  return {{"MaxIntensity", "", "maximum intensity", 0, +INF, -1},
101  {"DomainSize", "nm", "domain size", 0, +INF, -1},
102  {"Kappa", "", "?", 0, +INF, -1}};
103  }

References INF.

◆ peakDistribution()

double GaussFisherPeakShape::peakDistribution ( R3  q,
R3  q_lattice_point 
) const
overridevirtual

Peak shape at q from a reciprocal lattice point at q_lattice_point.

Implements IPeakShape.

Definition at line 154 of file IPeakShape.cpp.

155 {
156  const double q_r = q.mag();
157  const double q_lat_r = q_lattice_point.mag();
158  const double dq2 = (q_r - q_lat_r) * (q_r - q_lat_r);
159  if (q_lat_r == 0.0)
160  return m_max_intensity * Gauss3D(dq2, m_radial_size);
161  const double norm_factor = m_radial_size / std::sqrt(M_TWOPI);
162  const double radial_part = norm_factor * std::exp(-dq2 * m_radial_size * m_radial_size / 2.0);
163  double angular_part = 1.0;
164  if (q_r * q_lat_r > 0.0) {
165  const double dot_norm = q.dot(q_lattice_point) / q_r / q_lat_r;
166  angular_part = FisherDistribution(dot_norm, m_kappa) / (q_r * q_r);
167  }
168  return m_max_intensity * radial_part * angular_part;
169 }
#define M_TWOPI
Definition: Constants.h:54

References m_kappa, m_max_intensity, m_radial_size, and M_TWOPI.

◆ transferToCPP()

virtual void ICloneable::transferToCPP ( )
inlinevirtualinherited

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

Definition at line 32 of file ICloneable.h.

Member Data Documentation

◆ m_kappa

double GaussFisherPeakShape::m_kappa
private

Definition at line 112 of file IPeakShape.h.

Referenced by clone(), and peakDistribution().

◆ m_max_intensity

double GaussFisherPeakShape::m_max_intensity
private

Definition at line 110 of file IPeakShape.h.

Referenced by clone(), and peakDistribution().

◆ m_P

◆ m_radial_size

double GaussFisherPeakShape::m_radial_size
private

Definition at line 111 of file IPeakShape.h.

Referenced by clone(), and peakDistribution().


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