/RayTracer.h

Go to the documentation of this file.

#ifndef _TelescopeSimulatorKG2_RayTracer_h_
#define _TelescopeSimulatorKG2_RayTracer_h_


#include "RTResult.h"

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

#include <vector>
#include <string>

class TObjArray;

namespace fdet {
  class Telescope;
}

namespace utl {
  class RandomEngine;
  class Photon;
  class Point;
}

class TPolyLine3D;
class TObjArray;

namespace TelescopeSimulatorKG2 {

  class Lens;
  class Mirror;
  class Camera;
  class Filter;

  class RayTracer {

    RayTracer();
    RayTracer& operator=(const RayTracer&);
    RayTracer(const RayTracer&);

  public:

    // helper class to store photon track coordinates.
    // just for plotting (visual debugging)
    class Track {
    public:
      Track(const utl::CoordinateSystemPtr& cs) : fCS(cs), fColor(1) {}
      void AddPoint(const utl::Point& p);
      void SetColor(const int color) {fColor=color;}
      TPolyLine3D* GetLine();
    private:
      const utl::CoordinateSystemPtr& fCS;
      int fColor;
      std::vector<float> fX; // TPolyLine3D can only float ...
      std::vector<float> fY;
      std::vector<float> fZ;
    };

    RayTracer(const int Verbosity,
              const fdet::Telescope& tel,
              utl::RandomEngine& random,
              bool simulateShadow,
              bool simualteCameraSupportShadow,
              bool simulateMercedesStars,
              bool simulateFilterStructure,
              bool simulateHaloEffects,
              bool simulateGhostEffects,
              const int maxMirrorReflections,
              const double mirrorSize,
              const double mirrorSegmentSigma,
              const double mirrorRadiusSigma,
              const double mirrorAbsorptionTop,
              const double mirrorAbsorptionBot,
              utl::TabulatedFunction* mirrorDiffusionTop,
              utl::TabulatedFunction* mirrorDiffusionBot,
              const double filterIncreaseReflInside,
              const double filterIncreaseReflOutside,
              const double filterPosition,
              const double filterPositionVertical,
              const double filterPositionVertical2,
              const double filterPositionHorizontal,
              const double filterDustAbsorption,
              const double lensIncreaseRefl,
              const double lensPosition,
              const double minLensThickness,
              const double torusRadius,
              const double tubeRadius,
              const double torusZ0,
              const double pmt_n,
              bool plotPhotonTracks=false,
              bool plotPhotonTracksAtMercedes=false,
              double drawPhotonsProbabilty=0);

    ~RayTracer();

    RTResult Trace(const utl::Photon& photonIn,
                   utl::Photon& photonOut,
                   int& nreflections,
                   int& nbackscattered,
                   int& col,
                   int& row);

  private:
    int fVerbosity;
    const fdet::Telescope *fTel;
    const utl::RandomEngine *fRandom;
    bool fSimulateCameraShadow;
    bool fSimualteCameraSupportShadow;
    bool fSimulateMercedesStars;
    bool fSimulateFilterStructure;
    bool fSimulateHaloEffects;
    bool fSimulateGhostEffects;
    const int fMaxMirrorReflections;
    bool fHasCorrectorRing;

    // for 3d plotting
    bool fPlotPhotonTracks;
    double fDrawPhotonProbability;
    TObjArray* fObjectsPhotons;

    Lens *fLens;
    Camera *fCamera;
    Mirror *fMirror;
    Filter *fFilter;
    int fNTooManyReflections;
  };

} // TelescopeSimulatorKG2

#endif // _TelescopeSimulatorKG2_RayTracer_h_

// Configure (x)emacs for this file ...
// Local Variables:
// mode:c++
// compile-command: "make -C .. -k"
// End: