Author: [log in to unmask]
Date: Tue Nov 25 16:28:14 2014
New Revision: 1589

Log:
Add Drivers for plotting ECAL channel data.

Added:
    java/trunk/monitoring-drivers/src/main/java/org/hps/monitoring/ecal/plots/EcalChannelADCSumPlotsDriver.java   (with props)
    java/trunk/monitoring-drivers/src/main/java/org/hps/monitoring/ecal/plots/EcalChannelEnergyPlotsDriver.java   (with props)

Added: java/trunk/monitoring-drivers/src/main/java/org/hps/monitoring/ecal/plots/EcalChannelADCSumPlotsDriver.java
 =============================================================================
--- java/trunk/monitoring-drivers/src/main/java/org/hps/monitoring/ecal/plots/EcalChannelADCSumPlotsDriver.java	(added)
+++ java/trunk/monitoring-drivers/src/main/java/org/hps/monitoring/ecal/plots/EcalChannelADCSumPlotsDriver.java	Tue Nov 25 16:28:14 2014
@@ -0,0 +1,104 @@
+package org.hps.monitoring.ecal.plots;
+
+import hep.aida.IAnalysisFactory;
+import hep.aida.IHistogram1D;
+import hep.aida.IPlotter;
+import hep.aida.IPlotterFactory;
+import hep.aida.IPlotterStyle;
+
+import java.util.ArrayList;
+import java.util.HashSet;
+import java.util.List;
+import java.util.Set;
+
+import org.hps.conditions.database.TableConstants;
+import org.hps.conditions.ecal.EcalChannel;
+import org.hps.conditions.ecal.EcalChannel.EcalChannelCollection;
+import org.hps.conditions.ecal.EcalConditions;
+import org.hps.recon.ecal.EcalRawConverter;
+import org.lcsim.conditions.ConditionsManager;
+import org.lcsim.event.EventHeader;
+import org.lcsim.event.RawTrackerHit;
+import org.lcsim.geometry.Detector;
+import org.lcsim.util.Driver;
+import org.lcsim.util.aida.AIDA;
+
+/**
+ * This Driver will create a histogram for every crystal in the ECAL and plot the sum of
+ * its ADC values for the raw hit, minus the channel's pedestal value from the conditions system.
+ * @author Jeremy McCormick <[log in to unmask]>
+ */
+public class EcalChannelADCSumPlotsDriver extends Driver {
+
+    EcalConditions conditions = null;
+    EcalChannelCollection channels = null;
+    EcalRawConverter ecalRawConverter = null;
+
+    List<List<IPlotter>> plotterLists = new ArrayList<List<IPlotter>>();
+
+    AIDA aida = AIDA.defaultInstance();
+    IAnalysisFactory analysisFactory = aida.analysisFactory();
+
+    public void detectorChanged(Detector detector) {
+
+        conditions = ConditionsManager.defaultInstance().getCachedConditions(EcalConditions.class, TableConstants.ECAL_CONDITIONS).getCachedData();
+        
+        ecalRawConverter = new EcalRawConverter();
+        ecalRawConverter.setDetector(null);
+
+        channels = conditions.getChannelCollection();
+
+        Set<Integer> crates = new HashSet<Integer>();
+        Set<Integer> slots = new HashSet<Integer>();
+        Set<Integer> channels = new HashSet<Integer>();
+
+        for (EcalChannel channel : conditions.getChannelCollection()) {
+            if (channel == null) {
+                throw new RuntimeException("EcalChannel in collection is null.");
+            }
+            crates.add(channel.getCrate());
+            slots.add(channel.getSlot());
+            channels.add(channel.getChannel());
+        }
+                                
+        for (Integer crate : crates) {
+            IPlotterFactory plotterFactory = analysisFactory.createPlotterFactory("ECAL Sum ADC - Crate " + crate);
+            IPlotterStyle style = plotterFactory.createPlotterStyle();
+            style.dataStyle().lineStyle().setVisible(false);
+            style.legendBoxStyle().setVisible(false);
+            int plottersIndex = crate - 1;
+            plotterLists.add(new ArrayList<IPlotter>());
+            List<IPlotter> plotters = plotterLists.get(plottersIndex);
+            for (Integer slot : slots) {
+                IPlotter plotter = plotterFactory.create("Slot " + slot);
+                plotters.add(plotter);
+                plotter.createRegions(4, 4);
+                for (Integer channel : channels) {
+                    String adcSumPlotName = "Sum ADC Values : " + crate + " : " + slot + " : " + channel;
+                    IHistogram1D adcSumHistogram = aida.histogram1D(adcSumPlotName, 600, -6000., 6000.);
+                    plotter.region(channel).plot(adcSumHistogram, style);
+                }
+                plotter.show();
+            }
+        }
+    }
+
+    public void process(EventHeader event) {
+        if (event.hasCollection(RawTrackerHit.class, "EcalReadoutHits")) {
+            List<RawTrackerHit> hits = event.get(RawTrackerHit.class, "EcalReadoutHits");
+            for (RawTrackerHit hit : hits) {
+                EcalChannel channel = channels.findGeometric(hit.getCellID());
+                if (channel != null) {
+                    try {
+                        aida.histogram1D("Sum ADC Values : " + channel.getCrate() + " : " + channel.getSlot() + " : " + channel.getChannel())
+                            .fill(ecalRawConverter.sumADC(hit));
+                    } catch (IllegalArgumentException e) {
+                        e.printStackTrace();
+                    }
+                } else {
+                    System.err.println("EcalChannel not found for cell ID 0x" + String.format("%08d", Long.toHexString(hit.getCellID())));
+                }
+            }
+        }
+    }
+}

Added: java/trunk/monitoring-drivers/src/main/java/org/hps/monitoring/ecal/plots/EcalChannelEnergyPlotsDriver.java
 =============================================================================
--- java/trunk/monitoring-drivers/src/main/java/org/hps/monitoring/ecal/plots/EcalChannelEnergyPlotsDriver.java	(added)
+++ java/trunk/monitoring-drivers/src/main/java/org/hps/monitoring/ecal/plots/EcalChannelEnergyPlotsDriver.java	Tue Nov 25 16:28:14 2014
@@ -0,0 +1,92 @@
+package org.hps.monitoring.ecal.plots;
+
+import hep.aida.IAnalysisFactory;
+import hep.aida.IHistogram1D;
+import hep.aida.IPlotter;
+import hep.aida.IPlotterFactory;
+import hep.aida.IPlotterStyle;
+
+import java.util.ArrayList;
+import java.util.HashSet;
+import java.util.List;
+import java.util.Set;
+
+import org.hps.conditions.database.TableConstants;
+import org.hps.conditions.ecal.EcalChannel;
+import org.hps.conditions.ecal.EcalChannel.EcalChannelCollection;
+import org.hps.conditions.ecal.EcalConditions;
+import org.lcsim.conditions.ConditionsManager;
+import org.lcsim.event.CalorimeterHit;
+import org.lcsim.event.EventHeader;
+import org.lcsim.geometry.Detector;
+import org.lcsim.util.Driver;
+import org.lcsim.util.aida.AIDA;
+
+/**
+ * This Driver will create a histogram for every channel in the ECAL and plot its corrected energy.
+ * @author Jeremy McCormick <[log in to unmask]>
+ */
+public class EcalChannelEnergyPlotsDriver extends Driver {
+
+    EcalConditions conditions = null;
+    EcalChannelCollection channels = null;
+
+    List<List<IPlotter>> plotterLists = new ArrayList<List<IPlotter>>();
+
+    AIDA aida = AIDA.defaultInstance();
+    IAnalysisFactory analysisFactory = aida.analysisFactory();
+
+    public void detectorChanged(Detector detector) {
+        
+        conditions = ConditionsManager.defaultInstance().getCachedConditions(EcalConditions.class, TableConstants.ECAL_CONDITIONS).getCachedData();
+        
+        channels = conditions.getChannelCollection();
+
+        Set<Integer> crates = new HashSet<Integer>();
+        Set<Integer> slots = new HashSet<Integer>();
+        Set<Integer> channels = new HashSet<Integer>();
+
+        for (EcalChannel channel : conditions.getChannelCollection()) {
+            crates.add(channel.getCrate());
+            slots.add(channel.getSlot());
+            channels.add(channel.getChannel());
+        }
+                
+                        
+        for (Integer crate : crates) {
+            IPlotterFactory plotterFactory = analysisFactory.createPlotterFactory("ECAL Energy - Crate " + crate);
+            IPlotterStyle style = plotterFactory.createPlotterStyle();
+            style.dataStyle().lineStyle().setVisible(false);
+            style.legendBoxStyle().setVisible(false);
+            int plottersIndex = crate - 1;
+            plotterLists.add(new ArrayList<IPlotter>());
+            List<IPlotter> plotters = plotterLists.get(plottersIndex);
+            for (Integer slot : slots) {
+                IPlotter plotter = plotterFactory.create("Slot " + slot);
+                plotters.add(plotter);
+                plotter.createRegions(4, 4);
+                for (Integer channel : channels) {
+                                        
+                    String energyPlotName = "Crystal Energy : " + crate + " : " + slot + " : " + channel;
+                    IHistogram1D energyHistogram = aida.histogram1D(energyPlotName, 200, -1., 1.);
+                    
+                    plotter.region(channel).plot(energyHistogram, style);
+                }
+                plotter.show();
+            }
+        }
+    }
+
+    public void process(EventHeader event) {
+        if (event.hasCollection(CalorimeterHit.class, "EcalCalHits")) {
+            List<CalorimeterHit> hits = event.get(CalorimeterHit.class, "EcalCalHits");
+            for (CalorimeterHit hit : hits) {
+                EcalChannel channel = channels.findGeometric(hit.getCellID());
+                if (channel != null) {
+                    aida.histogram1D("Crystal Energy : " + channel.getCrate() + " : " + channel.getSlot() + " : " + channel.getChannel())
+                        .fill(hit.getCorrectedEnergy());
+                }                   
+            }
+        }
+    }
+}