Particles.hpp

Go to the documentation of this file.

#pragma once
 
#include <vector>
 
namespace kariba {
 
struct PlParams {
    double s;
    double n;
};
 
struct BknParams {
    double s1;
    double s2;
    double brk;
    double max;
    double m;
};
 
struct ThParams {
    double t;
    double n;
    double m;
};
 
struct KParams {
    double t;
    double k;
};
 
struct InjectionMixedParams {
    double s;
    double t;
    double nth;
    double npl;
    double m;
    double min;
    double max;
    double cutoff;
};
 
struct InjectionKappaParams {
    double t;
    double k;
    double n;
    double m;
};
 
struct InjectionPlParams {
    double s;
    double n;
    double m;
    double max;
};
 
struct InjectionBknParams {
    double s1;
    double s2;
    double brk;
    double max;
    double m;
    double n;
};
 
class Particles {
  protected:
    double mass_gr;     
    double mass_kev;    
 
    std::vector<double> p;        
    std::vector<double> ndens;    
    std::vector<double> gamma;    
    std::vector<double> gdens;    
    std::vector<double> gdens_diff;    
 
  public:
    Particles(size_t size);
 
    void set_mass(double m);
    void initialize_gdens();
    void initialize_pdens();
    void gdens_differentiate();
 
    const std::vector<double>& get_p() const { return p; }
 
    const std::vector<double>& get_pdens() const { return ndens; }
 
    const std::vector<double>& get_gamma() const { return gamma; }
 
    const std::vector<double>& get_gdens() const { return gdens; }
 
    const std::vector<double>& get_gdens_diff() const { return gdens_diff; }
 
    double count_particles();
    double count_particles_energy();
    double av_p();
    double av_gamma();
    double av_psq();
    double av_gammasq();
 
    void test_arrays();
};
 
}    // namespace kariba