import os
import sys
import numpy as np

dir = os.path.split(os.path.realpath(__file__))[0]
sys.path.append(dir + '/..')

import models
from mesh 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