1 #include "Fdimpar.h" //(DIMPAR) fluka include
2 #include "Ftrackr.h" //(TRACKR) fluka common
3 #include "Fiounit.h" //(IOUNIT) fluka common
5 #include "TGeoMaterial.h"
6 #include "TFlukaCerenkov.h"
9 # define rfrndx rfrndx_
11 # define rfrndx RFRNDX
14 Double_t rfrndx(Double_t& wvlngt, Double_t& /*omgpho*/, Int_t& mmat)
17 // Return refraction index for given photon energy and material
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));