Skip to content

Commit

Permalink
neutron shielding
Browse files Browse the repository at this point in the history
  • Loading branch information
@lobis
lobis committed Oct 6, 2024
1 parent 6d52bd5 commit 914740a
Showing 1 changed file with 96 additions and 89 deletions.
185 changes: 96 additions & 89 deletions generator/src/main/kotlin/BabyIAXO/ShieldingNeutrons.kt
View file
Original file line number Diff line number Diff line change
Expand Up @@ -5,14 +5,17 @@ import Geometry
import space.kscience.gdml.*

open class ShieldingNeutrons(
open val copperBox: Boolean = true,
open val innerVetoes: Boolean = false,
) : Geometry() {
companion object Parameters {

const val vetoOuterThickness: Double = 100.0
const val vetoInnerThickness: Double = 50.0
const val vetoInnerSpacing = 50.0
const val shieldingThicknessBase: Double = 200.0

const val shieldingThickness: Double = shieldingThicknessBase + vetoInnerThickness

const val shieldingThickness: Double = 200.0 + vetoInnerThickness
const val chamberHoleXY: Double = 200.0
const val chamberHoleZ: Double = 150.0

Expand All @@ -29,8 +32,7 @@ open class ShieldingNeutrons(
const val OffsetZ: Double =
DetectorToShieldingSeparation + Chamber.Height / 2 + Chamber.ReadoutKaptonThickness + Chamber.BackplateThickness

val vetoInnerSpacing = 50.0
val vetoInnerLength = shieldingThickness + chamberHoleZ + vetoInnerSpacing
const val vetoInnerLength = shieldingThickness + chamberHoleZ + vetoInnerSpacing
}

override fun generate(gdml: Gdml): GdmlRef<GdmlAssembly> {
Expand Down Expand Up @@ -108,59 +110,63 @@ open class ShieldingNeutrons(
position(z = SizeZ / 2 - ShaftLongSide / 2) { unit = LUnit.MM }
}

val leadShieldingVolumeHoleInnerVetoes1 =
gdml.solids.subtraction(
leadBoxWithShaftSolid,
vetoInnerBoxLateralSolid,
"leadShieldingVolumeHoleInnerVetoes1"
) {
position(
y = vetoInnerThickness / 2 + chamberHoleXY / 2 + vetoInnerSpacing,
x = vetoInnerThickness / 2,
z = -vetoInnerLength / 2 + SizeZ / 2
) { unit = LUnit.MM }
}
val leadShieldingVolumeHoleInnerVetoes2 =
gdml.solids.subtraction(
leadShieldingVolumeHoleInnerVetoes1,
vetoInnerBoxLateralSolid,
"leadShieldingVolumeHoleInnerVetoes2"
) {
position(
y = -vetoInnerThickness / 2 - chamberHoleXY / 2 - vetoInnerSpacing,
x = -vetoInnerThickness / 2,
z = -vetoInnerLength / 2 + SizeZ / 2
) { unit = LUnit.MM }
}
val leadShieldingVolumeHoleInnerVetoes3 =
gdml.solids.subtraction(
leadShieldingVolumeHoleInnerVetoes2,
vetoInnerBoxLateralSolid,
"leadShieldingVolumeHoleInnerVetoes3"
) {
rotation(z = 90.0) { unit = AUnit.DEGREE }
position(
x = -vetoInnerThickness / 2 - chamberHoleXY / 2 - vetoInnerSpacing,
y = vetoInnerThickness / 2,
z = -vetoInnerLength / 2 + SizeZ / 2
) { unit = LUnit.MM }
}
val leadShieldingVolumeHoleInnerVetoes4 =
gdml.solids.subtraction(
leadShieldingVolumeHoleInnerVetoes3,
vetoInnerBoxLateralSolid,
"leadShieldingVolumeHoleInnerVetoes4"
) {
rotation(z = 90.0) { unit = AUnit.DEGREE }
position(
x = vetoInnerThickness / 2 + chamberHoleXY / 2 + vetoInnerSpacing,
y = -vetoInnerThickness / 2,
z = -vetoInnerLength / 2 + SizeZ / 2
) { unit = LUnit.MM }
}

val leadShieldingVolume =
val leadShieldingVolume = if (innerVetoes) {
val leadShieldingVolumeHoleInnerVetoes1 =
gdml.solids.subtraction(
leadBoxWithShaftSolid,
vetoInnerBoxLateralSolid,
"leadShieldingVolumeHoleInnerVetoes1"
) {
position(
y = vetoInnerThickness / 2 + chamberHoleXY / 2 + vetoInnerSpacing,
x = vetoInnerThickness / 2,
z = -vetoInnerLength / 2 + SizeZ / 2
) { unit = LUnit.MM }
}
val leadShieldingVolumeHoleInnerVetoes2 =
gdml.solids.subtraction(
leadShieldingVolumeHoleInnerVetoes1,
vetoInnerBoxLateralSolid,
"leadShieldingVolumeHoleInnerVetoes2"
) {
position(
y = -vetoInnerThickness / 2 - chamberHoleXY / 2 - vetoInnerSpacing,
x = -vetoInnerThickness / 2,
z = -vetoInnerLength / 2 + SizeZ / 2
) { unit = LUnit.MM }
}
val leadShieldingVolumeHoleInnerVetoes3 =
gdml.solids.subtraction(
leadShieldingVolumeHoleInnerVetoes2,
vetoInnerBoxLateralSolid,
"leadShieldingVolumeHoleInnerVetoes3"
) {
rotation(z = 90.0) { unit = AUnit.DEGREE }
position(
x = -vetoInnerThickness / 2 - chamberHoleXY / 2 - vetoInnerSpacing,
y = vetoInnerThickness / 2,
z = -vetoInnerLength / 2 + SizeZ / 2
) { unit = LUnit.MM }
}
val leadShieldingVolumeHoleInnerVetoes4 =
gdml.solids.subtraction(
leadShieldingVolumeHoleInnerVetoes3,
vetoInnerBoxLateralSolid,
"leadShieldingVolumeHoleInnerVetoes4"
) {
rotation(z = 90.0) { unit = AUnit.DEGREE }
position(
x = vetoInnerThickness / 2 + chamberHoleXY / 2 + vetoInnerSpacing,
y = -vetoInnerThickness / 2,
z = -vetoInnerLength / 2 + SizeZ / 2
) { unit = LUnit.MM }
}

gdml.structure.volume(Materials.Lead.ref, leadShieldingVolumeHoleInnerVetoes4, "shieldingVolume")
} else {
gdml.structure.volume(Materials.Lead.ref, leadBoxWithShaftSolid, "shieldingVolume")
}

return@lazy gdml.structure.assembly {
name = "shielding"
Expand Down Expand Up @@ -213,43 +219,44 @@ open class ShieldingNeutrons(
}
}

// inner veto
physVolume(vetoInnerBoxLateralVolume, name = "vetoInnerTop") {
position(
y = vetoInnerThickness / 2 + chamberHoleXY / 2 + vetoInnerSpacing,
x = vetoInnerThickness / 2,
z = -OffsetZ - vetoInnerLength / 2 + SizeZ / 2
) {
unit = LUnit.MM
if (innerVetoes) {
physVolume(vetoInnerBoxLateralVolume, name = "vetoInnerTop") {
position(
y = vetoInnerThickness / 2 + chamberHoleXY / 2 + vetoInnerSpacing,
x = vetoInnerThickness / 2,
z = -OffsetZ - vetoInnerLength / 2 + SizeZ / 2
) {
unit = LUnit.MM
}
}
}
physVolume(vetoInnerBoxLateralVolume, name = "vetoInnerBottom") {
position(
y = -vetoInnerThickness / 2 - chamberHoleXY / 2 - vetoInnerSpacing,
x = -vetoInnerThickness / 2,
z = -OffsetZ - vetoInnerLength / 2 + SizeZ / 2
) {
unit = LUnit.MM
physVolume(vetoInnerBoxLateralVolume, name = "vetoInnerBottom") {
position(
y = -vetoInnerThickness / 2 - chamberHoleXY / 2 - vetoInnerSpacing,
x = -vetoInnerThickness / 2,
z = -OffsetZ - vetoInnerLength / 2 + SizeZ / 2
) {
unit = LUnit.MM
}
}
}
physVolume(vetoInnerBoxLateralVolume, name = "vetoInnerRight") {
rotation(z = 90.0) { unit = AUnit.DEGREE }
position(
x = -vetoInnerThickness / 2 - chamberHoleXY / 2 - vetoInnerSpacing,
y = vetoInnerThickness / 2,
z = -OffsetZ - vetoInnerLength / 2 + SizeZ / 2
) {
unit = LUnit.MM
physVolume(vetoInnerBoxLateralVolume, name = "vetoInnerRight") {
rotation(z = 90.0) { unit = AUnit.DEGREE }
position(
x = -vetoInnerThickness / 2 - chamberHoleXY / 2 - vetoInnerSpacing,
y = vetoInnerThickness / 2,
z = -OffsetZ - vetoInnerLength / 2 + SizeZ / 2
) {
unit = LUnit.MM
}
}
}
physVolume(vetoInnerBoxLateralVolume, name = "vetoInnerLeft") {
rotation(z = 90.0) { unit = AUnit.DEGREE }
position(
x = vetoInnerThickness / 2 + chamberHoleXY / 2 + vetoInnerSpacing,
y = -vetoInnerThickness / 2,
z = -OffsetZ - vetoInnerLength / 2 + SizeZ / 2
) {
unit = LUnit.MM
physVolume(vetoInnerBoxLateralVolume, name = "vetoInnerLeft") {
rotation(z = 90.0) { unit = AUnit.DEGREE }
position(
x = vetoInnerThickness / 2 + chamberHoleXY / 2 + vetoInnerSpacing,
y = -vetoInnerThickness / 2,
z = -OffsetZ - vetoInnerLength / 2 + SizeZ / 2
) {
unit = LUnit.MM
}
}
}
}
Expand Down

0 comments on commit 914740a

Please sign in to comment.