Changeset 2779 for trunk/MagicSoft/Mars/manalysis

Timestamp:

01/13/04 09:33:30 (21 years ago)

Author:

tbretz

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars/manalysis

Files:

• MEventRate.cc (modified) (1 diff)
• MEventRate.h (modified) (1 diff)
• MEventRateCalc.cc (modified) (1 diff)

trunk/MagicSoft/Mars/manalysis/MEventRate.cc

@@ -27 +27 @@
 // MEventRate
 //
-// Storage Container for the event rate
+// Storage Container for the event rate.
+//
+//    rate = Number of events / time in which the events were recorded
+//       r = N / t
+//       s = t / N  = 1 / r    mean time interval between consecutive events
+//
+//    for an exponential distribution of the time differences d between
+//    consecutive events:
+//
+//       s = <d>
+//       sigma(d)     = <d> = s
+//       delta(s)     = sigma(d) /sqrt(N) = s / sqrt(N)
+//       delta(s) / s = 1 / sqrt(N)
+//
+//       delta(r) / r = delta(s) / s = 1 / sqrt(N)
 //
 /////////////////////////////////////////////////////////////////////////////

trunk/MagicSoft/Mars/manalysis/MEventRate.h

@@ -17 +17 @@
     void SetRate(Double_t r, UInt_t n) { fRate = r; fNumEvents = n; }
     Double_t GetRate() const { return fRate; }
+    Double_t GetError() const { return TMath::Sqrt(1./fNumEvents); }
     UInt_t GetNumEvents() const { return fNumEvents; }
 

trunk/MagicSoft/Mars/manalysis/MEventRateCalc.cc

@@ -30 +30 @@
 //
 //  The algorithm is explained in Process()
+//
+//
+//    rate = Number of events / time in which the events were recorded
+//       r = N / t
+//       s = t / N  = 1 / r    mean time interval between consecutive events
+//
+//    for an exponential distribution of the time differences d between
+//    consecutive events:
+//
+//       s = <d>
+//       sigma(d)     = <d> = s
+//       delta(s)     = sigma(d) /sqrt(N) = s / sqrt(N)
+//       delta(s) / s = 1 / sqrt(N)
+//
+//       delta(r) / r = delta(s) / s = 1 / sqrt(N)
+//
 //
 //  Input Containers: