Add CsPbBr3 and lead materials

JuliaPhysics · Nov 11, 2024 · cd712ac · cd712ac
1 parent 739f2ce
commit cd712ac
Show file tree

Hide file tree

Showing 3 changed files with 41 additions and 1 deletion.
diff --git a/ext/Geant4/io_gdml.jl b/ext/Geant4/io_gdml.jl
@@ -463,6 +463,8 @@ end
     elseif material == "Polyethylene Naphthalate" return "G4_PEN"
     elseif material == "Polytetrafluorethylen" return "G4_PTFE"
     elseif material == "Cadmium zinc telluride" return "G4_CdZnTe"
+    elseif material == "CsPbBr3" return "G4_CsPbBr3"
+    elseif material == "Lead" return "G4_Pb"
     end
     throw("Material characteristics for \"$(material)\" not defined yet")
 end
diff --git a/ext/Geant4/materials.xml b/ext/Geant4/materials.xml
@@ -64,4 +64,26 @@
         <fraction n="0.0278" ref="G4_Zn" />
         <fraction n="0.5422" ref="G4_Te" />
     </material>
+
+    <element name="G4_Cs" Z="55">
+        <atom value="132.9054" />
+    </element>
+    <element name="G4_Pb" Z="82">
+        <atom value="207.2" />
+    </element>
+    <element name="G4_Br" Z="35">
+        <atom value="79.904" />
+    </element>
+
+    <material name="G4_CsPbBr3" formula="CsPbBr3">
+        <D value="4.73" unit="g/cm3" />
+        <fraction n="0.2292" ref="G4_Cs" />
+        <fraction n="0.3574" ref="G4_Pb" />
+        <fraction n="0.4134" ref="G4_Br" />
+    </material>
+
+    <material name="G4_Pb" Z="82">
+        <D value="11.34" />
+        <atom value="207.2" />
+    </material>
 </materials>
diff --git a/src/MaterialProperties/MaterialProperties.jl b/src/MaterialProperties/MaterialProperties.jl
@@ -89,6 +89,21 @@ material_properties[:CdZnTe] = (
     ϵ_r = 10.9,
     ρ = 5.78u"g*cm^-3"
 )
+
+abstract type CsPbBr3 <: AbstractDriftMaterial end
+Symbol(::Type{CsPbBr3}) = :CsPbBr3
+material_properties[:CsPbBr3] = (
+    name = "CsPbBr3",
+    E_ionisation = 6.9u"eV",
+    ϵ_r = 16.46,
+    ρ = 4.73u"g*cm^-3"
+)
+
+material_properties[:Pb] = (
+    name = "Lead",
+    ϵ_r = 1e6, # high value for lead as for conductor. Usage of lead is shielding.
+    ρ = 11.35u"g*cm^-3"
+)
 # Add new materials above this line
 # and just put different spellings into the dict `materials` below
 
@@ -101,7 +116,8 @@ materials = Dict{String, Symbol}(
     "Al"  => :Al,
     "LAr" => :LAr,
     "CZT" => :CdZnTe,
-    "Si" => :Si
+    "Si" => :Si,
+    "Lead" => :Pb
 )