/Lens.h

Go to the documentation of this file.

#ifndef _TelescopeSimulatorKG2_Lens_h_
#define _TelescopeSimulatorKG2_Lens_h_

#include "RTResult.h"
#include "RayTracer.h"
#include "RTFunctions.h"

#include <utl/CoordinateSystemPtr.h>
#include <utl/Point.h>

#include <string>


namespace utl {
  class Photon;
  class Vector;
  class TabulatedFunction;
  class RandomEngine;
}

namespace fdet {
  class Telescope;
}

class TObjArray;

namespace TelescopeSimulatorKG2 {

  class Lens {

  public:
    Lens(utl::RandomEngine& rndm,
         const fdet::Telescope& tel,
         const bool simulateHaloEffects,
         const double lensIncreaseRefl,
         const double lensPosition,
         const double minLensThickness,
         const double torusRadius,
         const double tubeRadius,
         const double torusZ0);
    ~Lens();

    RTNext Trace(const utl::Photon& photonIn, utl::Photon& photonOut);
    RTNext TraceWithTorus(const utl::Photon& photonIn, utl::Photon& photonOut);

    TObjArray* Draw();

    // this is just needed for photon plotting
    void SetPhotonTrack(RayTracer::Track* const track) { fPhotonTrack = track; }

  private:
    double ProfilePrototype(const double r) const;
    utl::Vector CurvaturePrototype(const double x, const double y) const;
    void ChoosePhysicalIntersection(unsigned int& index,
                                    const RTFunctions::IntersectionList& intersections) const;

  private:
    utl::RandomEngine& fRandom;
    const bool fSimulateHaloEffects;

    // to artificially increase reflection
    const double fIncreaseReflection;

    // geometry of the lens (and torus)
    const double fPosZ;
    const double fMinLensThickness;
    const double fTorusRadius;
    const double fTubeRadius;
    const double fTorusZ0;
    double fLXTorusRadius = 0;
    double fLXTubeRadius = 0;
    double fLXz0 = 0;

    utl::CoordinateSystemPtr fTelCS; // reference system
    const utl::Point fOrigin;
    const utl::Point fTorusOrigin;
    const utl::Point fPosPlane;
    //utl::Point fFilterCenter;

    const double fR1;  // geometry R1=.85m
    const double fR2;  // geometry R2=1.1m
    const double fIndexAir = 1.000277;  // reflection index of air, at STP
    const utl::TabulatedFunction& fRefractiveIndex;
    const utl::TabulatedFunction& fTransmittance;
    const double fSigRho;  // uncertaint in the definition of normal of the surface of the Lens
    unsigned int fNbug;

    // for visual plotting of ray tracing details
    RayTracer::Track* fPhotonTrack = nullptr;
  };

}


#endif