beam replaced by beamcm
[u/mrichter/AliRoot.git] / 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 }