Changeset 10054 for trunk/Mars/msimreflector

Timestamp:

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

Author:

tbretz

Message:

Fixed arrival time for a point source.

File:

• trunk/Mars/msimreflector/MSimRays.cc (modified) (3 diffs)

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);