Rev | Line | |
---|
[1199] | 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.