Changeset 19648


Ignore:
Timestamp:
09/15/19 22:23:08 (5 years ago)
Author:
tbretz
Message:
Minor fixes and improvements.
File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/Mars/hawc/fresnellens_psf.C

    r19645 r19648  
    1313
    1414******************************************************************/
    15 void fresnellens_psf2()
     15void fresnellens_psf()
    1616{
    1717    // ------------------- setup lens -----------------------
     
    2929    const double Rcam = 1.5*4.5; // [cm] radius of camera
    3030
    31     // const double lambda = 546; // [nm] Wavelength for lens surface definition
     31    const double lambda = 546; // [nm] Wavelength for simulated rays
    3232
    3333    double psf = 0;  // Do not change! It might have unpredicted side effectes
    3434
    35     double angle = 8;
     35    double angle = 6; // [cm] angle of incidence of simulates rays
    3636
    3737    MFresnelLens lens;
     
    5050
    5151    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);
    5353    TH2F h_cam("CAM", "Point Spread Function",     4*37*4/3, -14*4/3, 14*4/3, 4*37, -14, 14);
    5454
     
    8282        // ---------------------------------------------------------------
    8383
    84         int ret = lens.ExecuteOptics(pp, uu, 546);
     84        int ret = lens.ExecuteOptics(pp, uu, lambda);
    8585
    8686        // Skip photons which will not hit the focal plane within r<R
Note: See TracChangeset for help on using the changeset viewer.