---

Elecin.h

Go to the documentation of this file.

// -*- C++ -*-  $Id: Elecin.h,v 1.1 2000/09/01 20:20:15 burnett Exp $
//
// This file is part of Gismo 2



#ifndef __ELECIN_H
#define __ELECIN_H

#include "facilities/error.h"

#include <iostream>

class Elecin  // electron cross sections in this material
{
   friend class PEGSData;    // this only accessible from PEGSData

   float eold,   // last energy (in MeV)
         loge;   // log of E
   int   index;  // index

   float eke0,eke1;

   float xr0,teff0,blcc,xcc;

   int   meke ;  // actual dimension following arrays
   float  *esig0,  *esig1,
          *psig0,  *psig1,
          *ededx0, *ededx1,
          *pdedx0, *pdedx1,
          *ebr10,  *ebr11,
          *pbr10,  *pbr11,
          *pbr20,  *pbr21,
          *tmxs0,  *tmxs1;

   Elecin():meke(0){eold=-1.;} // constructor

   void updateIndex(float enew) const;
   // calculate new index if energy changes

   float  esig()const;
   float  psig()const;
   float ededx()const;
   float pdedx()const;
   float  ebr1()const;
   float  pbr1()const;
   float  pbr2()const;
   float  tmxs()const;
   // the interpolating functions (should have designed a macro)

   void readPEGS(std::istream&, int);
   void rescale(float);
   void fixtmx(float);
   // special routine that recalculates tmxs

public:
    ~Elecin();
};
/*
inline void
Elecin::updateIndex(float enew) const
{
   if (enew == eold) return;
   Elecin * self = (Elecin *)this;
   self->loge = log(enew);
   self->index= (int)(eke1*loge + eke0 -1); // note -1 for C arrays
   if (index>=meke)
   {  FATAL("Electin: energy too large for table");
   }
   self->eold = enew;
}
inline float Elecin::esig()const{return  esig1[index]*loge + esig0[index];}
inline float Elecin::psig()const{return  psig1[index]*loge + psig0[index];}
inline float Elecin::ededx()const{return ededx1[index]*loge +ededx0[index];}
inline float Elecin::pdedx()const{return pdedx1[index]*loge +pdedx0[index];}
inline float Elecin::ebr1()const{return  ebr11[index]*loge + ebr10[index];}
inline float Elecin::pbr1()const{return  pbr11[index]*loge + pbr10[index];}
inline float Elecin::pbr2()const{return  pbr21[index]*loge + pbr20[index];}
inline float Elecin::tmxs()const{return  tmxs1[index]*loge + tmxs0[index];}
 */
#endif

---