| 1 | /*----------------------------------------------------------------------
|
|---|
| 2 | jcio.c :
|
|---|
| 3 |
|
|---|
| 4 | Rutinas para la gestion de ficheros en CORSIKA 5.20
|
|---|
| 5 | Author J.C.Gonzalez
|
|---|
| 6 |
|
|---|
| 7 | Las posibilidades que contempla son:
|
|---|
| 8 |
|
|---|
| 9 | - Creacion de un fichero de inicio de run: run######
|
|---|
| 10 | - Creacion de un fichero de fin de run: end######
|
|---|
| 11 | - Creacion de fichero separados para cada cascada, tanto para
|
|---|
| 12 | particulas como para Cherenkov. Ademas, para la salida Cherenkov se
|
|---|
| 13 | genera en ficheros separados en caso de que se reutilice la cascada.
|
|---|
| 14 | Los fichero son por tanto: cer######.# y dat######, donde el indice
|
|---|
| 15 | del fichero cer va de 0 a 9.
|
|---|
| 16 | - La grabacion se realiza en C, por lo que no existen marcas de inicio
|
|---|
| 17 | ni final de bloque en los ficheros.
|
|---|
| 18 |
|
|---|
| 19 | Modified by C.Bigongiari 2000 Dec 19.
|
|---|
| 20 | No more RUN###### and END###### files created.
|
|---|
| 21 | Three files, CER######, DAT###### and STA###### files, per run instead
|
|---|
| 22 | of three files per shower.
|
|---|
| 23 | Now ###### is the RUN NUMBER not the shower number !!!
|
|---|
| 24 | ----------------------------------------------------------------------*/
|
|---|
| 25 |
|
|---|
| 26 | #include <stdio.h>
|
|---|
| 27 | #include <stdlib.h>
|
|---|
| 28 |
|
|---|
| 29 | #define MAXBUF 273 /*- 39*7 -*/
|
|---|
| 30 | #define NSUBBL 21
|
|---|
| 31 | #define NPLONG 1041
|
|---|
| 32 |
|
|---|
| 33 | static int nshow = 0;
|
|---|
| 34 | static float timefirst = 0.;
|
|---|
| 35 | static float timelast = 0.;
|
|---|
| 36 |
|
|---|
| 37 | FILE *cetape;
|
|---|
| 38 | FILE *patape;
|
|---|
| 39 | FILE *sttape;
|
|---|
| 40 | /* CBC
|
|---|
| 41 | FILE *runtape;
|
|---|
| 42 | FILE *endtape;
|
|---|
| 43 | CBC*/
|
|---|
| 44 | char tpl[100];
|
|---|
| 45 | char certpl[100];
|
|---|
| 46 | char dattpl[100];
|
|---|
| 47 | /* CBC char statpl[100]; */
|
|---|
| 48 | char cerfile[100];
|
|---|
| 49 | char datfile[100];
|
|---|
| 50 | /* CBC
|
|---|
| 51 | char stafile[100];
|
|---|
| 52 | char runfile[100];
|
|---|
| 53 | char endfile[100];
|
|---|
| 54 | CBC */
|
|---|
| 55 |
|
|---|
| 56 | /*--------------------------------------------------
|
|---|
| 57 | jcinitio:
|
|---|
| 58 |
|
|---|
| 59 | Define file names
|
|---|
| 60 | --------------------------------------------------*/
|
|---|
| 61 | #ifdef JC_UNDERSCORES
|
|---|
| 62 | void
|
|---|
| 63 | jcinitio_ (char *path, int *runnum)
|
|---|
| 64 | #else /* JC_NO_UNDERSCORES */
|
|---|
| 65 | void
|
|---|
| 66 | jcinitio (char *path, int *runnum)
|
|---|
| 67 | #endif /* JC_UNDERSCORES */
|
|---|
| 68 | {
|
|---|
| 69 | register int i = 0;
|
|---|
| 70 |
|
|---|
| 71 | while (*(path + (++i)) != ' ');
|
|---|
| 72 | strncpy (tpl, path, i);
|
|---|
| 73 | if (*(path + i - 1) != '/')
|
|---|
| 74 | strcat (tpl, "/");
|
|---|
| 75 |
|
|---|
| 76 | strcpy (certpl, tpl);
|
|---|
| 77 | strcat (certpl, "cer%06d");
|
|---|
| 78 | strcpy (dattpl, tpl);
|
|---|
| 79 | strcat (dattpl, "dat%06d");
|
|---|
| 80 | /* CBC
|
|---|
| 81 | strcpy (statpl, tpl);
|
|---|
| 82 | strcat (statpl, "sta%06d");
|
|---|
| 83 | */
|
|---|
| 84 | /* CBC
|
|---|
| 85 | sprintf (runfile, "%srun%06d", tpl, *runnum);
|
|---|
| 86 | sprintf (endfile, "%send%06d", tpl, *runnum);
|
|---|
| 87 | CBC */
|
|---|
| 88 |
|
|---|
| 89 | sprintf (cerfile, certpl, *runnum); /* Added by CB */
|
|---|
| 90 | sprintf (datfile, dattpl, *runnum); /* Added by CB */
|
|---|
| 91 | }
|
|---|
| 92 |
|
|---|
| 93 | /*--------------------------------------------------
|
|---|
| 94 | jcdatsave:
|
|---|
| 95 | graba un bloque al fichero de particulas
|
|---|
| 96 | --------------------------------------------------*/
|
|---|
| 97 | #ifdef JC_UNDERSCORES
|
|---|
| 98 | void
|
|---|
| 99 | jcdatsave_ (float *outbuf)
|
|---|
| 100 | #else /* JC_NO_UNDERSCORES */
|
|---|
| 101 | void
|
|---|
| 102 | jcdatsave (float *outbuf)
|
|---|
| 103 | #endif /* JC_UNDERSCORES */
|
|---|
| 104 | {
|
|---|
| 105 | /* puts("JCIO:: saving particle buffer..."); */
|
|---|
| 106 | fwrite (outbuf, sizeof (float) * MAXBUF * NSUBBL, 1, patape);
|
|---|
| 107 | }
|
|---|
| 108 |
|
|---|
| 109 | /*--------------------------------------------------
|
|---|
| 110 | jccersave:
|
|---|
| 111 | graba un bloque al fichero de particulas
|
|---|
| 112 | --------------------------------------------------*/
|
|---|
| 113 | #ifdef JC_UNDERSCORES
|
|---|
| 114 | void
|
|---|
| 115 | jccersave_ (float *outbuf)
|
|---|
| 116 | #else /* JC_NO_UNDERSCORES */
|
|---|
| 117 | void
|
|---|
| 118 | jccersave (float *outbuf)
|
|---|
| 119 | #endif /* JC_UNDERSCORES */
|
|---|
| 120 | {
|
|---|
| 121 | /* puts("JCIO:: saving cerenkov buffer..."); */
|
|---|
| 122 | fwrite (outbuf, sizeof (float) * MAXBUF * NSUBBL, 1, cetape);
|
|---|
| 123 | }
|
|---|
| 124 |
|
|---|
| 125 | /*--------------------------------------------------
|
|---|
| 126 | jcstartrun:
|
|---|
| 127 | abre el fichero run###### y graba el contenido
|
|---|
| 128 |
|
|---|
| 129 | NOW open CER and DAT files
|
|---|
| 130 | --------------------------------------------------*/
|
|---|
| 131 | #ifdef JC_UNDERSCORES
|
|---|
| 132 | void
|
|---|
| 133 | jcstartrun_ (void)
|
|---|
| 134 | #else /* JC_NO_UNDERSCORES */
|
|---|
| 135 | void
|
|---|
| 136 | jcstartrun (void)
|
|---|
| 137 | #endif /* JC_UNDERSCORES */
|
|---|
| 138 | {
|
|---|
| 139 | /* CBC
|
|---|
| 140 |
|
|---|
| 141 | if ((runtape = fopen (runfile, "wb")) == NULL) {
|
|---|
| 142 | printf ("JCIO:: Cannot open RUN file %s. Exiting.\n", runfile);
|
|---|
| 143 | exit (1);
|
|---|
| 144 | }
|
|---|
| 145 | CBC */
|
|---|
| 146 |
|
|---|
| 147 | /* puts("JCIO:: saving runheader buffer..."); */
|
|---|
| 148 |
|
|---|
| 149 | /* CBC
|
|---|
| 150 | fwrite (runh, sizeof (float) * MAXBUF, 1, runtape);
|
|---|
| 151 |
|
|---|
| 152 | CBC */
|
|---|
| 153 |
|
|---|
| 154 | /* Following lines up to CBC added by CB */
|
|---|
| 155 |
|
|---|
| 156 |
|
|---|
| 157 | if ((cetape = fopen (cerfile, "wb")) == NULL) {
|
|---|
| 158 | printf ("JCIO:: Cannot open CER file %s. Exiting.\n", cerfile);
|
|---|
| 159 | exit (1);
|
|---|
| 160 | }
|
|---|
| 161 | if ((patape = fopen (datfile, "wb")) == NULL) {
|
|---|
| 162 | printf ("JCIO:: Cannot open DAT file %s. Exiting.\n", datfile);
|
|---|
| 163 | exit (1);
|
|---|
| 164 | }
|
|---|
| 165 |
|
|---|
| 166 | /* CBC */
|
|---|
| 167 | }
|
|---|
| 168 |
|
|---|
| 169 | /* WHAT FOLLOWS IS NO MORE USED !!!!!! CB */
|
|---|
| 170 |
|
|---|
| 171 | /*--------------------------------------------------
|
|---|
| 172 | jcstartshower:
|
|---|
| 173 |
|
|---|
| 174 | Function added by CB
|
|---|
| 175 | Sets timefirst and timelast variables to their initial value at the beginnng
|
|---|
| 176 | of new shower
|
|---|
| 177 | --------------------------------------------------*/
|
|---|
| 178 | #ifdef JC_UNDERSCORES
|
|---|
| 179 | void
|
|---|
| 180 | jcstartshower_ (float *evth)
|
|---|
| 181 | #else /* JC_NO_UNDERSCORES */
|
|---|
| 182 | void
|
|---|
| 183 | jcstartshower (float *evth)
|
|---|
| 184 | #endif /* JC_UNDERSCORES */
|
|---|
| 185 | {
|
|---|
| 186 | timefirst = 9.0e10;
|
|---|
| 187 | timelast = -9.0e10;
|
|---|
| 188 | }
|
|---|
| 189 |
|
|---|
| 190 |
|
|---|
| 191 | /*--------------------------------------------------
|
|---|
| 192 | jctime:
|
|---|
| 193 | va echando cuentas del tiempo de los fotones
|
|---|
| 194 | --------------------------------------------------*/
|
|---|
| 195 | #ifdef JC_UNDERSCORES
|
|---|
| 196 | void
|
|---|
| 197 | jctime_ (float *cartim)
|
|---|
| 198 | #else /* JC_NO_UNDERSCORES */
|
|---|
| 199 | void
|
|---|
| 200 | jctime (float *cartim)
|
|---|
| 201 | #endif /* JC_UNDERSCORES */
|
|---|
| 202 | {
|
|---|
| 203 | if (*cartim>timelast) timelast = *cartim;
|
|---|
| 204 | if (*cartim<timefirst) timefirst = *cartim;
|
|---|
| 205 | }
|
|---|
| 206 |
|
|---|
| 207 | /*--------------------------------------------------
|
|---|
| 208 | jcstadata:
|
|---|
| 209 | graba las estadisticas de la cascada
|
|---|
| 210 | --------------------------------------------------*/
|
|---|
| 211 | #ifdef JC_UNDERSCORES
|
|---|
| 212 | void
|
|---|
| 213 | jcstadata_ (float *evth, float *evte,
|
|---|
| 214 | double *nproto, double *nprotb, double *nneutr,
|
|---|
| 215 | double *nneutb, double *nphoto, double *nelect,
|
|---|
| 216 | double *nposit, double *nnu, double *nmum,
|
|---|
| 217 | double *nmup, double *npi0, double *npim, double *npip,
|
|---|
| 218 | double *nk0l, double *nk0s, double *nkmi, double *nkpl,
|
|---|
| 219 | double *nhyp, double *ndeut, double *ntrit, double *nalpha,
|
|---|
| 220 | double *nother, int *ifinnu, int *ifinpi,
|
|---|
| 221 | int *ifinet, int *ifinka, int *ifinhy,
|
|---|
| 222 | double *cerele, double *cerhad,
|
|---|
| 223 | double *plong, int *lpct1, int *nstep, double *thstep)
|
|---|
| 224 | #else /* JC_NO_UNDERSCORES */
|
|---|
| 225 | void
|
|---|
| 226 | jcstadata (float *evth, float *evte,
|
|---|
| 227 | double *nproto, double *nprotb, double *nneutr,
|
|---|
| 228 | double *nneutb, double *nphoto, double *nelect,
|
|---|
| 229 | double *nposit, double *nnu, double *nmum,
|
|---|
| 230 | double *nmup, double *npi0, double *npim, double *npip,
|
|---|
| 231 | double *nk0l, double *nk0s, double *nkmi, double *nkpl,
|
|---|
| 232 | double *nhyp, double *ndeut, double *ntrit, double *nalpha,
|
|---|
| 233 | double *nother, int *ifinnu, int *ifinpi,
|
|---|
| 234 | int *ifinet, int *ifinka, int *ifinhy,
|
|---|
| 235 | double *cerele, double *cerhad,
|
|---|
| 236 | double *plong, int *lpct1, int *nstep, double *thstep)
|
|---|
| 237 | #endif /* JC_UNDERSCORES */
|
|---|
| 238 | {
|
|---|
| 239 | register int i, ii, k;
|
|---|
| 240 | int np[22];
|
|---|
| 241 | float f;
|
|---|
| 242 |
|
|---|
| 243 | fwrite (evth, sizeof (float) * MAXBUF, 1, sttape);
|
|---|
| 244 | fwrite (evte, sizeof (float) * MAXBUF, 1, sttape);
|
|---|
| 245 |
|
|---|
| 246 | fwrite (&timefirst, sizeof (float), 1, sttape);
|
|---|
| 247 | fwrite (&timelast, sizeof (float), 1, sttape);
|
|---|
| 248 |
|
|---|
| 249 | for (i = 0; i < 10; i++) {
|
|---|
| 250 | np[0] = (int) (*(nproto + i));
|
|---|
| 251 | np[1] = (int) (*(nprotb + i));
|
|---|
| 252 | np[2] = (int) (*(nneutr + i));
|
|---|
| 253 | np[3] = (int) (*(nneutb + i));
|
|---|
| 254 | np[4] = (int) (*(nphoto + i));
|
|---|
| 255 | np[5] = (int) (*(nelect + i));
|
|---|
| 256 | np[6] = (int) (*(nposit + i));
|
|---|
| 257 | np[7] = (int) (*(nnu + i));
|
|---|
| 258 | np[8] = (int) (*(nmum + i));
|
|---|
| 259 | np[9] = (int) (*(nmup + i));
|
|---|
| 260 | np[10] = (int) (*(npi0 + i));
|
|---|
| 261 | np[11] = (int) (*(npim + i));
|
|---|
| 262 | np[12] = (int) (*(npip + i));
|
|---|
| 263 | np[13] = (int) (*(nk0l + i));
|
|---|
| 264 | np[14] = (int) (*(nk0s + i));
|
|---|
| 265 | np[15] = (int) (*(nkmi + i));
|
|---|
| 266 | np[16] = (int) (*(nkpl + i));
|
|---|
| 267 | np[17] = (int) (*(nhyp + i));
|
|---|
| 268 | np[18] = (int) (*(ndeut + i));
|
|---|
| 269 | np[19] = (int) (*(ntrit + i));
|
|---|
| 270 | np[20] = (int) (*(nalpha + i));
|
|---|
| 271 | np[21] = (int) (*(nother + i));
|
|---|
| 272 | fwrite (np, sizeof (int) * 22, 1, sttape);
|
|---|
| 273 | }
|
|---|
| 274 |
|
|---|
| 275 | np[0] = (int) (*ifinnu);
|
|---|
| 276 | np[1] = (int) (*ifinpi);
|
|---|
| 277 | np[2] = (int) (*ifinet);
|
|---|
| 278 | np[3] = (int) (*ifinka);
|
|---|
| 279 | np[4] = (int) (*ifinhy);
|
|---|
| 280 | np[5] = (int) (*cerele);
|
|---|
| 281 | np[6] = (int) (*cerhad);
|
|---|
| 282 | fwrite (np, sizeof (int) * 7, 1, sttape);
|
|---|
| 283 |
|
|---|
| 284 | fwrite (lpct1, sizeof (int), 1, sttape);
|
|---|
| 285 | fwrite (nstep, sizeof (int), 1, sttape);
|
|---|
| 286 |
|
|---|
| 287 | f = (float) (*thstep);
|
|---|
| 288 | fwrite (&f, sizeof (float), 1, sttape);
|
|---|
| 289 |
|
|---|
| 290 | for (k=0;k<9;k++)
|
|---|
| 291 | for (i = 0; i < *nstep; i++) {
|
|---|
| 292 | f = (float)(*(plong + i + k*NPLONG));
|
|---|
| 293 | fwrite (&f, sizeof (float), 1, sttape);
|
|---|
| 294 | }
|
|---|
| 295 |
|
|---|
| 296 |
|
|---|
| 297 | }
|
|---|
| 298 |
|
|---|
| 299 | /*--------------------------------------------------
|
|---|
| 300 | jcendrun:
|
|---|
| 301 |
|
|---|
| 302 | Closes STA file
|
|---|
| 303 | --------------------------------------------------*/
|
|---|
| 304 | #ifdef JC_UNDERSCORES
|
|---|
| 305 | void
|
|---|
| 306 | jcendrun_ (float *rune)
|
|---|
| 307 | #else /* JC_NO_UNDERSCORES */
|
|---|
| 308 | void
|
|---|
| 309 | jcendrun (float *rune)
|
|---|
| 310 | #endif /* JC_UNDERSCORES */
|
|---|
| 311 | {
|
|---|
| 312 | /* fwrite (rune, sizeof (float) * MAXBUF, 1, sttape); */
|
|---|
| 313 | fclose (sttape);
|
|---|
| 314 | }
|
|---|
