---

Brempr.h

Go to the documentation of this file.

// -*- C++ -*-  $Id: Brempr.h,v 1.1 2000/09/01 20:20:11 burnett Exp $
//
// This file is part of Gismo 2
//
// Contents ----------------------------------------------------------------
//
//      class Brempr
//
// Description
//
//  This class  corresponds to egs common block /brempr/
//
// End --------------------------------------------------------------------

#ifndef __BREMPR_H
#define __BREMPR_H

#ifdef _MSC_VER
# pragma warning(disable:4305) //truncation from 'const double' to 'float'
#endif
#include <cmath>

#include <iostream>

class Brempr
{

    friend class PEGSData;

    float alphi[2];
    float bpar[2];
    float delpos[2];
    float delcm;
    float dl1[6],dl2[6],dl3[6],dl4[6],dl5[6],dl6[6];
    void readPEGS(std::istream&);


    inline float screening(float del, int lvx, int lvl)const
    { if (del > delpos[lvx]) return 0;
       //   a(delta) and b(delta) must always be positive
      float delta = del*delcm;
      return (delta<1.0)
        ? dl1[lvl] + delta*( dl2[lvl] + delta*dl3[lvl])
        : dl4[lvl] + dl5[lvl]*log(delta+dl6[lvl]);
    }

};

struct Kedge
{   friend class PEGSData;


    int   iedgfl; // flag that following data set up
    float ekalph; // energy of the k-alpha x-ray.
    float ekbeta; // energy of the k-beta x-ray.
    float bkr1;   // probability for a k-shell emission times the
                  //   probability for a fluorescent photon emission
                  //   due to an electron transition from the x -> k shell.
    float bkr2;   // bkr1 times the probability that the fluorescent photon
                  // is due to an electron transition from the l -> k shell.
    //------------------------------------------------------------------
    // ebinda is the energy of the k-edge (in class Photin)

    Kedge(){iedgfl=0;}  // make sure not used unless setup

    void edgeSetup(int z =0);
    // setup above constants from tables for given atomic number
};

#endif

---