lcsim-analysis/src/main/java/org/lcsim/analysis

diff -u -r1.3 -r1.4
--- SinglePhotonAnalysis.java	17 Jun 2010 18:06:10 -0000	1.3
+++ SinglePhotonAnalysis.java	17 Jun 2010 18:31:11 -0000	1.4
@@ -31,7 +31,7 @@

  *
  * @author Norman A. Graf
  *

- * @version $Id: SinglePhotonAnalysis.java,v 1.3 2010/06/17 18:06:10 ngraf Exp $

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

  */
 public class SinglePhotonAnalysis extends Driver
 {

@@ -51,7 +51,7 @@

 
     protected void process(EventHeader event)
     {

-         _detectorName = event.getDetectorName();

+        _detectorName = event.getDetectorName();

         _tree.mkdirs(_detectorName);
         _tree.cd(_detectorName);
         // organize the histogram tree by species and energy

@@ -265,7 +265,8 @@

                 style.statisticsBoxStyle().setVisibileStatistics("011");
                 style.legendBoxStyle().setVisible(true);
                 IFitResult fitLine = jminuit.fit(energyMeans, line);

-                System.out.println(" fit status: " + fitLine.fitStatus());

+                System.out.println("Linearity fit:");
+                printLineFitResult(fitLine);

                 double eMaxBin = eMax + 10.;
                 results.region(0).setXLimits(0., eMaxBin);
 //                results.region(0).setYLimits(0., eMaxBin);

@@ -292,7 +293,9 @@

                 style.statisticsBoxStyle().setVisibileStatistics("011");
                 style.legendBoxStyle().setVisible(false);
                 IFitResult resFitLine = jminuit.fit(resolutionFit, line);

-                System.out.println(" fit status: " + resFitLine.fitStatus());

+                System.out.println("Resolution fit:");
+                printLineFitResult(resFitLine);
+                

 //        resolution2.region(0).setXLimits(0., 1.05);
 //        resolution2.region(0).setYLimits(0., .2);
                 resolution2.region(0).plot(resolutionFit);

@@ -392,5 +395,18 @@

     {
         _defaultAidaFileName = aidaFileName;
     }

+
+    private void printLineFitResult(IFitResult fit)
+    {
+        System.out.println(" fit status: " + fit.fitStatus());
+        //String[] parNames = fit.fittedParameterNames();
+        String[] parNames = {"intercept","slope   "};
+        double[] parVals = fit.fittedParameters();
+        double[] parErrs = fit.errors();
+        for (int i = 0; i < parNames.length; ++i)
+        {
+            System.out.println(parNames[i] + " : " + parVals[i] + " +/- " + parErrs[i]);
+        }
+    }

 }