1 #ifndef __JSIRENE__JSIRENETOOLKIT__
2 #define __JSIRENE__JSIRENETOOLKIT__
25 namespace JPP {
using namespace JSIRENE; }
62 for (
int i = N; i != 0; --i) {
63 if (gRandom->Rndm() * NPE < npe) {
72 return gRandom->Poisson(npe);
89 static const int NPE = 20;
93 const int n = (int) (NPE * log10((
double) npe / (
double) NPE));
114 using namespace JAANET;
115 using namespace JPHYSICS;
335 std::vector<
JVertex>(1, vertex)
345 double getE(
const double z)
const
347 if (!this->empty()) {
351 if (pos != this->end() && pos != this->begin()) {
355 return pos->getE() - (z - pos->getZ()) *
gWater.
getA();
374 return sqrt((E - m) * (E + m));
JHitType_t
Enumeration of hit types based on km3 codes.
double getEs() const
Get shower energy.
JTrack(const JVertex &vertex)
Constructor.
Scattered light from primary shower.
double getT() const
Get time.
JVertex & step(const double ds)
Step.
JComparator< JResult_t T::*, JComparison::lt > make_comparator(JResult_t T::*member)
Helper method to create comparator between values of data member.
double getRange(double a) const
Get range of muon using precalculated a parameter This method applies only ionisation energy loss...
void applyEloss(const double Es)
Apply energy loss energy.
Direct light from delta-rays.
Scattered light from muon.
double getIndexOfRefraction()
Get average index of refraction of water corresponding to group velocity.
scattered light from EM shower
JHitType_t getHitType(const JPDFType_t pdf, const bool shower=false)
Get hit type corresponding to given PDF type.
static const double MASS_MUON
muon mass [GeV]
double getE() const
Get muon energy.
Scattered light from Bremsstrahlung.
static const JGeaneWater gWater
Function object for energy loss of muon in sea water.
Direct light from Bremsstrahlung.
Direct light from primary shower.
JVertex()
Default constructor.
JVertex(const double z, const double t, const double E)
Constructor.
Numbering scheme for PDF types.
JVertex & step(const JGeane &geane, const double ds)
Step.
scattered light from muon
double getZ() const
Get position.
double getE(const double z) const
Get muon energy at given position.
scattered light from delta-rays
direct light from EM shower
Interface for muon energy loss.
double getKineticEnergy(const double E, const double m)
Get kinetic energy of particle with given mass.
double getRange(const JGeane &geane) const
Get range of muon.
Vertex of energy loss of muon.
direct light from delta-rays
int getNumberOfPhotoElectrons(const double NPE, const int N, const double npe)
Get number of photo-electrons of a hit given the expectation values of the number of photo-electrons ...
const double getSpeedOfLight()
Get speed of light.
Scattered light from delta-rays.
alias put_queue eval echo n
double getRange() const
Get range of muon.
JVertex & step(const double ds, const double a)
then usage $script[input file[working directory[option]]] nWhere option can be N
then usage $script[input file[working directory[option]]] nWhere option can be E
virtual double getA() const override
Get energy loss constant.