28 using namespace JGEOMETRY3D;
33 const JVector3D l(sin(theta)*cos(phi),sin(theta)*sin(phi),cos(theta));
41 const double hitdist = (p0-crossing).getLength();
43 const double dist = sqrt(z*z+R*R);
45 if(hitdist<__Rdom&&dscat>=dist)
return true;
bool operator()(const double theta, const double phi, const double z, const double R, const double dscat)
Data structure for vector in three dimensions.
double getY() const
Get y position.
then usage $script[distance] fi case set_variable R
Class to check whether a photon emitted from a position along the z-axis crosses a DOM (approximated ...
then JMuonMCEvt f $INPUT_FILE o $INTERMEDIATE_FILE d
alias put_queue eval echo n
double getX() const
Get x position.
double getDot(const JVector3D &vector) const
Get dot product.
double getZ() const
Get z position.
DomCrossing(const double Rdom)