parse_trim_collision.c

   1 /*
   2  * parse srim 2003.26 collision data
   3  *
   4  * author: hackbard@hackdaworld.dyndns.org
   5  *
   6  * this may just work with srim 2003.26 and 180keV c+ -> si,
   7  * not quite sure though! :)
   8  *
   9  */
  10
  11 #define _GNU_SOURCE
  12 #include <stdio.h>
  13 #include <stdlib.h>
  14 #include <string.h>
  15 #include <unistd.h>
  16 #include <sys/types.h>
  17 #include <sys/stat.h>
  18 #include <fcntl.h>
  19
  20 #define MAXION 5000
  21 #define MAXZ 233
  22
  23 typedef struct s_z {
  24         int cell_hits_from_ions;
  25         int total_hits;
  26         double energy;
  27 } t_z;
  28
  29 int main(int argc,char **argv) {
  30
  31         int fd,ion,i,j;
  32         int z,skipped;
  33         int avg1,avg2;
  34         double en;
  35         double max_e;
  36         int max_chfi,max_th;
  37         t_z Z[MAXZ];
  38         unsigned char hit[MAXZ];
  39         char buf[256],*p,value[6],value2[10];
  40
  41         i=0;
  42
  43         if(argc!=2) return 0;
  44
  45         printf("opening file %s ...\n",argv[1]);
  46         if((fd=open(argv[1],O_RDONLY))<0) {
  47                 printf("unable to open file %s\n",argv[1]);
  48                 return 0;
  49         }
  50
  51         ion=1;
  52         skipped=0;
  53         avg1=0;
  54         avg2=0;
  55         for(i=0;i<MAXZ;i++) {
  56                 Z[i].cell_hits_from_ions=0;
  57                 Z[i].total_hits=0;
  58                 Z[i].energy=0;
  59         }
  60         memset(hit,0,MAXZ);
  61
  62         while(1) {
  63
  64                 /* read line */
  65                 i=0;
  66                 memset(buf,0,256);
  67                 while(1) {
  68                         j=read(fd,&buf[i],1);
  69                         if(j<=0) {
  70                                 close(fd);
  71                                 /* norm */
  72                                 max_e=Z[0].energy;
  73                                 max_th=Z[0].total_hits;
  74                                 max_chfi=Z[0].cell_hits_from_ions;
  75                                 for(i=0;i<MAXZ;i++) {
  76                                         if(Z[i].energy>max_e) max_e=Z[i].energy;
  77                                         if(Z[i].total_hits>max_th) max_th=Z[i].total_hits;
  78                                         if(Z[i].cell_hits_from_ions>max_chfi) max_chfi=Z[i].cell_hits_from_ions;
  79                                         avg1+=Z[i].total_hits;
  80                                         avg2+=Z[i].cell_hits_from_ions;
  81                                 }
  82                                 for(i=0;i<MAXZ;i++)
  83                                         printf("%d %f %f %f\n",i*3,1.0*Z[i].total_hits/max_th,
  84                                                                 Z[i].energy/max_e,
  85                                                                 1.0*Z[i].cell_hits_from_ions/max_chfi);
  86                                 printf("skipped = %d\n",skipped);
  87                                 printf("average hits per ion %d / %d\n",avg1/ion,avg2/ion);
  88                                 return 1;
  89                         }
  90                         if(buf[i]=='\n') break;
  91                         i++;
  92                 }
  93
  94                 /* parse line */
  95                 if((buf[0]=='³')&&(buf[1]!='=')) {
  96                         p=strtok(buf,"³");
  97                         value[0]=p[0];
  98                         value[1]=p[1];
  99                         value[2]=p[2];
 100                         value[3]=p[3];
 101                         value[4]=p[4];
 102                         value[5]='\0';
 103                         p=strtok(NULL,"³");
 104                         value2[0]=p[0];
 105                         value2[1]=p[1];
 106                         value2[2]='.';
 107                         value2[3]=p[3];
 108                         value2[4]=p[4];
 109                         value2[5]=p[5];
 110                         value2[6]=p[6];
 111                         value2[7]=p[7];
 112                         value2[8]=p[8];
 113                         value2[9]='\0';
 114                         en=atof(value2);
 115                         p=strtok(NULL,"³");
 116                         value2[0]=p[0];
 117                         value2[1]=p[1];
 118                         value2[2]=p[2];
 119                         value2[3]=p[3];
 120                         value2[4]=p[4];
 121                         value2[5]='.';
 122                         value2[6]=p[6];
 123                         value2[7]=p[7];
 124                         value2[8]=p[8];
 125                         value2[9]=p[9];
 126                         value2[10]='\0';
 127                         z=(int)(atof(value2)/30.);
 128                         if(z>232) skipped+=1;
 129                         else {
 130                                 Z[z].energy+=en;
 131                                 Z[z].total_hits+=1;
 132                                 hit[z]=1;
 133                         }
 134                         if(ion!=atoi(value)) {
 135                                 /* new ion */
 136                                 for(i=0;i<MAXZ;i++) Z[i].cell_hits_from_ions+=hit[i];
 137                                 memset(hit,0,MAXZ);
 138                                 ion=atoi(value);
 139                         }
 140                         // printf("%d %d %f %d\n",z*3,Z[z].total_hits,
 141                         //      Z[z].energy,Z[z].cell_hits_from_ions);
 142                         // return 0;
 143                 }
 144
 145
 146         }
 147         return 0;
 148 }