Prepare demo geometry

[fiedl]

scenario including:

• [x] a DOM
• [x] a cable
• [x] drill hole
• [x] bubble column
• [x] several photons

requirements:

• [x] kernel code ➝ #12
• [x] passing parameters ➝ #15
• [ ] cylinders with bottom and top
• [ ] improve memory allocation
• [ ] store more photon scatterings

[fiedl]

DOM Position

In 2014, for debugging purposes, I've noted down the position of DOM (1,1):

I3Position(-256.02301025390625, -521.281982421875, 500)

Can this position change, e.g. by position calibration, or by changing the coordinate system? Because the coordinates seem to be off here:

I haven't found a list of the DOMs and their positions. dataio-shovel does not show the geometry information. As a workaround, I'm inserting a debugging line into steamshovel in steamshovel/python/artists/DOMLabel.py:

print(geom.omgeo[omkey].position)
# => I3Position(-256.14,-521.08,496.03)

[fiedl]

Drill hole

Let's create a geometry where the DOM sits at one side of the drill hole.

The drill-hole radius is 30cm. The DOM radius is 0.16510m.

[fiedl]

DOM Radius

The real DOM radius is 0.16510m, which cannot be displayed in steamshovel. In steamshovel the DOM radius is set to 1m and can be changed in steps of 10cm. I need to modify steamshovel in order to be able to display the DOM with a radius of 0.16510m.

See also: https://github.com/fiedl/hole-ice-study/commit/a9e663cd0c48b5b8c343fcca94d035973eea30e7

[fiedl]

# Configure a drill hole in which the DOM sits at one side of the hole.
#
drill_hole_cylinder = I3CLSimMediumCylinder()
drill_hole_cylinder.x = -256.14 + (0.30 - 0.16510)
drill_hole_cylinder.y = -521.08
drill_hole_cylinder.radius = 0.30
drill_hole_cylinder.scattering_length = 1.0
drill_hole_cylinder.absorption_length = 100.0

# Configure bubble column with small scattering length in the center
# of the drill hole.
#
bubble_column_cylinder = I3CLSimMediumCylinder()
bubble_column_cylinder.x = -256.14 + (0.30 - 0.16510)
bubble_column_cylinder.y = -521.08
bubble_column_cylinder.radius = 0.10
bubble_column_cylinder.scattering_length = 0.1
bubble_column_cylinder.absorption_length = 100.0

# Configure opaque cable right beside the DOM.
#
cable_cylinder = I3CLSimMediumCylinder()
cable_cylinder.x = -256.14 + (0.16510 + 0.03)
cable_cylinder.y = -521.08
cable_cylinder.top = 496.03 + 0.5
cable_cylinder.bottom = 496.03 - 0.5
cable_cylinder.radius = 0.03
cable_cylinder.scattering_length = 100.0
cable_cylinder.absorption_length = 0.0

cylinders = I3CLSimMediumCylinderSeries([drill_hole_cylinder, bubble_column_cylinder, cable_cylinder])
frame.Put("I3CLSimMediumCylinders", cylinders)