diff options
Diffstat (limited to 'detectors/lbne.py')
| -rw-r--r-- | detectors/lbne.py | 10 |
1 files changed, 5 insertions, 5 deletions
diff --git a/detectors/lbne.py b/detectors/lbne.py index ca13c2d..97d09e9 100644 --- a/detectors/lbne.py +++ b/detectors/lbne.py @@ -6,18 +6,18 @@ from chroma.transform import rotate, make_rotation_matrix from itertools import product import chroma.make as make -def build_lbne(radius, height, nstrings, pmts_per_string, endcap_spacing, physical_model=True): +def build_lbne(radius, height, nstrings, pmts_per_string, endcap_spacing, detector_material=water_wcsim, physical_model=True): if physical_model: - pmt = solids.build_12inch_pmt_with_lc(outer_material=water_wcsim) + pmt = solids.build_12inch_pmt_with_lc(outer_material=detector_material) else: - pmt = solids.build_12inch_pmt_shell(outer_material=water_wcsim) + pmt = solids.build_12inch_pmt_shell(outer_material=detector_material) - lbne = Geometry() + lbne = Geometry(detector_material) # outer cylinder cylinder_mesh = make.segmented_cylinder(radius, height+height/(pmts_per_string-1), nsteps=16*nstrings, nsegments=1200) 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)) + lbne.add_solid(Solid(cylinder_mesh, detector_material, vacuum, black_surface, 0xff0000ff)) lbne.pmtids = [] |