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 abscff abscff_
12 # define abscff ABSCFF
15 Double_t abscff(Double_t& wvlngt, Double_t& /*omgpho*/, Int_t& mmat)
17 // printf("abscff%f %d\n", wvlngt, mmat);
20 // Return absorption length for given photon energy and material
23 TFluka* fluka = (TFluka*) gMC;
25 // Check if stopping has been required by user
27 if (fluka->GetStoppingCondition()) {
28 fluka->ResetStoppingCondition();
32 // Get absorption coefficient for current material
34 TGeoMaterial* material = (TGeoMaterial*) (fluka->GetFlukaMaterials())->At(fluka->GetMaterialIndex(mmat));
35 TFlukaCerenkov* cerenkov = dynamic_cast<TFlukaCerenkov*> (material->GetCerenkovProperties());
36 Double_t y = (cerenkov->GetAbsorptionCoefficientByWaveLength(wvlngt));