Changeset 19786

Timestamp:

10/24/19 15:15:03 (5 years ago)

Author:

tbretz

Message:

Implemneted the option of a (non ideal) point source.

File:

• trunk/Mars/hawc/fresnellens_psf.C (modified) (2 diffs)

trunk/Mars/hawc/fresnellens_psf.C

@@ -33 +33 @@
     double psf = 0;  // Do not change! It might have unpredicted side effectes
 
-    double angle = 6; // [cm] angle of incidence of simulates rays
+    double angle = 0; // [cm] angle of incidence of simulates rays
+
+    //TVector3 point_source(0, -F*atan(angle*TMath::DegToRad()), F);  // Point source at x=0, y=0, z=F
+    TVector3 point_source; // No point source
 
     MFresnelLens lens;
@@ -71 +74 @@
 
         TVector3 dir;
-        dir.SetMagThetaPhi(1, (180-angle)*TMath::DegToRad(), TMath::Pi()/2);
+        if (point_source.Mag()<1e-10)
+            dir.SetMagThetaPhi(1, (180-angle)*TMath::DegToRad(), TMath::Pi()/2);
+        else
+        {
+            // Note that this is not a perfect point source as the
+            // flux is homogeneous over the surface of the lens
+            TVector3 ps = pos - point_source;
+            dir.SetMagThetaPhi(1, ps.Theta(), ps.Phi());
+        }
 
         double v = 1./(TMath::C()*100/1e9); // cm/ns