NS_PER_TDC = FloatCol()
adc = FloatArrayCol(GPUDaqUBooNE.NTDC)
q = FloatCol()
t = FloatCol()
class OpMap(TreeModel):
opid = IntCol()
x = FloatCol()
y = FloatCol()
z = FloatCol()
lar1nd_wireplane = Surface('lar1nd_wireplane')
lar1nd_wireplane.nplanes = 3.0
lar1nd_wireplane.wire_pitch = 0.3
lar1nd_wireplane.wire_diameter = 0.015
lar1nd_wireplane.transmissive = 1
lar1nd_wireplane.model = Surface.SURFACE_WIREPLANE
def add_wireplane_surface(solid):
# function detector class will use to add a wireplane surface to the geometry
# LAr1ND has two drift regions, so we need two planes
# set surface for triangles on x=-2023.25 and x=2023.25 planes
for n, triangle in enumerate(solid.mesh.triangles):
#print [ solid.mesh.vertices[x] for x in triangle ] # for debug
nxplane = 0
for ivert in triangle:
if solid.mesh.vertices[
ivert, 0] == -2023.25 or solid.mesh.vertices[ivert,
surface_reemit = IntCol()
def reset(self):
self.reflect_diffuse = 0
self.reflect_specular = 0
self.bulk_scatter = 0
self.bulk_absorb = 0
self.surface_detect = 0
self.surface_absorb = 0
self.surface_reemit = 0
uboone_wireplane = Surface('uboone_wireplane')
uboone_wireplane.nplanes = 3.0
uboone_wireplane.wire_pitch = 0.3
uboone_wireplane.wire_diameter = 0.015
uboone_wireplane.transmissive = 1
uboone_wireplane.model = Surface.SURFACE_WIREPLANE
def add_wireplane_surface(solid):
# function detector class will use to add a wireplane surface to the geometry
# set surface for triangles on x=-1281.0 plane
for n, triangle in enumerate(solid.mesh.triangles):
nxplane = 0
for ivert in triangle:
if solid.mesh.vertices[ivert, 0] == -1281.0:
nxplane += 1
if nxplane == 3:
print[solid.mesh.vertices[x] for x in triangle]
solid.surface[n] = uboone_wireplane
import os, sys
from chroma.importgeo import UserVG4DEAGeo, load_hist_data
from chroma.geometry import Surface
import numpy as np
wireplane = Surface('wireplane')
wireplane.nplanes = 3.0
wireplane.wire_pitch = 0.3
wireplane.wire_diameter = 0.015
wireplane.transmissive = 1
wireplane.model = Surface.SURFACE_WIREPLANE
class SBND(UserVG4DEAGeo):
def __init__(self):
super(SBND,
self).__init__("SBND",
"dae/lar1nd_lightguides_nowires_chroma.dae")
def surfacesdict(self):
# Steel/LAr
steel_surface = Surface("steel_surface")
steel_surface.set('reflect_diffuse', 0.25)
steel_surface.set('reflect_specular', 0.0)
steel_surface.set('detect', 0.0)
steel_surface.set('absorb', 0.75)
steel_surface.set('reemit', 0.0)
steel_surface.transmissive = 0
# Titanium/LAr
titanium_surface = Surface("titanium_surface")
titanium_surface.set('reflect_diffuse', 0.125)