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());
+ }
+ }
+ }
+ }
+}
