# Macro file for geant4 simulation
                                  
# Set the random seeds 
/seed/setSeeds 1151113622152 1151115496606
                                  
# 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
               
# General Particle Source (GPS)
# Flat distribution
/gps/verbose 0
/gps/particle mu-
/gps/pos/type Point
/gps/pos/centre 0 0 0 cm
/gps/ang/type iso
/gps/ang/mintheta 165.0 deg
/gps/ang/maxtheta 180.0 deg
/gps/ang/minphi 0. deg
/gps/ang/maxphi 360. deg
/gps/energy 100 GeV
                    
# Output
/THGCal/hitwriter/writerTxtDataFileName sensors_d2_mum_E100_theta0-15_14.dat
               
/THGCal/hitwriter/writerPrimPartInfosFileName primaryParticles_d2_mum_E100_theta0-15_14.dat
               
/run/beamOn 10000