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
|
import os
import sys
import numpy as np
dir = os.path.split(os.path.realpath(__file__))[0]
sys.path.append(dir + '/..')
import models
from stl import mesh_from_stl
from geometry import *
from materials import *
r7081_outer_mesh = mesh_from_stl(models.dir + '/hamamatsu_12inch_outer.stl')
r7081_inner_mesh = mesh_from_stl(models.dir + '/hamamatsu_12inch_inner.stl')
photocathode_triangles = np.mean(r7081_inner_mesh[:], axis=1)[:,1] > 0
inner_color = np.empty(len(r7081_inner_mesh.triangles), np.uint32)
inner_color[photocathode_triangles] = 0xff0000
inner_color[~photocathode_triangles] = 0x00ff00
inner_surface = np.empty(len(r7081_inner_mesh.triangles), np.object)
inner_surface[photocathode_triangles] = black_surface
inner_surface[~photocathode_triangles] = shiny_surface
r7081_inner_solid = Solid(r7081_inner_mesh, vacuum, glass, inner_surface, color=inner_color)
r7081_outer_solid = Solid(r7081_outer_mesh, glass, lightwater_sno)
r7081 = r7081_inner_solid + r7081_outer_solid
if __name__ == '__main__':
from view import view
from copy import deepcopy
r7081_outer_mesh_cutaway = deepcopy(r7081_outer_mesh)
r7081_outer_mesh_cutaway.triangles = \
r7081_outer_mesh_cutaway.triangles[\
np.mean(r7081_outer_mesh_cutaway[:], axis=1)[:,0] > 0]
r7081_outer_solid_cutaway = Solid(r7081_outer_mesh_cutaway, glass, lightwater_sno)
r7081_cutaway = r7081_inner_solid + r7081_outer_solid_cutaway
geometry = Geometry()
geometry.add_solid(r7081_cutaway)
geometry.build(bits=8)
view(geometry, 'r7081_cutaway')
|
