Changeset 3005


Ignore:
Timestamp:
02/02/04 17:03:36 (21 years ago)
Author:
gaug
Message:
*** empty log message ***
Location:
trunk/MagicSoft/Mars
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • trunk/MagicSoft/Mars/Changelog

    r3004 r3005  
    44
    55                                                 -*-*- END OF LINE -*-*-
     6
     7 2004/02/02: Markus Gaug
     8
     9  * manalysis/MPedCalcPedRun.cc
     10    - corrected formulae in Class descriptino
     11
    612
    713 2004/02/02: Abelardo Moralejo
  • trunk/MagicSoft/Mars/manalysis/MPedCalcPedRun.cc

    r2997 r3005  
    3939// Actually, MPedCalcPedRun applies the following formula (1):
    4040//
    41 // PedRMS = Sqrt(  (sum(x_i^2) - sum(x_i)/n) / n-1 / 14 )
     41// PedRMS = Sqrt(  (sum(x_i^2) - sum(x_i)^2/n) / n-1 / 14 )
    4242//
    4343// where x_i is the sum of 14 FADC slices and sum means the sum over all
     
    4646// For a high number of events, this formula is equivalent to formula (2):
    4747//
    48 // PedRMS = Sqrt(  (<x_i*x_i> - <x_i>*<x_i>*n) / 14  )
     48// PedRMS = Sqrt(  (<x_i*x_i> - <x_i>*<x_i>) / 14  )
    4949//
    5050// where <> is the mean over all events and x_i again the sum over the 14
     
    5555// equivalent to (old formula) (3):
    5656//
    57 // PedRMS = Sqrt(  (<p_i*p_i> - <p_i>*<p_i>*m) / m  ) * Sqrt(14)
    58 //
    59 // which is the RMS of a single slice (p_i) with m being the total number of
    60 // measurements, i.e. m = n*14, later re-scaled to the number of used slices
    61 // (the factor sqrt(14)).
     57// PedRMS = Sqrt(  (<p_i*p_i> - <p_i>*<p_i>) )
     58//
     59// which is the RMS per slice of a single slice (p_i) and
     60// <> the mean over the total number of measurements, i.e. n*14.
    6261//
    6362// If we assume that at least our pairs fluctuate independently and Gaussian,
    64 // then we can use the actual formula (1) in order to get what you call
     63// then we can use the actual formula (1) in order to get
    6564// fluctuations of pairs by the transformation:
    6665//
    67 // PedRMS/pair = PedRMS (form. (3)) / Sqrt(7)
     66// PedRMS/pair = PedRMS (form. (3)) * Sqrt(2)
    6867//
    6968// (However, we know that our slice-to-slice fluctuations are not Gaussian
    7069// (and moreover asymmetric) and that they are also correlated.)
    7170//
    72 // We could still measure also the pair-to-pair fluctuations and add another
    73 // value to be investigated. What do you think?
    74 //
    7571//
    7672//  Input Containers:
Note: See TracChangeset for help on using the changeset viewer.