Warning, /geant4/examples/extended/exoticphysics/channeling/ch2/special_macros/run_Radiation.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 #THESE PARAMETERS REPRESENT GAMMA RADIATION SOURCE
0006 #SUCH AS AXIAL OR PLANAR CHANNELING RADIATION
0007 #AND COHERENT BREMSSTRAHLUNG:
0008 #B. Ferretti, Il Nuovo Cimento 7, 118–134 (1950)
0009 #M.L. Ter-Mikaelian, High-Energy Electromagnetic Processes in Condensed Media,
0010 # Wiley-Interscience, New York, 1972.
0011 #V. N. Baier, V. M. Katkov, V. M. Strakhovenko et al.,
0012 #Electromagnetic Processes At High Energies In Oriented Single
0013 #Crystals, World Scientific Publishing Company, 1998.
0014 #V.G. Baryshevskii, I.Y. Dubovskaya, Sov. Phys. Dokl. 21, 741–743 (1976).
0015
0016 /random/setSeeds 19577794 424238336
0017 #setting number of cores
0018 /run/numberOfThreads 2
0019
0020 #crystal geometry
0021 #change it as you want
0022 /crystal/setCrystalSize 20. 20. 0.05 mm
0023 /crystal/setCrystalAngleX 0. mrad #for coherent bremsstrahlung change the angle (usually few mrad would work fine)
0024 /crystal/setCrystalAngleY 0. mrad #default is 0
0025
0026 #crystal lattice parameters
0027 #change it according to the data available in G4CHANNELINGDATA
0028 /crystal/setCrystalMaterial G4_Diamond # check G4CHANNELINGDATA dataset;
0029 # G4_Si and G4_Ge are also available.
0030 # For less monochromatic but more intense
0031 # radiation use G4_W.
0032 /crystal/setCrystalLattice (110) # check G4CHANNELINGDATA dataset;
0033 # planes (110) and (100) as well as
0034 # axes <100>, <110> and <111> are also available
0035 # For G4_W check G4CHANNELINGDATA dataset
0036 # for the availability.
0037
0038 #detector parameters
0039 #change it as you want
0040 /crystal/setDetectorSize 10 10 0.03 cm
0041 /crystal/setFrontPositionZ 1000 mm
0042
0043 #G4ChannelingFastSimModel parameters
0044 #change it as you want
0045 /crystal/setChannelingModel true #switch on/off G4ChannelingFastSimModel
0046 /crystal/setRadiationModel true #switch on/off G4BaierKatkov
0047 #(G4ChannelingFastSimModel must be switched on)
0048
0049 #the following can be changed or commented:
0050 #the low energy threshold for particle to enter the G4ChannelingFastSimModel:
0051 /crystal/setParticleMinKinEnergy/e- 200 MeV # for energies of primaries below 200 MeV
0052 # you may decrease the threshold;
0053 # be careful with the energies
0054 # below 100 MeV where quantum effects
0055 # particle dynamics may become important.
0056 # More information -
0057 # Baier, Katkov, Strakhovenko book (above)
0058 /crystal/setParticleMinKinEnergy/e+ 200 MeV # the same as for e-.
0059
0060 #the following can be changed or commented:
0061 #high angular threshold for particle to enter the G4ChannelingFastSimModel expressed in Lindhard angles:
0062 /crystal/setLindhardAngles/e+ 100. #change if necessary to speed up the simulations
0063 /crystal/setLindhardAngles/e- 100. #change if necessary to speed up the simulations
0064
0065 #the following can be changed or commented:
0066 #G4BaierKatkov parameters
0067 #(see also comments in ConstructSDandField in DetectorConstruction):
0068 /crystal/setMinPhotonEnergy 1 MeV # G4BaierKatkov default
0069 /crystal/setSamplingPhotonsNumber 150 # G4BaierKatkov default
0070 /crystal/setNSmallTrajectorySteps 10000 # G4BaierKatkov default
0071 /crystal/setRadiationAngleFactor 4. # G4BaierKatkov default
0072
0073 #the following can be changed or commented:
0074 #For scoring purposes only: the maximal energy in the output spectrum from the Baier-Katkov method.
0075 #Note it does not modifies the simulations.
0076 #Note: the minimal energy in the spectrum = the minimal energy defined above.
0077 /crystal/MaxBKPhotonEnergyInSpectrum 0.5 GeV # must be > setMinPhotonEnergy
0078 /crystal/NBinsInSpectrum 50 # note these bins are not equidistant.
0079 #Virtual collimator sets the selection of photon angles to be written in the output spectrum
0080 /crystal/setVirtualCollimatorHalfSize 1. mrad #setting the angular radius of virtual collimator
0081
0082 /run/initialize
0083
0084 #settings
0085 /control/verbose 0
0086 /run/verbose 0
0087 /tracking/verbose 0
0088
0089 #beam
0090 /gps/particle e- #e+ also work fine
0091
0092 #coordinate distribution (radial Gauss)
0093 /gps/pos/centre 0. 0. -1. cm
0094 /gps/pos/type Beam
0095 /gps/pos/sigma_x 0.1 mm
0096 /gps/pos/sigma_y 0.1 mm
0097
0098 #angular distribution (radial Gauss)
0099 /gps/ang/type beam2d
0100 /gps/ang/rot1 1 0 0
0101 /gps/ang/rot2 0 -1 0
0102 /gps/ang/sigma_x 50.E-6 rad
0103 /gps/ang/sigma_y 50.E-6 rad
0104
0105 #energy distribution (constant)
0106 /gps/ene/mono 0.5 GeV
0107
0108 #statistics
0109 /run/printProgress 1000
0110 /run/beamOn 10000
