CERN Program Library Long Writeup W5013 - CERNLIB ...

More documents

Recommendations

Info

100 E+ dep − E− dep √ (%) σ+ 2 + σ− 2 Depth C Pb (X 0 units) No e ± diff e ± diff No e ± diff e ± diff 0.5 -11.7 -13.0 -0.8 -3.9 1.0 -5.3 -4.9 -1.0 -4.1 1.5 +7.3 +8.0 -1.4 -3.5 2.0 +7.1 +5.3 -0.7 -0.0 2.5 +4.9 +4.3 +1.7 +3.6 3.0 +4.8 +4.1 +1.1 +4.3 3.5 +3.3 +2.7 +2.7 +3.1 4.0 +3.6 +5.3 +2.9 +3.0 4.5 +1.7 +2.8 +0.5 +2.3 5.0 +3.4 +3.5 -1.9 +1.8 Table 1.5: Difference in the shower development for e − /e + in Carbon and Lead. No diff refers to the value without the correction for the difference e − /e + . We have defined weight factors F l and F σ for the positron continuous energy loss and discrete bremsstrahlung cross section: F l = 1 ɛ 0 ∫ ɛ0 0 f(ɛ)dɛ F σ = 1 ∫ 1 f(ɛ)dɛ (13) 1 − ɛ 0 ɛ 0 where ɛ 0 = kc T and k c is the photon cut BCUTE. In this scheme the positron energy loss and discrete bremsstrahlung can be calculated as: ( − dE dx ) + = F l ( − dE dx ) − σ + brems = F σσ − brems As in this approximation the photon spectra are identical, the same SUBROUTINE is used for generating e − /e + bremsstrahlung. The following relations hold: F σ = η(1 − ɛ 0 ) 1 η −1 η 1 − (1 − ɛ 0) ɛ 0 = η ⇒ { Fl >η F σ
Geant 3.16 GEANT User’s Guide PHYS341 Origin : G.N. Patrick, L.Urbán Submitted: 26.10.84 Revision : Revised: 16.12.93 Documentation : 1 Subroutines Simulation of discrete bremsstrahlung by electrons CALL GBREME GBREME generates a bremsstrahlung photon from an electron as a discrete process. The photon energy is sampled from a parameterisation of the bremsstrahlung cross-section of Seltzer and Berger [75] for electron energies below 10 GeV, and from the screened Bethe-Heitler cross-section above 10 GeV. Midgal corrections are applied in both cases. The angular distribution of the photon is calculated by the function GBTETH. Input : via common block /GCTRAK/ Output: via common block /GCKING/. GBREME is called from the tracking routine GTELEC when the parent electron reaches a radiation point during tracking. THETA = GBTETH (ENER,PARTM,EFRAC) THETA ENER PARTM EFRAC (REAL) angle of the radiated photon or e ± pair; (REAL) energy of the particle; (REAL) mass of the particle; (REAL) ratio between the energy of the photon and the energy of the particle. GBTETH calculates the angular distribution of the e + e − -pair in photon pair production and of the emitted photon in µ and e ± bremsstrahlung. GBTETH is called by GBREME. 2 Method The photon energy is sampled according to the Seltzer and Berger bremsstrahlung spectrum [75]. Seltzer and Berger have calculated the spectra for materials with atomic numbers Z = 6,13,29,47,74,92 in the electron (kinetic) energy range 1 keV - 10 GeV. Their tabulated results have been used as input in a parametrising-fitting procedure. The functional form of the parameterisation for the quantity: S(x) =Ck dσ dk can be written as S(x) = { (1 − ah ɛ)F 1 (δ)+b h ɛ 2 F 2 (δ) T ≥ 1MeV 1+a l x + b l x 2 T
Page 1 and 2:
CERN Program Library Long Writeup W
Page 3 and 4:
Chapter 1: Catalog of Geant section
Page 5 and 6:
IOPA500 KINE001 KINE100 KINE199 KIN
Page 7 and 8:
Geant 3.21 GEANT User’s Guide AAA
Page 9 and 10:
The routines made available for GEA
Page 11 and 12:
outines, then you can type the foll
Page 13 and 14:
AAAA Bibliography [1] H.J. Klein an
Page 15 and 16:
2 Event simulation framework The fr
Page 17 and 18:
• after each tracking step along
Page 19 and 20:
GUDIGI GUOUT GTRIGC UGLAST GLAST GU
Page 21 and 22:
Particles JPART Materials JMATE/JTM
Page 23 and 24:
3 Summary of GEANT data records 3.1
Page 25 and 26:
3.4 User applications KEY N I VAR S
Page 27 and 28:
Geant 3.16 GEANT User’s Guide BAS
Page 29 and 30:
Page 31 and 32:
SUBROUTINE UGEOM + (’TGT ’,2,
Page 33 and 34:
Page 35 and 36:
Page 37 and 38:
Page 39 and 40:
LQ(JHEAD-1) user link IQ(JHEAD+1) I
Page 41 and 42:
Page 43 and 44:
Page 45 and 46:
VALUE = GARNDM (DUMMY) DUMMY (REAL)
Page 47 and 48:
Geant 3.16 GEANT User’s Guide CON
Page 49 and 50:
Page 51 and 52:
Page 53 and 54:
Page 55 and 56:
Page 57 and 58:
1 Energy binning IDECAD Bin number
Page 59 and 60:
Page 61 and 62:
Page 63 and 64:
Page 65 and 66:
Particle No. Mass(GeV) Charge Life
Page 67 and 68:
Page 69 and 70:
¯Ω + ¯Ξ 0 π + 23.60 ¯ΛK + 67
Page 71 and 72:
Geant 3.16 GEANT User’s Guide DRA
Page 73 and 74:
The main drawing routines are: GDRA
Page 75 and 76:
Page 77 and 78:
Another feature which is available
Page 79 and 80:
z y x C CHARACTER*4 CHNAMS(5) INTEG
Page 81 and 82:
Page 83 and 84:
CALL GDRAWC(’OPAL’,2,0.,10.,10.
Page 85 and 86:
CALL GSATT(’HB’,’FILL’,3) C
Page 87 and 88:
CALL GDRAW(’CAVE’,40.,40.,0.,10
Page 89 and 90:
Page 91 and 92:
Page 93 and 94:
Figure 16: Example of use of GDSPEC
Page 95 and 96:
CALL GDOPEN(3) CALL GDRAWC(’ALEF
Page 97 and 98:
6 NKVIEW IVIEW JDRAW NKVIEW JV = IB
Page 99 and 100:
SHAD EDGE when HIDE is ON, selects
Page 101 and 102:
x1xxxx 1xxxxx add the text EVENT NR
Page 103 and 104:
DRAW Bibliography [1] R.Bock et al.
Page 105 and 106:
2 Volumes with contents A volume ca
Page 107 and 108:
The user can position a volume thro
Page 109 and 110:
0 1 2 1 0 1 0 1 2 3 4 5 6 7 Figure
Page 111 and 112:
2 Divisions along arbitrary axis As
Page 113 and 114:
9.BL2 half-length along x of the si
Page 115 and 116:
1.DZ half-length along the z axis;
Page 117 and 118:
ECAL BOX specifications 18/11/93 EC
Page 119 and 120:
ECAL SPHE specifications 18/11/93 E
Page 121 and 122:
Geant 3.16 GEANT User’s Guide GEO
Page 123 and 124:
• if the CHONLY flag is set to MA
Page 125 and 126:
Page 127 and 128:
y y x x y y y y z z x x z z y y x x
Page 129 and 130:
CAVE PAR1 specifications 20/10/94 y
Page 131 and 132:
Page 133 and 134:
Page 135 and 136:
Page 137 and 138:
Page 139 and 140:
Page 141 and 142:
Page 143 and 144:
Dynamic ordering The list of neighb
Page 145 and 146:
Page 147 and 148:
Page 149 and 150:
Page 151 and 152:
are represented by their values. A
Page 153 and 154:
Where • @302 is the block number
Page 155 and 156:
GEOM Bibliography [1] French Standa
Page 157 and 158:
Geant 3.16 GEANT User’s Guide HIT
Page 159 and 160:
Page 161 and 162:
Page 163 and 164:
Page 165 and 166:
Page 167 and 168:
ITRA NUMBV HITS NHITS (INTEGER) is
Page 169 and 170:
Page 171 and 172:
Page 173 and 174:
Page 175 and 176:
Page 177 and 178:
HITS Bibliography 175 HITS510 - 3
Page 179 and 180:
Geant 3.16 GEANT User’s Guide IOP
Page 181 and 182:
0 if only IER] structures read in o
Page 183 and 184:
Geant 3.16 GEANT User’s Guide IOP
Page 185 and 186:
IOPA Bibliography [1] J.Zoll. ZEBRA
Page 187 and 188:
Geant 3.16 GEANT User’s Guide KIN
Page 189 and 190:
Page 191 and 192:
NTRACK ITRA JKINE NTRACK JK 1 2 3 4
Page 193 and 194:
Page 195 and 196:
Geant 3.16 GEANT User’s Guide PHY
Page 197 and 198:
Below are listed the data record ke
Page 199 and 200:
Computation of total crosssection o
Page 201 and 202:
7. Update the number of interaction
Page 203 and 204:
1= generation of secondaries enable
Page 205 and 206:
DEEMAX The formula used by GEANT is
Page 207 and 208:
The second problem has been solved
Page 209 and 210:
Page 211 and 212:
GPHYSI Initialisation of physics pr
Page 213 and 214:
Page 215 and 216:
• the mean number of tries is not
Page 217 and 218:
as function of the variable u = Eθ
Page 219 and 220:
Page 221 and 222:
GEANT 3.16 GEANT User’s Guide PHY
Page 223 and 224:
where N Z(Z i ) A(A i ) ρ σ Avoga
Page 225 and 226:
Word # PHFN bank (data area, contin
Page 227 and 228:
2 Method If the energy of the radia
Page 229 and 230:
The Auger or Coster-Kronig transiti
Page 231 and 232:
Page 233 and 234: where: n number of energy bins q i
Page 235 and 236: • Case dielectric → metal. The
Page 237 and 238: | ⃗ k| = | ⃗ k ′′ | = k =
Page 239 and 240: Case dielectric → metal In this c
Page 241 and 242: 5 Processes involving Čerenkov pho
Page 243 and 244: Much of the difficulty in approxima
Page 245 and 246: CALL GMOLIE (OMEGA,BETA2,DIN*) OMEG
Page 247 and 248: where λ 0 = ¯h/p Compton waveleng
Page 249 and 250: 2.3 Sampling of the distribution fu
Page 251 and 252: X 0 is the radiation length. From (
Page 253 and 254: where we have set Θ=θ/χ α . To
Page 255 and 256: The variable DCUTE in common /GCCUT
Page 257 and 258: 2.2 Total cross-sections The integr
Page 259 and 260: For the other charged particles the
Page 261 and 262: POTI POTIL (REAL) average ionisatio
Page 263 and 264: where φ(λ) = 1 ∫ c+i∞ exp (u
Page 265 and 266: Counts 900 800 Landau 700 40 600 20
Page 267 and 268: Fast simulation for n 3 ≥ 16 If t
Page 269 and 270: Geant 3.16 GEANT User’s Guide PHY
Page 271 and 272: ξ/E Max (hereafter κ) relative im
Page 275 and 276: ln(∫ σ γ dE/norm.factor) -4 -6
Page 279 and 280: VALUE = GBFSIG (T,C) GBFSIG calcula
Page 281 and 282: ∆σ σ = { 12 − 15% for T ≤ 1
Page 283: T(MeV) C Pb ∆El 0 ∆E l ∆σ l
Page 287 and 288: 1. sample x from 1 1 ln 1 x c x set
Page 291 and 292: E bind , eV 10 5 10 20 30 40 50 60
Page 301 and 302: eV). At present the values recommen
Page 303 and 304: calculated dE dx − dE measured dx
Page 307 and 308: Figure 43: Stopping powers in Carbo
Page 309 and 310: esults can be seen in table 2 and 3
Page 311 and 312: VALUE = GBRSGM (Z,T,BCUT) Z T BCUT
Page 317 and 318: The functions F i (Z,X,Y) (i = σ,
Page 319 and 320: where, (a − 1) 1 f(ν) = ( ) 1 ν
Page 321 and 322: where: Γ = kα 2π [ ɛ = W E 1 1
Page 323 and 324: NLM DENS CORR IPART (INTEGER) numbe
Page 329 and 330: Hydrogen (bound) Sodium Copper Hydr
Page 331 and 332: [25] M. Gavrila. Relativistic l-she
Page 333 and 334: [78] R.W.Williams. Fundamental Form
Page 335 and 336:
[CONS]). Usually, this is done toge
Page 337 and 338:
Geant 3.16 GEANT User’s Guide TRA
Page 339 and 340:
Page 341 and 342:
VECT,SNEXT,SAFETY Compute distance
Page 343 and 344:
VECT,SNEXT,SAFETY Compute distance
Page 345 and 346:
VECT,SNEXT,SAFETY Compute SFIELD,SM
Page 347 and 348:
SUBROUTINE GUSTEP +SEQ,GCKING,GCVOL
Page 349 and 350:
Page 351 and 352:
Page 353 and 354:
CALL GRKUTA (CHARGE,STEP,VECT,VOUT*
Page 355 and 356:
Geant 3.16 GEANT User’s Guide XIN
Page 357 and 358:
Clicking then on the right button o
Page 359 and 360:
Page 361 and 362:
YMED R “ Center Y coordinate ”
Page 363 and 364:
HIDE is ON, the detector can be exp
Page 365 and 366:
e drawn. NAMNUM contain the arrays
Page 367 and 368:
following levels (red arrows) and o
Page 369 and 370:
Draw a polyline of ’npoint’ poi
Page 371 and 372:
a dynamical 3-D analysis of the sim
Page 373 and 374:
1. position the cursor at (uz0,vz0)
Page 375 and 376:
CHUSET C “User set identifier”
Page 377 and 378:
Set current attribute. 4.8 SSETVA [
Page 379 and 380:
CALL GDCOL(-abs(icol)) 4.12 LWID lw
Page 381 and 382:
5 GEANT/GEOMETRY Geometry commands.
Page 383 and 384:
Print matrixes’ specifications. 5
Page 385 and 386:
6 GEANT/CREATE It creates volumes o
Page 387 and 388:
YES NO 6.7 SCONS name numed inrdw o
Page 389 and 390:
7 GEANT/CONTROL Control commands. 7
Page 391 and 392:
To change lout in /GCUNIT/ Note: un
Page 393 and 394:
IPART I “Particle number” NAPAR
Page 395 and 396:
8.5 CDIR [ chpath chopt ] CHPATH C
Page 397 and 398:
9 GEANT/FZ ZEBRA/FZ commands 9.1 FZ
Page 399 and 400:
Check the structure of one or more
Page 401 and 402:
FACTL R “Scale factor for SL” D
Page 403 and 404:
12 GEANT/PHYSICS Commands to set ph
Page 405 and 406:
Possible IDRAY values are: 0 1 2 To
Page 407 and 408:
12.14 PAIR [ ipair ] IPAIR C “Fla
Page 409 and 410:
PPCUTM R “Cut for e+e- pairs by m
Page 411 and 412:
LGET_15 C “user word” LGET_16 C
Page 413 and 414:
13.6 GEOM [ lgeom˙1 lgeom˙2 lgeom
Page 415 and 416:
LSTAT_7 C “user word” LSTAT_8 C
Page 417 and 418:
Page 419 and 420:
XINT Bibliography [1] R.Brun and P.
Page 421 and 422:
This common contains the threshold
Page 423 and 424:
ITR3D track being scanned (used tog
Page 425 and 426:
NCLAS1 NCLAS2 NCLAS3 IHOLE CGXMIN C
Page 427 and 428:
* PARAMETER (MAXJMP=30) COMMON/GCJU
Page 429 and 430:
GYMAX GZMIN GZMAX GXXXX GYYYY GZZZZ
Page 431 and 432:
DENS density of current material in
Page 433 and 434:
RADDEG radiants to degrees conversi
Page 435 and 436:
PCUTNE parameterisation threshold f
Page 437 and 438:
SMULS SOMULS STMULS DPHYS3 ILABS SL
Page 439 and 440:
NTETA TETMIN TETMAX MODTET number o
Page 441 and 442:
S2 S3 SS1 SS2 SS3 LEP IPORLI ISUBLI
Page 443 and 444:
CFIELD constant for field step eval
Page 445 and 446:
NEXT 1 particle has reached the bou
Page 447 and 448:
C COMMON/GCVOL2/NLEVE2,NAMES2(15),N
Page 449 and 450:
STEP PLIN PLOG BE2 PLASM TRNSMA BOS
Page 451 and 452:
C ESHELL - Shells potentials in eV
Page 453 and 454:
Geant 3.11 GEANT User’s Guide ZZZ
Page 455 and 456:
GFKINE KINE001, KINE100, KINE199 GF
Page 457 and 458:
GPSTAT GEOM700 GPTMED CONS001, CONS
Page 459 and 460:
ZZZZ Bibliography [1] T.Sjöstrand
Page 461 and 462:
CDRSGA, 297 CGMULO, 221 CHANGEWK, 3
Page 463 and 464:
GEANT/DRAWING/KHITS, 372 GEANT/DRAW
Page 465 and 466:
GPART, 13, 46, 61, 293 GPCXYZ, 40,
Page 467:
PATCHY, 278 PCUTS, 399 PDIGI, 389 P
show all

CERN Program Library Long Writeup W5013 - CERNLIB ...

Create successful ePaper yourself

Delete template?

Save as template?