From 332ec7a527857cc8c5f92c16ee94326b4cbfcbe1 Mon Sep 17 00:00:00 2001 From: Stan Seibert Date: Fri, 12 Aug 2011 15:51:20 -0400 Subject: Segmented cylinder FOR MOAR SPEED! 1.8 Mphotons/sec! --- detectors/lbne.py | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) (limited to 'detectors') diff --git a/detectors/lbne.py b/detectors/lbne.py index 3e27907..d3c12da 100644 --- a/detectors/lbne.py +++ b/detectors/lbne.py @@ -10,7 +10,7 @@ from solids import build_12inch_pmt, build_12inch_pmt_with_lc, build_12inch_pmt_ from optics import * from transform import rotate, make_rotation_matrix from itertools import product -from make import cylinder +from make import cylinder, segmented_cylinder def build_lbne(radius, height, nstrings, pmts_per_string, endcap_spacing, physical_model=True): if physical_model: @@ -21,7 +21,7 @@ def build_lbne(radius, height, nstrings, pmts_per_string, endcap_spacing, physic lbne = Geometry() # outer cylinder - cylinder_mesh = cylinder(radius, radius, height+height/(pmts_per_string-1), theta=(2*np.pi/nstrings)/4) + cylinder_mesh = segmented_cylinder(radius, height+height/(pmts_per_string-1), theta=(2*np.pi/nstrings)/8, n=200) cylinder_mesh.vertices = rotate(cylinder_mesh.vertices, np.pi/2, (-1,0,0)) lbne.add_solid(Solid(cylinder_mesh, water_wcsim, vacuum, black_surface, 0xff0000ff)) -- cgit