Warning, /geant4/examples/extended/exoticphysics/channeling/ch2/special_macros/run_Positron_Source.mac is written in an unsupported language. File is not indexed.
0001 #FOR SIMPLICITY THIS MACRO USES MODEL DEFAULTS;
0002 #TO MODIFY PLEASE REFER run.mac
0003 #AND ADD THE CORRESPONDING OPTIONS TO THIS MACRO
0004
0005 #THIS IS A SIMPLIFIED SETUP FOR A CRYSTAL-BASED POSITRON SOURCE
0006 #FOR FUTURE LEPTON COLLIDERS.
0007 #PLEASE REFER THE FOLLOWING PAPERS FOR MORE INFORMATION:
0008 #L. Bandiera et al. Eur. Phys. J. C 82, 699 (2022).
0009 #F. Alharthi et al. NIM A 1075, 170412 (2025).
0010
0011 #CAUTION: NO COHERENT PAIR PRODUCTION IS CONSIDERED,
0012 #FOR THIS, PLEASE REFER TO THE EXAMPLE ch3.
0013 #THIS IS NOT AN ISSUE BELOW 10 GeV FOR W.
0014
0015 #CAUTION: ONLY A SINGLE CRYSTAL TARGET IS CONSIDERED.
0016 #NO DEPOSITED ENERGY IS SCORED. FOR MORE FUNCTIONALITY
0017 #PLEASE REFER TO THE DEDICATED EXAMPLE ch5.
0018
0019 /random/setSeeds 19577794 424238336
0020 #setting number of cores
0021 /run/numberOfThreads 2
0022
0023 #crystal geometry
0024 #change it as you want
0025 /crystal/setCrystalSize 20. 20. 0.5 mm #feel free to increase the crystal thickness
0026 /crystal/setCrystalAngleX 0 mrad #default is 0
0027 /crystal/setCrystalAngleY 0 mrad #default is 0
0028
0029 #crystal lattice parameters
0030 #change it according to the data available in G4CHANNELINGDATA
0031 /crystal/setCrystalMaterial G4_W # check G4CHANNELINGDATA dataset
0032 /crystal/setCrystalLattice <111> # check G4CHANNELINGDATA dataset
0033
0034 #detector parameters
0035 #change it as you want
0036 /crystal/setDetectorSize 10 10 0.03 cm
0037 /crystal/setFrontPositionZ 1000 mm
0038
0039 #G4ChannelingFastSimModel parameters
0040 #change it as you want
0041 /crystal/setChannelingModel true #switch on/off G4ChannelingFastSimModel
0042 /crystal/setRadiationModel true #switch on/off G4BaierKatkov
0043 #(G4ChannelingFastSimModel must be switched on)
0044
0045 #the pass to channeling data if different from the default:
0046 #/crystal/setChannelingDataPath your_path
0047
0048 #the following can be changed or commented:
0049 #the low energy threshold for particle to enter the G4ChannelingFastSimModel:
0050 /crystal/setParticleMinKinEnergy/e+ 500 MeV #to speed up the simulations, please evaluate to change
0051 /crystal/setParticleMinKinEnergy/e- 500 MeV #to speed up the simulations, please evaluate to change
0052
0053
0054 #the following can be changed or commented:
0055 #high angular threshold for particle to enter the G4ChannelingFastSimModel expressed in Lindhard angles:
0056 /crystal/setLindhardAngles/e+ 10. #to speed up the simulations, please evaluate to change
0057 /crystal/setLindhardAngles/e- 10. #to speed up the simulations, please evaluate to change
0058
0059 #the following can be changed or commented:
0060 #G4BaierKatkov parameters
0061 #(see also comments in ConstructSDandField in DetectorConstruction):
0062 /crystal/setMinPhotonEnergy 1 MeV # G4BaierKatkov default
0063 /crystal/setSamplingPhotonsNumber 150 # G4BaierKatkov default
0064 /crystal/setNSmallTrajectorySteps 10000 # G4BaierKatkov default
0065 /crystal/setRadiationAngleFactor 4. # G4BaierKatkov default
0066
0067 /run/initialize
0068
0069 #settings
0070 /control/verbose 0
0071 /run/verbose 0
0072 /tracking/verbose 0
0073
0074 #beam
0075 /gps/particle e-
0076
0077 #coordinate distribution (radial Gauss)
0078 /gps/pos/centre 0. 0. -1. cm
0079 /gps/pos/type Beam
0080 /gps/pos/sigma_x 1. mm
0081 /gps/pos/sigma_y 1. mm
0082
0083 #angular distribution (radial Gauss)
0084 /gps/ang/type beam2d
0085 /gps/ang/rot1 1 0 0
0086 /gps/ang/rot2 0 -1 0
0087 /gps/ang/sigma_x 10.E-6 rad
0088 /gps/ang/sigma_y 10.E-6 rad
0089
0090 #energy distribution (constant)
0091 /gps/ene/mono 6. GeV
0092
0093 #statistics
0094 /run/printProgress 10
0095 /run/beamOn 50