harness bug

icecube · Aug 26, 2024 · 7bf1528 · 7bf1528
1 parent ce023b2
commit 7bf1528
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/common/geometry_construction/src/OMSimMDOMHarness.cc b/common/geometry_construction/src/OMSimMDOMHarness.cc
@@ -255,7 +255,7 @@ void mDOMHarness::PCA()
     const G4double lWireRadius = 17.78 * cm;                                                                                                                                                                                                                                                                                                                                                                                                                                                       // minimum bending diameter of the cable (14 inches)
     const G4double lDistTubeToRope = ((mRopeStartingPoint) + mRopeRMax * (1 - cos(mRopeRotationAngleX)) - ((mOM->mGlassOutRad + lMushHeight - lTubeRadius) * cos(mPlugAngle) + lTubeLength * sin(mPlugAngle) + mOM->mCylHigh) * tan(mRopeRotationAngleX) - (mOM->mGlassOutRad + lMushHeight - lTubeRadius) * sin(mPlugAngle) + lTubeLength * cos(mPlugAngle)) / cos(mRopeRotationAngleX) - (lWireThickness / 2 + mRopeRMax) - mRopeDz * (1 - cos(mRopeRotationAngleX)) * sin(mRopeRotationAngleX); // Modified distance between start of first torus and end of second torus
 
-    const G4double lPCA2angle = acos((0.5) * (-lDistTubeToRope / lWireRadius + sin(mPlugAngle + mRopeRotationAngleX) + 1.));
+    const G4double lPCA2angle = acos((0.5) * (lDistTubeToRope / lWireRadius + sin(mPlugAngle + mRopeRotationAngleX) + 1.));
     const G4double lPCA1angle = 90. * deg + lPCA2angle - (mPlugAngle + mRopeRotationAngleX);
     const G4double lExtnLength = 30. * 2.54 * cm - lWireRadius * (lPCA1angle + lPCA2angle) - lTubeLength;