Changeset 9342 for trunk/MagicSoft/Mars/msim

Timestamp:

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

Author:

tbretz

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars/msim

Files:

• MPhotonEvent.cc (modified) (3 diffs)
• MPhotonEvent.h (modified) (4 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*/);
 }
 

trunk/MagicSoft/Mars/msim/MPhotonEvent.h

@@ -80 +80 @@
     Int_t GetNumExternal() const;
 
+    Float_t  GetTimeFirst() const;
+    Float_t  GetTimeLast() const;
+    Double_t GetTimeMedianDev() const;
+
     TClonesArray &GetArray() { return fData; }
     const TClonesArray &GetArray() const { return fData; }
@@ -131 +135 @@
 {
 private:
-    Float_t fTimeFirst;
-    Float_t fTimeLast;
+    Float_t fTimeFirst;  //! Start of (simulated) sampling window
+    Float_t fTimeLast;   //! Start of (simulated) sampling window
+
+    Float_t fLength;     // Time between first and last photon
+    Float_t fTimeMedDev; // Median deviation
 
 //    Float_t fOffset;
 //    Float_t fWindow;
 
-    Int_t fMaxIndex;
+    Int_t fMaxIndex;     //!
 
 public:
@@ -147 +154 @@
 
     void SetTime(Float_t first, Float_t last) { fTimeFirst=first; fTimeLast=last; }
+    void SetLength(Float_t len) { fLength=len; }
     void SetMaxIndex(UInt_t idx) { fMaxIndex=idx; }
+    void SetTimeMedDev(Float_t dev) { fTimeMedDev=dev; }
 
 //    Float_t GetRawTimeFirst() const { return fTimeFirst; }
@@ -154 +163 @@
     Float_t GetTimeFirst() const { return fTimeFirst; }
     Float_t GetTimeLast() const { return fTimeLast; }
+
+    Float_t GetLength() const { return fLength; }
+    Float_t GetTimeMedDev() const { return fTimeMedDev; }
 
     Int_t GetMaxIndex() const { return fMaxIndex; }