ditto for Icosahedron

Doc/UserManual/FormFactors.tex

@@ -731,7 +731,7 @@ with inversion symmetry~\cite{ba:ffp},
 
 \begin{figure}[H]
 \begin{center}
-\includegraphics[width=\textwidth]{fig/ff2/ff_Dodecahedron.pdf}
+\includegraphics[width=\textwidth]{fig/ff2/ff_Dodecahedron_sym.pdf}
 \end{center}
 \caption{Normalized intensity $|F|^2/V^2$,
 computed with $a=3.2$~nm,
@@ -741,6 +741,18 @@ vertex on the $x$ axis;
 edge in the $xy$ plane and perpendicular to the $x$ axis.}
 \end{figure}
 
+\begin{figure}[H]
+\begin{center}
+\includegraphics[width=\textwidth]{fig/ff2/ff_Dodecahedron_asy.pdf}
+\end{center}
+\caption{Normalized intensity $|F|^2/V^2$,
+computed with $a=3.2$~nm,
+for three orientations of decreasing symmetry:
+base pentagon in $xy$ plane and pointing in $x$ direction;
+rotated by $13^\circ$ around the $z$ axis;
+ditto, and tilted by $9^\circ$ around the $x$ axis.}
+\end{figure}
+
 \paragraph{History}\strut\\
 New in \BornAgain-1.6,
 based on the generic form factor of the polyhedron~\cite{ba:ffp}.
@@ -1053,11 +1065,26 @@ with inversion symmetry~\cite{ba:ffp},
 
 \begin{figure}[H]
 \begin{center}
-\includegraphics[width=\textwidth]{fig/ff2/ff_Icosahedron.pdf}
+\includegraphics[width=\textwidth]{fig/ff2/ff_Icosahedron_sym.pdf}
+\end{center}
+\caption{Normalized intensity $|F|^2/V^2$,
+computed with $a=4.8$~nm,
+for three orientations of high symmetry:
+$x$ axis perpendicular to a polygonal face;
+vertex on the $x$ axis;
+edge in the $xy$ plane and perpendicular to the $x$ axis.}
+\end{figure}
+
+\begin{figure}[H]
+\begin{center}
+\includegraphics[width=\textwidth]{fig/ff2/ff_Icosahedron_asy.pdf}
 \end{center}
 \caption{Normalized intensity $|F|^2/V^2$,
 computed with $a=4.8$~nm,
-for three different angles~$\omega$ of rotation around the $z$ axis.}
+for three orientations of decreasing symmetry:
+base pentagon in $xy$ plane and pointing in $x$ direction;
+rotated by $13^\circ$ around the $z$ axis;
+ditto, and tilted by $9^\circ$ around the $x$ axis.}
 \end{figure}
 
 \paragraph{History}\strut\\

Doc/UserManual/fig/ff2/sim_Icosahedron.py

-"""
-Plot form factors.
-"""
-import bornagain as ba
-from   bornagain import nanometer, degree
-import bornplot as bp
-
-det = bp.Detector( 200, -5, 5, -5, 5 )
-n    = 3
-results = []
-for i in range(n):
-    omega=30*i/(n-1)
-    title = r'$\omega=%d^\circ$' % omega
-    ff = ba.FormFactorIcosahedron(4.8*nanometer)
-    trafo = ba.RotationZ(omega*degree)
-    data = bp.run_simulation(det,ff,trafo)
-    results.append( bp.Result(i, data, title) )
-    
-bp.make_plot( results, det, "ff_Icosahedron" )

Doc/UserManual/fig/ff2/sim_Icosahedron_asy.py

+"""
+Plot form factors.
+"""
+import bornagain as ba
+from   bornagain import nanometer, degree
+import bornplot as bp
+import math
+
+det = bp.Detector( 200, -5, 5, -5, 5 )
+n    = 3
+results = []
+edge = 4.8
+
+title = 'default'
+ff = ba.FormFactorIcosahedron(edge*nanometer)
+data = bp.run_simulation(det,ff)
+results.append( bp.Result(0, data, title) )
+
+title = 'rotated'
+trafo = ba.RotationZ(13*degree)
+ff = ba.FormFactorIcosahedron(edge*nanometer)
+data = bp.run_simulation(det,ff,trafo)
+results.append( bp.Result(1, data, title) )
+
+title = 'rotated, tilted'
+trafo = ba.CreateProduct( ba.RotationX(9*degree), ba.RotationZ(13*degree))
+ff = ba.FormFactorIcosahedron(edge*nanometer)
+data = bp.run_simulation(det,ff,trafo)
+results.append( bp.Result(2, data, title) )
+
+bp.make_plot( results, det, "ff_Icosahedron_asy" )

Doc/UserManual/fig/ff2/sim_Icosahedron_sym.py

+"""
+Plot form factors.
+"""
+import bornagain as ba
+from   bornagain import nanometer, degree
+import bornplot as bp
+import math
+
+det = bp.Detector( 200, -5, 5, -5, 5 )
+n    = 3
+results = []
+edge = 4.8
+
+title = 'face normal'
+trafo = ba.RotationY(48.1897*degree)
+ff = ba.FormFactorIcosahedron(edge*nanometer)
+data = bp.run_simulation(det,ff,trafo)
+results.append( bp.Result(0, data, title) )
+
+title = 'vertex normal'
+trafo = ba.RotationY(-52.6226*degree)
+ff = ba.FormFactorIcosahedron(edge*nanometer)
+data = bp.run_simulation(det,ff,trafo)
+results.append( bp.Result(1, data, title) )
+
+title = 'edge normal'
+trafo = ba.RotationY(69.0948*degree)
+ff = ba.FormFactorIcosahedron(edge*nanometer)
+data = bp.run_simulation(det,ff,trafo)
+results.append( bp.Result(2, data, title) )
+
+bp.make_plot( results, det, "ff_Icosahedron_sym" )