+ TGeoPcon *hole31shape = new TGeoPcon(-holePhi/2., holePhi, 4);
+
+ hole31shape->Rmin(0) = kHole3RMin + kHole3DeltaR;
+ hole31shape->Rmax(0) = hole31shape->GetRmin(0);
+ hole31shape->Z(0) = ZFromRminpCone(coneshape,3,kConeTheta,
+ hole31shape->GetRmin(0));
+
+ hole31shape->Rmax(1) = hole31shape->GetRmax(0);
+ hole31shape->Z(1) = ZFromRminpCone(coneinsertshape,3,kConeTheta,
+ hole31shape->GetRmax(1));
+ hole31shape->Rmin(1) = RminFromZpCone(coneshape,3,kConeTheta,
+ hole31shape->GetZ(1));
+
+ hole31shape->Rmin(2) = kHole3RMin;
+ hole31shape->Z(2) = ZFromRminpCone(coneshape,3,kConeTheta,
+ hole31shape->GetRmin(2));
+ hole31shape->Rmax(2) = RminFromZpCone(coneinsertshape,3,kConeTheta,
+ hole31shape->GetZ(2));
+
+ hole31shape->Rmin(3) = hole31shape->GetRmin(2);
+ hole31shape->Rmax(3) = hole31shape->GetRmin(3);
+ hole31shape->Z(3) = ZFromRminpCone(coneinsertshape,3,kConeTheta,
+ hole31shape->GetRmax(3));
+
+ TGeoPcon *hole32shape = new TGeoPcon(-holePhi/2., holePhi, 4);
+
+ hole32shape->Rmin(0) = kHole3RMin + kHole3DeltaR;
+ hole32shape->Rmax(0) = hole32shape->GetRmin(0);
+ hole32shape->Z(0) = ZFromRmaxpCone(coneinsertshape,4,kConeTheta,
+ hole32shape->GetRmin(0));
+
+ hole32shape->Rmax(1) = hole32shape->GetRmax(0);
+ hole32shape->Z(1) = ZFromRmaxpCone(coneshape,4,kConeTheta,
+ hole32shape->GetRmax(1));
+ hole32shape->Rmin(1) = RmaxFromZpCone(coneinsertshape,4,kConeTheta,
+ hole32shape->GetZ(1));
+
+ hole32shape->Rmin(2) = kHole3RMin;
+ hole32shape->Z(2) = ZFromRmaxpCone(coneinsertshape,4,kConeTheta,
+ hole32shape->GetRmin(2));
+ hole32shape->Rmax(2) = RmaxFromZpCone(coneshape,4,kConeTheta,
+ hole32shape->GetZ(2));
+
+ hole32shape->Rmin(3) = hole32shape->GetRmin(2);
+ hole32shape->Rmax(3) = hole32shape->GetRmin(3);
+ hole32shape->Z(3) = ZFromRmaxpCone(coneshape,4,kConeTheta,
+ hole32shape->GetRmax(3));
+
+ //