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"
10 # define rflctv rflctv_
12 # define rflctv RFLCTV
15 Double_t rflctv(Double_t& wvlngt, Double_t& /*omgpho*/, Int_t& mmat)
18 // Return reflectivity (1-r) for given photon energy and material
20 TFluka* fluka = (TFluka*) gMC;
21 TGeoMaterial* material = (TGeoMaterial*) (fluka->GetFlukaMaterials())->At(fluka->GetMaterialIndex(mmat));
22 TFlukaCerenkov* cerenkov = dynamic_cast<TFlukaCerenkov*> (material->GetCerenkovProperties());
24 if (cerenkov->IsMetal()) y = (cerenkov->GetReflectivityByWaveLength(wvlngt));