lcsim/src/org/lcsim/contrib/LeiXia
diff -u -r1.1 -r1.2
--- rms90_order.java 29 Sep 2006 16:51:13 -0000 1.1
+++ rms90_order.java 29 Sep 2006 20:54:13 -0000 1.2
@@ -52,6 +52,7 @@
double CentralFraction = 0.90;
double rms_min = 1000.; // minimum rms for whatever CentralFraction range
+ double mean_min = 0.; // mean of the 90% entries, when rms_min is achieved
double rms = 0.;
int NTail = (int)(Entry*(1-CentralFraction));
@@ -72,10 +73,15 @@
rms = Math.sqrt(xm2/xn - (xm/xn)*(xm/xn));
test.fill(i*1.0, rms);
- if (rms < rms_min) rms_min = rms;
+ if (rms < rms_min)
+ {
+ rms_min = rms;
+ mean_min = xm/xn;
+ }
}
System.out.println("rms_90 : "+rms_min);
+ System.out.println("mean of the 90% event : "+mean_min);
}
}