DomCrossing.hh

Go to the documentation of this file.

#ifndef __DOMCROSSING__
#define __DOMCROSSING__
 
#include <cmath>
 
#include "JGeometry3D/JVector3D.hh"
 
/**
 * \author mdejong
 */
 
    /**
     * Class to check whether a photon emitted from a position along the
     * z-axis crosses a DOM (approximated by a disc)
     * first, inefficient version
     */
 
    class DomCrossing {
 
    private:
      double __Rdom;
    public:
 
      DomCrossing(const double Rdom):__Rdom(Rdom)
      {}
      ~DomCrossing()
      {}
 
      bool operator()(const double theta,
                      const double phi,
                      const double z,
                      const double R,
                      const double dscat)
      {
        using namespace JGEOMETRY3D;
 
        const JVector3D p0(R,0,0);
        const JVector3D n(sin(theta),0,cos(theta));
        const JVector3D l0(0,0,z);
        const JVector3D l(sin(theta)*cos(phi),sin(theta)*sin(phi),cos(theta));
 
        const double d = JVector3D(p0-l0).getDot(n)/l.getDot(n);
 
        if(d<0) return false;
 
        const JVector3D crossing = l0+JVector3D(l.getX()*d,l.getY()*d,l.getZ()*d);
 
        const double hitdist = (p0-crossing).getLength();
 
        const double dist = sqrt(z*z+R*R);
 
        if(hitdist<__Rdom&&dscat>=dist) return true;
        else return false;
 
      }
 
    };
 
    //  DomCrossing shadow;
#endif