SecondaryEmissionPhysics.h

Go to the documentation of this file.

#ifndef OPAL_SECONDARY_EMISSION_PHYSICS_HH
#define OPAL_SECONDARY_EMISSION_PHYSICS_HH

#include <hdf5.h>
#include <cmath>
#include <limits>
#include <sys/stat.h>
#include "Physics/Physics.h"
#include "Algorithms/PBunchDefs.h"
#include "Algorithms/PartBunchBase.h"

extern Inform *gmsg;

namespace myeps {
    const double EPS = std::numeric_limits<double>::epsilon();
    const double FPMIN = std::numeric_limits<double>::min() / EPS;
}

class SecondaryEmissionPhysics

{
public:

    SecondaryEmissionPhysics();

    ~SecondaryEmissionPhysics();

    void nSec(const double &incEnergy,
              const double &cosTheta,
              const int &matNumber,
              int &seNum,
              int &seType,
              const double &incQ,
              const Vector_t &TriNorm,
              const Vector_t &inteCoords,
              const Vector_t &localX,
              PartBunchBase<double, 3> *itsBunch,
              double &seyNum, const
              double &ppVw,
              const double &vVThermal,
              const bool nEmissionMode);

    void nSec(const double &incEnergy,
              const double &cosTheta,
              int &seNum,
              int &seType,
              const double &incQ,
              const Vector_t &TriNorm,
              const Vector_t &inteCoords,
              const Vector_t &localX,
              PartBunchBase<double, 3> *itsBunch,
              double &seyNum,
              const double &ppVw,
              const double &vSeyZero,
              const double &vEzero,
              const double &vSeyMax,
              const double &vEmax,
              const double &vKenergy,
              const double &vKtheta,
              const double &vVThermal,
              const bool nEmissionMode);

    double deltae_m;
    double deltar_m;
    double deltats_m;

    private:

    IpplTimings::TimerRef TPnSec_m;

    double seAlpha_m;


    /*==============================================================================================
      parameters for backscattered electrons
      ==============================================================================================*/



    double sePScat_m;
    double sePScatPeak_m;
    double seEScatPeak_m;
    double seW_m;
    double seP_m;
    double seDelta_m;
    double seEOne_m;
    double seETwo_m;



    /*==============================================================================================
      parameters for rediffused electrons
      ==============================================================================================*/



    double sePRed_m;
    double seERed_m;
    double seR_m;
    double seQ_m;
    double seROne_m;
    double seRTwo_m;


    /*==============================================================================================
      parameters for true secondary electrons
      ==============================================================================================*/

    double seYPeakTS_m;
    double seEPeakTS_m;
    double seSTS_m;
    double seTOneTS_m;
    double seTTwoTS_m;
    double seTThreeTS_m;
    double seTFourTS_m;
    double seEPeakTot_m;
    double seYPeakTot_m;

    /*==============================================================================================
      parameters in he TABLE II of Phys. Rev. ST Accel. Beams 5, 124404 (2002)
      ==============================================================================================*/
    double sePn_m[10];
    double seEpsn_m[10];


    void calcEmiNum(const double incEnergy,
                    const double cosTheta,
                    int &seNum,
                    const double &vSeyZero,
                    const double &vEzero,
                    const double &vSeyMax,
                    const double &vEmax,
                    const double &vKenergy,
                    const double &vKtheta,
                    double &seyNum);

    void calcEmiNum(const double incEnergy, const double cosTheta, const double *prob, int &seNum);

    double calcDeltats(const double incEnergy, const double cosTheta);

    double calcDeltar(const double incEnergy, const double cosTheta);

    double calcDeltae(const double incEnergy, const double cosTheta);

    double calcProb(const double incEnergy, const double cosTheta, double *prob);

    void setSeMaterial(int material_num);

    /*==========================================
     *http://www.taygeta.com/random/gaussian.html
     *return a gaussian distributed random number
     *==========================================*/

    double gaussRand();

    double gammp(const double a, const double x);

    double gser(const double a, const double x);
    double gcf(const double a, const double x);
    double gammpapprox(double a, double x, int psig);
    double gammln(const double xx);

    double invgammp(double p, double a);
    double min_Gamma(double x, double y);
    double max_Gamma(double x, double y);


    double betai(const double x, const double a, const double b);
    double betacf(const double a, const double b, const double x);
    double betaiapprox(double a, double b, double x);
    double invbetai(double p, double a, double b);
    void coordConverter (const Vector_t &/*TriNormal*/, Vector_t &/*x*/) {  }



};

inline
double SecondaryEmissionPhysics::min_Gamma(double x, double y) {
    if(x > y)
        return y;
    else
        return x;
}

inline
double SecondaryEmissionPhysics::max_Gamma(double x, double y) {
    if(x < y)
        return y;
    else
        return x;
}


/*  ==========================================================================*/







#endif //OPAL_SECONDARY_EMISSION_PHYSICS_HH