X-Git-Url: https://hackdaworld.org/gitweb/?p=physik%2Fnlsop.git;a=blobdiff_plain;f=parse_trim_collision.c;h=eed12c57e682779728f5869a31c490230ad76a19;hp=5469bfab49a9d9a5d8b5711269b82b2b33e252c9;hb=b5643e82d7e77ee5e324beb16e9f91008c4dafaf;hpb=f8eab7046ff397059600f9f25181b38c101e6f81

diff --git a/parse_trim_collision.c b/parse_trim_collision.c
index 5469bfa..eed12c5 100644
--- a/parse_trim_collision.c
+++ b/parse_trim_collision.c
@@ -30,6 +30,7 @@ int main(int argc,char **argv) {
 
 	int fd,ion,i,j;
 	int z,skipped;
+	int avg1,avg2;
 	double en;
 	double max_e;
 	int max_chfi,max_th;
@@ -49,6 +50,8 @@ int main(int argc,char **argv) {
 
 	ion=1;
 	skipped=0;
+	avg1=0;
+	avg2=0;
 	for(i=0;i<MAXZ;i++) {
 		Z[i].cell_hits_from_ions=0;
 		Z[i].total_hits=0;
@@ -73,12 +76,15 @@ int main(int argc,char **argv) {
 					if(Z[i].energy>max_e) max_e=Z[i].energy;
 					if(Z[i].total_hits>max_th) max_th=Z[i].total_hits;
 					if(Z[i].cell_hits_from_ions>max_chfi) max_chfi=Z[i].cell_hits_from_ions;
+					avg1+=Z[i].total_hits;
+					avg2+=Z[i].cell_hits_from_ions;
 				}
 				for(i=0;i<MAXZ;i++) 
 					printf("%d %f %f %f\n",i*3,1.0*Z[i].total_hits/max_th,
 								Z[i].energy/max_e,
 								1.0*Z[i].cell_hits_from_ions/max_chfi);
 				printf("skipped = %d\n",skipped);
+				printf("average hits per ion %d / %d\n",avg1/ion,avg2/ion);
 				return 1;
 			}
 			if(buf[i]=='\n') break;