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;