CRPacket.H

Go to the documentation of this file.

 
#ifndef _CRPACKET_H_
#define _CRPACKET_H_
 
#include "../Utils/PartTypes.H"
#include "../Utils/SR.H"
#include "../Utils/Types.H"
 
namespace criptic {
 
  class CRPacket {
 
  public:
 
    // Flag values for bitmasking sources. Binary representations are:
    // noSrc         = 111111...
    // secondaryFlag = 100000...
    static constexpr
    IdxType noSrc = maxIdx; 
    static constexpr
    IdxType secondaryFlag = ~(noSrc >> 1); 
    // Data
    partTypes::pType type;   
    IdxType src;             
    IdxType uniqueID;        
    Real tInj;               
    Real wInj;               
    Real p;                  
    Real w;                  
    Real gr;                 
#ifdef TRACK_PITCH_ANGLE
    Real mu;                 
#endif
    
    // Pass-throughs for particle properties
 
    constexpr Real Z() const {
      return partTypes::charge[type];
    }
 
    constexpr Real m() const {
      return partTypes::mass[type];
    }
 
    constexpr bool isSecondary() const {
      return src & secondaryFlag;
    }
    
    // Special relativistic transformations
 
    Real gamma() const {
      return sr::gamma_from_p(type, p);
    }
 
    Real T() const {
      return sr::T_from_p(type, p);
    }
 
    void set_T(const Real T) {
      p = sr::p_from_T(type, T);
    }
    
    Real E() const {
      return sr::E_from_p(type, p);
    }
    
    Real v() const {
      return sr::v_from_p(type, p);
    }
 
    Real R() const {
      return p / std::abs(Z());
    }
 
    Real dT_dp() const {
      return sr::dT_dp(type, p);
    }
 
    Real wG(const Real B) const {
#ifdef CRIPTIC_UNITS_CGS
      return constants::e * std::abs(Z()) * B / (gamma() * m() * constants::mp_c);
#else
      return constants::e * std::abs(Z()) * B / (gamma() * m() * constants::mp);
#endif
    }
 
    Real rG(const Real B) const {
      // Perpendicular momentum in *physical* instead of code units
#ifdef TRACK_PITCH_ANGLE
      // Pitch angle tracking case
      const Real p_perp = p * constants::mp_c * std::sqrt(1 - mu*mu);
#else
      // Pitch angle averaged case
      const Real p_perp = p * constants::mp_c * M_PI_4;
#endif
#ifdef CRIPTIC_UNITS_CGS
      return p_perp * constants::c / (std::abs(Z()) * constants::e * B);
#else
      return p_perp / (std::abs(Z()) * constants::e * B);
#endif
    }
  };
  
}
 
#endif
// _CRPACKET_H_