diff -u -r1.4 -r1.5
--- SinglePhotonAnalysis.java	17 Jun 2010 18:31:11 -0000	1.4
+++ SinglePhotonAnalysis.java	17 Jun 2010 19:55:48 -0000	1.5
@@ -31,7 +31,7 @@

  *
  * @author Norman A. Graf
  *

- * @version $Id: SinglePhotonAnalysis.java,v 1.4 2010/06/17 18:31:11 ngraf Exp $

+ * @version $Id: SinglePhotonAnalysis.java,v 1.5 2010/06/17 19:55:48 ngraf Exp $

  */
 public class SinglePhotonAnalysis extends Driver
 {

@@ -191,20 +191,24 @@

                         System.out.println("Energy " + energies[i]);
                         // +/- 5 sigma
                         int nSigma = 5;

-                        double lowE = energies[i] - (.18 * sqrt(energies[i]) * nSigma);
-                        double hiE = energies[i] + (.18 * sqrt(energies[i]) * nSigma);

+                        double expectedSigma = _expectedResolution * sqrt(energies[i]);
+                        double lowE = energies[i] - ( expectedSigma* nSigma);
+                        double hiE  = energies[i] + ( expectedSigma * nSigma);

 
                         ICloud1D e = (ICloud1D) _tree.find(objects[i] + "Cluster Energy");
                         if (!e.isConverted())
                         {

-                            int nbins = 50;

+                            int nbins = 100;

 //                            System.out.println(energies[i] + " - " + lowE + " + " + hiE);
                             e.convert(nbins, lowE, hiE);
                         }
                         IHistogram1D eHist = e.histogram();
                         gauss.setParameter("amplitude", eHist.maxBinHeight());

-                        gauss.setParameter("mean", eHist.mean());
-                        gauss.setParameter("sigma", eHist.rms());

+//                        gauss.setParameter("mean", eHist.mean());
+//                        gauss.setParameter("sigma", eHist.rms());
+                        // robustify...
+                        gauss.setParameter("mean", energies[i]);
+                        gauss.setParameter("sigma", expectedSigma);

                         style = plotter.region(i).style();
                         style.legendBoxStyle().setVisible(false);
                         style.statisticsBoxStyle().setVisible(false);