#
# G4SEE Single Event Effect simulation toolkit
# ============================================
# SPDX-FileCopyrightText: © 2022 CERN for the benefit of the G4SEE Collaboration <https://cern.ch/g4see>
# SPDX-License-Identifier: CC-BY-4.0
#
##########################################
### Example input macro file for G4SEE ###
##########################################

/run/numberOfThreads                4
/run/printProgress                  1000
/tracking/verbose                   0                   # use only for testing

#############
### Geometry
/SEE/geometry/setCylindricalGeometry   true

# BULK                              MATERIAL            WIDTH unit  THICK unit  BIAS
/SEE/geometry/Bulk                  G4_Si               8 mm        308 um      true
# SV                                POSITION unit       WIDTH unit  THICK unit  BIAS
/SEE/geometry/SV                    0 0 -0.5 um         7. 7. mm    307 um      true
# BEOL                              MATERIAL            WIDTH unit  THICK unit  BIAS    NAME
/SEE/geometry/BEOL/addLayer         G4_Al               8 mm        500 nm      false   Metal

##############
### Biasing
/SEE/biasing/biasParticle           neutron
/SEE/biasing/biasProcess            hadElastic
/SEE/biasing/biasProcess            neutronInelastic
/SEE/biasing/biasProcess            nCapture
/SEE/biasing/biasFactor             1000

##############
### Physics
/SEE/physics/addPhysics             G4EmStandardPhysics_option4         # Default EM physics
/SEE/physics/addPhysics             G4HadronElasticPhysicsHP            # hadron elastic physics w/ HP neutron models
/SEE/physics/addPhysics             G4HadronPhysicsFTFP_BERT_HP         # hadron inelastic physics w/ HP neutron models

# Particle production range cuts
/SEE/physics/setGammaCut            1 mm                                # only gamma
/SEE/physics/setElectronCut         1 um                                # only e-
/SEE/physics/setPositronCut         1 um                                # only e+
/SEE/physics/setHadronCut           1 um                                # applicable to protons and all other hadrons


##############
### Initialize
/run/initialize

##############
### Primaries (General Particle Source)
/gps/particle                       neutron
/gps/ene/type                       Gauss
/gps/ene/mono                       15 MeV
/gps/ene/sigma                      75 keV
/gps/direction                      0 0 -1
/gps/pos/centre                     0 0 1 mm
/gps/pos/type                       Point

##############
### Scoring
/SEE/scoring/addScoring             0       Edep
/SEE/scoring/setHistogram           lin     0 eV    10 MeV      200
/SEE/scoring/addScoring             1       Edep
/SEE/scoring/setHistogram           log     10 eV   10 MeV      200

/SEE/scoring/addScoring             2       Ekin    neutron
/SEE/scoring/setHistogram           log     10 eV   15 MeV      200
/SEE/scoring/addScoring             3       Ekin    gamma
/SEE/scoring/setHistogram           log     1 keV   10 MeV      200

/SEE/scoring/dumpHistogramsAfter    5000

##############
### Run
/run/beamOn                         10000