From cd712acb92e9f888f658c82a12c744361bd3689b Mon Sep 17 00:00:00 2001
From: Padniuk <ypadniuk@gmail.com>
Date: Mon, 11 Nov 2024 15:00:22 +0100
Subject: [PATCH] Add CsPbBr3 and lead materials

---
 ext/Geant4/io_gdml.jl                        |  2 ++
 ext/Geant4/materials.xml                     | 22 ++++++++++++++++++++
 src/MaterialProperties/MaterialProperties.jl | 18 +++++++++++++++-
 3 files changed, 41 insertions(+), 1 deletion(-)

diff --git a/ext/Geant4/io_gdml.jl b/ext/Geant4/io_gdml.jl
index f79b9410f..6d937a19b 100644
--- 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
index d996b6a8e..31cd127c0 100644
--- 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>
\ No newline at end of file
diff --git a/src/MaterialProperties/MaterialProperties.jl b/src/MaterialProperties/MaterialProperties.jl
index 444118fe4..5da89d6c6 100644
--- 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
 )