Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/examples/advanced/stim_pixe_tomography/scripts/Extract_Projection.C

Version: [ ReleaseNotes ] [ 1.0 ] [ 1.1 ] [ 2.0 ] [ 3.0 ] [ 3.1 ] [ 3.2 ] [ 4.0 ] [ 4.0.p1 ] [ 4.0.p2 ] [ 4.1 ] [ 4.1.p1 ] [ 5.0 ] [ 5.0.p1 ] [ 5.1 ] [ 5.1.p1 ] [ 5.2 ] [ 5.2.p1 ] [ 5.2.p2 ] [ 6.0 ] [ 6.0.p1 ] [ 6.1 ] [ 6.2 ] [ 6.2.p1 ] [ 6.2.p2 ] [ 7.0 ] [ 7.0.p1 ] [ 7.1 ] [ 7.1.p1 ] [ 8.0 ] [ 8.0.p1 ] [ 8.1 ] [ 8.1.p1 ] [ 8.1.p2 ] [ 8.2 ] [ 8.2.p1 ] [ 8.3 ] [ 8.3.p1 ] [ 8.3.p2 ] [ 9.0 ] [ 9.0.p1 ] [ 9.0.p2 ] [ 9.1 ] [ 9.1.p1 ] [ 9.1.p2 ] [ 9.1.p3 ] [ 9.2 ] [ 9.2.p1 ] [ 9.2.p2 ] [ 9.2.p3 ] [ 9.2.p4 ] [ 9.3 ] [ 9.3.p1 ] [ 9.3.p2 ] [ 9.4 ] [ 9.4.p1 ] [ 9.4.p2 ] [ 9.4.p3 ] [ 9.4.p4 ] [ 9.5 ] [ 9.5.p1 ] [ 9.5.p2 ] [ 9.6 ] [ 9.6.p1 ] [ 9.6.p2 ] [ 9.6.p3 ] [ 9.6.p4 ] [ 10.0 ] [ 10.0.p1 ] [ 10.0.p2 ] [ 10.0.p3 ] [ 10.0.p4 ] [ 10.1 ] [ 10.1.p1 ] [ 10.1.p2 ] [ 10.1.p3 ] [ 10.2 ] [ 10.2.p1 ] [ 10.2.p2 ] [ 10.2.p3 ] [ 10.3 ] [ 10.3.p1 ] [ 10.3.p2 ] [ 10.3.p3 ] [ 10.4 ] [ 10.4.p1 ] [ 10.4.p2 ] [ 10.4.p3 ] [ 10.5 ] [ 10.5.p1 ] [ 10.6 ] [ 10.6.p1 ] [ 10.6.p2 ] [ 10.6.p3 ] [ 10.7 ] [ 10.7.p1 ] [ 10.7.p2 ] [ 10.7.p3 ] [ 10.7.p4 ] [ 11.0 ] [ 11.0.p1 ] [ 11.0.p2 ] [ 11.0.p3, ] [ 11.0.p4 ] [ 11.1 ] [ 11.1.1 ] [ 11.1.2 ] [ 11.1.3 ] [ 11.2 ] [ 11.2.1 ] [ 11.2.2 ] [ 11.3.0 ] 

  1 //***********************************************************************************************************
  2 // Extract_Slice.C
  3 // Root command file
  4 // Use it by typing in the command line of Root terminal: root Extract_Slice.C
  5 //
  6 //
  7 // More information is available in UserGuide
  8 // Created by Z.LI LP2i Bordeaux 2022
  9 //***********************************************************************************************************
 10 
 11 #include <math.h>
 12 #include <stdint.h>
 13 #include <stdio.h>
 14 #include <string.h>
 15 
 16 #include <vector>
 17 // using namespace std;
 18 
 19 #define PI 3.14159265f
 20 
 21 // Define a structure to read and write each event in the required binary format
 22 struct PixeEvent
 23 {
 24   uint16_t energy_10eV;
 25   uint16_t pixelIndex;
 26   uint16_t sliceIndex;
 27   uint8_t projectionIndex;
 28 };
 29 
 30 // to extract a certain slice or slices
 31 
 32 void Extract_Projection()
 33 {
 34   // FILE *in =fopen("PixeEvent_std_AtCreation.DAT","rb");
 35   FILE* in = fopen("../build/PixeEvent_std_AtExit_Detector135_Aperture70.DAT", "rb");
 36 
 37   // FILE* out = fopen("PixeEvent_std_AtCreation_50Projections.DAT","wb");
 38   FILE* out = fopen("../build/PixeEvent_std_AtExit_Detector135_Aperture70_50Projections.DAT", "wb");
 39 
 40   if (in == NULL) {
 41     printf("error for opening the intput file\n");
 42     return;
 43   }
 44 
 45   PixeEvent p;
 46   PixeEvent pp;
 47   vector<int> valid_projections;
 48   for (int i = 0; i < 50; ++i) {
 49     int p = 2 * i;
 50     valid_projections.push_back(p);
 51   }
 52 
 53   while (fread(&p, 7, 1, in)) {
 54     int key = p.projectionIndex;
 55 
 56     if (std::find(valid_projections.begin(), valid_projections.end(), key)
 57         != valid_projections.end()) {
 58       pp.energy_10eV = p.energy_10eV;
 59       pp.projectionIndex = p.projectionIndex / 2;
 60       pp.sliceIndex = p.sliceIndex;  // index of slices should be reset, starting from 0
 61       pp.pixelIndex = p.pixelIndex;
 62       pp.pixelIndex = p.pixelIndex;
 63       // printf("__ProjectionIndex=%d, SliceIndex=%d, PixelIndex=%d, Energy_10eV=%d\n",
 64       // pp.projectionIndex, pp.sliceIndex, pp.pixelIndex, pp.energy_10eV);
 65       fwrite(&pp, 7, 1, out);
 66     }
 67   }
 68   fclose(in);
 69   fclose(out);
 70 }
 71