Fix: AHDC: use ahdcExtractor to process simulation (#607)

* Add ahdcExtractor in AHDCEngine to be able to process simulation.

* Add simulation flag in HitReader to bypass criteria for a hit in simulation.

Author: mathieuouillon

reconstruction/alert/src/main/java/org/jlab/rec/ahdc/Hit/HitReader.java

@@ -12,8 +12,10 @@ public class HitReader {
 
 	private ArrayList<Hit>     _AHDCHits;
 	private ArrayList<TrueHit> _TrueAHDCHits;
+	private boolean sim = false;
 
-        public HitReader(DataEvent event, AlertDCDetector detector, boolean simulation) {
+	public HitReader(DataEvent event, AlertDCDetector detector, boolean simulation) {
+		sim = simulation;
 		fetch_AHDCHits(event, detector);
 		if (simulation) fetch_TrueAHDCHits(event);
 	}
@@ -52,7 +54,7 @@ public final void fetch_AHDCHits(DataEvent event, AlertDCDetector detector) {
 				double ped_max = rawHitCuts[7];
 				//System.out.println("t_min : " + t_min + " t_max : " + t_max + " tot_min : " + tot_min + " tot_max : " + tot_max + " adc_min : " + adc_min + " adc_max : " + adc_max + " ped_min : " + ped_min + " ped_max : " + ped_max);
 				// Apply these cuts
-				if ((adc >= adc_min) && (adc <= adc_max) && (leadingEdgeTime >= t_min) && (leadingEdgeTime <= t_max) && (timeOverThreshold >= tot_min) && (timeOverThreshold <= tot_max) && (adcOffset >= ped_min) && (adcOffset <= ped_max)) {
+				if (((adc >= adc_min) && (adc <= adc_max) && (leadingEdgeTime >= t_min) && (leadingEdgeTime <= t_max) && (timeOverThreshold >= tot_min) && (timeOverThreshold <= tot_max) && (adcOffset >= ped_min) && (adcOffset <= ped_max)) || sim) {
 					// Retrieve others CCDB
 					double[] timeOffsets = CalibrationConstantsLoader.AHDC_TIME_OFFSETS.get( key_value );
 					double[] time2distance = CalibrationConstantsLoader.AHDC_TIME_TO_DISTANCE.get( 10101 ); // the time to distance table has only one row ! (10101 is its only key)

reconstruction/alert/src/main/java/org/jlab/service/ahdc/AHDCEngine.java

@@ -28,6 +28,7 @@
 import org.jlab.geom.detector.alert.AHDC.AlertDCDetector;
 import org.jlab.geom.detector.alert.AHDC.AlertDCFactory;
 import org.jlab.rec.constants.CalibrationConstantsLoader;
+import org.jlab.detector.pulse.ModeAHDC;
 
 /** AHDCEngine reconstruction service.
  *
@@ -61,7 +62,8 @@ public class AHDCEngine extends ReconstructionEngine {
     private Mode mode = Mode.CV_Track_Finding;
 
     private AlertDCDetector factory = null;
-    
+    private ModeAHDC ahdcExtractor = new ModeAHDC();
+
     public AHDCEngine() {
         super("ALERT", "ouillon", "1.0.1");
     }
@@ -124,6 +126,9 @@ public boolean processDataEvent(DataEvent event) {
 	    simulation = true;
 	}
 
+        // TODO: this code should perhaps be in the if statement MC::Particle
+        ahdcExtractor.update(30, null, event, "AHDC::wf", "AHDC::adc");
+
         if (event.hasBank("RUN::config")) {
             DataBank bank = event.getBank("RUN::config");
             runNo          = bank.getInt("run", 0);