1 | macro run que=513-2-asd |
---|
2 | * Start the G4PhotoNuclearCrossSectionTest (Is geant4/bin in the PATH? Did you rehash?) |
---|
3 | shell G4PhotoNuclearCrossSectionTest |
---|
4 | * Start the G4ElectroNuclearCrossSectionTest (Is geant4/bin in the PATH? Did you rehash?) |
---|
5 | shell G4ElectroNuclearCrossSectionTest |
---|
6 | pi/del * |
---|
7 | * pln=4 - portret, pln=5 - landscape |
---|
8 | pln = 4 |
---|
9 | * npx - subdivision of the screan in X direction, npy - in Y direction |
---|
10 | npx=2 |
---|
11 | npy=3 |
---|
12 | * define the range of the assembling pictures in memory |
---|
13 | first=1 |
---|
14 | last=6 |
---|
15 | * define the file name for the output |
---|
16 | for/fil 99 photonuclear.ps |
---|
17 | * define the metafile type to assemble the pictures |
---|
18 | ***************meta -99 -[npx][npy][pln] |
---|
19 | meta 99 -112 |
---|
20 | zone 2 3 |
---|
21 | * TITLE NAME ENAME NTH XMIN XMAX YMIN YMAX XTXT YTXT |
---|
22 | exec phot ^1!H h1 h1_e 227 1. 4.e7 1.e-5 1001. 1.e5 1.e-4 |
---|
23 | exec phot ^2!H h2 h2_e 227 1. 4.e7 1.e-5 1001. 1.e5 1.e-4 |
---|
24 | exec phot ^12!C c12 c12_e 227 1. 4.e7 1.e-5 1001. 1.e5 1.e-4 |
---|
25 | exec phot ^27!Al al27 al27_e 227 1. 4.e7 1.e-5 1001. 1.e5 1.e-4 |
---|
26 | exec phot Cu cu cu_e 227 1. 4.e7 1.e-5 1001. 1.e5 1.e-4 |
---|
27 | exec phot Pb pb pb_e 227 1. 4.e7 1.e-5 1001. 1.e5 1.e-4 |
---|
28 | print: |
---|
29 | meta 0 |
---|
30 | for/clo 99 |
---|
31 | mess '*** PhotoNucl.kumac is printing photonuclear.ps file on: '//[que] |
---|
32 | shell xprint -h -P [que] photonuclear.ps |
---|
33 | return |
---|
34 | |
---|
35 | macro phot targ=^12!C file=h1 efile=h1_e nth=227 xmin=1 _ |
---|
36 | xmax=4.e7 ymin=0.03 ymax=700 xt=8. yt=10. |
---|
37 | * |
---|
38 | ****** dir = ./ |
---|
39 | hi/del * |
---|
40 | ve/del * |
---|
41 | ve/cr de(999) R 999*0. |
---|
42 | opt utit |
---|
43 | * Switch to high quality picture |
---|
44 | gr/opt ndat |
---|
45 | gr/set pass 5 |
---|
46 | igset lwid 5. |
---|
47 | gr/set bwid 5. |
---|
48 | gr/set fwid 3. |
---|
49 | gr/set hwid 3. |
---|
50 | gr/set pwid 5. |
---|
51 | * |
---|
52 | title 'E?e! (MeV)' u |
---|
53 | title 'Nuclear electrooabsorbtion cross sections (in mb)' |
---|
54 | opt date |
---|
55 | opt logx |
---|
56 | opt logy |
---|
57 | *** opt liny |
---|
58 | null [xmin] [xmax] [ymin] [ymax] |
---|
59 | * |
---|
60 | ****atitle ' ' '[s]?abs! (mb)' |
---|
61 | title ' ' u |
---|
62 | ****title ' ' |
---|
63 | igset MTYP 20 |
---|
64 | SET TXFP -20 |
---|
65 | SET LTYP 1 |
---|
66 | SET CHHE 1. |
---|
67 | igset txci 1 |
---|
68 | itx [xt] [yt] [targ] |
---|
69 | ve/re p,ps [file].out |
---|
70 | igset PLCI 1 |
---|
71 | set DMOD 2 |
---|
72 | igset LWID 5 |
---|
73 | gra/prim/graph [nth] p ps C |
---|
74 | ve/re e,es [efile].out |
---|
75 | igset PLCI 1 |
---|
76 | set DMOD 1 |
---|
77 | igset LWID 5 |
---|
78 | gra/prim/graph [nth] e es C |
---|
79 | igset LWID 1 |
---|
80 | return |
---|
81 | |
---|
82 | |
---|
83 | |
---|
84 | |
---|
85 | |
---|
86 | |
---|