Line | |
---|
1 | #example: |
---|
2 | #setenv INPUT ~/g4/geant4alpha/prototype/particle+matter/processes/hadronic/test/AIX-AFS/KaonMinusAbsorptionAtRest.out |
---|
3 | #awk -fG4KaonMinusA_pionmomentum.awk $INPUT > piminus.vec |
---|
4 | # |
---|
5 | begin{i=0; pim-0; pip=0; piz=0; lam=0; sigm=0; sigp=0; sigz=0} |
---|
6 | { |
---|
7 | if($2="Name") i++ |
---|
8 | if($1!="!") { |
---|
9 | if($1=="pi-") pim++ |
---|
10 | if($1=="pi+") pip++ |
---|
11 | if($1=="pi0") piz++ |
---|
12 | if($1=="lambda") lam++ |
---|
13 | if($1=="sigma-") sigm++ |
---|
14 | if($1=="sigma+") sigp++ |
---|
15 | if($1=="sigma0") sigz++ |
---|
16 | |
---|
17 | px=$2 |
---|
18 | py=$3 |
---|
19 | pz=$4 |
---|
20 | e=$5 |
---|
21 | if($1=="pi-") printf("%6.1f\n", sqrt(px*px+py*py+pz*pz)) |
---|
22 | } |
---|
23 | } |
---|
24 | |
---|
25 | END { |
---|
26 | #printf(" pi- pi+ pi0 lam sigm sigp sigz \n") |
---|
27 | #printf("----------------------------------\n") |
---|
28 | #printf("%5d%5d%5d%5d%5d%5d%5d\n",pim,pip,piz,lam,sigm,sigp,sigz) |
---|
29 | #printf("%5.2f%5.2f%5.2f%5.2f%5.2f%5.2f%5.2f\n",pim/i,pip/i,piz/i,lam/i,sigm/i,sigp/i,sigz/i) |
---|
30 | #for(x=1;x<=j;x++) |
---|
31 | #printf("%6.1f ",ks3[x]) |
---|
32 | #printf("\n") |
---|
33 | #for(x=1;x<=j;x++) |
---|
34 | #printf("%6.1f ",ks4[x]) |
---|
35 | #printf("\n") |
---|
36 | } |
---|
37 | |
---|
Note: See
TracBrowser
for help on using the repository browser.