Changeset 9252 for trunk/MagicSoft/Mars/msim

Timestamp:

01/24/09 01:49:53 (16 years ago)

Author:

tbretz

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars/msim

Files:

• MPhotonData.cc (modified) (2 diffs)
• MPhotonData.h (modified) (3 diffs)

trunk/MagicSoft/Mars/msim/MPhotonData.cc

@@ -40 +40 @@
 #include <iostream>
 
+#include <TMath.h>
+
 #include "MLog.h"
 #include "MLogManip.h"
@@ -95 +97 @@
 // --------------------------------------------------------------------------
 //
+// return the cosine of the Theta-angle == sqrt(1-CosU^2-CosV^2)
+//
+Double_t MPhotonData::GetCosW() const
+{
+    return TMath::Sqrt(1 - fCosU*fCosU - fCosV*fCosV);
+}
+
+// --------------------------------------------------------------------------
+//
+// Return the theta angle in radians
+//
+Double_t MPhotonData::GetTheta() const
+{
+    return TMath::ACos(GetCosW());
+}
+
+// --------------------------------------------------------------------------
+//
+// Return a TQuaternion with the first three components x, y, and z
+// and the fourth component the time.
+//
+TQuaternion MPhotonData::GetPosQ() const
+{
+    return TQuaternion(GetPos3(), fTime);
+}
+
+// --------------------------------------------------------------------------
+//
+// return a TQuaternion with the first three components the direction
+// moving in space (GetDir3()) and the fourth component is the
+// one devided by the speed of light (converted to cm/ns)
+//
+// FIXME: v in air!
+//
+TQuaternion MPhotonData::GetDirQ() const
+{
+    return TQuaternion(GetDir3(), 1./(TMath::C()*100/1e9));
+}
+
+// --------------------------------------------------------------------------
+//
 // Set the data member according to the 8 floats read from a reflector-file.
 // This function MUST reset all data-members, no matter whether these are

trunk/MagicSoft/Mars/msim/MPhotonData.h

@@ -21 +21 @@
 //class ifstream;
 #include <iosfwd>
-/*
-class TH1;
-class TH2;
-class TH3;
-*/
+
 class MCorsikaRunHeader;
 
@@ -67 +63 @@
     Float_t  GetCosU()  const { return fCosU; }
     Float_t  GetCosV()  const { return fCosV; }
-    Double_t GetCosW()  const { return TMath::Sqrt(1 - fCosU*fCosU - fCosV*fCosV); } // cos(Theta)
-    Double_t GetTheta() const { return TMath::ACos(GetCosW()); }
+    Double_t GetCosW()  const;
+    Double_t GetTheta() const;
 
     Double_t GetTime()  const { return fTime; }
@@ -81 +77 @@
     TVector3 GetDir3()  const { return TVector3(fCosU, fCosV, -GetCosW()); }
 
-    // Getter 4D                                      // cm   // ns
-    TQuaternion GetPosQ() const { return TQuaternion(GetPos3(), fTime); }
-    // FIXME: v in air!                                       // m/s  cm? / ns
-    TQuaternion GetDirQ() const { return TQuaternion(GetDir3(), 1./(TMath::C()*100/1e9)); }
+    // Getter 4D
+    TQuaternion GetPosQ() const;
+    TQuaternion GetDirQ() const;
 
     // Getter Others