2 * parse srim 2003.26 collision data
4 * author: hackbard@hackdaworld.dyndns.org
6 * this may just work with srim 2003.26 and 180keV c+ -> si,
7 * not quite sure though! :)
16 #include <sys/types.h>
24 int cell_hits_from_ions;
29 int main(int argc,char **argv) {
37 unsigned char hit[MAXZ];
38 char buf[256],*p,value[6],value2[10];
44 printf("opening file %s ...\n",argv[1]);
45 if((fd=open(argv[1],O_RDONLY))<0) {
46 printf("unable to open file %s\n",argv[1]);
53 Z[i].cell_hits_from_ions=0;
70 max_th=Z[0].total_hits;
71 max_chfi=Z[0].cell_hits_from_ions;
73 if(Z[i].energy>max_e) max_e=Z[i].energy;
74 if(Z[i].total_hits>max_th) max_th=Z[i].total_hits;
75 if(Z[i].cell_hits_from_ions>max_chfi) max_chfi=Z[i].cell_hits_from_ions;
78 printf("%d %f %f %f\n",i*3,1.0*Z[i].total_hits/max_th,
80 1.0*Z[i].cell_hits_from_ions/max_chfi);
81 printf("skipped = %d\n",skipped);
84 if(buf[i]=='\n') break;
89 if((buf[0]=='³')&&(buf[1]!='=')) {
121 z=(int)(atof(value2)/30.);
122 if(z>232) skipped+=1;
128 if(ion!=atoi(value)) {
130 for(i=0;i<MAXZ;i++) Z[i].cell_hits_from_ions+=hit[i];
134 // printf("%d %d %f %d\n",z*3,Z[z].total_hits,
135 // Z[z].energy,Z[z].cell_hits_from_ions);