//
//
if (fMedium > -1) {
- fprintf(fgFile,"*\n*Material specific process and cut settings for #%8d \n", fMedium);
+ TFluka* fluka = (TFluka*) gMC;
+ TObjArray *matList = fluka->GetFlukaMaterials();
+ Int_t nmaterial = matList->GetEntriesFast();
+ TGeoMaterial* material = 0;
+ for (Int_t im = 0; im < nmaterial; im++)
+ {
+ material = dynamic_cast<TGeoMaterial*> (matList->At(im));
+ Int_t idmat = material->GetIndex();
+ if (idmat == fMedium) break;
+ }
+
+
+//
+// Check if global option
+
+ fprintf(fgFile,"*\n*Material specific process and cut settings for #%8d %s\n", fMedium, material->GetName());
fCMatMin = fMedium;
fCMatMax = fMedium;
} else {
Float_t neutronCut = cut;
if (neutronCut < 0.0196) {
neutronCut = 0.0196;
- printf("Cut on neutron lower than upper limit if first energy group.\n");
+ printf("Cut on neutron lower than upper limit of first energy group.\n");
printf("Cut reset to 19.6 MeV !\n");
}
fprintf(fgFile,"PART-THR %10.4g%10.1f%10.1f\n", -neutronCut, 8.0, 9.0);
protected:
Double_t fCutValue[11]; // User cut
Int_t fProcessFlag[15]; // User flag assigned to processes
- Int_t fMedium; // Materials assigned to user settings
+ Int_t fMedium; // Material assigned to user settings
Float_t fCMatMin; // Minimum material number used for current card
Float_t fCMatMax; // Maximum material number used for current card
}
Float_t hasfield = (vol->GetMedium()->GetParam(1) > 0) ? flagfield : 0.;
-
+ out << "* Assigning material: " << vol->GetMedium()->GetMaterial()->GetName() << " to Volume: " << vol->GetName();
+ out << endl;
+
out << setw(10) << "ASSIGNMAT ";
out.setf(static_cast<std::ios::fmtflags>(0),std::ios::floatfield);
out << setw(10) << setiosflags(ios::fixed) << Double_t(idmat);