From e7c893a6df1a5ba033e482fb5fce7be90e0ceb6c Mon Sep 17 00:00:00 2001
From: Clenain <lenain@lpccaen.in2p3.fr>
Date: Wed, 2 Dec 2020 11:26:39 +0100
Subject: [PATCH] * debugging NPCrossTalk
---
NPLib/Neutron/NPCrossTalk.cxx | 5 +++--
1 file changed, 3 insertions(+), 2 deletions(-)
diff --git a/NPLib/Neutron/NPCrossTalk.cxx b/NPLib/Neutron/NPCrossTalk.cxx
index 7dc77bf11..f0aedf6ba 100644
--- a/NPLib/Neutron/NPCrossTalk.cxx
+++ b/NPLib/Neutron/NPCrossTalk.cxx
@@ -80,10 +80,10 @@ vector<int> CrossTalk::ComputeCrossTalk(){
Dist = sqrt((x2-x1)*(x2-x1) + (y2-y1)*(y2-y1) + (z2-z1)*(z2-z1));
if(Dist < coef*dR1){
if(ID_ClustHit[jj] < ID_ClustHit[j]){
- ID_ClustHit[j] == ID_ClustHit[jj];
+ ID_ClustHit[j] = ID_ClustHit[jj];
}
else if(ID_ClustHit[jj] > ID_ClustHit[j]){
- ID_ClustHit[jj] == ID_ClustHit[j];
+ ID_ClustHit[jj] = ID_ClustHit[j];
}
}
}
@@ -134,6 +134,7 @@ vector<int> CrossTalk::ComputeCrossTalk(){
x2 = (*HitX)[mapOfFirstN[i]], y2 = (*HitY)[mapOfFirstN[i]], z2 = (*HitZ)[mapOfFirstN[i]], t2 = (*HitT)[mapOfFirstN[i]];
Dist = sqrt((x2-x1)*(x2-x1) + (y2-y1)*(y2-y1) + (z2-z1)*(z2-z1));
Dmax = (t2-t1)*v_n;
+ cout << Dmax << endl;
if(Dist < Dmax){
if(t1 < t2){
CrossTalk.push_back(j);
--
GitLab