JSIRENE::JTrack Struct Reference

Muon trajectory. More...

#include <JSireneToolkit.hh>

Inheritance diagram for JSIRENE::JTrack:

Public Member Functions
	JTrack (const JPoint &point)
	Constructor. More...
double	getE (const double z) const
	Get muon energy at given position along trajectory. More...
JVector2D	getPosition (const double z) const
	Get muon position at given position along trajectory. More...

Detailed Description

Muon trajectory.

Definition at line 247 of file JSireneToolkit.hh.

Constructor & Destructor Documentation

JTrack()

JSIRENE::JTrack::JTrack ( const JPoint &  point )

inline

Constructor.

Parameters

point

muon starting point

Definition at line 255 of file JSireneToolkit.hh.

                                :
      std::vector<JPoint>(1, point)
    {}

Member Function Documentation

getE()

double JSIRENE::JTrack::getE ( const double  z ) const

inline

Get muon energy at given position along trajectory.

Parameters

z	position [m]

Returns

energy [GeV]

Definition at line 266 of file JSireneToolkit.hh.

    {
       if (!this->empty()) {
        
        const_iterator p = lower_bound(this->begin(), this->end(), z, JLANG::make_comparator(&JPoint::getZ));
 
        if (p != this->end() && p != this->begin()) {
          
          --p;
        
          return p->getE() -  (z - p->getZ()) * gWater.getA();
        }
      }
 
      return MASS_MUON;
    }

getPosition()

JVector2D JSIRENE::JTrack::getPosition ( const double  z ) const

inline

Get muon position at given position along trajectory.

Parameters

z	position [m]

Returns

position

Definition at line 290 of file JSireneToolkit.hh.

    {
      using namespace JPP;
 
      const double precision = 1.0e-2;
 
      if (!this->empty()) {
        
        const_iterator p = lower_bound(this->begin(), this->end(), z, JLANG::make_comparator(&JPoint::getZ));
 
        if (p == this->end()) { 
          --p;
        }
 
        if (p == this->begin()) {
 
          return JVector2D(p->getX(),
                           p->getY());
 
        } else {
 
          JVector3D pos;
 
          pos  = p->getPosition();
 
          --p;
 
          pos -= p->getPosition();
 
          const double u = (pos.getZ() > precision ? (z - p->getZ()) / pos.getZ() : 0.0);
 
          return JVector2D(p->getX() + u * pos.getX(),
                           p->getY() + u * pos.getY());
        }
      }
 
      return JVector2D(0.0, 0.0);
    }

---

The documentation for this struct was generated from the following file:

• JSireneToolkit.hh