[bug] reset replay.tof when tof exceeds tend, like fangq/mcxcl#57

fangq · Jun 11, 2024 · 1e7d0f1 · 1e7d0f1
1 parent 026eebf
commit 1e7d0f1
Showing 1 changed file with 2 additions and 0 deletions.
diff --git a/src/mcx_utils.c b/src/mcx_utils.c
@@ -3653,6 +3653,7 @@ void mcx_replayprep(int* detid, float* ppath, History* his, Config* cfg) {
             }
 
             cfg->replay.weight[cfg->nphoton] = 1.f;
+            cfg->replay.tof[cfg->nphoton] = 0.f;
             cfg->replay.detid[cfg->nphoton] = (detid != NULL) ? detid[i] : 1;
 
             for (j = 0; j < his->maxmedia; j++) {
@@ -3976,6 +3977,7 @@ void mcx_loadseedfile(Config* cfg) {
                 }
 
                 cfg->replay.weight[cfg->nphoton] = 1.f;
+                cfg->replay.tof[cfg->nphoton] = 0.f;
                 cfg->replay.detid[cfg->nphoton] = (hasdetid) ? (int)(ppath[i * his.colcount]) : 1;
 
                 for (j = hasdetid; j < his.maxmedia + hasdetid; j++) {