1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
|
// ************************************************************************************************
//
// BornAgain: simulate and fit reflection and scattering
//
//! @file Device/Resolution/ConvolutionDetectorResolution.h
//! @brief Defines class ConvolutionDetectorResolution.
//!
//! @homepage http://www.bornagainproject.org
//! @license GNU General Public License v3 or higher (see COPYING)
//! @copyright Forschungszentrum Jülich GmbH 2018
//! @authors Scientific Computing Group at MLZ (see CITATION, AUTHORS)
//
// ************************************************************************************************
#ifdef SWIG
#error no need to expose this header to Swig
#endif // SWIG
#ifndef BORNAGAIN_DEVICE_RESOLUTION_CONVOLUTIONDETECTORRESOLUTION_H
#define BORNAGAIN_DEVICE_RESOLUTION_CONVOLUTIONDETECTORRESOLUTION_H
#include "Device/Resolution/IDetectorResolution.h"
#include "Device/Resolution/IResolutionFunction2D.h"
//! Convolutes the intensity in 1 or 2 dimensions with a resolution function.
//! Limitation: this class assumes that the data points are evenly distributed on each axis
class ConvolutionDetectorResolution : public IDetectorResolution {
public:
using cumulative_DF_1d = double (*)(double);
//! Constructor taking a 1 dimensional resolution function as argument.
ConvolutionDetectorResolution(cumulative_DF_1d res_function_1d);
//! Constructor taking a 2 dimensional resolution function as argument.
ConvolutionDetectorResolution(const IResolutionFunction2D& p_res_function_2d);
~ConvolutionDetectorResolution() override;
#ifndef SWIG
ConvolutionDetectorResolution* clone() const override;
std::vector<const INode*> nodeChildren() const override;
#endif // SWIG
std::string className() const final { return "ConvolutionDetectorResolution"; }
//! Convolve given intensities with the encapsulated resolution.
void execDetectorResolution(Datafield* df) const override;
const IResolutionFunction2D* getResolutionFunction2D() const;
protected:
ConvolutionDetectorResolution(const ConvolutionDetectorResolution& other);
private:
void setResolutionFunction(const IResolutionFunction2D& resFunc);
void apply1dConvolution(Datafield* df) const;
void apply2dConvolution(Datafield* df) const;
double getIntegratedPDF1d(double x, double step) const;
double getIntegratedPDF2d(double x, double step_x, double y, double step_y) const;
size_t m_rank;
cumulative_DF_1d m_res_function_1d;
std::unique_ptr<IResolutionFunction2D> m_res_function_2d;
};
inline const IResolutionFunction2D* ConvolutionDetectorResolution::getResolutionFunction2D() const
{
return m_res_function_2d.get();
}
#endif // BORNAGAIN_DEVICE_RESOLUTION_CONVOLUTIONDETECTORRESOLUTION_H
|
