Skip to content

Place the L3 coil turns explicitly instead of dividing the polyhedra #15767

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Open
sawenzel wants to merge 1 commit into
base: dev
Choose a base branch
from
+13 −2
Open

Place the L3 coil turns explicitly instead of dividing the polyhedra #15767

Changes from all commits
Commits
Show all changes
1 commit
Select commit
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
15 changes: 13 additions & 2 deletions Detectors/Passive/src/Magnet.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -225,8 +225,19 @@ void Magnet::ConstructGeometry()
// Coils
TGeoVolume* voCoilMother = new TGeoVolume("L3CM", shCoilMother, medAir);
voBMother->AddNode(voCoilMother, 1, new TGeoTranslation(0., 0., 0.));
// Divide into the 168 turns
TGeoVolume* voCoilTurn = voCoilMother->Divide("L3CD", 3, 168, 0., 0.);
// The 168 turns, placed explicitly rather than as a TGeoPgon division.
// Geant4's G4ParameterisationPolyhedraZ rebuilds the shared master solid while
// navigating, so a divided polyhedra crashes a multithreaded native-Geant4 run.
const Int_t kNCoilTurns = 168;
const Float_t kDzCoilTurn = kLCoil / kNCoilTurns;
TGeoPgon* shCoilTurn = new TGeoPgon(kStartAngle, kFullAngle, kNSides, 2);
shCoilTurn->DefineSection(0, -kDzCoilTurn, kRCoilInner - 2. * kRCoolingOuter, kRCoilOuter + 2. * kRCoolingOuter);
shCoilTurn->DefineSection(1, kDzCoilTurn, kRCoilInner - 2. * kRCoolingOuter, kRCoilOuter + 2. * kRCoolingOuter);
TGeoVolume* voCoilTurn = new TGeoVolume("L3CD", shCoilTurn, medAir);
for (Int_t iTurn = 0; iTurn < kNCoilTurns; ++iTurn) {
voCoilMother->AddNode(voCoilTurn, iTurn + 1,
new TGeoTranslation(0., 0., -kLCoil + (2 * iTurn + 1) * kDzCoilTurn));
}
TGeoPgon* shCoils = new TGeoPgon(kStartAngle, kFullAngle, kNSides, 2);
shCoils->DefineSection(0, -3., kRCoilInner, kRCoilOuter);
shCoils->DefineSection(1, 3., kRCoilInner, kRCoilOuter);
Expand Down