BornAgain  1.18.0
Simulate and fit neutron and x-ray scattering at grazing incidence
OffSpecularConverter Class Reference
Inheritance diagram for OffSpecularConverter:
Collaboration diagram for OffSpecularConverter:

Public Member Functions

 OffSpecularConverter (const IDetector2D &detector, const Beam &beam, const IAxis &alpha_axis)
 
 ~OffSpecularConverter () final
 
OffSpecularConverterclone () const final
 
Axes::Units defaultUnits () const final
 
virtual size_t dimension () const final
 
double calculateMin (size_t i_axis, Axes::Units units_type) const final
 
double calculateMax (size_t i_axis, Axes::Units units_type) const final
 
size_t axisSize (size_t i_axis) const final
 
std::vector< Axes::UnitsavailableUnits () const override
 
std::unique_ptr< IAxiscreateConvertedAxis (size_t i_axis, Axes::Units units) const final
 
std::string axisName (size_t i_axis, Axes::Units units_type=Axes::Units::DEFAULT) const
 
Axes::Units substituteDefaultUnits (Axes::Units units) const
 
virtual std::unique_ptr< OutputData< double > > createConvertedData (const OutputData< double > &data, Axes::Units units) const
 
virtual void transferToCPP ()
 

Protected Member Functions

void addDetectorAxis (const IDetector &detector, size_t i_axis)
 
void addAxisData (std::string name, double min, double max, Axes::Units default_units, size_t nbins)
 
void checkIndex (size_t i_axis) const
 
void throwUnitsError (std::string method, std::vector< Axes::Units > available) const
 

Protected Attributes

std::vector< AxisDatam_axis_data_table
 
double m_wavelength
 
double m_alpha_i
 
double m_phi_i
 

Private Member Functions

 OffSpecularConverter (const OffSpecularConverter &other)
 
double calculateValue (size_t i_axis, Axes::Units units_type, double value) const final
 
std::vector< std::map< Axes::Units, std::string > > createNameMaps () const final
 
void addDetectorYAxis (const IDetector2D &detector)
 

Detailed Description

IUnitConverter class that handles the unit translations for off-specular simulations with a spherical detector Its default units are radians for both axes.

Definition at line 127 of file SimpleUnitConverters.h.

Constructor & Destructor Documentation

◆ OffSpecularConverter() [1/2]

OffSpecularConverter::OffSpecularConverter ( const IDetector2D detector,
const Beam beam,
const IAxis alpha_axis 
)

Definition at line 311 of file SimpleUnitConverters.cpp.

313  : UnitConverterSimple(beam)
314 {
315  if (detector.dimension() != 2)
316  throw std::runtime_error("Error in OffSpecularConverter constructor: "
317  "detector has wrong dimension: "
318  + std::to_string(static_cast<int>(detector.dimension())));
319  addAxisData(axisName(0), alpha_axis.getMin(), alpha_axis.getMax(), defaultUnits(),
320  alpha_axis.size());
321  addDetectorYAxis(detector);
322 }
virtual double getMin() const =0
Returns value of first point of axis.
virtual size_t size() const =0
retrieve the number of bins
virtual double getMax() const =0
Returns value of last point of axis.
size_t dimension() const
Returns actual dimensionality of the detector (number of defined axes)
Definition: IDetector.cpp:44
std::string axisName(size_t i_axis, Axes::Units units_type=Axes::Units::DEFAULT) const
Axes::Units defaultUnits() const final
void addDetectorYAxis(const IDetector2D &detector)
void addAxisData(std::string name, double min, double max, Axes::Units default_units, size_t nbins)
UnitConverterSimple(const Beam &beam)

References UnitConverterSimple::addAxisData(), addDetectorYAxis(), IUnitConverter::axisName(), defaultUnits(), IDetector::dimension(), IAxis::getMax(), IAxis::getMin(), and IAxis::size().

Referenced by clone().

Here is the call graph for this function:

◆ ~OffSpecularConverter()

OffSpecularConverter::~OffSpecularConverter ( )
finaldefault

◆ OffSpecularConverter() [2/2]

OffSpecularConverter::OffSpecularConverter ( const OffSpecularConverter other)
private

Definition at line 336 of file SimpleUnitConverters.cpp.

337  : UnitConverterSimple(other)
338 {
339 }

Member Function Documentation

◆ clone()

OffSpecularConverter * OffSpecularConverter::clone ( ) const
finalvirtual

Implements IUnitConverter.

Definition at line 326 of file SimpleUnitConverters.cpp.

327 {
328  return new OffSpecularConverter(*this);
329 }
OffSpecularConverter(const IDetector2D &detector, const Beam &beam, const IAxis &alpha_axis)

References OffSpecularConverter().

Here is the call graph for this function:

◆ defaultUnits()

Axes::Units OffSpecularConverter::defaultUnits ( ) const
finalvirtual

Implements IUnitConverter.

Definition at line 331 of file SimpleUnitConverters.cpp.

332 {
333  return Axes::Units::DEGREES;
334 }

Referenced by addDetectorYAxis(), and OffSpecularConverter().

◆ calculateValue()

double OffSpecularConverter::calculateValue ( size_t  i_axis,
Axes::Units  units_type,
double  value 
) const
finalprivatevirtual

Implements UnitConverterSimple.

Definition at line 341 of file SimpleUnitConverters.cpp.

342 {
343  switch (units_type) {
344  case Axes::Units::RADIANS:
345  return value;
346  case Axes::Units::DEGREES:
347  return Units::rad2deg(value);
348  default:
349  throwUnitsError("OffSpecularConverter::calculateValue", availableUnits());
350  }
351 }
void throwUnitsError(std::string method, std::vector< Axes::Units > available) const
std::vector< Axes::Units > availableUnits() const override
Returns the list of all available units.
double rad2deg(double angle)
Definition: Units.h:43

References UnitConverterSimple::availableUnits(), Units::rad2deg(), and IUnitConverter::throwUnitsError().

Here is the call graph for this function:

◆ createNameMaps()

std::vector< std::map< Axes::Units, std::string > > OffSpecularConverter::createNameMaps ( ) const
finalprivatevirtual

Implements IUnitConverter.

Definition at line 353 of file SimpleUnitConverters.cpp.

354 {
355  std::vector<std::map<Axes::Units, std::string>> result;
356  result.push_back(AxisNames::InitOffSpecAxis0());
357  result.push_back(AxisNames::InitOffSpecAxis1());
358  return result;
359 }
std::map< Axes::Units, std::string > InitOffSpecAxis1()
Definition: AxisNames.cpp:74
std::map< Axes::Units, std::string > InitOffSpecAxis0()
Definition: AxisNames.cpp:66

References AxisNames::InitOffSpecAxis0(), and AxisNames::InitOffSpecAxis1().

Here is the call graph for this function:

◆ addDetectorYAxis()

void OffSpecularConverter::addDetectorYAxis ( const IDetector2D detector)
private

Definition at line 361 of file SimpleUnitConverters.cpp.

362 {
363  const auto& axis = detector.getAxis(1);
364  const auto* p_roi = detector.regionOfInterest();
365  const auto& axis_name = axisName(1);
366  std::unique_ptr<IAxis> P_new_axis;
367  if (p_roi) {
368  P_new_axis = p_roi->clipAxisToRoi(1, axis);
369  } else {
370  P_new_axis.reset(axis.clone());
371  }
372  if (!P_new_axis)
373  throw std::runtime_error("Error in OffSpecularConverter::addDetectorYAxis: "
374  "could not retrieve the y-axis of the detector");
375  if (const auto* P_rect_det = dynamic_cast<const RectangularDetector*>(&detector)) {
376  std::unique_ptr<RectangularPixel> P_det_pixel(P_rect_det->regionOfInterestPixel());
377  const auto k00 = P_det_pixel->getPosition(0.0, 0.0);
378  const auto k01 = P_det_pixel->getPosition(0.0, 1.0);
379  const double alpha_f_min = M_PI_2 - k00.theta();
380  const double alpha_f_max = M_PI_2 - k01.theta();
381  addAxisData(axis_name, alpha_f_min, alpha_f_max, defaultUnits(), P_new_axis->size());
382  } else if (dynamic_cast<const SphericalDetector*>(&detector)) {
383  const double alpha_f_min = P_new_axis->getMin();
384  const double alpha_f_max = P_new_axis->getMax();
385  addAxisData(axis_name, alpha_f_min, alpha_f_max, defaultUnits(), P_new_axis->size());
386  } else {
387  throw std::runtime_error("Error in OffSpecularConverter::addDetectorYAxis: "
388  "wrong detector type");
389  }
390 }
#define M_PI_2
Definition: MathConstants.h:40
const RegionOfInterest * regionOfInterest() const override
Returns region of interest if exists.
Definition: IDetector2D.cpp:43
const IAxis & getAxis(size_t index) const
Definition: IDetector.cpp:54
A flat rectangular detector with axes and resolution function.
A spherical detector with axes and resolution function.

References UnitConverterSimple::addAxisData(), IUnitConverter::axisName(), defaultUnits(), IDetector::getAxis(), M_PI_2, and IDetector2D::regionOfInterest().

Referenced by OffSpecularConverter().

Here is the call graph for this function:

◆ dimension()

size_t UnitConverterSimple::dimension ( ) const
finalvirtualinherited

Implements IUnitConverter.

Definition at line 41 of file SimpleUnitConverters.cpp.

42 {
43  return m_axis_data_table.size();
44 }
std::vector< AxisData > m_axis_data_table

References UnitConverterSimple::m_axis_data_table.

◆ calculateMin()

double UnitConverterSimple::calculateMin ( size_t  i_axis,
Axes::Units  units_type 
) const
finalvirtualinherited

Implements IUnitConverter.

Definition at line 53 of file SimpleUnitConverters.cpp.

54 {
55  checkIndex(i_axis);
56  units_type = substituteDefaultUnits(units_type);
57  const auto& axis_data = m_axis_data_table[i_axis];
58  if (units_type == Axes::Units::NBINS) {
59  return 0.0;
60  }
61  return calculateValue(i_axis, units_type, axis_data.min);
62 }
Axes::Units substituteDefaultUnits(Axes::Units units) const
void checkIndex(size_t i_axis) const
virtual double calculateValue(size_t i_axis, Axes::Units units_type, double value) const =0

References UnitConverterSimple::calculateValue(), IUnitConverter::checkIndex(), UnitConverterSimple::m_axis_data_table, and IUnitConverter::substituteDefaultUnits().

Referenced by UnitConverterSimple::createConvertedAxis().

Here is the call graph for this function:

◆ calculateMax()

double UnitConverterSimple::calculateMax ( size_t  i_axis,
Axes::Units  units_type 
) const
finalvirtualinherited

Implements IUnitConverter.

Definition at line 64 of file SimpleUnitConverters.cpp.

65 {
66  checkIndex(i_axis);
67  units_type = substituteDefaultUnits(units_type);
68  const auto& axis_data = m_axis_data_table[i_axis];
69  if (units_type == Axes::Units::NBINS) {
70  return static_cast<double>(axis_data.nbins);
71  }
72  return calculateValue(i_axis, units_type, axis_data.max);
73 }

References UnitConverterSimple::calculateValue(), IUnitConverter::checkIndex(), UnitConverterSimple::m_axis_data_table, and IUnitConverter::substituteDefaultUnits().

Referenced by UnitConverterSimple::createConvertedAxis().

Here is the call graph for this function:

◆ axisSize()

size_t UnitConverterSimple::axisSize ( size_t  i_axis) const
finalvirtualinherited

Implements IUnitConverter.

Definition at line 75 of file SimpleUnitConverters.cpp.

76 {
77  checkIndex(i_axis);
78  return m_axis_data_table[i_axis].nbins;
79 }

References IUnitConverter::checkIndex(), and UnitConverterSimple::m_axis_data_table.

Referenced by UnitConverterSimple::createConvertedAxis().

Here is the call graph for this function:

◆ availableUnits()

std::vector< Axes::Units > UnitConverterSimple::availableUnits ( ) const
overridevirtualinherited

Returns the list of all available units.

Implements IUnitConverter.

Reimplemented in DepthProbeConverter, RectangularConverter, and SphericalConverter.

Definition at line 81 of file SimpleUnitConverters.cpp.

82 {
83  return {Axes::Units::NBINS, Axes::Units::RADIANS, Axes::Units::DEGREES};
84 }

Referenced by SphericalConverter::availableUnits(), RectangularConverter::availableUnits(), DepthProbeConverter::availableUnits(), and calculateValue().

◆ createConvertedAxis()

std::unique_ptr< IAxis > UnitConverterSimple::createConvertedAxis ( size_t  i_axis,
Axes::Units  units 
) const
finalvirtualinherited

Implements IUnitConverter.

Definition at line 86 of file SimpleUnitConverters.cpp.

88 {
89  const double min = calculateMin(i_axis, units);
90  const double max = calculateMax(i_axis, units);
91  const auto& axis_name = axisName(i_axis, units);
92  const auto axis_size = axisSize(i_axis);
93  return std::make_unique<FixedBinAxis>(axis_name, axis_size, min, max);
94 }
double calculateMax(size_t i_axis, Axes::Units units_type) const final
double calculateMin(size_t i_axis, Axes::Units units_type) const final
size_t axisSize(size_t i_axis) const final

References IUnitConverter::axisName(), UnitConverterSimple::axisSize(), UnitConverterSimple::calculateMax(), and UnitConverterSimple::calculateMin().

Here is the call graph for this function:

◆ addDetectorAxis()

void UnitConverterSimple::addDetectorAxis ( const IDetector detector,
size_t  i_axis 
)
protectedinherited

Definition at line 102 of file SimpleUnitConverters.cpp.

103 {
104  const auto& axis = detector.getAxis(i_axis);
105  const auto* p_roi = detector.regionOfInterest();
106  const auto& axis_name = axisName(i_axis);
107  if (p_roi) {
108  auto P_roi_axis = p_roi->clipAxisToRoi(i_axis, axis);
109  addAxisData(axis_name, P_roi_axis->getMin(), P_roi_axis->getMax(), defaultUnits(),
110  P_roi_axis->size());
111  } else {
112  addAxisData(axis_name, axis.getMin(), axis.getMax(), defaultUnits(), axis.size());
113  }
114 }
virtual const RegionOfInterest * regionOfInterest() const =0
Returns region of interest if exists.
virtual Axes::Units defaultUnits() const =0

References UnitConverterSimple::addAxisData(), IUnitConverter::axisName(), IUnitConverter::defaultUnits(), IDetector::getAxis(), and IDetector::regionOfInterest().

Referenced by RectangularConverter::RectangularConverter(), and SphericalConverter::SphericalConverter().

Here is the call graph for this function:

◆ addAxisData()

void UnitConverterSimple::addAxisData ( std::string  name,
double  min,
double  max,
Axes::Units  default_units,
size_t  nbins 
)
protectedinherited

Definition at line 46 of file SimpleUnitConverters.cpp.

48 {
49  AxisData axis_data{name, min, max, default_units, nbins};
50  m_axis_data_table.push_back(axis_data);
51 }

References UnitConverterSimple::m_axis_data_table.

Referenced by UnitConverterSimple::addDetectorAxis(), addDetectorYAxis(), DepthProbeConverter::DepthProbeConverter(), and OffSpecularConverter().

◆ axisName()

std::string IUnitConverter::axisName ( size_t  i_axis,
Axes::Units  units_type = Axes::Units::DEFAULT 
) const
inherited

Definition at line 20 of file IUnitConverter.cpp.

21 {
22  const auto& name_maps = createNameMaps();
23  if (name_maps.size() <= i_axis)
24  throw std::runtime_error("Error in IUnitConverter::axisName: the size of name map vector "
25  "is smaller or equal to the axis index"
26  + std::to_string(static_cast<int>(i_axis)));
27  const auto& name_map = name_maps[i_axis];
28  units_type = substituteDefaultUnits(units_type);
29  const auto& it = name_map.find(units_type);
30  if (it == name_map.cend())
31  throwUnitsError("IUnitConverter::axisName", availableUnits());
32  return it->second;
33 }
virtual std::vector< Axes::Units > availableUnits() const =0
virtual std::vector< std::map< Axes::Units, std::string > > createNameMaps() const =0

References IUnitConverter::availableUnits(), IUnitConverter::createNameMaps(), IUnitConverter::substituteDefaultUnits(), and IUnitConverter::throwUnitsError().

Referenced by UnitConverterSimple::addDetectorAxis(), addDetectorYAxis(), UnitConverterSimple::createConvertedAxis(), UnitConverter1D::createConvertedAxis(), DepthProbeConverter::DepthProbeConverter(), OffSpecularConverter(), and UnitConverterConvSpec::UnitConverterConvSpec().

Here is the call graph for this function:

◆ substituteDefaultUnits()

Axes::Units IUnitConverter::substituteDefaultUnits ( Axes::Units  units) const
inherited

Definition at line 65 of file IUnitConverter.cpp.

66 {
67  return units == Axes::Units::DEFAULT ? defaultUnits() : units;
68 }

References IUnitConverter::defaultUnits().

Referenced by IUnitConverter::axisName(), UnitConverterSimple::calculateMax(), UnitConverter1D::calculateMax(), UnitConverterSimple::calculateMin(), UnitConverter1D::calculateMin(), and UnitConverter1D::createConvertedAxis().

Here is the call graph for this function:

◆ createConvertedData()

std::unique_ptr< OutputData< double > > IUnitConverter::createConvertedData ( const OutputData< double > &  data,
Axes::Units  units 
) const
virtualinherited

Creates OutputData array in converter units.

Reimplemented in UnitConverter1D.

Definition at line 36 of file IUnitConverter.cpp.

37 {
38  const size_t dim = data.getRank();
39  std::unique_ptr<OutputData<double>> result(new OutputData<double>);
40  for (size_t i = 0; i < dim; ++i)
41  result->addAxis(*createConvertedAxis(i, units));
42  result->setRawDataVector(data.getRawDataVector());
43  return result;
44 }
virtual std::unique_ptr< IAxis > createConvertedAxis(size_t i_axis, Axes::Units units) const =0
std::vector< T > getRawDataVector() const
Returns copy of raw data vector.
Definition: OutputData.h:335
size_t getRank() const
Returns number of dimensions.
Definition: OutputData.h:59

References IUnitConverter::createConvertedAxis(), OutputData< T >::getRank(), and OutputData< T >::getRawDataVector().

Here is the call graph for this function:

◆ checkIndex()

void IUnitConverter::checkIndex ( size_t  i_axis) const
protectedinherited

Definition at line 46 of file IUnitConverter.cpp.

47 {
48  if (i_axis < dimension())
49  return;
50  throw std::runtime_error("Error in IUnitConverter::checkIndex: passed axis index too big: "
51  + std::to_string(static_cast<int>(i_axis)));
52 }
virtual size_t dimension() const =0

References IUnitConverter::dimension().

Referenced by UnitConverterSimple::axisSize(), UnitConverterConvSpec::axisSize(), UnitConverterQSpec::axisSize(), UnitConverterSimple::calculateMax(), UnitConverter1D::calculateMax(), UnitConverterSimple::calculateMin(), UnitConverter1D::calculateMin(), and UnitConverter1D::createConvertedAxis().

Here is the call graph for this function:

◆ throwUnitsError()

void IUnitConverter::throwUnitsError ( std::string  method,
std::vector< Axes::Units available 
) const
protectedinherited

Definition at line 54 of file IUnitConverter.cpp.

55 {
56  std::stringstream ss;
57  ss << "Unit type error in " << method
58  << ": unknown or unsupported unit type. Available units "
59  "are:\n";
60  for (auto unit : available)
61  ss << axisUnitName.at(unit) << "\n";
62  throw std::runtime_error(ss.str());
63 }
const std::map< Axes::Units, const char * > axisUnitName

References axisUnitName.

Referenced by IUnitConverter::axisName(), SphericalConverter::calculateValue(), RectangularConverter::calculateValue(), calculateValue(), DepthProbeConverter::checkUnits(), UnitConverterConvSpec::getTraslatorFrom(), UnitConverterConvSpec::getTraslatorTo(), and UnitConverterQSpec::getTraslatorTo().

◆ transferToCPP()

virtual void ICloneable::transferToCPP ( )
inlinevirtualinherited

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

Definition at line 34 of file ICloneable.h.

Member Data Documentation

◆ m_axis_data_table

◆ m_wavelength

◆ m_alpha_i

double UnitConverterSimple::m_alpha_i
protectedinherited

◆ m_phi_i

double UnitConverterSimple::m_phi_i
protectedinherited

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