Changeset 9342 for trunk/MagicSoft/Mars/msim/MPhotonEvent.cc

Timestamp:

02/15/09 15:52:47 (16 years ago)

Author:

tbretz

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/msim/MPhotonEvent.cc (modified) (3 diffs)

trunk/MagicSoft/Mars/msim/MPhotonEvent.cc

@@ -113 +113 @@
 #include <TMarker.h>
 
+#include <MMath.h>
+
+#include "MArrayF.h"
+
 #include "MLog.h"
 #include "MLogManip.h"
@@ -213 +217 @@
 }
 
+// --------------------------------------------------------------------------
+//
+// Return the number of "external" photons, i.e. which are not NightSky
+//
 Int_t MPhotonEvent::GetNumExternal() const
 {
@@ -222 +230 @@
 
     return n;
+}
+
+// --------------------------------------------------------------------------
+//
+// Return time of first photon, 0 if none in array.
+// Note: If you want this to be the earliest make sure that the array
+// is properly sorted.
+//
+Float_t MPhotonEvent::GetTimeFirst() const
+{
+    const MPhotonData *dat=GetFirst();
+    return dat ? dat->GetTime() : 0;
+}
+
+// --------------------------------------------------------------------------
+//
+// Return time of first photon, 0 if none in array.
+// Note: If you want this to be the latest make sure that the array
+// is properly sorted.
+//
+Float_t MPhotonEvent::GetTimeLast() const
+{
+    const MPhotonData *dat=GetLast();
+    return dat ? dat->GetTime() : 0;
+}
+
+// --------------------------------------------------------------------------
+//
+// Return the median devian from the median of all arrival times.
+// The median deviation is calculated using MMath::MedianDev.
+// It is the half width in which one sigma (~68%) of all times are
+// contained around the median.
+//
+Double_t MPhotonEvent::GetTimeMedianDev() const
+{
+    const UInt_t n = GetNumPhotons();
+
+    MArrayF arr(n);
+    for (UInt_t i=0; i<n; i++)
+        arr[i] = operator[](i).GetTime();
+
+    return MMath::MedianDev(n, arr.GetArray()/*, Double_t &med*/);
 }