Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/examples/extended/electromagnetic/TestEm2/egs4/egs4.kumac

Version: [ ReleaseNotes ] [ 1.0 ] [ 1.1 ] [ 2.0 ] [ 3.0 ] [ 3.1 ] [ 3.2 ] [ 4.0 ] [ 4.0.p1 ] [ 4.0.p2 ] [ 4.1 ] [ 4.1.p1 ] [ 5.0 ] [ 5.0.p1 ] [ 5.1 ] [ 5.1.p1 ] [ 5.2 ] [ 5.2.p1 ] [ 5.2.p2 ] [ 6.0 ] [ 6.0.p1 ] [ 6.1 ] [ 6.2 ] [ 6.2.p1 ] [ 6.2.p2 ] [ 7.0 ] [ 7.0.p1 ] [ 7.1 ] [ 7.1.p1 ] [ 8.0 ] [ 8.0.p1 ] [ 8.1 ] [ 8.1.p1 ] [ 8.1.p2 ] [ 8.2 ] [ 8.2.p1 ] [ 8.3 ] [ 8.3.p1 ] [ 8.3.p2 ] [ 9.0 ] [ 9.0.p1 ] [ 9.0.p2 ] [ 9.1 ] [ 9.1.p1 ] [ 9.1.p2 ] [ 9.1.p3 ] [ 9.2 ] [ 9.2.p1 ] [ 9.2.p2 ] [ 9.2.p3 ] [ 9.2.p4 ] [ 9.3 ] [ 9.3.p1 ] [ 9.3.p2 ] [ 9.4 ] [ 9.4.p1 ] [ 9.4.p2 ] [ 9.4.p3 ] [ 9.4.p4 ] [ 9.5 ] [ 9.5.p1 ] [ 9.5.p2 ] [ 9.6 ] [ 9.6.p1 ] [ 9.6.p2 ] [ 9.6.p3 ] [ 9.6.p4 ] [ 10.0 ] [ 10.0.p1 ] [ 10.0.p2 ] [ 10.0.p3 ] [ 10.0.p4 ] [ 10.1 ] [ 10.1.p1 ] [ 10.1.p2 ] [ 10.1.p3 ] [ 10.2 ] [ 10.2.p1 ] [ 10.2.p2 ] [ 10.2.p3 ] [ 10.3 ] [ 10.3.p1 ] [ 10.3.p2 ] [ 10.3.p3 ] [ 10.4 ] [ 10.4.p1 ] [ 10.4.p2 ] [ 10.4.p3 ] [ 10.5 ] [ 10.5.p1 ] [ 10.6 ] [ 10.6.p1 ] [ 10.6.p2 ] [ 10.6.p3 ] [ 10.7 ] [ 10.7.p1 ] [ 10.7.p2 ] [ 10.7.p3 ] [ 10.7.p4 ] [ 11.0 ] [ 11.0.p1 ] [ 11.0.p2 ] [ 11.0.p3, ] [ 11.0.p4 ] [ 11.1 ] [ 11.1.1 ] [ 11.1.2 ] [ 11.1.3 ] [ 11.2 ] [ 11.2.1 ] [ 11.2.2 ] [ 11.3.0 ] 

  1 macro egs4.kumac
  2 *
  3 set * ; opt * 
  4 *
  5 filecase keep
  6 
  7 set *fon -40
  8 set mscf 1.5
  9 
 10 set ygti 0.5  ; set gsiz 0.35 ; set ymgu 1.5
 11 set vsiz 0.25 ; set yval 0.3  ; set ylab 0.55
 12 set asiz 0.25 ; set yhti 0.9  ; set xlab 1.2
 13 
 14 title '30 GeV e- on 20 X0 iron'
 15 h/file 1 9.0p1.hbook
 16 
 17 set dmod 1 ; set hcol 4 ; set plci 4
 18 h/pl //lun1/4
 19 key 14 9.0  1 'Geant4-9.0p1' ! l
 20 
 21 exec data
 22  
 23 set dmod 2 ; set hcol 2 ; set plci 2
 24 h/pl 1 s 
 25 key 14 8.0  2 'EGS4' ! l
 26 
 27 atitle t  '(100/E0) (dE/dt)'
 28 exec Redraw-Axis
 29 
 30 close 1
 31 
 32 return
 33 
 34 macro data
 35 *
 36 * e- 30 GeV on 20 X0 Iron
 37 * Particle Data Group. Physics Review D 50-3 (1994)
 38 * longitudinal energy deposition
 39 * (from 0 to 20 X0 by bin of 0.5 X0 --> 40 bins) 
 40 *
 41 * *** read data
 42 nbin=40
 43 vector/create vect(2,[nbin])
 44 vector/read vect pdg.ascii
 45 **vector/write vect ! (1x,f6.1,3x,e12.6) 
 46 *
 47 * *** create and fill histogram       
 48 Emin= 0.
 49 Emax= 20.        
 50 create/1dhisto 1 '20 X0 iron.  e- GeV' [nbin] [Emin] [Emax]
 51 put/contents 1 vect(2)
 52    
 53 **h/pl 1
 54 
 55 * *** save histo
 56 *
 57 **h/file 111 pdg.hbook ! n
 58 **hrout 0
 59  
 60 return 
 61        
 62 Macro Redraw-Axis
 63 set plci 1
 64    X1 = $GRAFINFO('WNXMIN')
 65    X2 = $GRAFINFO('WNXMAX')
 66    Y1 = $GRAFINFO('WNYMIN')
 67    Y2 = $GRAFINFO('WNYMAX')
 68    IS = $GRAFINFO('?FAIS')
 69    IC = $GRAFINFO('?FACI')
 70    Set FAIS 0 ; Set FACI 1
 71    Box [X1] [X2] [Y1] [Y2]
 72    Set FAIS [IS] ; Set FACI [IC]
 73    Tick XL [X1] !
 74 *   Tick YB ! [Y1]
 75    Return
 76