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
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
|
// ************************************************************************************************
//
// BornAgain: simulate and fit reflection and scattering
//
//! @file Img3D/Model/ParticleFromFF.cpp
//! @brief Implements namespace GUI::View::TransformTo3D.
//!
//! @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)
//
// ************************************************************************************************
#include "Img3D/Model/ParticleFromFF.h"
#include "Img3D/Model/Particles.h"
#include "Sample/HardParticle/HardParticles.h"
std::unique_ptr<Img3D::PlotParticle> Img3D::particle3DfromFF(const IFormfactor* ff)
{
if (const auto* f = dynamic_cast<const Pyramid2*>(ff)) {
double length = f->length();
double width = f->width();
double height = f->height();
double alpha = f->alpha();
return std::make_unique<Img3D::Particles::Pyramid2>(length, width, height, alpha);
}
if (const auto* f = dynamic_cast<const BarGauss*>(ff)) {
double length = f->length();
double width = f->width();
double height = f->height();
return std::make_unique<Img3D::Particles::BarGauss>(length, width, height);
}
if (const auto* f = dynamic_cast<const BarLorentz*>(ff)) {
double length = f->length();
double width = f->width();
double height = f->height();
return std::make_unique<Img3D::Particles::BarLorentz>(length, width, height);
}
if (const auto* f = dynamic_cast<const Box*>(ff)) {
double length = f->length();
double width = f->width();
double height = f->height();
return std::make_unique<Img3D::Particles::Box>(length, width, height);
}
if (const auto* f = dynamic_cast<const Cone*>(ff)) {
double radius = f->radius();
double height = f->height();
double alpha = f->alpha();
return std::make_unique<Img3D::Particles::Cone>(radius, height, alpha);
}
if (const auto* f = dynamic_cast<const Pyramid6*>(ff)) {
double baseedge = f->baseEdge();
double height = f->height();
double alpha = f->alpha();
return std::make_unique<Img3D::Particles::Pyramid6>(baseedge, height, alpha);
}
if (const auto* f = dynamic_cast<const Bipyramid4*>(ff)) {
double length = f->length();
double base_height = f->base_height();
double height_ratio = f->heightRatio();
double alpha = f->alpha();
return std::make_unique<Img3D::Particles::Bipyramid4>(length, base_height, height_ratio,
alpha);
}
if (const auto* f = dynamic_cast<const Cylinder*>(ff)) {
double radius = f->radius();
double height = f->height();
return std::make_unique<Img3D::Particles::Cylinder>(radius, height);
}
if (const auto* f = dynamic_cast<const Dodecahedron*>(ff)) {
double edge = f->edge();
return std::make_unique<Img3D::Particles::Dodecahedron>(edge);
}
if (const auto* f = dynamic_cast<const EllipsoidalCylinder*>(ff)) {
double radius_x = f->radiusX();
double radius_y = f->radiusY();
double height = f->height();
return std::make_unique<Img3D::Particles::EllipsoidalCylinder>(radius_x, radius_y, height);
}
if (const auto* f = dynamic_cast<const Sphere*>(ff)) {
double radius = f->radius();
return std::make_unique<Img3D::Particles::Sphere>(radius);
}
if (const auto* f = dynamic_cast<const Spheroid*>(ff)) {
double radius_xy = f->radiusXY();
double radius_z = f->radiusZ();
return std::make_unique<Img3D::Particles::Spheroid>(radius_xy, radius_z);
}
if (const auto* f = dynamic_cast<const HemiEllipsoid*>(ff)) {
double radius_xy = f->radiusX();
double radius_z = f->radiusY();
double height = f->radiusZ();
return std::make_unique<Img3D::Particles::HemiEllipsoid>(radius_xy, radius_z, height);
}
if (const auto* f = dynamic_cast<const Icosahedron*>(ff)) {
double edge = f->edge();
return std::make_unique<Img3D::Particles::Icosahedron>(edge);
}
if (const auto* f = dynamic_cast<const Prism3*>(ff)) {
double baseedge = f->baseEdge();
double height = f->height();
return std::make_unique<Img3D::Particles::Prism3>(baseedge, height);
}
if (const auto* f = dynamic_cast<const Prism6*>(ff)) {
double baseedge = f->baseEdge();
double height = f->height();
return std::make_unique<Img3D::Particles::Prism6>(baseedge, height);
}
if (const auto* f = dynamic_cast<const Pyramid4*>(ff)) {
double baseedge = f->baseEdge();
double height = f->height();
double alpha = f->alpha();
return std::make_unique<Img3D::Particles::Pyramid4>(baseedge, height, alpha);
}
if (const auto* f = dynamic_cast<const CosineRippleBox*>(ff)) {
double length = f->length();
double width = f->width();
double height = f->height();
return std::make_unique<Img3D::Particles::CosineRippleBox>(length, width, height);
}
if (const auto* f = dynamic_cast<const CosineRippleGauss*>(ff)) {
double length = f->length();
double width = f->width();
double height = f->height();
return std::make_unique<Img3D::Particles::CosineRippleGauss>(length, width, height);
}
if (const auto* f = dynamic_cast<const CosineRippleLorentz*>(ff)) {
double length = f->length();
double width = f->width();
double height = f->height();
return std::make_unique<Img3D::Particles::CosineRippleLorentz>(length, width, height);
}
if (const auto* f = dynamic_cast<const SawtoothRippleBox*>(ff)) {
double length = f->length();
double width = f->width();
double height = f->height();
return std::make_unique<Img3D::Particles::SawtoothRippleBox>(length, width, height);
}
if (const auto* f = dynamic_cast<const SawtoothRippleGauss*>(ff)) {
double length = f->length();
double width = f->width();
double height = f->height();
return std::make_unique<Img3D::Particles::SawtoothRippleGauss>(length, width, height);
}
if (const auto* f = dynamic_cast<const SawtoothRippleLorentz*>(ff)) {
double length = f->length();
double width = f->width();
double height = f->height();
return std::make_unique<Img3D::Particles::SawtoothRippleLorentz>(length, width, height);
}
if (const auto* f = dynamic_cast<const Pyramid3*>(ff)) {
double baseedge = f->baseEdge();
double height = f->height();
double alpha = f->alpha();
return std::make_unique<Img3D::Particles::Pyramid3>(baseedge, height, alpha);
}
if (const auto* f = dynamic_cast<const TruncatedCube*>(ff)) {
double length = f->length();
double removed_length = f->removedLength();
return std::make_unique<Img3D::Particles::TruncatedCube>(length, removed_length);
}
if (const auto* f = dynamic_cast<const SphericalSegment*>(ff)) {
double radius = f->radius();
double untruncated_height = 2 * radius - f->cutFromBottom();
double cutFromTop = f->cutFromTop();
return std::make_unique<Img3D::Particles::SphericalSegment>(radius, untruncated_height,
cutFromTop);
}
if (const auto* f = dynamic_cast<const SpheroidalSegment*>(ff)) {
double radius = f->radiusXY();
double untruncated_height = 2 * f->radiusZ() - f->cutFromBottom();
double hfc = f->radiusZ() / f->radiusXY();
double deltaH = f->cutFromTop();
return std::make_unique<Img3D::Particles::SpheroidalSegment>(radius, untruncated_height,
hfc, deltaH);
}
if (const auto* f = dynamic_cast<const CantellatedCube*>(ff)) {
double length = f->length();
double removed_length = f->removedLength();
return std::make_unique<Img3D::Particles::CantellatedCube>(length, removed_length);
}
if (const auto* f = dynamic_cast<const HorizontalCylinder*>(ff)) {
double radius = f->radius();
double length = f->length();
double slice_bottom = f->slice_bottom();
double slice_top = f->slice_top();
return std::make_unique<Img3D::Particles::HorizontalCylinder>(radius, length, slice_bottom,
slice_top);
}
if (const auto* f = dynamic_cast<const PlatonicOctahedron*>(ff)) {
double edge = f->edge();
return std::make_unique<Img3D::Particles::PlatonicOctahedron>(edge);
}
if (const auto* f = dynamic_cast<const PlatonicTetrahedron*>(ff)) {
double edge = f->edge();
return std::make_unique<Img3D::Particles::PlatonicTetrahedron>(edge);
}
return {};
}
|
