Changeset 10054

Timestamp:

11/22/10 16:10:06 (14 years ago)

Author:

tbretz

Message:

Fixed arrival time for a point source.

Location:

trunk/Mars

Files:

• Changelog (modified) (1 diff)
• msimreflector/MSimRays.cc (modified) (3 diffs)

trunk/Mars/Changelog

@@ -23 +23 @@
    * Makefile.rules:
      - use -C option of make instead of cd
+
+   * msimreflector/MSimRays.cc:
+     - added some comments
+     - fixed timing for a point source
 
 

trunk/Mars/msimreflector/MSimRays.cc

@@ -164 +164 @@
         MPhotonData &dat = *static_cast<MPhotonData*>(arr.UncheckedAt(idx));
 
+        // Get radom incident point on the mirror plane.
         const Double_t x = gRandom->Uniform(-maxr, maxr);
         const Double_t y = gRandom->Uniform(-maxr, maxr);
@@ -170 +171 @@
             continue;
 
+        // The is the incident direction of the photon
         // h==0 means infinitiy
         const TVector3 u = fHeight>0 ? TVector3(x, y, -h).Unit() :  TVector3(0, 0, -1);
@@ -176 +178 @@
         // CORSIKA-orig: x(north), y(west),  z(up), t(time)
         // NOW:          x(east),  y(north), z(up), t(time)
-        MQuaternion p(TVector3(x, y, 0));
+        MQuaternion p(TVector3(x, y, 0), fHeight>0 ? TMath::Sqrt(x*x + y*y + h*h): 0);
         MQuaternion w(u, conv);