SlicDiagnostics/src/org/lcsim/slic/diagnostics

diff -u -r1.26 -r1.27
--- MCParticlePlots.java	1 Apr 2006 00:30:02 -0000	1.26
+++ MCParticlePlots.java	5 Apr 2006 00:20:35 -0000	1.27
@@ -22,7 +22,7 @@

  * A set of MCParticle plots, including plots for each particle type by PDGID.
  * 
  * @author jeremym

- * @version $Id: MCParticlePlots.java,v 1.26 2006/04/01 00:30:02 jeremy Exp $

+ * @version $Id: MCParticlePlots.java,v 1.27 2006/04/05 00:20:35 jeremy Exp $

  */
 class MCParticlePlots extends AbstractPlots
 {

@@ -39,6 +39,7 @@

     ICloud1D _mcpNegativeCount;
     ICloud1D _mcpChargedCount;
     ICloud1D _endpointDau;

+    ICloud1D _prodTimelog10;

     ICloud1D _prodTime;
 
     // Map of ParticleTypeInfo

@@ -103,8 +104,12 @@

         _mcpChargedCount.annotation().addItem("xAxisLabel", "Number of Charged Particles");
         _mcpChargedCount.annotation().addItem("yAxisLabel", "Number of Events");
 

+        _prodTimelog10 = aida().cloud1D("log10(Production Time)");
+        _prodTimelog10.annotation().addItem("xAxisLabel", "Production Time (log10(nanoseconds))");
+        _prodTimelog10.annotation().addItem("yAxisLabel", "Number of Particles");
+        

         _prodTime = aida().cloud1D("Production Time");

-        _prodTime.annotation().addItem("xAxisLabel", "Production Time (log10(seconds))");

+        _prodTime.annotation().addItem("xAxisLabel", "Production Time (nanoseconds)");

         _prodTime.annotation().addItem("yAxisLabel", "Number of Particles");
                 
         for (MCParticleFlag flag : MCParticleFlags.allFlags())

@@ -210,9 +215,19 @@

 
             // Time
             double time = particle.getProductionTime();

-
-            // Fill production time

+            
+            // Fill prod time

             _prodTime.fill(time);

+            
+            // Fill log10 prod time
+            if (time != 0)
+            {
+                _prodTimelog10.fill(log10(time));
+            }
+            else 
+            {
+                _prodTimelog10.fill(0);
+            }

 
             // End Point
             boolean hasEndpoint = true;