add utilities for rebuilding scalers

baltzell · baltzell · commit 054a212d811c · 2024-09-24T18:13:06.000+02:00
diff --git a/common-tools/clas-detector/src/main/java/org/jlab/detector/scalers/DaqScalers.java b/common-tools/clas-detector/src/main/java/org/jlab/detector/scalers/DaqScalers.java
@@ -50,7 +50,10 @@ public class DaqScalers {
     public StruckScalers struck=null;
 
     private long timestamp=0;
-    public void setTimestamp(long timestamp) { this.timestamp=timestamp; }
+    public DaqScalers setTimestamp(long timestamp) {
+        this.timestamp=timestamp;
+        return this;
+    }
     public long getTimestamp(){ return this.timestamp; }
 
     /**
@@ -137,6 +140,37 @@ public static DaqScalers create(Bank rawScalerBank,IndexedTable fcupTable,Indexe
         return DaqScalers.create(rawScalerBank,fcupTable,slmTable,helTable,dsc2.getGatedClockSeconds());
     }
 
+    /**
+     * 
+     * @param conman
+     * @param runConfig
+     * @param rawScaler
+     * @return 
+     */
+    public static DaqScalers create(ConstantsManager conman, Bank runConfig, Bank rawScaler) {
+        int run = runConfig.getInt("run", 0);
+        IndexedTable fcup = conman.getConstants(run, "/runcontrol/fcup");
+        IndexedTable slm = conman.getConstants(run, "/runcontrol/slm");
+        IndexedTable hel = conman.getConstants(run, "/runcontrol/helicity");
+        IndexedTable dsc = conman.getConstants(run, "/daq/config/scalers/dsc1");
+        if (fcup!=null && dsc!=null) {
+            if (dsc.getIntValue("frequency",0,0,0) < Dsc2Scaler.MAX_DSC2_CLOCK_FREQ) {
+                return DaqScalers.create(rawScaler, fcup, slm, hel, dsc)
+                    .setTimestamp(runConfig.getLong("timestamp", 0));
+            }
+            else {
+                try {
+                    Time rst = (Time)conman.getRcdbConstant(run,"run_start_time").getValue();
+                    Date uet = new Date(runConfig.getInt("unixtime",0)*1000L);
+                    return DaqScalers.create(rawScaler, fcup, slm, hel, rst, uet)
+                        .setTimestamp(runConfig.getLong("timestamp", 0));
+                }
+                catch (Exception e) {}
+            }
+        }
+        return null;
+    }
+
     /**
      * @param schema bank schema
      * @return RUN::scaler banks
@@ -253,7 +287,7 @@ public static List<Bank> createBanks(int runnumber, SchemaFactory schema, Event
         IndexedTable hel = conman.getConstants(runnumber, "/runcontrol/helicity");
         IndexedTable dsc = conman.getConstants(runnumber, "/daq/config/scalers/dsc1");
 
-        if (dsc.getIntValue("frequency", 0,0,0) < 2e5) {
+        if (dsc.getIntValue("frequency", 0,0,0) < Dsc2Scaler.MAX_DSC2_CLOCK_FREQ) {
             ret.addAll(createBanks(schema,rawScaler,fcup, slm, hel, dsc));
         }
         else {
@@ -274,5 +308,6 @@ public static List<Bank> createBanks(int runnumber, SchemaFactory schema, Event
 
         return ret;
     }
+
 }
 
diff --git a/common-tools/clas-detector/src/main/java/org/jlab/detector/scalers/DaqScalersSequence.java b/common-tools/clas-detector/src/main/java/org/jlab/detector/scalers/DaqScalersSequence.java
@@ -5,6 +5,8 @@
 import java.util.Arrays;
 import java.util.Collections;
 import java.util.Comparator;
+import java.util.logging.Logger;
+import org.jlab.detector.calib.utils.ConstantsManager;
 
 import org.jlab.jnp.hipo4.io.HipoReader;
 import org.jlab.jnp.hipo4.data.Event;
@@ -27,6 +29,8 @@ public class DaqScalersSequence implements Comparator<DaqScalers> {
     
     private Bank rcfgBank=null;
   
+    static final Logger logger = Logger.getLogger(DaqScalersSequence.class.getName());
+    
     public static class Interval {
         private DaqScalers previous = null;
         private DaqScalers next = null;
@@ -187,6 +191,7 @@ public Interval getInterval(Event event1, Event event2) {
      * @return  sequence
      */
     public static DaqScalersSequence readSequence(List<String> filenames) {
+        logger.info("DaqScalersSequence::  Reading scaler sequence from "+String.join(",", filenames));
        
         DaqScalersSequence seq=new DaqScalersSequence();
 
@@ -229,6 +234,38 @@ public static DaqScalersSequence readSequence(List<String> filenames) {
         
         return seq;
     }
+   
+    /**
+     * 
+     * @param tags
+     * @param conman
+     * @param filenames
+     * @return 
+     */
+    public static DaqScalersSequence rebuildSequence(int tags, ConstantsManager conman, List<String> filenames) {
+        logger.info("DaqScalersSequence::  Rebuilding scaler sequence from "+String.join(",", filenames));
+        DaqScalersSequence seq=new DaqScalersSequence();
+        for (String filename : filenames) {
+            HipoReader reader = new HipoReader();
+            reader.setTags(tags);
+            reader.open(filename);
+            SchemaFactory schema = reader.getSchemaFactory();
+            if (seq.rcfgBank==null)
+                seq.rcfgBank = new Bank(schema.getSchema("RUN::config"));
+            while (reader.hasNext()) {
+                Event event=new Event();
+                Bank scaler=new Bank(schema.getSchema("RAW::scaler"));
+                Bank config=new Bank(schema.getSchema("RUN::config"));
+                reader.nextEvent(event);
+                event.read(scaler);
+                event.read(config);
+                if (scaler.getRows()<1 || config.getRows()<1) continue;
+                seq.add(DaqScalers.create(conman, config, scaler));
+            }
+            reader.close();
+        }
+        return seq;
+    }
     
     public static void main(String[] args) {
         
diff --git a/common-tools/clas-detector/src/main/java/org/jlab/detector/scalers/Dsc2Scaler.java b/common-tools/clas-detector/src/main/java/org/jlab/detector/scalers/Dsc2Scaler.java
@@ -16,6 +16,8 @@
  */
 public class Dsc2Scaler extends DaqScaler{
 
+    public static final long MAX_DSC2_CLOCK_FREQ = (long)2e5;
+
     private static final boolean GATEINVERTED=true;
 
     // DSC has TRG and TDC thresholds, we use only TDC here:
