Changeset 19648
- Timestamp:
- 09/15/19 22:23:08 (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/Mars/hawc/fresnellens_psf.C
r19645 r19648 13 13 14 14 ******************************************************************/ 15 void fresnellens_psf 2()15 void fresnellens_psf() 16 16 { 17 17 // ------------------- setup lens ----------------------- … … 29 29 const double Rcam = 1.5*4.5; // [cm] radius of camera 30 30 31 // const double lambda = 546; // [nm] Wavelength for lens surface definition31 const double lambda = 546; // [nm] Wavelength for simulated rays 32 32 33 33 double psf = 0; // Do not change! It might have unpredicted side effectes 34 34 35 double angle = 8;35 double angle = 6; // [cm] angle of incidence of simulates rays 36 36 37 37 MFresnelLens lens; … … 50 50 51 51 TH2F h_out("OUT", "Photon Distribution (out)", NN*4/3, -R*4/3, R*4/3, NN, -R, R); 52 TH2F h_psf("PSF", "Point Spread Function", 1000*4/3, -R*4/3, R*4/3, 1000, -R, R);52 TH2F h_psf("PSF", "Point Spread Function", 2000*4/3, -R*4/3, R*4/3, 2000, -R, R); 53 53 TH2F h_cam("CAM", "Point Spread Function", 4*37*4/3, -14*4/3, 14*4/3, 4*37, -14, 14); 54 54 … … 82 82 // --------------------------------------------------------------- 83 83 84 int ret = lens.ExecuteOptics(pp, uu, 546);84 int ret = lens.ExecuteOptics(pp, uu, lambda); 85 85 86 86 // Skip photons which will not hit the focal plane within r<R
Note:
See TracChangeset
for help on using the changeset viewer.