1 ! constants used
for pair/brem in Air
4 ! we obtain probability /radation length.
6 real*8 X0g ! radiation length. in
g/cm^2
8 real*8
cScrC1 !
const which appears in the complete screening
14 real*8 CompScrE ! Energy above which we can use complete screening
15 ! cross-sections. evaluate at Eg/Ee=
x= 0.99.
17 ! This can be applied
to pair creation, too.
18 ! 10.0*( 0.05*0.99/0.01/media.Zeff**0.333)
27 !
for brems. 5*rho_Pb/rho GeV.
38 !
for brems in the region A.: 20
48 ! ---------------------------------------------
53 ! B.H original xsection
is used
as xsec=
55 ! Norm * B.H(10MeV) = Pair(10MeV)
56 real*8 PairNonSc ! see above.
59 real*8
PrScrE ! below
this, screened cross-section
is used.
70 !
for pair in the region A.: 20
83 !
for Seltzer cross-section
! constants used for pair brem in Air real *mbtoPX0 ! mbtoPgrm x X0g If multiplied to sigma in ! we obtain probability radation length real *X0g ! radiation length in g cm real *cScrC1 ! const which appears in the complete screening ! crossection real *cScrC2 ! the other such one real too **media Zeff **Z **real *BrScrE ! below screened cross section is used real *BremEgmin ! Brems gamma min energy in Gev(100 keV) real *8 BremEemin ! Brems by e+/e- is considered at E > BremEemin !(Ek ~ 300 keV) real *8 BremLEemin ! log10 of BremEemin real *8 BremEemaxL ! Energy of e+/e- below which LPM is negelected ! for brems. 5 *rho_Pb/rho GeV. integer BremTXTL ! Size of the Brems total x-section table. ! in the energy region BremEemin ~ BremEemaxL B PairTXTL real *PairUminLA !min of uniform random number in the region A ! at energies PairEgmin * BremUmaxLB
! constants used for pair brem in Air real *mbtoPX0 ! mbtoPgrm x X0g If multiplied to sigma in ! we obtain probability radation length real *X0g ! radiation length in g cm real *cScrC1 ! const which appears in the complete screening ! crossection real *cScrC2 ! the other such one real too **media Zeff **Z **real *BrScrE ! below screened cross section is used real *BremEgmin ! Brems gamma min energy in Gev(100 keV) real *8 BremEemin ! Brems by e+/e- is considered at E > BremEemin !(Ek ~ 300 keV) real *8 BremLEemin ! log10 of BremEemin real *8 BremEemaxL ! Energy of e+/e- below which LPM is negelected ! for brems. 5 *rho_Pb/rho GeV. integer BremTXTL ! Size of the Brems total x-section table. ! in the energy region BremEemin ~ BremEemaxL B PairTXTL real *PairUminLA !min of uniform random number in the region A ! at energies PairEgmin BrdES
! constants used for pair brem in Air real *mbtoPX0 ! mbtoPgrm x X0g If multiplied to sigma in ! we obtain probability radation length real *X0g ! radiation length in g cm real *cScrC1 ! const which appears in the complete screening ! crossection real *cScrC2 ! the other such one real too **media Zeff **Z **real *BrScrE ! below screened cross section is used real *BremEgmin ! Brems gamma min energy in Gev(100 keV) real *8 BremEemin ! Brems by e+/e- is considered at E > BremEemin !(Ek ~ 300 keV) real *8 BremLEemin ! log10 of BremEemin real *8 BremEemaxL ! Energy of e+/e- below which LPM is negelected ! for brems. 5 *rho_Pb/rho GeV. integer BremTXTL ! Size of the Brems total x-section table. ! in the energy region BremEemin ~ BremEemaxL B PairTXTL real *PairUminLA !min of uniform random number in the region A ! at energies PairEgmin BrEeminS
! constants used for pair brem in Air real *mbtoPX0 ! mbtoPgrm x X0g If multiplied to sigma in ! we obtain probability radation length real *X0g ! radiation length in g cm real *cScrC1 ! const which appears in the complete screening ! crossection real *cScrC2 ! the other such one real too **media Zeff **Z **real *BrScrE ! below screened cross section is used real *BremEgmin ! Brems gamma min energy in Gev(100 keV) real *8 BremEemin ! Brems by e+/e- is considered at E > BremEemin !(Ek ~ 300 keV) real *8 BremLEemin ! log10 of BremEemin real *8 BremEemaxL ! Energy of e+/e- below which LPM is negelected ! for brems. 5 *rho_Pb/rho GeV. integer BremTXTL ! Size of the Brems total x-section table. ! in the energy region BremEemin ~ BremEemaxL B PairTXTL real *PairUminLA !min of uniform random number in the region A ! at energies PairEgmin X0g muPrdE common Zicnst BrUszSA
! constants used for pair brem in Air real *mbtoPX0 ! mbtoPgrm x X0g If multiplied to sigma in ! we obtain probability radation length real *X0g ! radiation length in g cm real *cScrC1 ! const which appears in the complete screening ! crossection real *cScrC2 ! the other such one real too **media Zeff **Z **real *BrScrE ! below screened cross section is used real *BremEgmin ! Brems gamma min energy in Gev(100 keV) real *8 BremEemin ! Brems by e+/e- is considered at E > BremEemin !(Ek ~ 300 keV) real *8 BremLEemin ! log10 of BremEemin real *8 BremEemaxL ! Energy of e+/e- below which LPM is negelected ! for brems. 5 *rho_Pb/rho GeV. integer BremTXTL ! Size of the Brems total x-section table. ! in the energy region BremEemin ~ BremEemaxL B PairTXTL real *PairUminLA !min of uniform random number in the region A ! at energies PairEgmin X0g muNEmax
! constants used for pair brem in Air real *mbtoPX0 ! mbtoPgrm x X0g If multiplied to sigma in ! we obtain probability radation length real *X0g ! radiation length in g cm real *cScrC1 ! const which appears in the complete screening ! crossection real *cScrC2 ! the other such one real too **media Zeff **Z **real *BrScrE ! below screened cross section is used real *BremEgmin ! Brems gamma min energy in Gev(100 keV) real *8 BremEemin ! Brems by e+/e- is considered at E > BremEemin !(Ek ~ 300 keV) real *8 BremLEemin ! log10 of BremEemin real *8 BremEemaxL ! Energy of e+/e- below which LPM is negelected ! for brems. 5 *rho_Pb/rho GeV. integer BremTXTL ! Size of the Brems total x-section table. ! in the energy region BremEemin ~ BremEemaxL B PairTXTL real *PairUminLA !min of uniform random number in the region A ! at energies PairEgmin X0g muNpwdEdxt
! constants used for pair brem in Air real *mbtoPX0 ! mbtoPgrm x X0g If multiplied to sigma in ! we obtain probability radation length real *X0g ! radiation length in g cm real *cScrC1 ! const which appears in the complete screening ! crossection real *cScrC2 ! the other such one real too **media Zeff **Z **real *BrScrE ! below screened cross section is used real *BremEgmin ! Brems gamma min energy in Gev(100 keV) real *8 BremEemin ! Brems by e+/e- is considered at E > BremEemin !(Ek ~ 300 keV) real *8 BremLEemin ! log10 of BremEemin real *8 BremEemaxL ! Energy of e+/e- below which LPM is negelected ! for brems. 5 *rho_Pb/rho GeV. integer BremTXTL ! Size of the Brems total x-section table. ! in the energy region BremEemin ~ BremEemaxL B PairTXTL real *PairUminLA !min of uniform random number in the region A ! at energies PairEgmin X0g muBrdE
! constants used for pair brem in Air real *mbtoPX0 ! mbtoPgrm x X0g If multiplied to sigma in ! we obtain probability radation length real *X0g ! radiation length in g cm real *cScrC1 ! const which appears in the complete screening ! crossection real *cScrC2 ! the other such one real too **media Zeff **Z **real *BrScrE ! below screened cross section is used real *BremEgmin ! Brems gamma min energy in Gev(100 keV) real *8 BremEemin ! Brems by e+/e- is considered at E > BremEemin !(Ek ~ 300 keV) real *8 BremLEemin ! log10 of BremEemin real *8 BremEemaxL ! Energy of e+/e- below which LPM is negelected ! for brems. 5 *rho_Pb/rho GeV. integer BremTXTL ! Size of the Brems total x-section table. ! in the energy region BremEemin ~ BremEemaxL B PairTXTL real *PairUminLA !min of uniform random number in the region A ! at energies PairEgmin BremdEL
! constants used for pair brem in Air real *mbtoPX0 ! mbtoPgrm x X0g If multiplied to sigma in ! we obtain probability radation length real *X0g ! radiation length in g cm real *cScrC1 ! const which appears in the complete screening ! crossection real *cScrC2 ! the other such one real too **media Zeff **Z **real *BrScrE ! below screened cross section is used real *BremEgmin ! Brems gamma min energy in Gev(100 keV) real *8 BremEemin ! Brems by e+/e- is considered at E > BremEemin !(Ek ~ 300 keV) real *8 BremLEemin ! log10 of BremEemin real *8 BremEemaxL ! Energy of e+/e- below which LPM is negelected ! for brems. 5 *rho_Pb/rho GeV. integer BremTXTL ! Size of the Brems total x-section table. ! in the energy region BremEemin ~ BremEemaxL B PairTXTL real *PairUminLA !min of uniform random number in the region A ! at energies PairEgmin X0g muPrdE common Zicnst PairEsize
! constants used for pair brem in Air real *mbtoPX0 ! mbtoPgrm x X0g If multiplied to sigma in ! we obtain probability radation length real *X0g ! radiation length in g cm real *cScrC1 ! const which appears in the complete screening ! crossection real *cScrC2 ! the other such one real too **media Zeff **Z **real *BrScrE ! below screened cross section is used real *BremEgmin ! Brems gamma min energy in Gev(100 keV) real *8 BremEemin ! Brems by e+/e- is considered at E > BremEemin !(Ek ~ 300 keV) real *8 BremLEemin ! log10 of BremEemin real *8 BremEemaxL ! Energy of e+/e- below which LPM is negelected ! for brems. 5 *rho_Pb/rho GeV. integer BremTXTL ! Size of the Brems total x-section table. ! in the energy region BremEemin ~ BremEemaxL B PairTXTL real *PairUminLA !min of uniform random number in the region A ! at energies PairEgmin X0g muPrdE common Zicnst PairUszLB
! constants used for pair brem in Air real *mbtoPX0 ! mbtoPgrm x X0g If multiplied to sigma in ! we obtain probability radation length real *X0g ! radiation length in g cm real *cScrC1 ! const which appears in the complete screening ! crossection real *cScrC2 ! the other such one real too **media Zeff **Z **real *BrScrE ! below screened cross section is used real *BremEgmin ! Brems gamma min energy in Gev(100 keV) real *8 BremEemin ! Brems by e+/e- is considered at E > BremEemin !(Ek ~ 300 keV) real *8 BremLEemin ! log10 of BremEemin real *8 BremEemaxL ! Energy of e+/e- below which LPM is negelected ! for brems. 5 *rho_Pb/rho GeV. integer BremTXTL ! Size of the Brems total x-section table. ! in the energy region BremEemin ~ BremEemaxL B PairTXTL real *PairUminLA !min of uniform random number in the region A ! at energies PairEgmin PairNonSc
! constants used for pair brem in Air real *mbtoPX0 ! mbtoPgrm x X0g If multiplied to sigma in ! we obtain probability radation length real *X0g ! radiation length in g cm real *cScrC1 ! const which appears in the complete screening ! crossection real *cScrC2 ! the other such one real too **media Zeff **Z **real *BrScrE ! below screened cross section is used real *BremEgmin ! Brems gamma min energy in Gev(100 keV) real *8 BremEemin ! Brems by e+/e- is considered at E > BremEemin !(Ek ~ 300 keV) real *8 BremLEemin ! log10 of BremEemin real *8 BremEemaxL ! Energy of e+/e- below which LPM is negelected ! for brems. 5 *rho_Pb/rho GeV. integer BremTXTL ! Size of the Brems total x-section table. ! in the energy region BremEemin ~ BremEemaxL B PairTXTL real *PairUminLA !min of uniform random number in the region A ! at energies PairEgmin X0g muBrEmin
! constants used for pair brem in Air real *mbtoPX0 ! mbtoPgrm x X0g If multiplied to sigma in ! we obtain probability radation length real *X0g ! radiation length in g cm real *cScrC1 ! const which appears in the complete screening ! crossection real *cScrC2 ! the other such one real too **media Zeff **Z **real *BrScrE ! below screened cross section is used real *BremEgmin ! Brems gamma min energy in Gev(100 keV) real *8 BremEemin ! Brems by e+/e- is considered at E > BremEemin !(Ek ~ 300 keV) real *8 BremLEemin ! log10 of BremEemin real *8 BremEemaxL ! Energy of e+/e- below which LPM is negelected ! for brems. 5 *rho_Pb/rho GeV. integer BremTXTL ! Size of the Brems total x-section table. ! in the energy region BremEemin ~ BremEemaxL B PairTXTL real *PairUminLA !min of uniform random number in the region A ! at energies PairEgmin BrUminSB
! constants used for pair brem in Air real *mbtoPX0 ! mbtoPgrm x X0g If multiplied to sigma in ! we obtain probability radation length real *X0g ! radiation length in g cm real *cScrC1 ! const which appears in the complete screening ! crossection real *cScrC2 ! the other such one real too **media Zeff **Z **real *BrScrE ! below screened cross section is used real *BremEgmin ! Brems gamma min energy in Gev(100 keV) real *8 BremEemin ! Brems by e+/e- is considered at E > BremEemin !(Ek ~ 300 keV) real *8 BremLEemin ! log10 of BremEemin real *8 BremEemaxL ! Energy of e+/e- below which LPM is negelected ! for brems. 5 *rho_Pb/rho GeV. integer BremTXTL ! Size of the Brems total x-section table. ! in the energy region BremEemin ~ BremEemaxL B PairTXTL real *PairUminLA !min of uniform random number in the region A ! at energies PairEgmin BremdULA
dE dx *! Nuc Int sampling table e
! constants used for pair brem in Air real *mbtoPX0 ! mbtoPgrm x X0g If multiplied to sigma in ! we obtain probability radation length real *X0g ! radiation length in g cm real *cScrC1 ! const which appears in the complete screening ! crossection real *cScrC2 ! the other such one real too **media Zeff **Z **real *BrScrE ! below screened cross section is used real *BremEgmin ! Brems gamma min energy in Gev(100 keV) real *8 BremEemin ! Brems by e+/e- is considered at E > BremEemin !(Ek ~ 300 keV) real *8 BremLEemin ! log10 of BremEemin real *8 BremEemaxL ! Energy of e+/e- below which LPM is negelected ! for brems. 5 *rho_Pb/rho GeV. integer BremTXTL ! Size of the Brems total x-section table. ! in the energy region BremEemin ~ BremEemaxL B PairTXTL real *PairUminLA !min of uniform random number in the region A ! at energies PairEgmin PairEgmin
! constants used for pair brem in Air real *mbtoPX0 ! mbtoPgrm x X0g If multiplied to sigma in ! we obtain probability radation length real *X0g ! radiation length in g cm real *cScrC1 ! const which appears in the complete screening ! crossection real *cScrC2 ! the other such one real too **media Zeff **Z **real *BrScrE ! below screened cross section is used real *BremEgmin ! Brems gamma min energy in Gev(100 keV) real *8 BremEemin ! Brems by e+/e- is considered at E > BremEemin !(Ek ~ 300 keV) real *8 BremLEemin ! log10 of BremEemin real *8 BremEemaxL ! Energy of e+/e- below which LPM is negelected ! for brems. 5 *rho_Pb/rho GeV. integer BremTXTL ! Size of the Brems total x-section table. ! in the energy region BremEemin ~ BremEemaxL B PairTXTL real *PairUminLA !min of uniform random number in the region A ! at energies PairEgmin * BremUminLA
! constants used for pair brem in Air real *mbtoPX0 ! mbtoPgrm x X0g If multiplied to sigma in ! we obtain probability radation length real *X0g ! radiation length in g cm real *cScrC1 ! const which appears in the complete screening ! crossection real *cScrC2 ! the other such one real too **media Zeff **Z **real *BrScrE ! below screened cross section is used real *BremEgmin ! Brems gamma min energy in Gev(100 keV) real *8 BremEemin ! Brems by e+/e- is considered at E > BremEemin !(Ek ~ 300 keV) real *8 BremLEemin ! log10 of BremEemin real *8 BremEemaxL ! Energy of e+/e- below which LPM is negelected ! for brems. 5 *rho_Pb/rho GeV. integer BremTXTL ! Size of the Brems total x-section table. ! in the energy region BremEemin ~ BremEemaxL B PairTXTL real *PairUminLA !min of uniform random number in the region A ! at energies PairEgmin X0g muNVmin
! constants used for pair brem in Air real *mbtoPX0 ! mbtoPgrm x X0g If multiplied to sigma in ! we obtain probability radation length real *X0g ! radiation length in g cm real *cScrC1 ! const which appears in the complete screening ! crossection real *cScrC2 ! the other such one real too **media Zeff **Z **real *BrScrE ! below screened cross section is used real *BremEgmin ! Brems gamma min energy in Gev(100 keV) real *8 BremEemin ! Brems by e+/e- is considered at E > BremEemin !(Ek ~ 300 keV) real *8 BremLEemin ! log10 of BremEemin real *8 BremEemaxL ! Energy of e+/e- below which LPM is negelected ! for brems. 5 *rho_Pb/rho GeV. integer BremTXTL ! Size of the Brems total x-section table. ! in the energy region BremEemin ~ BremEemaxL B PairTXTL real *PairUminLA !min of uniform random number in the region A ! at energies PairEgmin X0g muBrLEmin
! constants used for pair brem in Air real *mbtoPX0 ! mbtoPgrm x X0g If multiplied to sigma in ! we obtain probability radation length real *X0g ! radiation length in g cm real *cScrC1 ! const which appears in the complete screening ! crossection real *cScrC2 ! the other such one real too **media Zeff **Z **real *BrScrE ! below screened cross section is used real *BremEgmin ! Brems gamma min energy in Gev(100 keV) real *8 BremEemin ! Brems by e+/e- is considered at E > BremEemin !(Ek ~ 300 keV) real *8 BremLEemin ! log10 of BremEemin real *8 BremEemaxL ! Energy of e+/e- below which LPM is negelected ! for brems. 5 *rho_Pb/rho GeV. integer BremTXTL ! Size of the Brems total x-section table. ! in the energy region BremEemin ~ BremEemaxL B PairTXTL real *PairUminLA !min of uniform random number in the region A ! at energies PairEgmin ~ PairEgmaxL
! constants used for pair brem in Air real *mbtoPX0 ! mbtoPgrm x X0g If multiplied to sigma in ! we obtain probability radation length real *X0g ! radiation length in g cm real *cScrC1 ! const which appears in the complete screening ! crossection real *cScrC2 ! the other such one real too **media Zeff **Z **real *BrScrE ! below screened cross section is used real *BremEgmin ! Brems gamma min energy in Gev(100 keV) real *8 BremEemin ! Brems by e+/e- is considered at E > BremEemin !(Ek ~ 300 keV) real *8 BremLEemin ! log10 of BremEemin real *8 BremEemaxL ! Energy of e+/e- below which LPM is negelected ! for brems. 5 *rho_Pb/rho GeV. integer BremTXTL ! Size of the Brems total x-section table. ! in the energy region BremEemin ~ BremEemaxL B PairTXTL real *PairUminLA !min of uniform random number in the region A ! at energies PairEgmin * PairUminLB
! constants used for pair brem in Air real *mbtoPX0 ! mbtoPgrm x X0g If multiplied to sigma in ! we obtain probability radation length real *X0g ! radiation length in g cm real *cScrC1 ! const which appears in the complete screening ! crossection real *cScrC2 ! the other such one real too **media Zeff **Z **real *BrScrE ! below screened cross section is used real *BremEgmin ! Brems gamma min energy in Gev(100 keV) real *8 BremEemin ! Brems by e+/e- is considered at E > BremEemin !(Ek ~ 300 keV) real *8 BremLEemin ! log10 of BremEemin real *8 BremEemaxL ! Energy of e+/e- below which LPM is negelected ! for brems. 5 *rho_Pb/rho GeV. integer BremTXTL ! Size of the Brems total x-section table. ! in the energy region BremEemin ~ BremEemaxL B PairTXTL real *PairUminLA !min of uniform random number in the region A ! at energies PairEgmin BrEemaxS
! constants used for pair brem in Air real *mbtoPX0 ! mbtoPgrm x X0g If multiplied to sigma in ! we obtain probability radation length real *X0g ! radiation length in g cm real *cScrC1 ! const which appears in the complete screening ! crossection real *cScrC2 ! the other such one real too **media Zeff **Z **real *BrScrE ! below screened cross section is used real *BremEgmin ! Brems gamma min energy in Gev(100 keV) real *8 BremEemin ! Brems by e+/e- is considered at E > BremEemin !(Ek ~ 300 keV) real *8 BremLEemin ! log10 of BremEemin real *8 BremEemaxL ! Energy of e+/e- below which LPM is negelected ! for brems. 5 *rho_Pb/rho GeV. integer BremTXTL ! Size of the Brems total x-section table. ! in the energy region BremEemin ~ BremEemaxL B PairTXTL real *PairUminLA !min of uniform random number in the region A ! at energies PairEgmin X0g * muBrVmin
! constants used for pair brem in Air real *mbtoPX0 ! mbtoPgrm x X0g If multiplied to sigma in ! we obtain probability radation length real *X0g ! radiation length in g cm real *cScrC1 ! const which appears in the complete screening ! crossection real *cScrC2 ! the other such one real too **media Zeff **Z **real *BrScrE ! below screened cross section is used real *BremEgmin ! Brems gamma min energy in Gev(100 keV) real *8 BremEemin ! Brems by e+/e- is considered at E > BremEemin !(Ek ~ 300 keV) real *8 BremLEemin ! log10 of BremEemin real *8 BremEemaxL ! Energy of e+/e- below which LPM is negelected ! for brems. 5 *rho_Pb/rho GeV. integer BremTXTL ! Size of the Brems total x-section table. ! in the energy region BremEemin ~ BremEemaxL B PairTXTL real *PairUminLA !min of uniform random number in the region A ! at energies PairEgmin * cScrC1
! constants used for pair brem in Air real *mbtoPX0 ! mbtoPgrm x X0g If multiplied to sigma in ! we obtain probability radation length real *X0g ! radiation length in g cm real *cScrC1 ! const which appears in the complete screening ! crossection real *cScrC2 ! the other such one real too **media Zeff **Z **real *BrScrE ! below screened cross section is used real *BremEgmin ! Brems gamma min energy in Gev(100 keV) real *8 BremEemin ! Brems by e+/e- is considered at E > BremEemin !(Ek ~ 300 keV) real *8 BremLEemin ! log10 of BremEemin real *8 BremEemaxL ! Energy of e+/e- below which LPM is negelected ! for brems. 5 *rho_Pb/rho GeV. integer BremTXTL ! Size of the Brems total x-section table. ! in the energy region BremEemin ~ BremEemaxL B PairTXTL real *PairUminLA !min of uniform random number in the region A ! at energies PairEgmin X0g muBrEmax1
! constants used for pair brem in Air real *mbtoPX0 ! mbtoPgrm x X0g If multiplied to sigma in ! we obtain probability radation length real *X0g ! radiation length in g cm real *cScrC1 ! const which appears in the complete screening ! crossection real *cScrC2 ! the other such one real too **media Zeff **Z **real *BrScrE ! below screened cross section is used real *BremEgmin ! Brems gamma min energy in Gev(100 keV) real *8 BremEemin ! Brems by e+/e- is considered at E > BremEemin !(Ek ~ 300 keV) real *8 BremLEemin ! log10 of BremEemin real *8 BremEemaxL ! Energy of e+/e- below which LPM is negelected ! for brems. 5 *rho_Pb/rho GeV. integer BremTXTL ! Size of the Brems total x-section table. ! in the energy region BremEemin ~ BremEemaxL B PairTXTL real *PairUminLA !min of uniform random number in the region A ! at energies PairEgmin X0g muPrEmin
! constants used for pair brem in Air real *mbtoPX0 ! mbtoPgrm x X0g If multiplied to sigma in ! we obtain probability radation length real *X0g ! radiation length in g cm real *cScrC1 ! const which appears in the complete screening ! crossection real *cScrC2 ! the other such one real too **media Zeff **Z **real *BrScrE ! below screened cross section is used real *BremEgmin ! Brems gamma min energy in Gev(100 keV) real *8 BremEemin ! Brems by e+/e- is considered at E > BremEemin !(Ek ~ 300 keV) real *8 BremLEemin ! log10 of BremEemin real *8 BremEemaxL ! Energy of e+/e- below which LPM is negelected ! for brems. 5 *rho_Pb/rho GeV. integer BremTXTL ! Size of the Brems total x-section table. ! in the energy region BremEemin ~ BremEemaxL B PairTXTL real *PairUminLA !min of uniform random number in the region A ! at energies PairEgmin BrLEeminS
! constants used for pair brem in Air real *mbtoPX0 ! mbtoPgrm x X0g If multiplied to sigma in ! we obtain probability radation length real *X0g ! radiation length in g cm real *cScrC1 ! const which appears in the complete screening ! crossection real *cScrC2 ! the other such one real too **media Zeff **Z **real *BrScrE ! below screened cross section is used real *BremEgmin ! Brems gamma min energy in Gev(100 keV) real *8 BremEemin ! Brems by e+/e- is considered at E > BremEemin !(Ek ~ 300 keV) real *8 BremLEemin ! log10 of BremEemin real *8 BremEemaxL ! Energy of e+/e- below which LPM is negelected ! for brems. 5 *rho_Pb/rho GeV. integer BremTXTL ! Size of the Brems total x-section table. ! in the energy region BremEemin ~ BremEemaxL B PairTXTL real *PairUminLA !min of uniform random number in the region A ! at energies PairEgmin X0g muPrdE common Zicnst BremEsize
! constants used for pair brem in Air real *mbtoPX0 ! mbtoPgrm x X0g If multiplied to sigma in ! we obtain probability radation length real *X0g ! radiation length in g cm real *cScrC1 ! const which appears in the complete screening ! crossection real *cScrC2 ! the other such one real too **media Zeff **Z **real *BrScrE ! below screened cross section is used real *BremEgmin ! Brems gamma min energy in Gev(100 keV) real *8 BremEemin ! Brems by e+/e- is considered at E > BremEemin !(Ek ~ 300 keV) real *8 BremLEemin ! log10 of BremEemin real *8 BremEemaxL ! Energy of e+/e- below which LPM is negelected ! for brems. 5 *rho_Pb/rho GeV. integer BremTXTL ! Size of the Brems total x-section table. ! in the energy region BremEemin ~ BremEemaxL B PairTXTL real *PairUminLA !min of uniform random number in the region A ! at energies PairEgmin PairdELA
! constants used for pair brem in Air real *mbtoPX0 ! mbtoPgrm x X0g If multiplied to sigma in ! we obtain probability radation length real *X0g ! radiation length in g cm real *cScrC1 ! const which appears in the complete screening ! crossection real *cScrC2 ! the other such one real too **media Zeff **Z **real *BrScrE ! below screened cross section is used real *BremEgmin ! Brems gamma min energy in Gev(100 keV) real *8 BremEemin ! Brems by e+/e- is considered at E > BremEemin !(Ek ~ 300 keV) real *8 BremLEemin ! log10 of BremEemin real *8 BremEemaxL ! Energy of e+/e- below which LPM is negelected ! for brems. 5 *rho_Pb/rho GeV. integer BremTXTL ! Size of the Brems total x-section table. ! in the energy region BremEemin ~ BremEemaxL B PairTXTL real *PairUminLA !min of uniform random number in the region A ! at energies PairEgmin BremEemaxL
! constants used for pair brem in Air real *mbtoPX0 ! mbtoPgrm x X0g If multiplied to sigma in ! we obtain probability radation length real *X0g ! radiation length in g cm real *cScrC1 ! const which appears in the complete screening ! crossection real *cScrC2 ! the other such one real too **media Zeff **Z **real *BrScrE ! below screened cross section is used real *BremEgmin ! Brems gamma min energy in Gev(100 keV) real *8 BremEemin ! Brems by e+/e- is considered at E > BremEemin !(Ek ~ 300 keV) real *8 BremLEemin ! log10 of BremEemin real *8 BremEemaxL ! Energy of e+/e- below which LPM is negelected ! for brems. 5 *rho_Pb/rho GeV. integer BremTXTL ! Size of the Brems total x-section table. ! in the energy region BremEemin ~ BremEemaxL B PairTXTL real *PairUminLA !min of uniform random number in the region A ! at energies PairEgmin X0g muPrdE common Zicnst PairTXTL
! constants used for pair brem in Air real *mbtoPX0 ! mbtoPgrm x X0g If multiplied to sigma in ! we obtain probability radation length real *X0g ! radiation length in g cm real *cScrC1 ! const which appears in the complete screening ! crossection real *cScrC2 ! the other such one real too **media Zeff **Z **real *BrScrE ! below screened cross section is used real *BremEgmin ! Brems gamma min energy in Gev(100 keV) real *8 BremEemin ! Brems by e+/e- is considered at E > BremEemin !(Ek ~ 300 keV) real *8 BremLEemin ! log10 of BremEemin real *8 BremEemaxL ! Energy of e+/e- below which LPM is negelected ! for brems. 5 *rho_Pb/rho GeV. integer BremTXTL ! Size of the Brems total x-section table. ! in the energy region BremEemin ~ BremEemaxL B PairTXTL real *PairUminLA !min of uniform random number in the region A ! at energies PairEgmin X0g muPrdE common Zicnst BrUszSB muNUsize muBrTXT
! constants used for pair brem in Air real *mbtoPX0 ! mbtoPgrm x X0g If multiplied to sigma in ! we obtain probability radation length real *X0g ! radiation length in g cm real *cScrC1 ! const which appears in the complete screening ! crossection real *cScrC2 ! the other such one real too **media Zeff **Z **real *BrScrE ! below screened cross section is used real *BremEgmin ! Brems gamma min energy in Gev(100 keV) real *8 BremEemin ! Brems by e+/e- is considered at E > BremEemin !(Ek ~ 300 keV) real *8 BremLEemin ! log10 of BremEemin real *8 BremEemaxL ! Energy of e+/e- below which LPM is negelected ! for brems. 5 *rho_Pb/rho GeV. integer BremTXTL ! Size of the Brems total x-section table. ! in the energy region BremEemin ~ BremEemaxL B PairTXTL real *PairUminLA !min of uniform random number in the region A ! at energies PairEgmin X0g muBrEmax
! constants used for pair brem in Air real *mbtoPX0 ! mbtoPgrm x X0g If multiplied to sigma in ! we obtain probability radation length real *X0g ! radiation length in g cm real *cScrC1 ! const which appears in the complete screening ! crossection real *cScrC2 ! the other such one real too **media Zeff **Z **real *BrScrE ! below screened cross section is used real *BremEgmin ! Brems gamma min energy in Gev(100 keV) real *8 BremEemin ! Brems by e+/e- is considered at E > BremEemin !(Ek ~ 300 keV) real *8 BremLEemin ! log10 of BremEemin real *8 BremEemaxL ! Energy of e+/e- below which LPM is negelected ! for brems. 5 *rho_Pb/rho GeV. integer BremTXTL ! Size of the Brems total x-section table. ! in the energy region BremEemin ~ BremEemaxL B PairTXTL real *PairUminLA !min of uniform random number in the region A ! at energies PairEgmin X0g muPrdE common Zicnst BrUszSB muNTXT
! constants used for pair brem in Air real *mbtoPX0 ! mbtoPgrm x X0g If multiplied to sigma in ! we obtain probability radation length real *X0g ! radiation length in g cm real *cScrC1 ! const which appears in the complete screening ! crossection real *cScrC2 ! the other such one real too **media Zeff **Z **real *BrScrE ! below screened cross section is used real *BremEgmin ! Brems gamma min energy in Gev(100 keV) real *8 BremEemin ! Brems by e+/e- is considered at E > BremEemin !(Ek ~ 300 keV) real *8 BremLEemin ! log10 of BremEemin real *8 BremEemaxL ! Energy of e+/e- below which LPM is negelected ! for brems. 5 *rho_Pb/rho GeV. integer BremTXTL ! Size of the Brems total x-section table. ! in the energy region BremEemin ~ BremEemaxL B PairTXTL real *PairUminLA !min of uniform random number in the region A ! at energies PairEgmin * BrEgminS
! constants used for pair brem in Air real *mbtoPX0 ! mbtoPgrm x X0g If multiplied to sigma in ! we obtain probability radation length real *X0g ! radiation length in g cm real *cScrC1 ! const which appears in the complete screening ! crossection real *cScrC2 ! the other such one real too **media Zeff **Z **real *BrScrE ! below screened cross section is used real *BremEgmin ! Brems gamma min energy in Gev(100 keV) real *8 BremEemin ! Brems by e+/e- is considered at E > BremEemin !(Ek ~ 300 keV) real *8 BremLEemin ! log10 of BremEemin real *8 BremEemaxL ! Energy of e+/e- below which LPM is negelected ! for brems. 5 *rho_Pb/rho GeV. integer BremTXTL ! Size of the Brems total x-section table. ! in the energy region BremEemin ~ BremEemaxL B PairTXTL real *PairUminLA !min of uniform random number in the region A ! at energies PairEgmin * PairEgmaxL
! constants used for pair brem in Air real *mbtoPX0 ! mbtoPgrm x X0g If multiplied to sigma in mb
! constants used for pair brem in Air real *mbtoPX0 ! mbtoPgrm x X0g If multiplied to sigma in ! we obtain probability radation length real *X0g ! radiation length in g cm real *cScrC1 ! const which appears in the complete screening ! crossection real *cScrC2 ! the other such one real too **media Zeff **Z **real *BrScrE ! below screened cross section is used real *BremEgmin ! Brems gamma min energy in Gev(100 keV) real *8 BremEemin ! Brems by e+/e- is considered at E > BremEemin !(Ek ~ 300 keV) real *8 BremLEemin ! log10 of BremEemin real *8 BremEemaxL ! Energy of e+/e- below which LPM is negelected ! for brems. 5 *rho_Pb/rho GeV. integer BremTXTL ! Size of the Brems total x-section table. ! in the energy region BremEemin ~ BremEemaxL B PairTXTL real *PairUminLA !min of uniform random number in the region A ! at energies PairEgmin PairdELB
! constants used for pair brem in Air real *mbtoPX0 ! mbtoPgrm x X0g If multiplied to sigma in ! we obtain probability radation length real *X0g ! radiation length in g cm real *cScrC1 ! const which appears in the complete screening ! crossection real *cScrC2 ! the other such one real too **media Zeff **Z **real *BrScrE ! below screened cross section is used real *BremEgmin ! Brems gamma min energy in Gev(100 keV) real *8 BremEemin ! Brems by e+/e- is considered at E > BremEemin !(Ek ~ 300 keV) real *8 BremLEemin ! log10 of BremEemin real *8 BremEemaxL ! Energy of e+/e- below which LPM is negelected ! for brems. 5 *rho_Pb/rho GeV. integer BremTXTL ! Size of the Brems total x-section table. ! in the energy region BremEemin ~ BremEemaxL B PairTXTL real *PairUminLA !min of uniform random number in the region A ! at energies PairEgmin X0g muPrdE common Zicnst BremUszLA
! constants used for pair brem in Air real *mbtoPX0 ! mbtoPgrm x X0g If multiplied to sigma in ! we obtain probability radation length real *X0g ! radiation length in g cm real *cScrC1 ! const which appears in the complete screening ! crossection real *cScrC2 ! the other such one real too **media Zeff **Z **real *BrScrE ! below screened cross section is used real *BremEgmin ! Brems gamma min energy in Gev(100 keV) real *8 BremEemin ! Brems by e+/e- is considered at E > BremEemin !(Ek ~ 300 keV) real *8 BremLEemin ! log10 of BremEemin real *8 BremEemaxL ! Energy of e+/e- below which LPM is negelected ! for brems. 5 *rho_Pb/rho GeV. integer BremTXTL ! Size of the Brems total x-section table. ! in the energy region BremEemin ~ BremEemaxL B PairTXTL real *PairUminLA !min of uniform random number in the region A ! at energies PairEgmin BrUmaxSB
! constants used for pair brem in Air real *mbtoPX0 ! mbtoPgrm x X0g If multiplied to sigma in ! we obtain probability radation length real *X0g ! radiation length in g cm real *cScrC1 ! const which appears in the complete screening ! crossection real *cScrC2 ! the other such one real too **media Zeff **Z **real *BrScrE ! below screened cross section is used real *BremEgmin ! Brems gamma min energy in Gev(100 keV) real *8 BremEemin ! Brems by e+/e- is considered at E > BremEemin !(Ek ~ 300 keV) real *8 BremLEemin ! log10 of BremEemin real *8 BremEemaxL ! Energy of e+/e- below which LPM is negelected ! for brems. 5 *rho_Pb/rho GeV. integer BremTXTL ! Size of the Brems total x-section table. ! in the energy region BremEemin ~ BremEemaxL B PairTXTL real *PairUminLA !min of uniform random number in the region A ! at energies PairEgmin X0g muNpwdEdx0
! constants used for pair brem in Air real *mbtoPX0 ! mbtoPgrm x X0g If multiplied to sigma in ! we obtain probability radation length real *X0g ! radiation length in g cm real *cScrC1 ! const which appears in the complete screening ! crossection real *cScrC2 ! the other such one real too **media Zeff **Z **real *BrScrE ! below screened cross section is used real *BremEgmin ! Brems gamma min energy in Gev(100 keV) real *8 BremEemin ! Brems by e+/e- is considered at E > BremEemin !(Ek ~ 300 keV) real *8 BremLEemin ! log10 of BremEemin real *8 BremEemaxL ! Energy of e+/e- below which LPM is negelected ! for brems. 5 *rho_Pb/rho GeV. integer BremTXTL ! Size of the Brems total x-section table. ! in the energy region BremEemin ~ BremEemaxL B PairTXTL real *PairUminLA !min of uniform random number in the region A ! at energies PairEgmin PrScrE
! Parameters used for hadronic cascade shower is generated newline ! For you may give as as or em as(order/case/separator insensitive) is to generate EM-cascade and AS. \newline ! Generate
! constants used for pair brem in Air real *mbtoPX0 ! mbtoPgrm x X0g If multiplied to sigma in ! we obtain probability radation length real *X0g ! radiation length in g cm real *cScrC1 ! const which appears in the complete screening ! crossection real *cScrC2 ! the other such one real too **media Zeff **Z **real *BrScrE ! below screened cross section is used real *BremEgmin ! Brems gamma min energy in Gev(100 keV) real *8 BremEemin ! Brems by e+/e- is considered at E > BremEemin !(Ek ~ 300 keV) real *8 BremLEemin ! log10 of BremEemin real *8 BremEemaxL ! Energy of e+/e- below which LPM is negelected ! for brems. 5 *rho_Pb/rho GeV. integer BremTXTL ! Size of the Brems total x-section table. ! in the energy region BremEemin ~ BremEemaxL B PairTXTL real *PairUminLA !min of uniform random number in the region A ! at energies PairEgmin PairdULA
! constants used for pair brem in Air real *mbtoPX0 ! mbtoPgrm x X0g If multiplied to sigma in ! we obtain probability radation length real *X0g ! radiation length in g cm real *cScrC1 ! const which appears in the complete screening ! crossection real *cScrC2 ! the other such one real too **media Zeff **Z **real *BrScrE ! below screened cross section is used real *BremEgmin ! Brems gamma min energy in Gev(100 keV) real *8 BremEemin ! Brems by e+/e- is considered at E > BremEemin !(Ek ~ 300 keV) real *8 BremLEemin ! log10 of BremEemin real *8 BremEemaxL ! Energy of e+/e- below which LPM is negelected ! for brems. 5 *rho_Pb/rho GeV. integer BremTXTL ! Size of the Brems total x-section table. ! in the energy region BremEemin ~ BremEemaxL B PairTXTL real *PairUminLA !min of uniform random number in the region A ! at energies PairEgmin X0g * muPrEmax
! constants used for pair brem in Air real *mbtoPX0 ! mbtoPgrm x X0g If multiplied to sigma in ! we obtain probability radation length real *X0g ! radiation length in g cm real *cScrC1 ! const which appears in the complete screening ! crossection real *cScrC2 ! the other such one real too **media Zeff **Z **real *BrScrE ! below screened cross section is used real *BremEgmin ! Brems gamma min energy in Gev(100 keV) real *8 BremEemin ! Brems by e+/e- is considered at E > BremEemin !(Ek ~ 300 keV) real *8 BremLEemin ! log10 of BremEemin real *8 BremEemaxL ! Energy of e+/e- below which LPM is negelected ! for brems. 5 *rho_Pb/rho GeV. integer BremTXTL ! Size of the Brems total x-section table. ! in the energy region BremEemin ~ BremEemaxL B PairTXTL real *PairUminLA !min of uniform random number in the region A ! at energies PairEgmin X0g * muNEmax1
! constants used for pair brem in Air real *mbtoPX0 ! mbtoPgrm x X0g If multiplied to sigma in ! we obtain probability radation length real *X0g ! radiation length in g cm real *cScrC1 ! const which appears in the complete screening ! crossection real *cScrC2 ! the other such one real too **media Zeff **Z **real *BrScrE ! below screened cross section is used real *BremEgmin ! Brems gamma min energy in Gev(100 keV) real *8 BremEemin ! Brems by e+/e- is considered at E > BremEemin !(Ek ~ 300 keV) real *8 BremLEemin ! log10 of BremEemin real *8 BremEemaxL ! Energy of e+/e- below which LPM is negelected ! for brems. 5 *rho_Pb/rho GeV. integer BremTXTL ! Size of the Brems total x-section table. ! in the energy region BremEemin ~ BremEemaxL B PairTXTL real *PairUminLA !min of uniform random number in the region A ! at energies PairEgmin PairdULB
! constants used for pair brem in Air real *mbtoPX0 ! mbtoPgrm x X0g If multiplied to sigma in ! we obtain probability radation length real *X0g ! radiation length in g cm real *cScrC1 ! const which appears in the complete screening ! crossection real *cScrC2 ! the other such one real too **media Zeff **Z **real *BrScrE ! below screened cross section is used real *BremEgmin ! Brems gamma min energy in Gev(100 keV) real *8 BremEemin ! Brems by e+/e- is considered at E > BremEemin !(Ek ~ 300 keV) real *8 BremLEemin ! log10 of BremEemin real *8 BremEemaxL ! Energy of e+/e- below which LPM is negelected ! for brems. 5 *rho_Pb/rho GeV. integer BremTXTL ! Size of the Brems total x-section table. ! in the energy region BremEemin ~ BremEemaxL B PairTXTL real *PairUminLA !min of uniform random number in the region A ! at energies PairEgmin BrUmaxSA
! constants used for pair brem in Air real *mbtoPX0 ! mbtoPgrm x X0g If multiplied to sigma in ! we obtain probability radation length real *X0g ! radiation length in g cm real *cScrC1 ! const which appears in the complete screening ! crossection real *cScrC2 ! the other such one real too **media Zeff **Z **real *BrScrE ! below screened cross section is used real *BremEgmin ! Brems gamma min energy in Gev(100 keV) real *8 BremEemin ! Brems by e+/e- is considered at E > BremEemin !(Ek ~ 300 keV) real *8 BremLEemin ! log10 of BremEemin real *8 BremEemaxL ! Energy of e+/e- below which LPM is negelected ! for brems. 5 *rho_Pb/rho GeV. integer BremTXTL ! Size of the Brems total x-section table. ! in the energy region BremEemin ~ BremEemaxL B PairTXTL real *PairUminLA !min of uniform random number in the region A ! at energies PairEgmin X0g muPrdE common Zicnst BremUszLB
integer Charge2heavyG charge of heavy rightarrow heavy group index conversion array integer HeavyG2massN heavy group index rightarrow mass number conversion array integer HeavyG2charge heavy group index rightarrow charge of heavy conversion array integer HeavyG2code heavy group index rightarrow particle code conversion array integer Code2massN particle code rightarrow mass number conversion array integer Code2heavyG particle code rightarrow heavy group index conversion array real *FragmentTbl the number of interacting nucleons among a projectile heavy nucleus is ! determined as the number of first collision of each interacting nucleon inside ! the nucleus If
! constants used for pair brem in Air real *mbtoPX0 ! mbtoPgrm x X0g If multiplied to sigma in ! we obtain probability radation length real *X0g ! radiation length in g cm real *cScrC1 ! const which appears in the complete screening ! crossection real *cScrC2 ! the other such one real too **media Zeff **Z **real *BrScrE ! below screened cross section is used real *BremEgmin ! Brems gamma min energy in Gev(100 keV) real *8 BremEemin ! Brems by e+/e- is considered at E > BremEemin !(Ek ~ 300 keV) real *8 BremLEemin ! log10 of BremEemin real *8 BremEemaxL ! Energy of e+/e- below which LPM is negelected ! for brems. 5 *rho_Pb/rho GeV. integer BremTXTL ! Size of the Brems total x-section table. ! in the energy region BremEemin ~ BremEemaxL B PairTXTL real *PairUminLA !min of uniform random number in the region A ! at energies PairEgmin X0g muNEmin
integer maxbin nregion c minsave drx2 ! drawing region real maxdep integer maxpos integer to
block data include Zlatfit h c fitting region data data data data data d0 data data d0 data data h g *is for param c g data low(1, 1)/1.d-5/
! constants used for pair brem in Air real *mbtoPX0 ! mbtoPgrm x X0g If multiplied to sigma in ! we obtain probability radation length real *X0g ! radiation length in g cm real *cScrC1 ! const which appears in the complete screening ! crossection real *cScrC2 ! the other such one real too **media Zeff **Z **real *BrScrE ! below screened cross section is used real *BremEgmin ! Brems gamma min energy in Gev(100 keV) real *8 BremEemin ! Brems by e+/e- is considered at E > BremEemin !(Ek ~ 300 keV) real *8 BremLEemin ! log10 of BremEemin real *8 BremEemaxL ! Energy of e+/e- below which LPM is negelected ! for brems. 5 *rho_Pb/rho GeV. integer BremTXTL ! Size of the Brems total x-section table. ! in the energy region BremEemin ~ BremEemaxL B PairTXTL real *PairUminLA !min of uniform random number in the region A ! at energies PairEgmin X0g muPrdE common Zicnst BrUszSB muNUsize muPrTXT
! constants used for pair brem in Air real *mbtoPX0 ! mbtoPgrm x X0g If multiplied to sigma in ! we obtain probability radation length real *X0g ! radiation length in g cm real *cScrC1 ! const which appears in the complete screening ! crossection real *cScrC2 ! the other such one real too **media Zeff **Z **real *BrScrE ! below screened cross section is used real *BremEgmin ! Brems gamma min energy in Gev(100 keV) real *8 BremEemin ! Brems by e+/e- is considered at E > BremEemin !(Ek ~ 300 keV) real *8 BremLEemin ! log10 of BremEemin real *8 BremEemaxL ! Energy of e+/e- below which LPM is negelected ! for brems. 5 *rho_Pb/rho GeV. integer BremTXTL ! Size of the Brems total x-section table. ! in the energy region BremEemin ~ BremEemaxL B PairTXTL real *PairUminLA !min of uniform random number in the region A ! at energies PairEgmin BrScrE
! constants used for pair brem in Air real *mbtoPX0 ! mbtoPgrm x X0g If multiplied to sigma in ! we obtain probability radation length real *X0g ! radiation length in g cm real *cScrC1 ! const which appears in the complete screening ! crossection real *cScrC2 ! the other such one real too **media Zeff **Z **real *BrScrE ! below screened cross section is used real *BremEgmin ! Brems gamma min energy in Gev(100 keV) real *8 BremEemin ! Brems by e+/e- is considered at E > BremEemin !(Ek ~ 300 keV) real *8 BremLEemin ! log10 of BremEemin real *8 BremEemaxL ! Energy of e+/e- below which LPM is negelected ! for brems. 5 *rho_Pb/rho GeV. integer BremTXTL ! Size of the Brems total x-section table. ! in the energy region BremEemin ~ BremEemaxL B PairTXTL real *PairUminLA !min of uniform random number in the region A ! at energies PairEgmin X0g muNdE
! constants used for pair brem in Air real *mbtoPX0 ! mbtoPgrm x X0g If multiplied to sigma in ! we obtain probability radation length real *X0g ! radiation length in g cm real *cScrC1 ! const which appears in the complete screening ! crossection real *cScrC2 ! the other such one real too **media Zeff **Z **real *BrScrE ! below screened cross section is used real *BremEgmin ! Brems gamma min energy in Gev(100 keV) real *8 BremEemin ! Brems by e+/e- is considered at E > BremEemin !(Ek ~ 300 keV) real *8 BremLEemin ! log10 of BremEemin real *8 BremEemaxL ! Energy of e+/e- below which LPM is negelected ! for brems. 5 *rho_Pb/rho GeV. integer BremTXTL ! Size of the Brems total x-section table. ! in the energy region BremEemin ~ BremEemaxL B PairTXTL real *PairUminLA !min of uniform random number in the region A ! at energies PairEgmin X0g muNdETX
! constants used for pair brem in Air real *mbtoPX0 ! mbtoPgrm x X0g If multiplied to sigma in ! we obtain probability radation length real *X0g ! radiation length in g cm real *cScrC1 ! const which appears in the complete screening ! crossection real *cScrC2 ! the other such one real too **media Zeff **Z **real *BrScrE ! below screened cross section is used real *BremEgmin ! Brems gamma min energy in Gev(100 keV) real *8 BremEemin ! Brems by e+/e- is considered at E > BremEemin !(Ek ~ 300 keV) real *8 BremLEemin ! log10 of BremEemin real *8 BremEemaxL ! Energy of e+/e- below which LPM is negelected ! for brems. 5 *rho_Pb/rho GeV. integer BremTXTL ! Size of the Brems total x-section table. ! in the energy region BremEemin ~ BremEemaxL B PairTXTL real *PairUminLA !min of uniform random number in the region A ! at energies PairEgmin * BrdETXS
! Parameters used for hadronic cascade shower is generated newline ! For you may give as as or em quick generation of AS for heavy primaries is tried See chookASbyH f character *Generate2 don t touch this for skeleton flesh use integer MagBrem no magnetic bremsstrahlung is considered newline ! if and Ee energy loss due to magnetic brems is considered newline ! if and Ee real sampling of gamma is performed WaitRatio ! must be made small so that WaitRatio *E0 sim MagBremEmin integer MagPair no magnetic pair creation is considered newline ! if and Eg real sampling is the LPM effect is considered when Ee LpmBremEmin for electrons and ! Eg LpmPairEmin for gamma rays real *MagBremEmin E magnetic bremsstrahlung by electrons may be considered However
! constants used for pair brem in Air real *mbtoPX0 ! mbtoPgrm x X0g If multiplied to sigma in ! we obtain probability radation length real *X0g ! radiation length in g cm real *cScrC1 ! const which appears in the complete screening ! crossection real *cScrC2 ! the other such one real too **media Zeff **Z **real *BrScrE ! below screened cross section is used real *BremEgmin ! Brems gamma min energy in Gev(100 keV) real *8 BremEemin ! Brems by e+/e- is considered at E > BremEemin !(Ek ~ 300 keV) real *8 BremLEemin ! log10 of BremEemin real *8 BremEemaxL ! Energy of e+/e- below which LPM is negelected ! for brems. 5 *rho_Pb/rho GeV. integer BremTXTL ! Size of the Brems total x-section table. ! in the energy region BremEemin ~ BremEemaxL B PairTXTL real *PairUminLA !min of uniform random number in the region A ! at energies PairEgmin mbtoPX0
! constants used for pair brem in Air real *mbtoPX0 ! mbtoPgrm x X0g If multiplied to sigma in ! we obtain probability radation length real *X0g ! radiation length in g cm real *cScrC1 ! const which appears in the complete screening ! crossection real *cScrC2 ! the other such one real too **media Zeff **Z **real *BrScrE ! below screened cross section is used real *BremEgmin ! Brems gamma min energy in Gev(100 keV) real *8 BremEemin ! Brems by e+/e- is considered at E > BremEemin !(Ek ~ 300 keV) real *8 BremLEemin ! log10 of BremEemin real *8 BremEemaxL ! Energy of e+/e- below which LPM is negelected ! for brems. 5 *rho_Pb/rho GeV. integer BremTXTL ! Size of the Brems total x-section table. ! in the energy region BremEemin ~ BremEemaxL B PairTXTL real *PairUminLA !min of uniform random number in the region A ! at energies PairEgmin X0g muPrEmax1
! constants used for pair brem in Air real *mbtoPX0 ! mbtoPgrm x X0g If multiplied to sigma in ! we obtain probability radation length real *X0g ! radiation length in g cm real *cScrC1 ! const which appears in the complete screening ! crossection real *cScrC2 ! the other such one real too **media Zeff **Z **real *BrScrE ! below screened cross section is used real *BremEgmin ! Brems gamma min energy in Gev(100 keV) real *8 BremEemin ! Brems by e+/e- is considered at E > BremEemin !(Ek ~ 300 keV) real *8 BremLEemin ! log10 of BremEemin real *8 BremEemaxL ! Energy of e+/e- below which LPM is negelected ! for brems. 5 *rho_Pb/rho GeV. integer BremTXTL ! Size of the Brems total x-section table. ! in the energy region BremEemin ~ BremEemaxL B PairTXTL real *PairUminLA !min of uniform random number in the region A ! at energies PairEgmin X0g muNpwtx
! constants used for pair brem in Air real *mbtoPX0 ! mbtoPgrm x X0g If multiplied to sigma in ! we obtain probability radation length real *X0g ! radiation length in g cm real *cScrC1 ! const which appears in the complete screening ! crossection real *cScrC2 ! the other such one real too **media Zeff **Z **real *BrScrE ! below screened cross section is used real *BremEgmin ! Brems gamma min energy in Gev(100 keV) real *8 BremEemin ! Brems by e+/e- is considered at E > BremEemin !(Ek ~ 300 keV) real *8 BremLEemin ! log10 of BremEemin real *8 BremEemaxL ! Energy of e+/e- below which LPM is negelected ! for brems. 5 *rho_Pb/rho GeV. integer BremTXTL ! Size of the Brems total x-section table. ! in the energy region BremEemin ~ BremEemaxL B PairTXTL real *PairUminLA !min of uniform random number in the region A ! at energies PairEgmin PairUmaxLA
! muon interaction related variables ! real Zeff3 ! Zeff **real muPrEmax1 integer muPrEsize
! constants used for pair brem in Air real *mbtoPX0 ! mbtoPgrm x X0g If multiplied to sigma in ! we obtain probability radation length real *X0g ! radiation length in g cm real *cScrC1 ! const which appears in the complete screening ! crossection real *cScrC2 ! the other such one real too **media Zeff **Z **real *BrScrE ! below screened cross section is used real *BremEgmin ! Brems gamma min energy in Gev(100 keV) real *8 BremEemin ! Brems by e+/e- is considered at E > BremEemin !(Ek ~ 300 keV) real *8 BremLEemin ! log10 of BremEemin real *8 BremEemaxL ! Energy of e+/e- below which LPM is negelected ! for brems. 5 *rho_Pb/rho GeV. integer BremTXTL ! Size of the Brems total x-section table. ! in the energy region BremEemin ~ BremEemaxL B PairTXTL real *PairUminLA !min of uniform random number in the region A ! at energies PairEgmin BrUminSA
! constants used for pair brem in Air real *mbtoPX0 ! mbtoPgrm x X0g If multiplied to sigma in ! we obtain probability radation length real *X0g ! radiation length in g cm real *cScrC1 ! const which appears in the complete screening ! crossection real *cScrC2 ! the other such one real too **media Zeff **Z **real *BrScrE ! below screened cross section is used real *BremEgmin ! Brems gamma min energy in Gev(100 keV) real *8 BremEemin ! Brems by e+/e- is considered at E > BremEemin !(Ek ~ 300 keV) real *8 BremLEemin ! log10 of BremEemin real *8 BremEemaxL ! Energy of e+/e- below which LPM is negelected ! for brems. 5 *rho_Pb/rho GeV. integer BremTXTL ! Size of the Brems total x-section table. ! in the energy region BremEemin ~ BremEemaxL B PairTXTL real *PairUminLA !min of uniform random number in the region A ! at energies PairEgmin PairUminLA
! muon interaction related variables ! real Zeff3 ! Zeff **real muPrEmax1 integer muPrTXT ! real muBrEmax1 integer muBrTXT ! real muNEmax1 integer muNUsize
! constants used for pair brem in Air real *mbtoPX0 ! mbtoPgrm x X0g If multiplied to sigma in ! we obtain probability radation length real *X0g ! radiation length in g cm real *cScrC1 ! const which appears in the complete screening ! crossection real *cScrC2 ! the other such one real too **media Zeff **Z **real *BrScrE ! below screened cross section is used real *BremEgmin ! Brems gamma min energy in Gev(100 keV) real *8 BremEemin ! Brems by e+/e- is considered at E > BremEemin !(Ek ~ 300 keV) real *8 BremLEemin ! log10 of BremEemin real *8 BremEemaxL ! Energy of e+/e- below which LPM is negelected ! for brems. 5 *rho_Pb/rho GeV. integer BremTXTL ! Size of the Brems total x-section table. ! in the energy region BremEemin ~ BremEemaxL B PairTXTL real *PairUminLA !min of uniform random number in the region A ! at energies PairEgmin PairUmaxLB
! constants used for pair brem in Air real *mbtoPX0 ! mbtoPgrm x X0g If multiplied to sigma in ! we obtain probability radation length real *X0g ! radiation length in g cm real *cScrC1 ! const which appears in the complete screening ! crossection real *cScrC2 ! the other such one real too **media Zeff **Z **real *BrScrE ! below screened cross section is used real *BremEgmin ! Brems gamma min energy in Gev(100 keV) real *8 BremEemin ! Brems by e+/e- is considered at E > BremEemin !(Ek ~ 300 keV) real *8 BremLEemin ! log10 of BremEemin real *8 BremEemaxL ! Energy of e+/e- below which LPM is negelected ! for brems. 5 *rho_Pb/rho GeV. integer BremTXTL ! Size of the Brems total x-section table. ! in the energy region BremEemin ~ BremEemaxL B PairTXTL real *PairUminLA !min of uniform random number in the region A ! at energies PairEgmin X0g muPrdE common Zicnst BrUszSB muNUsize muBrEsize
! constants used for pair brem in Air real *mbtoPX0 ! mbtoPgrm x X0g If multiplied to sigma in ! we obtain probability radation length real *X0g ! radiation length in g cm real *cScrC1 ! const which appears in the complete screening ! crossection real *cScrC2 ! the other such one real too **media Zeff **Z **real *BrScrE ! below screened cross section is used real *BremEgmin ! Brems gamma min energy in Gev(100 keV) real *8 BremEemin ! Brems by e+/e- is considered at E > BremEemin !(Ek ~ 300 keV) real *8 BremLEemin ! log10 of BremEemin real *8 BremEemaxL ! Energy of e+/e- below which LPM is negelected ! for brems. 5 *rho_Pb/rho GeV. integer BremTXTL ! Size of the Brems total x-section table. ! in the energy region BremEemin ~ BremEemaxL B PairTXTL real *PairUminLA !min of uniform random number in the region A ! at energies PairEgmin X0g muPrdE common Zicnst BrUszSB muNUsize muNEsize
! constants used for pair brem in Air real *mbtoPX0 ! mbtoPgrm x X0g If multiplied to sigma in ! we obtain probability radation length real *X0g ! radiation length in g cm real *cScrC1 ! const which appears in the complete screening ! crossection real *cScrC2 ! the other such one real too **media Zeff **Z **real *BrScrE ! below screened cross section is used real *BremEgmin ! Brems gamma min energy in Gev(100 keV) real *8 BremEemin ! Brems by e+/e- is considered at E > BremEemin !(Ek ~ 300 keV) real *8 BremLEemin ! log10 of BremEemin real *8 BremEemaxL ! Energy of e+/e- below which LPM is negelected ! for brems. 5 *rho_Pb/rho GeV. integer BremTXTL ! Size of the Brems total x-section table. ! in the energy region BremEemin ~ BremEemaxL B PairTXTL real *PairUminLA !min of uniform random number in the region A ! at energies PairEgmin X0g muNdU
! constants used for pair brem in Air real *mbtoPX0 ! mbtoPgrm x X0g If multiplied to sigma in ! we obtain probability radation length real *X0g ! radiation length in g cm real *cScrC1 ! const which appears in the complete screening ! crossection real *cScrC2 ! the other such one real too **media Zeff **Z **real *BrScrE ! below screened cross section is used real *BremEgmin ! Brems gamma min energy in Gev(100 keV) real *8 BremEemin ! Brems by e+/e- is considered at E > BremEemin !(Ek ~ 300 keV) real *8 BremLEemin ! log10 of BremEemin real *8 BremEemaxL ! Energy of e+/e- below which LPM is negelected ! for brems. 5 *rho_Pb/rho GeV. integer BremTXTL ! Size of the Brems total x-section table. ! in the energy region BremEemin ~ BremEemaxL B PairTXTL real *PairUminLA !min of uniform random number in the region A ! at energies PairEgmin PairLEgmin
! constants used for pair brem in Air real *mbtoPX0 ! mbtoPgrm x X0g If multiplied to sigma in ! we obtain probability radation length real *X0g ! radiation length in g cm real *cScrC1 ! const which appears in the complete screening ! crossection real *cScrC2 ! the other such one real too **media Zeff **Z **real *BrScrE ! below screened cross section is used real *BremEgmin ! Brems gamma min energy in Gev(100 keV) real *8 BremEemin ! Brems by e+/e- is considered at E > BremEemin !(Ek ~ 300 keV) real *8 BremLEemin ! log10 of BremEemin real *8 BremEemaxL ! Energy of e+/e- below which LPM is negelected ! for brems. 5 *rho_Pb/rho GeV. integer BremTXTL ! Size of the Brems total x-section table. ! in the energy region BremEemin ~ BremEemaxL B PairTXTL real *PairUminLA !min of uniform random number in the region A ! at energies PairEgmin BrdUSA
! constants used for pair brem in Air real *mbtoPX0 ! mbtoPgrm x X0g If multiplied to sigma in ! we obtain probability radation length real *X0g ! radiation length in g cm real *cScrC1 ! const which appears in the complete screening ! crossection real *cScrC2 ! the other such one real too **media Zeff **Z **real *BrScrE ! below screened cross section is used real *BremEgmin ! Brems gamma min energy in Gev(100 keV) real *8 BremEemin ! Brems by e+/e- is considered at E > BremEemin !(Ek ~ 300 keV) real *8 BremLEemin ! log10 of BremEemin real *8 BremEemaxL ! Energy of e+/e- below which LPM is negelected ! for brems. 5 *rho_Pb/rho GeV. integer BremTXTL ! Size of the Brems total x-section table. ! in the energy region BremEemin ~ BremEemaxL B PairTXTL real *PairUminLA !min of uniform random number in the region A ! at energies PairEgmin BrdUSB
! constants used for pair brem in Air real *mbtoPX0 ! mbtoPgrm x X0g If multiplied to sigma in ! we obtain probability radation length real *X0g ! radiation length in g cm real *cScrC1 ! const which appears in the complete screening ! crossection real *cScrC2 ! the other such one real too **media Zeff **Z **real *BrScrE ! below screened cross section is used real *BremEgmin ! Brems gamma min energy in Gev(100 keV) real *8 BremEemin ! Brems by e+/e- is considered at E > BremEemin !(Ek ~ 300 keV) real *8 BremLEemin ! log10 of BremEemin real *8 BremEemaxL ! Energy of e+/e- below which LPM is negelected ! for brems. 5 *rho_Pb/rho GeV. integer BremTXTL ! Size of the Brems total x-section table. ! in the energy region BremEemin ~ BremEemaxL B PairTXTL real *PairUminLA !min of uniform random number in the region A ! at energies PairEgmin X0g muPrLEmin
! constants used for pair brem in Air real *mbtoPX0 ! mbtoPgrm x X0g If multiplied to sigma in ! we obtain probability radation length real *X0g ! radiation length in g cm real *cScrC1 ! const which appears in the complete screening ! crossection real *cScrC2 ! the other such one real too **media Zeff **Z **real *BrScrE ! below screened cross section is used real *BremEgmin ! Brems gamma min energy in Gev(100 keV) real *8 BremEemin ! Brems by e+/e- is considered at E > BremEemin !(Ek ~ 300 keV) real *8 BremLEemin ! log10 of BremEemin real *8 BremEemaxL ! Energy of e+/e- below which LPM is negelected ! for brems. 5 *rho_Pb/rho GeV. integer BremTXTL ! Size of the Brems total x-section table. ! in the energy region BremEemin ~ BremEemaxL B PairTXTL real *PairUminLA !min of uniform random number in the region A ! at energies PairEgmin X0g muPrdU
! constants used for pair brem in Air real *mbtoPX0 ! mbtoPgrm x X0g If multiplied to sigma in ! we obtain probability radation length real *X0g ! radiation length in g cm real *cScrC1 ! const which appears in the complete screening ! crossection real *cScrC2 ! the other such one real too **media Zeff **Z **real *BrScrE ! below screened cross section is used real *BremEgmin ! Brems gamma min energy in Gev(100 keV) real *8 BremEemin ! Brems by e+/e- is considered at E > BremEemin !(Ek ~ 300 keV) real *8 BremLEemin ! log10 of BremEemin real *8 BremEemaxL ! Energy of e+/e- below which LPM is negelected ! for brems. 5 *rho_Pb/rho GeV. integer BremTXTL ! Size of the Brems total x-section table. ! in the energy region BremEemin ~ BremEemaxL B PairTXTL real *PairUminLA !min of uniform random number in the region A ! at energies PairEgmin cScrMain
! constants used for pair brem in Air real *mbtoPX0 ! mbtoPgrm x X0g If multiplied to sigma in ! we obtain probability radation length real *X0g ! radiation length in g cm real *cScrC1 ! const which appears in the complete screening ! crossection real *cScrC2 ! the other such one real too **media Zeff **Z **real *BrScrE ! below screened cross section is used real *BremEgmin ! Brems gamma min energy in Gev(100 keV) real *8 BremEemin ! Brems by e+/e- is considered at E > BremEemin !(Ek ~ 300 keV) real *8 BremLEemin ! log10 of BremEemin real *8 BremEemaxL ! Energy of e+/e- below which LPM is negelected ! for brems. 5 *rho_Pb/rho GeV. integer BremTXTL ! Size of the Brems total x-section table. ! in the energy region BremEemin ~ BremEemaxL B PairTXTL real *PairUminLA !min of uniform random number in the region A ! at energies PairEgmin BremLEemin
! constants used for pair brem in Air real *mbtoPX0 ! mbtoPgrm x X0g If multiplied to sigma in ! we obtain probability radation length real *X0g ! radiation length in g cm real *cScrC1 ! const which appears in the complete screening ! crossection real *cScrC2 ! the other such one real too **media Zeff **Z **real *BrScrE ! below screened cross section is used real *BremEgmin ! Brems gamma min energy in Gev(100 keV) real *8 BremEemin ! Brems by e+/e- is considered at E > BremEemin !(Ek ~ 300 keV) real *8 BremLEemin ! log10 of BremEemin real *8 BremEemaxL ! Energy of e+/e- below which LPM is negelected ! for brems. 5 *rho_Pb/rho GeV. integer BremTXTL ! Size of the Brems total x-section table. ! in the energy region BremEemin ~ BremEemaxL B PairTXTL real *PairUminLA !min of uniform random number in the region A ! at energies PairEgmin BremEemin
! constants used for pair brem in Air real *mbtoPX0 ! mbtoPgrm x X0g If multiplied to sigma in ! we obtain probability radation length real *X0g ! radiation length in g cm real *cScrC1 ! const which appears in the complete screening ! crossection real *cScrC2 ! the other such one real too **media Zeff **Z **real *BrScrE ! below screened cross section is used real *BremEgmin ! Brems gamma min energy in Gev(100 keV) real *8 BremEemin ! Brems by e+/e- is considered at E > BremEemin !(Ek ~ 300 keV) real *8 BremLEemin ! log10 of BremEemin real *8 BremEemaxL ! Energy of e+/e- below which LPM is negelected ! for brems. 5 *rho_Pb/rho GeV. integer BremTXTL ! Size of the Brems total x-section table. ! in the energy region BremEemin ~ BremEemaxL B PairTXTL real *PairUminLA !min of uniform random number in the region A ! at energies PairEgmin X0g muPrdE common Zicnst BrUszSB muNUsize muPrUsize
! constants used for pair brem in Air real *mbtoPX0 ! mbtoPgrm x X0g If multiplied to sigma in ! we obtain probability radation length real *X0g ! radiation length in g cm real *cScrC1 ! const which appears in the complete screening ! crossection real *cScrC2 ! the other such one real too **media Zeff **Z **real *BrScrE ! below screened cross section is used real *BremEgmin ! Brems gamma min energy in Gev(100 keV) real *8 BremEemin ! Brems by e+/e- is considered at E > BremEemin !(Ek ~ 300 keV) real *8 BremLEemin ! log10 of BremEemin real *8 BremEemaxL ! Energy of e+/e- below which LPM is negelected ! for brems. 5 *rho_Pb/rho GeV. integer BremTXTL ! Size of the Brems total x-section table. ! in the energy region BremEemin ~ BremEemaxL B PairTXTL real *PairUminLA !min of uniform random number in the region A ! at energies PairEgmin X0g muPrVmin
! constants used for pair brem in Air real *mbtoPX0 ! mbtoPgrm x X0g If multiplied to sigma in ! we obtain probability radation length real *X0g ! radiation length in g cm real *cScrC1 ! const which appears in the complete screening ! crossection real *cScrC2 ! the other such one real too **media Zeff **Z **real *BrScrE ! below screened cross section is used real *BremEgmin ! Brems gamma min energy in Gev(100 keV) real *8 BremEemin ! Brems by e+/e- is considered at E > BremEemin !(Ek ~ 300 keV) real *8 BremLEemin ! log10 of BremEemin real *8 BremEemaxL ! Energy of e+/e- below which LPM is negelected ! for brems. 5 *rho_Pb/rho GeV. integer BremTXTL ! Size of the Brems total x-section table. ! in the energy region BremEemin ~ BremEemaxL B PairTXTL real *PairUminLA !min of uniform random number in the region A ! at energies PairEgmin BremUmaxLA
! constants used for pair brem in Air real *mbtoPX0 ! mbtoPgrm x X0g If multiplied to sigma in ! we obtain probability radation length real *X0g ! radiation length in g cm real *cScrC1 ! const which appears in the complete screening ! crossection real *cScrC2 ! the other such one real too **media Zeff **Z **real *BrScrE ! below screened cross section is used real *BremEgmin ! Brems gamma min energy in Gev(100 keV) real *8 BremEemin ! Brems by e+/e- is considered at E > BremEemin !(Ek ~ 300 keV) real *8 BremLEemin ! log10 of BremEemin real *8 BremEemaxL ! Energy of e+/e- below which LPM is negelected ! for brems. 5 *rho_Pb/rho GeV. integer BremTXTL ! Size of the Brems total x-section table. ! in the energy region BremEemin ~ BremEemaxL B PairTXTL real *PairUminLA !min of uniform random number in the region A ! at energies PairEgmin X0g muPrdE common Zicnst * PairUszLA
! constants used for pair brem in Air real *mbtoPX0 ! mbtoPgrm x X0g If multiplied to sigma in ! we obtain probability radation length real *X0g ! radiation length in g cm real *cScrC1 ! const which appears in the complete screening ! crossection real *cScrC2 ! the other such one real too **media Zeff **Z **real *BrScrE ! below screened cross section is used real *BremEgmin ! Brems gamma min energy in Gev(100 keV) real *8 BremEemin ! Brems by e+/e- is considered at E > BremEemin !(Ek ~ 300 keV) real *8 BremLEemin ! log10 of BremEemin real *8 BremEemaxL ! Energy of e+/e- below which LPM is negelected ! for brems. 5 *rho_Pb/rho GeV. integer BremTXTL ! Size of the Brems total x-section table. ! in the energy region BremEemin ~ BremEemaxL B PairTXTL real *PairUminLA !min of uniform random number in the region A ! at energies PairEgmin X0g muNLEmin
dE dx *! Nuc Int sampling table g
! constants used for pair brem in Air real *mbtoPX0 ! mbtoPgrm x X0g If multiplied to sigma in ! we obtain probability radation length real *X0g ! radiation length in g cm real *cScrC1 ! const which appears in the complete screening ! crossection real *cScrC2 ! the other such one real too **media Zeff **Z **real *BrScrE ! below screened cross section is used real *BremEgmin ! Brems gamma min energy in Gev(100 keV) real *8 BremEemin ! Brems by e+/e- is considered at E > BremEemin !(Ek ~ 300 keV) real *8 BremLEemin ! log10 of BremEemin real *8 BremEemaxL ! Energy of e+/e- below which LPM is negelected ! for brems. 5 *rho_Pb/rho GeV. integer BremTXTL ! Size of the Brems total x-section table. ! in the energy region BremEemin ~ BremEemaxL B PairTXTL real *PairUminLA !min of uniform random number in the region A ! at energies PairEgmin PairdETXL
! constants used for pair brem in Air real *mbtoPX0 ! mbtoPgrm x X0g If multiplied to sigma in ! we obtain probability radation length real *X0g ! radiation length in g cm real *cScrC1 ! const which appears in the complete screening ! crossection real *cScrC2 ! the other such one real too **media Zeff **Z **real *BrScrE ! below screened cross section is used real *BremEgmin ! Brems gamma min energy in Gev(100 keV) real *8 BremEemin ! Brems by e+/e- is considered at E > BremEemin !(Ek ~ 300 keV) real *8 BremLEemin ! log10 of BremEemin real *8 BremEemaxL ! Energy of e+/e- below which LPM is negelected ! for brems. 5 *rho_Pb/rho GeV. integer BremTXTL ! Size of the Brems total x-section table. ! in the energy region BremEemin ~ BremEemaxL B PairTXTL real *PairUminLA !min of uniform random number in the region A ! at energies PairEgmin BremEgmin
! constants used for pair brem in Air real *mbtoPX0 ! mbtoPgrm x X0g If multiplied to sigma in ! we obtain probability radation length real *X0g ! radiation length in g cm real *cScrC1 ! const which appears in the complete screening ! crossection real *cScrC2 ! the other such one real too **media Zeff **Z **real *BrScrE ! below screened cross section is used real *BremEgmin ! Brems gamma min energy in Gev(100 keV) real *8 BremEemin ! Brems by e+/e- is considered at E > BremEemin !(Ek ~ 300 keV) real *8 BremLEemin ! log10 of BremEemin real *8 BremEemaxL ! Energy of e+/e- below which LPM is negelected ! for brems. 5 *rho_Pb/rho GeV. integer BremTXTL ! Size of the Brems total x-section table. ! in the energy region BremEemin ~ BremEemaxL B PairTXTL real *PairUminLA !min of uniform random number in the region A ! at energies PairEgmin BremUminLB
! constants used for pair brem in Air real *mbtoPX0 ! mbtoPgrm x X0g If multiplied to sigma in ! we obtain probability radation length real *X0g ! radiation length in g cm real *cScrC1 ! const which appears in the complete screening ! crossection real *cScrC2 ! the other such one real too **media Zeff **Z **real *BrScrE ! below screened cross section is used real *BremEgmin ! Brems gamma min energy in Gev(100 keV) real *8 BremEemin ! Brems by e+/e- is considered at E > BremEemin !(Ek ~ 300 keV) real *8 BremLEemin ! log10 of BremEemin real *8 BremEemaxL ! Energy of e+/e- below which LPM is negelected ! for brems. 5 *rho_Pb/rho GeV. integer BremTXTL ! Size of the Brems total x-section table. ! in the energy region BremEemin ~ BremEemaxL B PairTXTL real *PairUminLA !min of uniform random number in the region A ! at energies PairEgmin X0g muBrdU
*Zfirst p fm *Zfirst p Zfirst p Zfirst p *Zfirst p *Zfirst pos xyz Zfirst pos xyz *Zfirst pos *Zfirst pos Zfirst pos *Zfirst pos *Zfirst Zfirst vec w *Zfirst vec w Zfirst vec *Zfirst Zfirst where
! constants used for pair brem in Air real *mbtoPX0 ! mbtoPgrm x X0g If multiplied to sigma in ! we obtain probability radation length real *X0g ! radiation length in g cm real *cScrC1 ! const which appears in the complete screening ! crossection real *cScrC2 ! the other such one real too **media Zeff **Z **real *BrScrE ! below screened cross section is used real *BremEgmin ! Brems gamma min energy in Gev(100 keV) real *8 BremEemin ! Brems by e+/e- is considered at E > BremEemin !(Ek ~ 300 keV) real *8 BremLEemin ! log10 of BremEemin real *8 BremEemaxL ! Energy of e+/e- below which LPM is negelected ! for brems. 5 *rho_Pb/rho GeV. integer BremTXTL ! Size of the Brems total x-section table. ! in the energy region BremEemin ~ BremEemaxL B PairTXTL real *PairUminLA !min of uniform random number in the region A ! at energies PairEgmin cScrC2
! constants used for pair brem in Air real *mbtoPX0 ! mbtoPgrm x X0g If multiplied to sigma in ! we obtain probability radation length real *X0g ! radiation length in g cm real *cScrC1 ! const which appears in the complete screening ! crossection real *cScrC2 ! the other such one real too **media Zeff **Z **real *BrScrE ! below screened cross section is used real *BremEgmin ! Brems gamma min energy in Gev(100 keV) real *8 BremEemin ! Brems by e+/e- is considered at E > BremEemin !(Ek ~ 300 keV) real *8 BremLEemin ! log10 of BremEemin real *8 BremEemaxL ! Energy of e+/e- below which LPM is negelected ! for brems. 5 *rho_Pb/rho GeV. integer BremTXTL ! Size of the Brems total x-section table. ! in the energy region BremEemin ~ BremEemaxL B PairTXTL real *PairUminLA !min of uniform random number in the region A ! at energies PairEgmin BremdULB
! constants used for pair brem in Air real *mbtoPX0 ! mbtoPgrm x X0g If multiplied to sigma in ! we obtain probability radation length real *X0g ! radiation length in g cm real *cScrC1 ! const which appears in the complete screening ! crossection real *cScrC2 ! the other such one real too **media Zeff **Z **real *BrScrE ! below screened cross section is used real *BremEgmin ! Brems gamma min energy in Gev(100 keV) real *8 BremEemin ! Brems by e+/e- is considered at E > BremEemin !(Ek ~ 300 keV) real *8 BremLEemin ! log10 of BremEemin real *8 BremEemaxL ! Energy of e+/e- below which LPM is negelected ! for brems. 5 *rho_Pb/rho GeV. integer BremTXTL ! Size of the Brems total x-section table. ! in the energy region BremEemin ~ BremEemaxL B PairTXTL real *PairUminLA !min of uniform random number in the region A ! at energies PairEgmin BremdETXL
integer maxbin nregion c minsave drx2 ! drawing region real maxdep integer maxpos integer from
! constants used for pair brem in Air real *mbtoPX0 ! mbtoPgrm x X0g If multiplied to sigma in ! we obtain probability radation length real *X0g ! radiation length in g cm real *cScrC1 ! const which appears in the complete screening ! crossection real *cScrC2 ! the other such one real too **media Zeff **Z **real *BrScrE ! below screened cross section is used real *BremEgmin ! Brems gamma min energy in Gev(100 keV) real *8 BremEemin ! Brems by e+/e- is considered at E > BremEemin !(Ek ~ 300 keV) real *8 BremLEemin ! log10 of BremEemin real *8 BremEemaxL ! Energy of e+/e- below which LPM is negelected ! for brems. 5 *rho_Pb/rho GeV. integer BremTXTL ! Size of the Brems total x-section table. ! in the energy region BremEemin ~ BremEemaxL B PairTXTL real *PairUminLA !min of uniform random number in the region A ! at energies PairEgmin X0g muPrdETX
! muon interaction related variables ! real Zeff3 ! Zeff **real muPrdE
block data cblkIncident data *Za1ry is
! constants used for pair brem in Air real *mbtoPX0 ! mbtoPgrm x X0g If multiplied to sigma in ! we obtain probability radation length real *X0g ! radiation length in g cm real *cScrC1 ! const which appears in the complete screening ! crossection real *cScrC2 ! the other such one real too **media Zeff **Z **real *BrScrE ! below screened cross section is used real *BremEgmin ! Brems gamma min energy in Gev(100 keV) real *8 BremEemin ! Brems by e+/e- is considered at E > BremEemin !(Ek ~ 300 keV) real *8 BremLEemin ! log10 of BremEemin real *8 BremEemaxL ! Energy of e+/e- below which LPM is negelected ! for brems. 5 *rho_Pb/rho GeV. integer BremTXTL ! Size of the Brems total x-section table. ! in the energy region BremEemin ~ BremEemaxL B PairTXTL real *PairUminLA !min of uniform random number in the region A ! at energies PairEgmin X0g muPrdE common Zicnst BrTXTS
! constants used for pair brem in Air real *mbtoPX0 ! mbtoPgrm x X0g If multiplied to sigma in ! we obtain probability radation length real *X0g ! radiation length in g cm real *cScrC1 ! const which appears in the complete screening ! crossection real *cScrC2 ! the other such one real too **media Zeff **Z **real *BrScrE ! below screened cross section is used real *BremEgmin ! Brems gamma min energy in Gev(100 keV) real *8 BremEemin ! Brems by e+/e- is considered at E > BremEemin !(Ek ~ 300 keV) real *8 BremLEemin ! log10 of BremEemin real *8 BremEemaxL ! Energy of e+/e- below which LPM is negelected ! for brems. 5 *rho_Pb/rho GeV. integer BremTXTL ! Size of the Brems total x-section table. ! in the energy region BremEemin ~ BremEemaxL B PairTXTL real *PairUminLA !min of uniform random number in the region A ! at energies PairEgmin X0g muPrdE common Zicnst BrUszSB muNUsize muBrUsize
! structure defining a particle at production ! Basic idea of what is to be contained in ! the particle structue is that dynamical ones should be included those derivable from the particle code ! is not included ******************************************************type fmom momentum sequence union map real e endmap map real * x
! constants used for pair brem in Air real *mbtoPX0 ! mbtoPgrm x X0g If multiplied to sigma in ! we obtain probability radation length real *X0g ! radiation length in g cm real *cScrC1 ! const which appears in the complete screening ! crossection real *cScrC2 ! the other such one real too **media Zeff **Z **real *BrScrE ! below screened cross section is used real *BremEgmin ! Brems gamma min energy in Gev(100 keV) real *8 BremEemin ! Brems by e+/e- is considered at E > BremEemin !(Ek ~ 300 keV) real *8 BremLEemin ! log10 of BremEemin real *8 BremEemaxL ! Energy of e+/e- below which LPM is negelected ! for brems. 5 *rho_Pb/rho GeV. integer BremTXTL ! Size of the Brems total x-section table. ! in the energy region BremEemin ~ BremEemaxL B PairTXTL real *PairUminLA !min of uniform random number in the region A ! at energies PairEgmin X0g muPrdE common Zicnst * BremTXTL
! constants used for pair brem in Air real *mbtoPX0 ! mbtoPgrm x X0g If multiplied to sigma in ! we obtain probability radation length real *X0g ! radiation length in g cm real *cScrC1 ! const which appears in the complete screening ! crossection real *cScrC2 ! the other such one real too **media Zeff **Z **real *BrScrE ! below screened cross section is used real *BremEgmin ! Brems gamma min energy in Gev(100 keV) real *8 BremEemin ! Brems by e+/e- is considered at E > BremEemin !(Ek ~ 300 keV) real *8 BremLEemin ! log10 of BremEemin real *8 BremEemaxL ! Energy of e+/e- below which LPM is negelected ! for brems. 5 *rho_Pb/rho GeV. integer BremTXTL ! Size of the Brems total x-section table. ! in the energy region BremEemin ~ BremEemaxL B PairTXTL real *PairUminLA !min of uniform random number in the region A ! at energies PairEgmin X0g muPrdE common Zicnst BrES
! constants used for pair brem in Air real *mbtoPX0 ! mbtoPgrm x X0g If multiplied to sigma in ! we obtain probability radation length real *X0g ! radiation length in g cm real *cScrC1 ! const which appears in the complete screening ! crossection real *cScrC2 ! the other such one real too **media Zeff **Z **real *BrScrE ! below screened cross section is used real *BremEgmin ! Brems gamma min energy in Gev(100 keV) real *8 BremEemin ! Brems by e+/e- is considered at E > BremEemin !(Ek ~ 300 keV) real *8 BremLEemin ! log10 of BremEemin real *8 BremEemaxL ! Energy of e+/e- below which LPM is negelected ! for brems. 5 *rho_Pb/rho GeV. integer BremTXTL ! Size of the Brems total x-section table. ! in the energy region BremEemin ~ BremEemaxL B PairTXTL real *PairUminLA !min of uniform random number in the region A ! at energies PairEgmin X0g * muBrdETX