# Macro file for geant4 simulation
                                  
# Set the random seeds 
/seed/setSeeds 1706416031822 1706417861852
                                  
# Set the geometry visualisation 
# if setGeomVisu=true, /run/beamOn should not be used 
/THGCal/det/setGeomVisu false
                                  
# Set the flag to write the geometry 
/THGCal/det/setWriteGeom false
# Set the name of the geometry file 
/THGCal/det/setGeomFileName geom_d2.dat
                                  
# ECAL PART  
## General information 
/THGCal/det/setEcalNumberOfLayers 26
/THGCal/det/setEcalRingsNumberPerLayer 35 35 35 35 35 35 35 35 35 35 35 35 35 35 35 35 35 35 35 35 35 35 35 35 35 35
                                        
## The ecal cells dimensions are in mm
/THGCal/det/setEcalCellSizePerLayer 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10
                                        
## The ecal absorber dimensions are in mm
/THGCal/det/setEcalAbsThicknessPerLayer 6.05 6.05 6.05 6.05 6.05 6.05 6.05 6.05 6.05 6.05 6.05 6.05 6.05 6.05 6.05 6.05 6.05 6.05 6.05 6.05 6.05 6.05 6.05 6.05 6.05 6.05
/THGCal/det/setEcalAbsMaterialNamePerLayer G4_Pb G4_Pb G4_Pb G4_Pb G4_Pb G4_Pb G4_Pb G4_Pb G4_Pb G4_Pb G4_Pb G4_Pb G4_Pb G4_Pb G4_Pb G4_Pb G4_Pb G4_Pb G4_Pb G4_Pb G4_Pb G4_Pb G4_Pb G4_Pb G4_Pb G4_Pb
               
# HADRONIC PART 
## General information 
/THGCal/det/setHadNumberOfLayers 24
/THGCal/det/setHadRingsNumberPerLayer 25 25 25 25 25 25 25 25 25 25 25 25 25 25 25 25 25 25 25 25 25 25 25 25
                                        
## The had cells dimensions are in mm
/THGCal/det/setHadCellSizePerLayer 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20
                                        
## The hadronic absorber dimensions are in mm
/THGCal/det/setHadAbsThicknessPerLayer 45.0 45.0 45.0 45.0 45.0 45.0 45.0 45.0 45.0 45.0 45.0 45.0 80.0 80.0 80.0 80.0 80.0 80.0 80.0 80.0 80.0 80.0 80.0 80.0
/THGCal/det/setHadAbsMaterialNamePerLayer G4_STAINLESS-STEEL G4_STAINLESS-STEEL G4_STAINLESS-STEEL G4_STAINLESS-STEEL G4_STAINLESS-STEEL G4_STAINLESS-STEEL G4_STAINLESS-STEEL G4_STAINLESS-STEEL G4_STAINLESS-STEEL G4_STAINLESS-STEEL G4_STAINLESS-STEEL G4_STAINLESS-STEEL G4_STAINLESS-STEEL G4_STAINLESS-STEEL G4_STAINLESS-STEEL G4_STAINLESS-STEEL G4_STAINLESS-STEEL G4_STAINLESS-STEEL G4_STAINLESS-STEEL G4_STAINLESS-STEEL G4_STAINLESS-STEEL G4_STAINLESS-STEEL G4_STAINLESS-STEEL G4_STAINLESS-STEEL
                                        
# Detector position, distance between the world volume and the center of the forward plan of the first absorber
/THGCal/det/setDetPosition 0.0 0.0 1000.0 mm
/THGCal/det/setLayerSpacing 3.0 mm
               
# Set the number of threads
/run/numberOfThreads 8
               
# Initialize kernel
/run/initialize
               
# Gun particle
/gun/particle pi+
# Momentum direction is a 3Vector
/gun/direction 0 0 1
# Particle energy, default unit is GeV
/gun/energy 100.0
#Set starting position of the particle, 3Vector, default unit cm
/gun/position 0 0 0
# Number of particles to be generated
/gun/number 1
                    
# Output
/THGCal/hitwriter/writerTxtDataFileName sensors_d2_pip_E100_theta0_98.dat
               
/THGCal/hitwriter/writerPrimPartInfosFileName primaryParticles_d2_pip_E100_theta0_98.dat
               
/run/beamOn 10000