Back to home page

EIC code displayed by LXR

 
 

    


Warning, /geant4/examples/extended/electromagnetic/TestEm11/sandia/sandia.kumac is written in an unsupported language. File is not indexed.

0001 macro sandia.kumac
0002 *
0003 set * ; opt * 
0004 *
0005 filecase keep
0006 
0007 set *fon -40
0008 set mscf 1.5
0009 
0010 set ygti 0.5  ; set gsiz 0.35 ; set ymgu 1.5
0011 set vsiz 0.25 ; set yval 0.3  ; set ylab 0.55
0012 set asiz 0.25 ; set yhti 0.9  ; set xlab 1.2
0013 * 
0014 title 'Depth dose distribution of 1033 keV e- in Al'
0015 *
0016 h/file 1 al.1033kev.hbook
0017 *
0018 h/list
0019 *
0020 opt liny
0021 *
0022 set dmod 1 ; set hcol 1 ; set plci 1
0023 h/pl //lun1/8
0024 key 0.6 3.0 1 'G4 9.1-cand-01' ! l
0025 atitle 'x/r0' 'MeV*cm**2/g'
0026 
0027 * *** data
0028 * G.J.Lockwood et al.
0029 *     Sandia report SAND79-0414.UC-34a, February 1987
0030 * O.Kadri et al. NIM B 258 (2007) 381
0031 *
0032 fname = 'data/Al_1033keV.ascii'
0033 vector/create nb(1)
0034 vector/read nb [fname] ! 'O'
0035 nbin = nb(1)
0036 vector/create vect(2,[nbin])
0037 vector/read vect [fname] ! 'C'
0038 ***vector/write vect ! (1x,f6.4,3x,f5.3)
0039 
0040 set mtyp 22
0041 set pmci 2
0042 graph [nbin] vect(1) vect(2)  P
0043 
0044 vector/create xmark(1)
0045 vector/create mark(1)
0046 vector/input xmark 0.57
0047 vector/input mark  2.85
0048 graph 1 xmark mark P
0049 itx   0.64 2.8 'Sandia data'
0050                                    
0051 exec Redraw-Axis
0052 
0053 return
0054 
0055 Macro Redraw-Axis
0056 set plci 1
0057    X1 = $GRAFINFO('WNXMIN')
0058    X2 = $GRAFINFO('WNXMAX')
0059    Y1 = $GRAFINFO('WNYMIN')
0060    Y2 = $GRAFINFO('WNYMAX')
0061    IS = $GRAFINFO('?FAIS')
0062    IC = $GRAFINFO('?FACI')
0063    Set FAIS 0 ; Set FACI 1
0064    Box [X1] [X2] [Y1] [Y2]
0065    Set FAIS [IS] ; Set FACI [IC]
0066    Tick XL [X1] !
0067 *   Tick YB ! [Y1]
0068    Return
0069 
0070