Changeset 19649 for trunk/Mars/hawc


Ignore:
Timestamp:
09/15/19 22:24:12 (5 years ago)
Author:
tbretz
Message:
Defining correct focal distance, disentengle wavelength for lens definition from photon wavelength
File:
1 edited

Legend:

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

    r19645 r19649  
    3232    const double D = 54.92;  // [cm] Diameter of the refractive surface
    3333    const double F = 50.21;  // [cm] Nominal focal length of the lens
    34     const double Z = F+H/2;  // [cm] camera position w.r.t. lens exit (see also MGeomCamFAMOUS!)
    3534
    36     const double R = D/2;    // [cm] radius of lens
     35    const double H = 0.25;   // [cm] Width of a single groove
    3736    const double w = 0.01;   // [cm] Width of a single groove
    38     const double H = 0.25;   // [cm] Width of a single groove
    3937    //const double w = 1;        // [cm] Width of a single groove
    4038    //const double H = 2.5;    // [cm] Thickness of lens
    4139
    42     const double lambda = 546; // [nm] Wavelength for lens surface definition
     40    const double Z = F+H;    // [cm] camera position w.r.t. lens exit (see also MGeomCamFAMOUS!)
     41    const double R = D/2;    // [cm] radius of lens
    4342
    44     double angle = 0;  // Check ix=0 at angle=-55 and N=10
     43    const double lambda0 = 546; // [nm] Wavelength for lens surface definition
     44    const double lambda  = 546; // [nm] Wavelength for the simulated photons
    4545
    46     double Z0 = 3;
     46    double angle = 6; // Angle of incidence of the simulated rays
    4747
    48     const double n0 = MFresnelLens::RefractiveIndex(lambda);
     48    double Z0 = 3;    // Starting Z-position of rays
    4949
    5050    MFresnelLens lens;
    51     lens.DefineLens(F, D, w, H, lambda);
     51    lens.DefineLens(F, D, w, H, lambda0);
    5252    //lens.SetPSF(psf);
    53     //lens.DefineLens(F, D, w, H, lambda);
    5453    //lens.EnableSlopeAbsorption();
    5554    //lens.EnableDraftAbsorption();
Note: See TracChangeset for help on using the changeset viewer.