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

Description

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

Definition at line 118 of file IPeakShape.h.

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

Public Member Functions

 LorentzFisherPeakShape (double max_intensity, double radial_size, double kappa)
 
 ~LorentzFisherPeakShape () 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...
 
LorentzFisherPeakShapeclone () 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

◆ LorentzFisherPeakShape()

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

Definition at line 175 of file IPeakShape.cpp.

177  : m_max_intensity(max_intensity)
178  , m_radial_size(radial_size)
179  , m_kappa(kappa)
180 {
181 }

Referenced by clone().

◆ ~LorentzFisherPeakShape()

LorentzFisherPeakShape::~LorentzFisherPeakShape ( )
overridedefault

Member Function Documentation

◆ angularDisorder()

bool LorentzFisherPeakShape::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 134 of file IPeakShape.h.

134 { 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 LorentzFisherPeakShape::className ( ) const
inlinefinalvirtual

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

Implements INode.

Definition at line 124 of file IPeakShape.h.

124 { return "LorentzFisherPeakShape"; }

◆ clone()

LorentzFisherPeakShape * LorentzFisherPeakShape::clone ( ) const
overridevirtual

Implements IPeakShape.

Definition at line 185 of file IPeakShape.cpp.

186 {
188 }
LorentzFisherPeakShape(double max_intensity, double radial_size, double kappa)
Definition: IPeakShape.cpp:175

References LorentzFisherPeakShape(), 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> LorentzFisherPeakShape::parDefs ( ) const
inlinefinalvirtual

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

Reimplemented from INode.

Definition at line 125 of file IPeakShape.h.

126  {
127  return {{"MaxIntensity", "", "maximum intensity", 0, +INF, -1},
128  {"DomainSize", "nm", "domain size", 0, +INF, -1},
129  {"Kappa", "", "?", 0, +INF, -1}};
130  }

References INF.

◆ peakDistribution()

double LorentzFisherPeakShape::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 190 of file IPeakShape.cpp.

191 {
192  const double q_r = q.mag();
193  const double q_lat_r = q_lattice_point.mag();
194  const double dq2 = (q_r - q_lat_r) * (q_r - q_lat_r);
195  if (q_lat_r == 0.0)
196  return m_max_intensity * Cauchy3D(dq2, m_radial_size);
197  const double radial_part = m_radial_size / (1.0 + dq2 * m_radial_size * m_radial_size) / M_PI;
198  double angular_part = 1.0;
199  if (q_r * q_lat_r > 0.0) {
200  const double dot_norm = q.dot(q_lattice_point) / q_r / q_lat_r;
201  angular_part = FisherDistribution(dot_norm, m_kappa) / (q_r * q_r);
202  }
203  return m_max_intensity * radial_part * angular_part;
204 }
#define M_PI
Definition: Constants.h:44

References m_kappa, m_max_intensity, M_PI, and m_radial_size.

◆ 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 LorentzFisherPeakShape::m_kappa
private

Definition at line 139 of file IPeakShape.h.

Referenced by clone(), and peakDistribution().

◆ m_max_intensity

double LorentzFisherPeakShape::m_max_intensity
private

Definition at line 137 of file IPeakShape.h.

Referenced by clone(), and peakDistribution().

◆ m_P

◆ m_radial_size

double LorentzFisherPeakShape::m_radial_size
private

Definition at line 138 of file IPeakShape.h.

Referenced by clone(), and peakDistribution().


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