lcsim-analysis/src/main/java/org/lcsim/analysis/dbd/evtgen
diff -N CalculateZSmearing.java
--- /dev/null 1 Jan 1970 00:00:00 -0000
+++ CalculateZSmearing.java 24 Aug 2012 20:25:47 -0000 1.1
@@ -0,0 +1,52 @@
+package org.lcsim.analysis.dbd.evtgen;
+
+import hep.aida.IAnalysisFactory;
+import hep.aida.IHistogram1D;
+import hep.aida.IHistogramFactory;
+import hep.aida.IPlotter;
+import hep.aida.ITree;
+import java.util.Random;
+
+/**
+ *
+ * @author Norman A. Graf
+ *
+ * @version $Id:
+ */
+public class CalculateZSmearing
+{
+
+ public static void main(String[] args) throws Exception
+ {
+ Random ran = new Random(12345);
+ IAnalysisFactory af = IAnalysisFactory.create();
+ ITree tree = af.createTreeFactory().create();
+ IHistogramFactory hf = af.createHistogramFactory(tree);
+
+ IHistogram1D h1 = hf.createHistogram1D("Physics Z", 100, -10, 10);
+ IHistogram1D h2 = hf.createHistogram1D("Physics Z - 1 bkgnd Z", 100, -10, 10);
+ IHistogram1D h3 = hf.createHistogram1D("Physics Z - 2 bkgnd Z", 100, -10, 10);
+ IHistogram1D h4 = hf.createHistogram1D("Physics Z - 3 bkgnd Z", 100, -10, 10);
+
+ IPlotter plotter = af.createPlotterFactory().create("Z vertex smearing plot");
+ plotter.createRegions(2, 2);
+
+ plotter.region(0).plot(h1);
+ plotter.region(1).plot(h2);
+ plotter.region(2).plot(h3);
+ plotter.region(3).plot(h4);
+
+ Random r = new Random();
+
+ for (int i = 0; i < 100000; i++)
+ {
+ double p = r.nextGaussian();
+ h1.fill(p);
+ for (int j = 0; j < 1; ++j) h2.fill(p - r.nextGaussian());
+ for (int k = 0; k < 2; ++k) h3.fill(p - r.nextGaussian());
+ for (int l = 0; l < 3; ++l) h4.fill(p-r.nextGaussian());
+ }
+ plotter.show();
+ plotter.writeToFile("energyPlots.png", "png");
+ }
+}
\ No newline at end of file