Modules/FdSimulation/TelescopeSimulatorKG/Mirror.cc

Go to the documentation of this file.

#include "Mirror.h"
//#include "Photon.h"
#include "RTFunctions.h"
#include "RayTracer.h"

#include <iostream>

#include <utl/MathConstants.h>
#include <utl/AugerUnits.h>
#include <utl/ErrorLogger.h>
#include <utl/Photon.h>
#include <utl/Point.h>
#include <utl/Vector.h>
#include <utl/CoordinateSystemPtr.h>

#include <fdet/Telescope.h>
#include <fdet/Mirror.h>

#include <TPolyLine3D.h>
#include <TObjArray.h>

using namespace TelescopeSimulatorKG;
using namespace utl;
using namespace std;

Mirror::Mirror (utl::RandomEngine& rndm, const fdet::Telescope& tel ) {

  fRandom = &rndm;

  fTelCS = tel.GetTelescopeCoordinateSystem();
  fOrigin = Point (0, 0, 0, fTelCS);
  
  const fdet::Mirror& mirror = tel.GetMirror();

  //const TabulatedFunction& fReflectivity = mirror.GetReflectivity();
  fRCurv = mirror.GetRadiusOfCurvature();
  fSigma = mirror.GetSigmaNormal();
  
} // End of Mirror::Mirror


Mirror::~Mirror() {
}


RTResult Mirror::Trace(const utl::Photon& photonIn, utl::Photon& photonOut ) {
  
  RTFunctions::IntersectionList intersections;
  const int sphereStatus = RTFunctions::Sphere(fOrigin, fRCurv, photonIn, intersections);
  if(sphereStatus<=0) {
    ostringstream err;
    err << "Missed Mirror : " << sphereStatus << " !!!";
    ERROR(err);
    return eAbsorbed;
  }
  
  // always chose latter intersection, we want to hit inside of sphere
  photonOut = intersections.back().first;
  Vector normal(intersections.back().second);
    
  //RTFunctions::Normal(normal, fSigma, &fISeed, normAux);   // todo not yet
  
  if(!RTFunctions::Reflection (photonOut, normal, photonOut)) {
    INFO("Not reflected at mirror");
    return eAbsorbed;
  }

  return eOK;
  
} // End of Mirror::Trace


TObjArray* Mirror::Draw() {

  TObjArray* objs = new TObjArray();
  
  int color = 1;
  int size = 2;

  double width = 2.4*m;
  int nSeg = 10;
  
  double segWidth = width/nSeg;
  
  for (int i=0; i<nSeg; i++) {
    
    for (int j=0; j<nSeg; j++) {
      
      double x_center = -(width-segWidth)/2 + (width-segWidth)/(nSeg-1)*i;
      double y_center = -(width-segWidth)/2 + (width-segWidth)/(nSeg-1)*j;
      
      double x1 = x_center-segWidth/2;
      double x2 = x_center+segWidth/2;
      double y1 = y_center-segWidth/2;
      double y2 = y_center+segWidth/2;

      double r1 = sqrt (x1*x1 + y1*y1);
      double r2 = sqrt (x2*x2 + y1*y1);
      double r3 = sqrt (x1*x1 + y2*y2);
      double r4 = sqrt (x2*x2 + y2*y2);
      
      double z1 = - sqrt (fRCurv*fRCurv - r1*r1);
      double z2 = - sqrt (fRCurv*fRCurv - r2*r2);
      double z3 = - sqrt (fRCurv*fRCurv - r3*r3);
      double z4 = - sqrt (fRCurv*fRCurv - r4*r4);
      
      TPolyLine3D *l1 = new TPolyLine3D (4);
      l1->SetLineColor (color);
      l1->SetLineWidth (size);
      l1->SetPoint (0, x1,y1,z1);
      l1->SetPoint (1, x2,y1,z2);
      l1->SetPoint (2, x2,y2,z4);
      l1->SetPoint (3, x1,y2,z3);
      l1->SetPoint (4, x1,y1,z1);
      
      l1->Draw();
      objs->AddLast(l1);

    }
  }
  return objs;
}



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