Files
geant4/examples/extended/optical/OpNovice2/scint_by_particle.mac
2023-12-08 10:43:34 +01:00

106 lines
4.3 KiB
Plaintext

/control/verbose 2
/tracking/verbose 0
/process/optical/verbose 1
/run/verbose 1
/control/cout/ignoreThreadsExcept 0
/opnovice2/boxMaterial G4_BGO
/opnovice2/worldMaterial G4_AIR
/opnovice2/boxProperty RINDEX 0.000002 1.3 0.000008 1.4
/opnovice2/boxProperty ABSLENGTH 0.000002 1 0.000005 2 0.000008 3
/opnovice2/boxProperty SCINTILLATIONCOMPONENT1 0.000002 1.0 0.000005 1.1 0.000008 1.3
/opnovice2/boxProperty SCINTILLATIONCOMPONENT2 0.000002 0.8 0.000005 0.6 0.000008 0.3
/opnovice2/boxConstProperty SCINTILLATIONTIMECONSTANT1 20 ## ns
/opnovice2/boxConstProperty SCINTILLATIONTIMECONSTANT2 200 ## ns
/opnovice2/boxConstProperty SCINTILLATIONRISETIME1 5 ## ns
/opnovice2/boxConstProperty SCINTILLATIONRISETIME2 10 ## ns
/opnovice2/boxConstProperty RESOLUTIONSCALE 1
#/opnovice2/boxProperty SCINTILLATIONYIELD 5000.
/opnovice2/boxProperty PROTONSCINTILLATIONYIELD 0 50 10 5000
/opnovice2/boxConstProperty PROTONSCINTILLATIONYIELD1 0.3
/opnovice2/boxConstProperty PROTONSCINTILLATIONYIELD2 0.5
/opnovice2/boxConstProperty PROTONSCINTILLATIONYIELD3 0.2
/opnovice2/boxConstProperty PROTONSCINTILLATIONTIMECONSTANT1 50
/opnovice2/boxConstProperty PROTONSCINTILLATIONTIMECONSTANT2 200
/opnovice2/boxConstProperty PROTONSCINTILLATIONTIMECONSTANT3 500
/opnovice2/boxProperty DEUTERONSCINTILLATIONYIELD 0 50 10 5000
/opnovice2/boxConstProperty DEUTERONSCINTILLATIONYIELD1 0.6
/opnovice2/boxConstProperty DEUTERONSCINTILLATIONYIELD2 0.4
# if time constants aren't specified for a particle, SCINTILLATIONTIMECONSTANTs
# will be used for that particle
#/opnovice2/boxConstProperty DEUTERONSCINTILLATIONTIMECONSTANT1 50
#/opnovice2/boxConstProperty DEUTERONSCINTILLATIONTIMECONSTANT2 100
/opnovice2/boxProperty TRITONSCINTILLATIONYIELD 0 50 10 5000
/opnovice2/boxConstProperty TRITONSCINTILLATIONYIELD1 0.5
/opnovice2/boxConstProperty TRITONSCINTILLATIONYIELD2 0.5
/opnovice2/boxConstProperty TRITONSCINTILLATIONTIMECONSTANT1 100
/opnovice2/boxConstProperty TRITONSCINTILLATIONTIMECONSTANT2 500
/opnovice2/boxProperty ALPHASCINTILLATIONYIELD 0 50 10 50000
/opnovice2/boxConstProperty ALPHASCINTILLATIONYIELD1 0.7
/opnovice2/boxConstProperty ALPHASCINTILLATIONYIELD2 0.3
/opnovice2/boxConstProperty ALPHASCINTILLATIONTIMECONSTANT1 30
/opnovice2/boxConstProperty ALPHASCINTILLATIONTIMECONSTANT2 60
/opnovice2/boxProperty IONSCINTILLATIONYIELD 0 50 10 5000
/opnovice2/boxConstProperty IONSCINTILLATIONYIELD1 1.0
/opnovice2/boxConstProperty IONSCINTILLATIONYIELD2 0.0
/opnovice2/boxConstProperty IONSCINTILLATIONTIMECONSTANT1 20
/opnovice2/boxConstProperty IONSCINTILLATIONTIMECONSTANT2 100
/opnovice2/boxProperty ELECTRONSCINTILLATIONYIELD 0 5000 10 500000
/opnovice2/boxConstProperty ELECTRONSCINTILLATIONYIELD1 0.5
/opnovice2/boxConstProperty ELECTRONSCINTILLATIONYIELD2 0.5
/opnovice2/boxConstProperty ELECTRONSCINTILLATIONTIMECONSTANT1 50
/opnovice2/boxConstProperty ELECTRONSCINTILLATIONTIMECONSTANT2 500
/opnovice2/worldProperty RINDEX 0.000002 1.01 0.000008 1.01
/opnovice2/worldProperty ABSLENGTH 0.000002 100 0.000005 100 0.000008 100
## no scintillation in the world, but we want to exercise
## this macro command
/opnovice2/worldConstProperty SCINTILLATIONTIMECONSTANT1 1
/process/optical/processActivation Cerenkov false
/process/optical/processActivation OpAbsorption true
/process/optical/processActivation OpBoundary false
/process/optical/processActivation Scintillation true
/process/optical/processActivation OpRayleigh false
/process/optical/processActivation OpMieHG false
/process/optical/processActivation OpWLS false
/process/optical/scintillation/setByParticleType true
/process/optical/scintillation/setTrackInfo false
/process/optical/scintillation/setFiniteRiseTime false
/process/optical/scintillation/setStackPhotons true
/process/optical/scintillation/setTrackSecondariesFirst true
/run/initialize
/analysis/setFileName scint_by_particle_e
/analysis/h1/set 2 100 0 10 # eV
/analysis/h1/set 3 400 0 200 # ns
#
/gun/particle e-
/gun/energy 1 MeV
/gun/position -1 0 0 m
/gun/direction 1 0 0
/run/beamOn 100
/gun/particle proton
/analysis/setFileName scint_by_particle_p
/run/beamOn 100
/gun/particle alpha
/analysis/setFileName scint_by_particle_a
/run/beamOn 100
/gun/particle deuteron
/analysis/setFileName scint_by_particle_d
/run/beamOn 100
/gun/particle ion
/analysis/setFileName scint_by_particle_i
/gun/ion 6 12 6
/run/beamOn 100