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