Print

Print


Commit in lcsim/src/org/lcsim/contrib/LeiXia on MAIN
rms90_order.java+7-11.1 -> 1.2


lcsim/src/org/lcsim/contrib/LeiXia
rms90_order.java 1.1 -> 1.2
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);
 	 
    }
 }
CVSspam 0.2.8