- printf("rflctv called %e %e %d \n", wvlngt, omgpho, mmat);
- return (0.);
+//
+// Return reflectivity (1-r) for given photon energy and material
+//
+ TFluka* fluka = (TFluka*) gMC;
+ TGeoMaterial* material = (TGeoMaterial*) (fluka->GetFlukaMaterials())->At(fluka->GetMaterialIndex(mmat));
+ TFlukaCerenkov* cerenkov = dynamic_cast<TFlukaCerenkov*> (material->GetCerenkovProperties());
+ Double_t y = 0.;
+ if (cerenkov->IsMetal()) y = (cerenkov->GetReflectivityByWaveLength(wvlngt));
+ return (y);