---

Photin.cxx

Go to the documentation of this file.

// -*- C++ -*-  $Id: Photin.cxx,v 1.1 2000/09/01 20:20:16 burnett Exp $
//
// This file is part of Gismo 2
//
#include "Photin.h"

void Photin::readPEGS(std::istream& fin, int m, int iray1)
{
   mge = m;

   fin >> ebinda >> ge0 >> ge1;

   gmfp0 = new float[mge];
   gmfp1 = new float[mge];
   gbr10 = new float[mge];
   gbr11 = new float[mge];
   gbr20 = new float[mge];
   gbr21 = new float[mge];

   int i;
   for (i=0; i< mge; i++)
   {
      fin >> gmfp0[i] >> gmfp1[i]
          >> gbr10[i] >> gbr11[i]
          >> gbr20[i] >> gbr21[i] ;
   }

   // optional Rayleigh scattering input

   if (iray1)
   {
      fin >> ngr >> rco0 >> rco1;
      cohe0 = new float[mge];
      cohe1 = new float[mge];
      rsct0 = new float[ngr];
      rsct1 = new float[ngr];

      for (i=0; i<ngr; i++)  fin >> rsct0[i] >> rsct1[i];
      for (i=0; i<mge; i++)  fin >> cohe0[i] >> cohe1[i];
   }else
          ngr = 0;

    
}

// Following is equivalent to the EGS4 macros
//        $set interval xxx,rco;
//        $evaluate x2 using rsct(xxx);

float Photin::rsct(float xxx)const
{
    int lxxx=static_cast<int>(rco1*xxx+rco0 -1);
    return rsct1[lxxx]*xxx+rsct0[lxxx];
}

void Photin::rescale(float dfact)
{
   for (int i=0; i<mge; i++)
   {
      gmfp0[i] *= dfact;
      gmfp1[i] *= dfact;
   }
}

Photin::~Photin()
{
    if( mge==0) return;
    delete [] gmfp0;
    delete [] gmfp1;
    delete [] gbr10;
    delete [] gbr11;
    delete [] gbr20;
    delete [] gbr21;
    if( ngr ) {
        delete [] cohe0;
        delete [] cohe1;
        delete [] rsct0;
        delete [] rsct1;
    }

}

---