From 60cf235eec76e883d7b4c3afbf6c1fea65d59655 Mon Sep 17 00:00:00 2001
From: Morfouace <pierre.morfouace@gmail.com>
Date: Fri, 30 Sep 2022 10:37:10 +0200
Subject: [PATCH] Updating GladFieldMap.cxx: fixing cm->mm issue
---
NPLib/Detectors/Sofia/GladFieldMap.cxx | 11 +++++++++--
1 file changed, 9 insertions(+), 2 deletions(-)
diff --git a/NPLib/Detectors/Sofia/GladFieldMap.cxx b/NPLib/Detectors/Sofia/GladFieldMap.cxx
index 0126bebb8..5da17e50f 100644
--- a/NPLib/Detectors/Sofia/GladFieldMap.cxx
+++ b/NPLib/Detectors/Sofia/GladFieldMap.cxx
@@ -332,6 +332,13 @@ void GladFieldMap::LoadMap(string filename) {
ifile >> m_y_min >> m_y_max >> m_Ny;
ifile >> m_z_min >> m_z_max >> m_Nz;
+ m_x_min = m_x_min*10;
+ m_x_max = m_x_max*10;
+ m_y_min = m_y_min*10;
+ m_y_max = m_y_max*10;
+ m_z_min = m_z_min*10;
+ m_z_max = m_z_max*10;
+
unsigned int count=0;
int index = 0;
//while(!ifile.eof()){
@@ -573,8 +580,8 @@ bool GladFieldMap::IsInside(TVector3 localpoint, int& ix, int& iy, int& iz, doub
bool isIn=false;
double x = localpoint.X();
- double y = localpoint.X();
- double z = localpoint.X();
+ double y = localpoint.Y();
+ double z = localpoint.Z();
if(!(x>=m_x_min && x<=m_x_max && y>=m_y_min && y<=m_y_max && z>=m_z_min && z<=m_z_max)){
ix = iy = iz = 0;
--
GitLab