Print

Print


Commit in SlicDiagnostics on MAIN
TODO+11.4 -> 1.5
src/org/lcsim/slic/diagnostics/CalorimeterHitPlots.java+12-11.29 -> 1.30
+13-1
2 modified files
JM: Add total corrected energy plot.

SlicDiagnostics
TODO 1.4 -> 1.5
diff -u -r1.4 -r1.5
--- TODO	30 Mar 2006 01:09:43 -0000	1.4
+++ TODO	6 Jul 2006 01:09:22 -0000	1.5
@@ -1,5 +1,6 @@
 -overlay 
   -takes input AIDA files and makes output AIDA file with all matching plots overlayed
+-create data points of means for all clouds, h1ds
 -comparison
   -do difference plots (Max's stuff)
 -multiple AIDA files

SlicDiagnostics/src/org/lcsim/slic/diagnostics
CalorimeterHitPlots.java 1.29 -> 1.30
diff -u -r1.29 -r1.30
--- CalorimeterHitPlots.java	30 Mar 2006 01:09:44 -0000	1.29
+++ CalorimeterHitPlots.java	6 Jul 2006 01:09:22 -0000	1.30
@@ -27,7 +27,7 @@
  * collections of SimCalorimeterHits are done by GenericHitPlots.
  * 
  * @author jeremym
- * @version $Id: CalorimeterHitPlots.java,v 1.29 2006/03/30 01:09:44 jeremy Exp $
+ * @version $Id: CalorimeterHitPlots.java,v 1.30 2006/07/06 01:09:22 jeremy Exp $
  */
 class CalorimeterHitPlots extends AbstractPlots
 {
@@ -39,6 +39,7 @@
     ICloud1D _hitCorrElog10;
     ICloud2D _hitRawVsCorrE;
     ICloud1D _mcpMissing;
+    ICloud1D _totCorrE;
 
     // Type of segmentation
     SegmentationType _segmentationType = SegmentationType.NONE;
@@ -95,6 +96,10 @@
         _mcpMissing = aida().cloud1D("MCParticle: Contribution Missing");
         _mcpMissing.annotation().addItem("xAxisLabel", "Number of Hits Missing an MCParticle Contribution per Event");
         _mcpMissing.annotation().addItem("yAxisLabel", "Number of Events");
+
+        _totCorrE = aida().cloud1D("Hit: Total Corrected Energy");
+        _totCorrE.annotation().addItem("xAxisLabel", "Event Energy (GeV)");
+        _totCorrE.annotation().addItem("yAxisLabel", "Number of Events");
     }
 
     /** Set the segmentation type based on the LCMetaData. */
@@ -161,6 +166,7 @@
     {
         // Get corrected energy
         boolean haveSF = true;
+        double totCorrE = 0;
         try
         {
             hits.get(0).getCorrectedEnergy();
@@ -201,6 +207,9 @@
 
                     // Plot raw vs corrected E
                     _hitRawVsCorrE.fill(hit.getRawEnergy(), hit.getCorrectedEnergy());
+
+                    // Add to corrected energy total.
+                    totCorrE += hit.getCorrectedEnergy();
                 }
 
                 // MCParticle Plots
@@ -226,6 +235,8 @@
             // Plot missing MCParticle contribution count
             _mcpMissing.fill(mcpMissingCount);
         }
+
+        _totCorrE.fill(totCorrE);
     }
 
     public abstract class SegmentationPlots
CVSspam 0.2.8