81 if (fRadius < newRadius) fRadius = newRadius;
83 else if (join.
mag () + newRadius <= fRadius) {
88 G4Point3D oldExtremity1 = fCentre - fRadius * unitJoin;
89 G4Point3D newExtremity1 = newCentre - newRadius * unitJoin;
90 G4Point3D oldExtremity2 = fCentre + fRadius * unitJoin;
91 G4Point3D newExtremity2 = newCentre + newRadius * unitJoin;
93 if (oldExtremity1 * unitJoin < newExtremity1 * unitJoin) {
94 extremity1 = oldExtremity1;
97 extremity1 = newExtremity1;
100 if (oldExtremity2 * unitJoin > newExtremity2 * unitJoin) {
101 extremity2 = oldExtremity2;
104 extremity2 = newExtremity2;
106 fCentre = 0.5 * (extremity2 + extremity1);
107 fRadius = 0.5 * (extremity2 - extremity1).mag ();