| 1 | macro run que=513-2-asd |
|---|
| 2 | * Start the G4PhotoNuclearCrossSectionTest (Is geant4/bin in the PATH? Did you rehash?) |
|---|
| 3 | shell G4PhotoNuclearCrossSectionTest |
|---|
| 4 | pi/del * |
|---|
| 5 | * pln=4 - portret, pln=5 - landscape |
|---|
| 6 | pln = 4 |
|---|
| 7 | * npx - subdivision of the screan in X direction, npy - in Y direction |
|---|
| 8 | npx=2 |
|---|
| 9 | npy=3 |
|---|
| 10 | * define the range of the assembling pictures in memory |
|---|
| 11 | first=1 |
|---|
| 12 | last=6 |
|---|
| 13 | * define the file name for the output |
|---|
| 14 | ** for/fil 99 photonuclear.ps |
|---|
| 15 | for/fil 99 photonuclear.eps |
|---|
| 16 | * define the metafile type to assemble the pictures |
|---|
| 17 | ***************meta -99 -[npx][npy][pln] |
|---|
| 18 | * this is for the Lanscape ps file |
|---|
| 19 | ** meta 99 -112 |
|---|
| 20 | * this is for the eps file |
|---|
| 21 | meta 99 -113 |
|---|
| 22 | zone 2 3 |
|---|
| 23 | * TITLE NAME NTH NEX XMIN XMAX YMIN YMAX XTXT YTXT |
|---|
| 24 | exec phot ^1!H h1 227 326 150. 4.e7 0.03 0.60 1.e4 0.15 |
|---|
| 25 | exec phot ^2!H h2 227 310 23. 4.e7 0.04 4.00 1.e3 0.5 |
|---|
| 26 | exec phot ^12!C c12 227 396 13. 4.e7 0.10 30.0 2.e3 2.0 |
|---|
| 27 | exec phot ^27!Al al27 227 303 9. 4.e7 0.30 90.0 800. 6. |
|---|
| 28 | exec phot Cu cu 227 144 7. 4.e7 2.00 120. 2.e3 20. |
|---|
| 29 | exec phot Pb pb 227 185 7.4 4.e7 7.00 700. 2.e3 100. |
|---|
| 30 | print: |
|---|
| 31 | meta 0 |
|---|
| 32 | for/clo 99 |
|---|
| 33 | ** mess '*** PhotoNucl.kumac is printing photonuclear.ps file on: '//[que] |
|---|
| 34 | ** shell xprint -h -P [que] photonuclear.ps |
|---|
| 35 | mess '*** PhotoNucl.kumac is printing photonuclear.eps file on: '//[que] |
|---|
| 36 | shell xprint -h -P [que] photonuclear.eps |
|---|
| 37 | return |
|---|
| 38 | |
|---|
| 39 | macro phot targ=^12!C file=h1 nth=227 nex=326 xmin=1 _ |
|---|
| 40 | xmax=4.e7 ymin=0.03 ymax=700 xt=8. yt=10. |
|---|
| 41 | * |
|---|
| 42 | ****** dir = ./ |
|---|
| 43 | hi/del * |
|---|
| 44 | ve/del * |
|---|
| 45 | ve/cr de(999) R 999*0. |
|---|
| 46 | ***** clr |
|---|
| 47 | ***** opt zfl |
|---|
| 48 | ***** swi gz |
|---|
| 49 | opt utit |
|---|
| 50 | * Switch to high quality picture |
|---|
| 51 | gr/opt ndat |
|---|
| 52 | gr/set pass 5 |
|---|
| 53 | igset lwid 5. |
|---|
| 54 | gr/set bwid 5. |
|---|
| 55 | gr/set fwid 3. |
|---|
| 56 | gr/set hwid 3. |
|---|
| 57 | gr/set pwid 5. |
|---|
| 58 | * |
|---|
| 59 | title 'E?[g]! (MeV)' u |
|---|
| 60 | title 'Nuclear photoabsorbtion cross sections (in mb)' |
|---|
| 61 | opt date |
|---|
| 62 | opt logx |
|---|
| 63 | opt logy |
|---|
| 64 | null [xmin] [xmax] [ymin] [ymax] |
|---|
| 65 | * Switch to standard quality picture |
|---|
| 66 | * set pass 1 |
|---|
| 67 | * igset lwid 1. |
|---|
| 68 | * set bwid 1. |
|---|
| 69 | * set fwid 1. |
|---|
| 70 | * set hwid 1. |
|---|
| 71 | * |
|---|
| 72 | ****atitle ' ' '[s]?abs! (mb)' |
|---|
| 73 | title ' ' u |
|---|
| 74 | ****title ' ' |
|---|
| 75 | set dmod 3 |
|---|
| 76 | igset MTYP 20 |
|---|
| 77 | SET TXFP -20 |
|---|
| 78 | SET LTYP 1 |
|---|
| 79 | ve/re eg,ex,er [file].dat |
|---|
| 80 | igset PMCI 1 |
|---|
| 81 | igset PLCI 1 |
|---|
| 82 | gra/hp/errors eg ex de er [nex] 20 .1 1 |
|---|
| 83 | SET CHHE 1. |
|---|
| 84 | igset txci 1 |
|---|
| 85 | itx [xt] [yt] [targ] |
|---|
| 86 | ve/re e,ap [file].out |
|---|
| 87 | igset PLCI 1 |
|---|
| 88 | igset LWID 1 |
|---|
| 89 | gra/prim/graph [nth] e ap C |
|---|
| 90 | igset LWID 1 |
|---|
| 91 | return |
|---|
| 92 | |
|---|
| 93 | |
|---|
| 94 | |
|---|
| 95 | |
|---|
| 96 | |
|---|
| 97 | |
|---|
