list-particles.c

Go to the documentation of this file.

#include <stdlib.h>
#include <stdio.h>

#include <math.h>

#include "corsika-tape.h"
#include "corsika-particle.h"


#define SKIP 4                  /* bytes to skip at beginning of file */

char trash[SKIP];


int
main(int argc, char** argv)
{
  FILE* tape;
  struct {
    char skip1[SKIP];
    block bl;
    char skip2[SKIP];
  } tb;
  int i, j;
  RWORD desc, A, Z;

  if ((tape = fopen(argv[1], "r")) == NULL) {
    fprintf(stderr, "Cannot open file %s. Aborting...\n", argv[1]);
    exit(1);
  }

  while (fread(&tb, sizeof(tb), 1, tape) == 1) {
    for (i = 0; i < SUB_BLOCKS_IN_BLOCK; ++i) {
      if (is_control(tb.bl[i].rh.id)) {
        printf("%4.4s\n", tb.bl[i].rh.id);
      } else {
        for (j = 0; j < PARTICLES_IN_SUB_BLOCK; ++j) {
          desc = tb.bl[i].pd.particle[j].description;
          printf("%8g: ", desc);
          if (is_particle_r(desc)) {
            printf("%-10s", particle[(int)desc/1000].name);
          } else if (is_nucleus_r(desc)) {
            A = (int)desc / 100000;
            Z = (int)desc / 1000 - 100*A;
            printf("N, Z=%g, A=%g", Z, A);
          } else {
            printf("OTHER: %g\n", desc);
          }
          if (desc != 0) {
            printf(", Distance %7.2g (m)", sqrt(tb.bl[i].pd.particle[j].x*tb.bl[i].pd.particle[j].x + tb.bl[i].pd.particle[j].y * tb.bl[i].pd.particle[j].y) / 100.0);
#ifdef THINNING
            printf(", Weight: %g", tb.bl[i].pd.particle[j].weight);
#endif
            printf("\n");
          }
        }
      }
    }
  }

  return 0;
}