different electron pid for zero-torus

baltzell · baltzell · commit 51abf8e50101 · 2026-02-02T11:20:56.000-05:00
diff --git a/reconstruction/eb/src/main/java/org/jlab/rec/eb/EBUtil.java b/reconstruction/eb/src/main/java/org/jlab/rec/eb/EBUtil.java
@@ -99,6 +99,27 @@ public static boolean isSimpleElectron(DetectorParticle p,EBCCDBConstants ccdb)
        
         return true;
     }
+     
+    /**
+     * Perform a basic true/false identification for electrons.
+     */
+    public static boolean isZeroFieldElectron(DetectorParticle p,EBCCDBConstants ccdb) {
+
+        // require ECAL:
+        final int sector = p.getSector(DetectorType.ECAL);
+        if (sector<1) return false;
+       
+        // requre HTCC photoelectrons:
+        final double nphe = p.getNphe(DetectorType.HTCC);
+        if (nphe < ccdb.getDouble(EBCCDBEnum.HTCC_NPHE_CUT)) return false;
+
+        // require PCAL minimum energy:
+        final double minPcalEnergy = ccdb.getSectorDouble(EBCCDBEnum.ELEC_PCAL_min_energy,sector);
+        final double pcalEnergy = p.getEnergy(DetectorType.ECAL,DetectorLayer.PCAL);
+        if (pcalEnergy < minPcalEnergy) return false;
+      
+        return true;
+    }
 
     /**
      * Calculate timing resolution from EventBuilder constants:
diff --git a/reconstruction/eb/src/main/java/org/jlab/service/eb/EBAnalyzer.java b/reconstruction/eb/src/main/java/org/jlab/service/eb/EBAnalyzer.java
@@ -350,7 +350,7 @@ public void assignBetas(DetectorEvent event,final boolean useStartTimeFromFT){
     }
     public void assignPids(DetectorEvent event,final boolean useStartTimeFromFT) {
 
-        PIDHypothesis pidHyp = new PIDHypothesis();
+        PIDHypothesis pidHyp = new PIDHypothesis(event.getEventHeader().getTorus());
         pidHyp.setEvent(event);
         pidHyp.setUseStartTimeFromFT(useStartTimeFromFT);
 
@@ -382,12 +382,13 @@ else if (p.getCharge()<0) {
 
     public class PIDHypothesis {
 
+        private float torus = 0;
         private int theoryPID = -1;
         private double PIDquality = 0.0;
         private DetectorEvent event;
         private boolean useStartTimeFromFT = false;
 
-        public PIDHypothesis() {}
+        public PIDHypothesis(float t) {torus=t;}
 
         public void setEvent(DetectorEvent e) {event = e;}
 
@@ -405,7 +406,9 @@ public void PIDMatch(DetectorParticle p, int pid) {
             
             final boolean pidFromTimingCheck = pid==pidFromTiming && p.getTheoryBeta(pid)>0;
            
-            final boolean isElectron = EBUtil.isSimpleElectron(p,ccdb);
+            final boolean isElectron = Math.abs(torus)>1e-8 ? 
+                EBUtil.isSimpleElectron(p,ccdb) :
+                EBUtil.isZeroFieldElectron(p, ccdb);
             
             final boolean htccSignalCheck = p.getNphe(DetectorType.HTCC)>ccdb.getDouble(EBCCDBEnum.HTCC_NPHE_CUT);
             final boolean ltccSignalCheck = p.getNphe(DetectorType.LTCC)>ccdb.getDouble(EBCCDBEnum.LTCC_NPHE_CUT);
