From a20f3d46e8ed0d3a1637c2a55c0550675e1cc126 Mon Sep 17 00:00:00 2001 From: Ryan Tang Date: Wed, 31 Jan 2024 19:54:25 -0500 Subject: [PATCH] fixed the CalTrack --- ClassAnasen.h | 6 +++--- 1 file changed, 3 insertions(+), 3 deletions(-) diff --git a/ClassAnasen.h b/ClassAnasen.h index 7140e72..7e313fe 100644 --- a/ClassAnasen.h +++ b/ClassAnasen.h @@ -472,11 +472,11 @@ inline void ANASEN::CalTrack(TVector3 sx3Pos, int anodeID, int cathodeID, bool v trackPos = sx3Pos; - TVector3 n1 = (P1[anodeID].first - P1[anodeID].second).Cross((sx3Pos - P1[anodeID].second)); - TVector3 n2 = (Q1[anodeID].first - Q1[anodeID].second).Cross((sx3Pos - Q1[anodeID].second)); + TVector3 n1 = (P1[anodeID].first - P1[anodeID].second).Cross((sx3Pos - P1[anodeID].second)).Unit(); + TVector3 n2 = (Q1[cathodeID].first - Q1[cathodeID].second).Cross((sx3Pos - Q1[cathodeID].second)).Unit(); // if the handiness of anode and cathode revered, it should be n2 cross n1 - trackVec = (n1.Cross(n2)).Unit(); + trackVec = (n2.Cross(n1)).Unit(); if( verbose ) printf("Theta, Phi = %f, %f \n", trackVec.Theta() *TMath::RadToDeg(), trackVec.Phi()*TMath::RadToDeg());