Gabriele Cosmo
111 lines
4.6 KiB
Plaintext
111 lines
4.6 KiB
Plaintext
#FOR SIMPLICITY THIS MACRO USES MODEL DEFAULTS;
|
||
#TO MODIFY PLEASE REFER run.mac
|
||
#AND ADD THE CORRESPONDING OPTIONS TO THIS MACRO
|
||
|
||
#THESE PARAMETERS REPRESENT GAMMA RADIATION SOURCE
|
||
#SUCH AS AXIAL OR PLANAR CHANNELING RADIATION
|
||
#AND COHERENT BREMSSTRAHLUNG:
|
||
#B. Ferretti, Il Nuovo Cimento 7, 118–134 (1950)
|
||
#M.L. Ter-Mikaelian, High-Energy Electromagnetic Processes in Condensed Media,
|
||
# Wiley-Interscience, New York, 1972.
|
||
#V. N. Baier, V. M. Katkov, V. M. Strakhovenko et al.,
|
||
#Electromagnetic Processes At High Energies In Oriented Single
|
||
#Crystals, World Scientific Publishing Company, 1998.
|
||
#V.G. Baryshevskii, I.Y. Dubovskaya, Sov. Phys. Dokl. 21, 741–743 (1976).
|
||
|
||
/random/setSeeds 19577794 424238336
|
||
#setting number of cores
|
||
/run/numberOfThreads 2
|
||
|
||
#crystal geometry
|
||
#change it as you want
|
||
/crystal/setCrystalSize 20. 20. 0.05 mm
|
||
/crystal/setCrystalAngleX 0. mrad #for coherent bremsstrahlung change the angle (usually few mrad would work fine)
|
||
/crystal/setCrystalAngleY 0. mrad #default is 0
|
||
|
||
#crystal lattice parameters
|
||
#change it according to the data available in G4CHANNELINGDATA
|
||
/crystal/setCrystalMaterial G4_Diamond # check G4CHANNELINGDATA dataset;
|
||
# G4_Si and G4_Ge are also available.
|
||
# For less monochromatic but more intense
|
||
# radiation use G4_W.
|
||
/crystal/setCrystalLattice (110) # check G4CHANNELINGDATA dataset;
|
||
# planes (110) and (100) as well as
|
||
# axes <100>, <110> and <111> are also available
|
||
# For G4_W check G4CHANNELINGDATA dataset
|
||
# for the availability.
|
||
|
||
#detector parameters
|
||
#change it as you want
|
||
/crystal/setDetectorSize 10 10 0.03 cm
|
||
/crystal/setFrontPositionZ 1000 mm
|
||
|
||
#G4ChannelingFastSimModel parameters
|
||
#change it as you want
|
||
/crystal/setChannelingModel true #switch on/off G4ChannelingFastSimModel
|
||
/crystal/setRadiationModel true #switch on/off G4BaierKatkov
|
||
#(G4ChannelingFastSimModel must be switched on)
|
||
|
||
#the following can be changed or commented:
|
||
#the low energy threshold for particle to enter the G4ChannelingFastSimModel:
|
||
/crystal/setParticleMinKinEnergy/e- 200 MeV # for energies of primaries below 200 MeV
|
||
# you may decrease the threshold;
|
||
# be careful with the energies
|
||
# below 100 MeV where quantum effects
|
||
# particle dynamics may become important.
|
||
# More information -
|
||
# Baier, Katkov, Strakhovenko book (above)
|
||
/crystal/setParticleMinKinEnergy/e+ 200 MeV # the same as for e-.
|
||
|
||
#the following can be changed or commented:
|
||
#high angular threshold for particle to enter the G4ChannelingFastSimModel expressed in Lindhard angles:
|
||
/crystal/setLindhardAngles/e+ 100. #change if necessary to speed up the simulations
|
||
/crystal/setLindhardAngles/e- 100. #change if necessary to speed up the simulations
|
||
|
||
#the following can be changed or commented:
|
||
#G4BaierKatkov parameters
|
||
#(see also comments in ConstructSDandField in DetectorConstruction):
|
||
/crystal/setMinPhotonEnergy 1 MeV # G4BaierKatkov default
|
||
/crystal/setSamplingPhotonsNumber 150 # G4BaierKatkov default
|
||
/crystal/setNSmallTrajectorySteps 10000 # G4BaierKatkov default
|
||
/crystal/setRadiationAngleFactor 4. # G4BaierKatkov default
|
||
|
||
#the following can be changed or commented:
|
||
#For scoring purposes only: the maximal energy in the output spectrum from the Baier-Katkov method.
|
||
#Note it does not modifies the simulations.
|
||
#Note: the minimal energy in the spectrum = the minimal energy defined above.
|
||
/crystal/MaxBKPhotonEnergyInSpectrum 0.5 GeV # must be > setMinPhotonEnergy
|
||
/crystal/NBinsInSpectrum 50 # note these bins are not equidistant.
|
||
#Virtual collimator sets the selection of photon angles to be written in the output spectrum
|
||
/crystal/setVirtualCollimatorHalfSize 1. mrad #setting the angular radius of virtual collimator
|
||
|
||
/run/initialize
|
||
|
||
#settings
|
||
/control/verbose 0
|
||
/run/verbose 0
|
||
/tracking/verbose 0
|
||
|
||
#beam
|
||
/gps/particle e- #e+ also work fine
|
||
|
||
#coordinate distribution (radial Gauss)
|
||
/gps/pos/centre 0. 0. -1. cm
|
||
/gps/pos/type Beam
|
||
/gps/pos/sigma_x 0.1 mm
|
||
/gps/pos/sigma_y 0.1 mm
|
||
|
||
#angular distribution (radial Gauss)
|
||
/gps/ang/type beam2d
|
||
/gps/ang/rot1 1 0 0
|
||
/gps/ang/rot2 0 -1 0
|
||
/gps/ang/sigma_x 50.E-6 rad
|
||
/gps/ang/sigma_y 50.E-6 rad
|
||
|
||
#energy distribution (constant)
|
||
/gps/ene/mono 0.5 GeV
|
||
|
||
#statistics
|
||
/run/printProgress 1000
|
||
/run/beamOn 10000
|