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]);
+ }
+ }
}