TFluka/rfrndx.cxx

   1 #include "Fdimpar.h"  //(DIMPAR) fluka include
   2 #include "Ftrackr.h"  //(TRACKR) fluka common
   3 #include "Fiounit.h"  //(IOUNIT) fluka common
   4 #include "TFluka.h"
   5 #include "TGeoMaterial.h"
   6 #include "TFlukaCerenkov.h"
   7
   8 #ifndef WIN32
   9 # define rfrndx rfrndx_
  10 #else
  11 # define rfrndx RFRNDX
  12 #endif
  13 extern "C" {
  14 Double_t rfrndx(Double_t& wvlngt, Double_t& /*omgpho*/, Int_t& mmat)
  15 {
  16 //
  17 //  Return refraction index for given photon energy and material
  18 //
  19     TFluka* fluka =  (TFluka*) gMC;
  20     TGeoMaterial*    material =  (TGeoMaterial*) (fluka->GetFlukaMaterials())->At(fluka->GetMaterialIndex(mmat));
  21     TFlukaCerenkov*  cerenkov = dynamic_cast<TFlukaCerenkov*> (material->GetCerenkovProperties());
  22     Double_t y =  (cerenkov->GetRefractionIndexByWaveLength(wvlngt));
  23     return (y);
  24 }
  25 }