[819] | 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 | * define the metafile type to assemble the pictures |
---|
| 16 | ***************meta -99 -[npx][npy][pln] |
---|
| 17 | meta 99 -112 |
---|
| 18 | **** zone 2 2 |
---|
| 19 | zone 2 2 |
---|
| 20 | * TITLE NAME NTH NEX XMIN XMAX YMIN YMAX XTXT YTXT |
---|
| 21 | exec phot e+^1!H(1GeV) h1_epl 100 326 0. 1.e3 1. 1.e7 5.e2 4.e5 |
---|
| 22 | exec phot e+^1!H(10TeV) h1_eph 100 310 0. 1.e7 1. 1.e7 5.e6 4.e5 |
---|
| 23 | exec phot e+Al(1GeV) cu_epl 100 396 0. 1.e3 1. 1.e7 5.e2 4.e5 |
---|
| 24 | exec phot e+Al(10GeV) cu_eph 100 303 0. 1.e4 1. 1.e7 5.e3 4.e5 |
---|
| 25 | print: |
---|
| 26 | meta 0 |
---|
| 27 | for/clo 99 |
---|
| 28 | mess '*** PhotoNucl.kumac is printing photonuclear.ps file on: '//[que] |
---|
| 29 | shell xprint -h -P [que] photonuclear.ps |
---|
| 30 | return |
---|
| 31 | |
---|
| 32 | macro phot targ=^12!C file=h1 nth=227 nex=326 xmin=1 _ |
---|
| 33 | xmax=4.e7 ymin=0.03 ymax=700 xt=8. yt=10. |
---|
| 34 | * |
---|
| 35 | ****** dir = ./ |
---|
| 36 | hi/del * |
---|
| 37 | ve/del * |
---|
| 38 | ve/cr de(999) R 999*0. |
---|
| 39 | ***** clr |
---|
| 40 | ***** opt zfl |
---|
| 41 | ***** swi gz |
---|
| 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?[g]! (MeV)' u |
---|
| 53 | title 'Energy distribution of the interacting effective photons' |
---|
| 54 | opt date |
---|
| 55 | ** opt logx |
---|
| 56 | opt linx |
---|
| 57 | opt logy |
---|
| 58 | *** opt liny |
---|
| 59 | null [xmin] [xmax] [ymin] [ymax] |
---|
| 60 | * Switch to standard quality picture |
---|
| 61 | * set pass 1 |
---|
| 62 | * igset lwid 1. |
---|
| 63 | * set bwid 1. |
---|
| 64 | * set fwid 1. |
---|
| 65 | * set hwid 1. |
---|
| 66 | * |
---|
| 67 | ****atitle ' ' '[s]?abs! (mb)' |
---|
| 68 | title ' ' u |
---|
| 69 | ****title ' ' |
---|
| 70 | set dmod 3 |
---|
| 71 | igset MTYP 20 |
---|
| 72 | SET TXFP -20 |
---|
| 73 | SET LTYP 1 |
---|
| 74 | * ve/re eg,ex,er [file].dat |
---|
| 75 | * igset PMCI 1 |
---|
| 76 | * igset PLCI 1 |
---|
| 77 | * gra/hp/errors eg ex de er [nex] 20 .1 1 |
---|
| 78 | SET CHHE .5 |
---|
| 79 | igset txci 1 |
---|
| 80 | itx [xt] [yt] [targ] |
---|
| 81 | ve/re e,ap [file].out |
---|
| 82 | igset PLCI 1 |
---|
| 83 | igset LWID 1 |
---|
| 84 | gra/prim/graph [nth] e ap C |
---|
| 85 | igset LWID 1 |
---|
| 86 | return |
---|
| 87 | |
---|
| 88 | |
---|
| 89 | |
---|
| 90 | |
---|
| 91 | |
---|
| 92 | |
---|