JeffersonLab · N-Plx · Jun 9, 2026 · Jun 4, 2026 · Jun 4, 2026 · Jun 4, 2026
@@ -28,6 +28,7 @@ public class ATOFHit {
     private int associatedClusterIndex;
     int idTDC;
     private IndexedTable atofTimeOffsetsTable;
+    private int Run;
 
 
     public int getSector() {
@@ -199,25 +200,13 @@ public final int convertTdcToTime() {
         //Time offsets
         if (atofTimeOffsetsTable == null) return 0;
         int order0 = 0;
+        //global t0
         double t0  = atofTimeOffsetsTable.getDoubleValue("t0",                  this.sector, this.layer, this.component, order0);
+        //time difference between bars up and down channels
         double tud = atofTimeOffsetsTable.getDoubleValue("upstream_downstream", this.sector, this.layer, this.component, order0);
-        //The rest of the constants are not used for now
-        /*double twb = timeOffsets[2];
-        double xtra1 = timeOffsets[3];
-        double xtra2 = timeOffsets[4];*/
-
-        //TW corrections TO BE IMPLEMENTED
-        /*
-        double[] timeWalks = CalibrationConstantsLoader.ATOF_TIME_WALK.get(key);
-	double tw0 = timeWalks[0];
-        double tw1 = timeWalks[1];
-        double tw2 = timeWalks[2];
-        double tw3 = timeWalks[3];
-        double dtw0 = timeWalks[4];
-        double dtw1 = timeWalks[5];
-        double dtw2 = timeWalks[6];
-        double dtw3 = timeWalks[7];
-        double chi2ndf = timeWalks[8];*/
+        //run-dependent correction for the radiation damage
+        double slope     = atofTimeOffsetsTable.getDoubleValue("extra1", this.sector, this.layer, this.component, order0);
+        double intersect = atofTimeOffsetsTable.getDoubleValue("extra2", this.sector, this.layer, this.component, order0);
 
         double veff, distance_to_sipm, timeOffset;
         if (null == this.type) {
@@ -229,6 +218,7 @@ public final int convertTdcToTime() {
                     veff = Parameters.VEFF;
                     //Wedge hits are placed at the center of wedges and sipm at their top
                     distance_to_sipm = Parameters.WEDGE_THICKNESS / 2.;
+                    //Run-dependent slope correction
                     timeOffset = - t0;
                 }
                 case "bar up" -> {
@@ -254,8 +244,10 @@ public final int convertTdcToTime() {
                 }
             }
         }
+        //Run-dependent slope correction
+        double dt = (this.Run - intersect) * slope;
         //Hit time.
-        this.time = this.time - distance_to_sipm / veff + timeOffset;
+        this.time = this.time - distance_to_sipm / veff + timeOffset - dt;
         return 0;
     }
 
@@ -396,11 +388,14 @@ public double getPhi() {
      * @param order Order of the hit.
      * @param tdc TDC value.
      * @param tot ToT value.
+     * @param startTime event start time
      * @param atof Detector object representing the atof, used to calculate
      * spatial coordinates.
+     * @param atofTimeOffsetsTable time offset constants
+     * @param Run for run-dependent slope correction
      */
     public ATOFHit(int sector, int layer, int component, int order, int tdc, int tot, Float startTime, Detector atof,
-                   IndexedTable atofTimeOffsetsTable) {
+                   IndexedTable atofTimeOffsetsTable, int Run) {
         this.sector = sector;
         this.layer = layer;
         this.component = component;
@@ -410,6 +405,8 @@ public ATOFHit(int sector, int layer, int component, int order, int tdc, int tot
         this.startTime = startTime;
         this.atofTimeOffsetsTable = atofTimeOffsetsTable;
         this.isInACluster = false;
+
+        this.Run = Run;
 
         this.makeType();
         this.convertTdcToTime();

@@ -1,7 +1,5 @@
 package org.jlab.rec.atof.hit;
 
-import org.jlab.io.base.DataEvent;
-import org.jlab.io.hipo.HipoDataSource;
 import org.jlab.rec.atof.constants.Parameters;
 import org.jlab.utils.groups.IndexedTable;
 

@@ -63,10 +63,15 @@ public void setWedgeHits(ArrayList<ATOFHit> wedge_hits) {
      * @param event the {@link DataEvent} containing hits.
      * @param atof the {@link Detector} representing the atof geometry to match
      * the sector/layer/component to x/y/z.
+     * @param startTime the event start time
+     * @param atofTimeOffsetsTable the constant table with time offsets
+     * @param atofEffectiveVelocityTable the constant table with effective velocity
+     * @param Run the run number corresponding to the event for run-dependent slope correction
      */
     public void findHits(DataEvent event, Detector atof, Float startTime,
                          IndexedTable atofTimeOffsetsTable,
-                         IndexedTable atofEffectiveVelocityTable) {
+                         IndexedTable atofEffectiveVelocityTable,
+                         int Run) {
         //For each event a list of bar hits and a list of wedge hits are filled
         this.barHits.clear();
         this.wedgeHits.clear();
@@ -92,7 +97,7 @@ public void findHits(DataEvent event, Detector atof, Float startTime,
             int tot = bank.getInt("ToT", i);
 
             //Building a Hit
-            ATOFHit hit = new ATOFHit(sector, layer, component, order, tdc, tot, startTime, atof, atofTimeOffsetsTable);
+            ATOFHit hit = new ATOFHit(sector, layer, component, order, tdc, tot, startTime, atof, atofTimeOffsetsTable, Run);
             if (hit.getEnergy() < 0.01) {
                 continue; //energy threshold
             }

@@ -1,6 +1,5 @@
 package org.jlab.service.alert;
 
-import ai.djl.repository.zoo.ZooModel;
 import ai.djl.translate.TranslateException;
 import java.io.File;
 import java.util.ArrayList;
@@ -396,35 +395,33 @@ public boolean processDataEventUser(DataEvent event) {
                 continue;
             }
 
-            ATOFHit hit_pred = new ATOFHit(sector_pred, layer_pred, wedge_pred, 0, 0, 0, 0f, ATOF, null);
+            ATOFHit hit_pred = new ATOFHit(sector_pred, layer_pred, wedge_pred, 0, 0, 0, 0f, ATOF, null, 0);
             double pred_x = hit_pred.getX();
             double pred_y = hit_pred.getY();
             double pred_z = hit_pred.getZ();
 
             double threshold = 20.0;
             double minDistanceSquared = threshold * threshold;
 
-            ATOFHit matchAtofHit = null; // Could be used later
             int matchHitId = -1;
 
             for (int k = 0; k < bank_ATOFHits.rows(); k++) {
                 int component = bank_ATOFHits.getInt("component", k);
                 if (component == 10) continue;
 
-                int sector = bank_ATOFHits.getInt("sector", k);
-                int layer = bank_ATOFHits.getInt("layer", k);
+                double hitX = bank_ATOFHits.getFloat("x", k);
+                double hitY = bank_ATOFHits.getFloat("y", k);
+                double hitZ = bank_ATOFHits.getFloat("z", k);
 
-                ATOFHit hit = new ATOFHit(sector, layer, component, 0, 0, 0, 0f, ATOF, null);
-
-                double dx = pred_x - hit.getX();
-                double dy = pred_y - hit.getY();
-                double dz = pred_z - hit.getZ();
+                double dx = pred_x - hitX;
+                double dy = pred_y - hitY;
+                double dz = pred_z - hitZ;
 
                 double distanceSquared = dx * dx + dy * dy + dz * dz;
 
                 if (distanceSquared < minDistanceSquared) {
                     minDistanceSquared = distanceSquared;
-                    matchAtofHit = hit;
+                    //matchAtofHit = hit;
                     matchHitId = bank_ATOFHits.getInt("id", k);
                 }
             }

@@ -12,6 +12,7 @@
 import org.jlab.detector.calib.utils.DatabaseConstantProvider;
 import org.jlab.geom.base.Detector;
 import org.jlab.geom.detector.alert.ATOF.AlertTOFFactory;
+import org.jlab.io.hipo.HipoDataSource;
 import org.jlab.rec.atof.banks.RecoBankWriter;
 import org.jlab.rec.atof.cluster.ATOFCluster;
 import org.jlab.rec.atof.cluster.ClusterFinder;
@@ -94,25 +95,10 @@ public boolean processDataEventUser(DataEvent event) {
             atofTimeOffsetsTable       = this.getConstantsManager().getConstants(newRun, "/calibration/alert/atof/time_offsets");
             Run = newRun;
         }
-
-        ////Do we need to read the event vx,vy,vz?
-        ////If not, this part can be moved in the initialization of the engine.
-        //double eventVx=0,eventVy=0,eventVz=0; //They should be in CM
-        ////Track Projector Initialisation with b field
-        //Swim swim = new Swim();
-        //float magField[] = new float[3];
-        //swim.BfieldLab(eventVx, eventVy, eventVz, magField); 
-        //this.b = Math.sqrt(Math.pow(magField[0],2) + Math.pow(magField[1],2) + Math.pow(magField[2],2));
-
-        ///// \todo move this to ALERTEngine
-        //TrackProjector projector = new TrackProjector();
-        //projector.setB(this.b);
-        //projector.projectTracks(event);
-        //rbc.appendMatchBanks(event, projector.getProjections());
 
         //Hit finder init
         HitFinder hitfinder = new HitFinder();
-        hitfinder.findHits(event, ATOF, startTime, atofTimeOffsetsTable, atofEffectiveVelocityTable);
+        hitfinder.findHits(event, ATOF, startTime, atofTimeOffsetsTable, atofEffectiveVelocityTable, Run);
         ArrayList<ATOFHit> WedgeHits = hitfinder.getWedgeHits();
         ArrayList<BarHit> BarHits = hitfinder.getBarHits();
         //Exit if hit lists are empty
@@ -165,5 +151,8 @@ else if ("false".equals(this.getEngineConfigString(useStartTimeString)))
     }
 
     public static void main(String arg[]) {
+
     }
 }
+
+