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
|
#!/usr/bin/env python3
import bornagain as ba
from bornagain import ba_plot as bp, deg, nm, std_samples, std_simulations
material_particle = ba.RefractiveMaterial("Particle", 1e-6, 0)
ff = ba.Box(12*nm, 4.1*nm, 6*nm)
def get_sample(omega):
particle = ba.Particle(material_particle, ff)
particle.rotate(ba.RotationZ(omega*deg))
return std_samples.sas_sample_with_particle(particle)
def get_simulation(sample):
n = 11
return std_simulations.sas(sample, n)
if __name__ == '__main__':
results = []
for omega in [0, 30, 45]:
title = r'$\omega=%d^\circ$' % omega
sample = get_sample(omega)
simulation = get_simulation(sample)
result = simulation.simulate()
result.setTitle(title)
results.append(result)
from bornagain import ba_check
ba_check.persistence_test(results)
|
