Changeset 3005 for trunk/MagicSoft/Mars/manalysis/MPedCalcPedRun.cc

Timestamp:

02/02/04 17:03:36 (21 years ago)

Author:

gaug

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/manalysis/MPedCalcPedRun.cc (modified) (3 diffs)

trunk/MagicSoft/Mars/manalysis/MPedCalcPedRun.cc

@@ -39 +39 @@
 // Actually, MPedCalcPedRun applies the following formula (1):
 // 
-// PedRMS = Sqrt(  (sum(x_i^2) - sum(x_i)/n) / n-1 / 14 )
+// PedRMS = Sqrt(  (sum(x_i^2) - sum(x_i)^2/n) / n-1 / 14 )
 // 
 // where x_i is the sum of 14 FADC slices and sum means the sum over all
@@ -46 +46 @@
 // For a high number of events, this formula is equivalent to formula (2):
 // 
-// PedRMS = Sqrt(  (<x_i*x_i> - <x_i>*<x_i>*n) / 14  )
+// PedRMS = Sqrt(  (<x_i*x_i> - <x_i>*<x_i>) / 14  )
 // 
 // where <> is the mean over all events and x_i again the sum over the 14
@@ -55 +55 @@
 // equivalent to (old formula) (3):
 // 
-// PedRMS = Sqrt(  (<p_i*p_i> - <p_i>*<p_i>*m) / m  ) * Sqrt(14)
-// 
-// which is the RMS of a single slice (p_i) with m being the total number of
-// measurements, i.e. m = n*14, later re-scaled to the number of used slices
-// (the factor sqrt(14)).
+// PedRMS = Sqrt(  (<p_i*p_i> - <p_i>*<p_i>) ) 
+// 
+// which is the RMS per slice of a single slice (p_i) and 
+// <> the mean over the total number of measurements, i.e. n*14.
 // 
 // If we assume that at least our pairs fluctuate independently and Gaussian,
-// then we can use the actual formula (1) in order to get what you call
+// then we can use the actual formula (1) in order to get 
 // fluctuations of pairs by the transformation:
 // 
-// PedRMS/pair = PedRMS (form. (3)) / Sqrt(7)
+// PedRMS/pair = PedRMS (form. (3)) * Sqrt(2)
 // 
 // (However, we know that our slice-to-slice fluctuations are not Gaussian
 // (and moreover asymmetric) and that they are also correlated.)
 // 
-// We could still measure also the pair-to-pair fluctuations and add another
-// value to be investigated. What do you think?
-//
 // 
 //  Input Containers: