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 "TGeoManager.h"
7 #include "TFlukaCerenkov.h"
10 # define rfrndx rfrndx_
12 # define rfrndx RFRNDX
15 Double_t rfrndx(Double_t& wvlngt, Double_t& /*omgpho*/, Int_t& mmat)
18 // Return refraction index for given photon energy and material
20 TFluka* fluka = (TFluka*) gMC;
21 TGeoMaterial* material = (TGeoMaterial*) (gGeoManager->GetListOfMaterials())->At(fluka->GetMaterialIndex(mmat));
22 TFlukaCerenkov* cerenkov = dynamic_cast<TFlukaCerenkov*> (material->GetCerenkovProperties());
23 Double_t y = (cerenkov->GetRefractionIndexByWaveLength(wvlngt));