Commit cb1acc60 authored by mphipps2's avatar mphipps2
Browse files

Materials is updated but something seems broken in RPD (simulation hangs when...

Materials is updated but something seems broken in RPD (simulation hangs when RPD included). appears further back than my changes today. most likely in ModTypeRPD
parent 955aa332
...@@ -20,7 +20,7 @@ ...@@ -20,7 +20,7 @@
/Detector/ZDC/SteelAbsorberHeight 290. mm /Detector/ZDC/SteelAbsorberHeight 290. mm
/Detector/ZDC/Optical false /Detector/ZDC/Optical false
/Detector/ZDC/CheckOverlaps false /Detector/ZDC/CheckOverlaps false
/Detector/ZDC/ReducedTree false /Detector/ZDC/ReducedTree true
/Detector/ZDC/HousingMaterial aluminum /Detector/ZDC/HousingMaterial aluminum
/Detector/ZDC/AbsorberMaterial pure /Detector/ZDC/AbsorberMaterial pure
...@@ -45,4 +45,4 @@ ...@@ -45,4 +45,4 @@
/Detector/RPD/RPDtype panflute /Detector/RPD/RPDtype panflute
/Detector/RPD/Optical true /Detector/RPD/Optical true
/Detector/RPD/CheckOverlaps false /Detector/RPD/CheckOverlaps false
/Detector/RPD/ReducedTree false /Detector/RPD/ReducedTree true
...@@ -132,17 +132,17 @@ void Materials::DefineOpticalProperties(void){ ...@@ -132,17 +132,17 @@ void Materials::DefineOpticalProperties(void){
G4double kapton_RIND_UI[nEntries_UI] = {1.8, 1.82, 1.86, 1.95, 2.18}; G4double kapton_RIND_UI[nEntries_UI] = {1.8, 1.82, 1.86, 1.95, 2.18};
G4double kapton_ABSL_UI[nEntries_UI] = {0.011*cm, 0.0022*cm, 0.0056*cm, 0.04*cm, 0.063*cm}; G4double kapton_ABSL_UI[nEntries_UI] = {0.011*cm, 0.0022*cm, 0.0056*cm, 0.04*cm, 0.063*cm};
G4double kapton_REFL_UI[nEntries_UI] = {0.35, 0.2, 0.02, 0.03, 0.04}; G4double kapton_REFL_UI[nEntries_UI] = {0.35, 0.2, 0.02, 0.03, 0.04};
MPT_Array.push_back(new G4MaterialPropertiesTable()); MPT_Array.push_back(new G4MaterialPropertiesTable());
MPT_Array.back()->AddProperty("RINDEX",photonEnergy_UI,kapton_RIND_UI,nEntries_UI); MPT_Array.back()->AddProperty("RINDEX",photonEnergy_UI,kapton_RIND_UI,nEntries_UI);
MPT_Array.back()->AddProperty("ABSLENGTH",photonEnergy_UI,kapton_ABSL_UI,nEntries_UI); MPT_Array.back()->AddProperty("ABSLENGTH",photonEnergy_UI,kapton_ABSL_UI,nEntries_UI);
MPT_Array.back()->AddProperty("REFLECTIVITY",photonEnergy_UI,kapton_REFL_UI,nEntries_UI); MPT_Array.back()->AddProperty("REFLECTIVITY",photonEnergy_UI,kapton_REFL_UI,nEntries_UI);
Kapton_UI->SetMaterialPropertiesTable(MPT_Array.back()); Kapton_UI->SetMaterialPropertiesTable(MPT_Array.back());
// https://refractiveindex.info/?shelf=main&book=W&page=Werner // https://refractiveindex.info/?shelf=main&book=W&page=Werner
G4double NiW_REFL_UI[nEntries_UI] = {0.62, 0.62, 0.62, 0.62, 0.62}; G4double NiW_REFL_UI[nEntries_UI] = {0.62, 0.62, 0.62, 0.62, 0.62};
MPT_Array.push_back(new G4MaterialPropertiesTable()); MPT_Array.push_back(new G4MaterialPropertiesTable());
MPT_Array.back()->AddProperty("REFLECTIVITY",photonEnergy_UI,NiW_REFL_UI,nEntries_UI); MPT_Array.back()->AddProperty("REFLECTIVITY",photonEnergy_UI,NiW_REFL_UI,nEntries_UI);
NiW->SetMaterialPropertiesTable(MPT_Array.back()); NiW->SetMaterialPropertiesTable(MPT_Array.back());
//Air optical properties //Air optical properties
...@@ -150,13 +150,13 @@ void Materials::DefineOpticalProperties(void){ ...@@ -150,13 +150,13 @@ void Materials::DefineOpticalProperties(void){
MPT_Array.push_back(new G4MaterialPropertiesTable()); MPT_Array.push_back(new G4MaterialPropertiesTable());
MPT_Array.back()->AddProperty("RINDEX", photonEnergy_UI, RefractiveIndexAir, nEntries_UI); MPT_Array.back()->AddProperty("RINDEX", photonEnergy_UI, RefractiveIndexAir, nEntries_UI);
Air->SetMaterialPropertiesTable(MPT_Array.back()); Air->SetMaterialPropertiesTable(MPT_Array.back());
//Aluminum optical properties //Aluminum optical properties
G4double Al_refl[nEntries_UI] = {0.89, 0.89, 0.89, 0.89, 0.89};//"Hard cut" G4double Al_refl[nEntries_UI] = {0.89, 0.89, 0.89, 0.89, 0.89};//"Hard cut"
MPT_Array.push_back(new G4MaterialPropertiesTable()); MPT_Array.push_back(new G4MaterialPropertiesTable());
MPT_Array.back()->AddProperty("REFLECTIVITY",photonEnergy_UI,Al_refl,nEntries_UI); MPT_Array.back()->AddProperty("REFLECTIVITY",photonEnergy_UI,Al_refl,nEntries_UI);
Al->SetMaterialPropertiesTable(MPT_Array.back()); Al->SetMaterialPropertiesTable(MPT_Array.back());
// Maryland RPD Materials // Maryland RPD Materials
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment