// Standard GASP Spy version

#include "/euroball/src/neo/include/neo++.h"

#include "../../online/event_params.h"

//////////
// main //
//////////

void        enable_read(int);
void        enable_other(int);
int          read_file(void);
int         read_calib(int);
int         read_shift(int);
int         read_banana(int);

float        ge_egain,recoil_vc;
float        ge_ecal[40][4],ge_tshift[40];


VectLocGate<4>    Ge_Ener_Gate(4096);    // Ge Gates
LocGate        Ge_Time_Gate(4096);
int           howManySets = 0;
int          low,high;

LocBanana      BGO_Ban(256,256);  // BGO Banana
  
int          do_read_file = 0;

void  enable_read(int signal)
{
  do_read_file = 1;
  printf("Read File at Next event !!\n");
}

void  enable_other(int signal)
{
  exit(0);
}

eb_main {
  if(sigset(SIGUSR1, enable_read) == SIG_ERR) {
    printf(" Error setting SIGUSR1\n");
  }
  if(sigset(SIGUSR2, enable_other) == SIG_ERR) {
    printf(" Error setting SIGUSR2\n");
  }
    Event event(FORMATTED,Spy(1,32768));    ->event type FORMATTED and spy on channel 1 
                                         ->(see spy def or spy use in  gasp_std.cc )
  event.setVerboseOff();

  Hist1D             classif(256,"a_stat_tot");

  Hist1D        Bgo_Mult(4096,"a_bgo_count");
  Hist1D         Bgo_Esum(4096,"a_bgo_ener");
  Hist2D         Mat_Bgo(256,256,"mat_bgo");
  Hist2D         Mat_Bgo_Ban(256,256,"mat_bgo_ban");
  Hist1D        Ger_Mult(4096,"a_ger_count");
  Hist1D         Ger_Esum(4096,"a_ger_ener");

  VectHist1D<40>    Ge_Ener(4096,"a_ge_ener");
  VectHist1D<40>    Ge_Time(4096,"a_ge_time");
  VectHist1D<7>    Ge_Rings(4096,"a_ge_rings");
  
  Hist1D        Ge_E_All(4096,"a_ge_ener_all");
  Hist1D        Ge_T_All(4096,"a_ge_time_all");
  
  VectHist1D<4>    Ge_E_Gate(4096,"a_ge_ener_gate");
  VectHist1D<4>    Ge_E_SGated(4096,"a_ge_gated1");
  VectHist1D<4>    Ge_E_DGated(4096,"a_ge_gated2");
  Hist1D        Ge_T_Gate(4096,"a_ge_time_gate");
  
  Hist1D        BGO_K_Gate(4096,"a_bgo_mult_ban");
  Hist1D        BGO_H_Gate(4096,"a_bgo_ener_ban");
  Hist1D        Ge_BGO_Gated(4096,"a_ge_bgo");
  
  
  int          howManyB, howManyG;
  int          rnd_step,rnd_mod,rnd_val;
  float        rnd_rnd,rnd_num;
  float        doppler[40],costheta[40];
  int          ring[40];
  float        rtmp,geval;
  long        kbval,hbval;
  long        tval,gtval;
  int         i,j,k,m,detKey,npar;
  int          ii,jj,kk;

// Initialisations
  rnd_step   = 359;      // random number increment
  rnd_val    = rnd_step;    // initial value
    rnd_mod    = 1024;         // "random" number modulo.
    rnd_rnd    = 1./1024.;     // 1/rnd_mod..

  for (i=0;i<40;i++) {      // Perform Ge Calibration                     
      for (j=0;j<4;j++)
        ge_ecal[i][j]  = 0.0;  
      ge_ecal[i][1]    = 1.0;  
      ge_tshift[i]    = 0;  
      doppler[i]      = 1.0;  
  }
  ge_egain        = 2.0;
    
  recoil_vc        = 0.0;

  costheta[ 0]  =  0.80902;
  costheta[ 1]  =  0.80902;
  costheta[ 2]  =  0.85065;
  costheta[ 3]  =  0.85065;
  costheta[ 4]  =  0.80902;
  costheta[ 5]  =  0.80902;
  costheta[ 6]  =  0.50000;
  costheta[ 7]  =  0.50000;
  costheta[ 8]  =  0.50000;
  costheta[ 9]  =  0.50000;
  costheta[10]  =  0.52573;
  costheta[11]  =  0.52573;
  costheta[12]  =  0.30902;
  costheta[13]  =  0.30902;
  costheta[14]  =  0.30902;
  costheta[15]  =  0.30902;
  costheta[16]  =  0.00000;
  costheta[17]  =  0.00000;
  costheta[18]  =  0.00000;
  costheta[19]  =  0.00000;
  costheta[20]  =  0.00000;
  costheta[21]  =  0.00000;
  costheta[22]  =  0.00000;
  costheta[23]  =  0.00000;
  costheta[24]  = -0.30902;
  costheta[25]  = -0.30902;
  costheta[26]  = -0.30902;
  costheta[27]  = -0.30902;
  costheta[28]  = -0.50000;
  costheta[29]  = -0.50000;
  costheta[30]  = -0.50000;
  costheta[31]  = -0.50000;
  costheta[32]  = -0.52573;
  costheta[33]  = -0.52573;
  costheta[34]  = -0.80902;
  costheta[35]  = -0.80902;
  costheta[36]  = -0.85065;
  costheta[37]  = -0.85065;
  costheta[38]  = -0.80902;
  costheta[39]  = -0.80902;

  ring[0]  = 0;            // 30 degrees
  ring[1]  = 0;      
  ring[2]  = 0;
  ring[3]  = 0;
  ring[4]  = 0;
  ring[5]  = 0;
  ring[6]  = 1;      // 60 degrees
  ring[7]  = 1;
  ring[8]  = 1;
  ring[9]  = 1;
  ring[10] = 1;
  ring[11] = 1;
  ring[12] = 2;      // 72 degrees
  ring[13] = 2;
  ring[14] = 2;
  ring[15] = 2;
  ring[16] = 3;      // 90 degrees
  ring[17] = 3;
  ring[18] = 3;
  ring[19] = 3;
  ring[20] = 3;
  ring[21] = 3;
  ring[22] = 3;
  ring[23] = 3;
  ring[24] = 4;      // 108 degrees
  ring[25] = 4;
  ring[26] = 4;
  ring[27] = 4;
  ring[28] = 5;      // 120 degrees
  ring[29] = 5;
  ring[30] = 5;
  ring[31] = 5;
  ring[32] = 5;
  ring[33] = 5;
  ring[34] = 6;      //  150 degrees
  ring[35] = 6;
  ring[36] = 6;
  ring[37] = 6;
  ring[38] = 6;
  ring[39] = 6;
  
  for(i=0;i<4;i++)
    Ge_Ener_Gate[i].set_Gate(0,4095); // Gate open to Gamma-Gamma
  Ge_Time_Gate.set_Gate(0,4095);      // Time Gate Open
  
  BGO_Ban.set_Zero();          // Gate Closed No BGO K-H Selection
  
  read_file();      // Read Initial Values from file if it exists
          for(ii=0;ii<256;ii++)
            for(jj=0;jj<256;jj++)
              if(BGO_Ban.Banana(ii,jj))
                for(kk=0;kk<1000;kk++)
                  Mat_Bgo_Ban.inc(ii,jj);
// End Initialisations

  int evNum = 0;
  int flush = 0;
  while(1) {
    try {
      if(do_read_file)  // if read_file flag is set read file with params.
        read_file();
      event.next();
      evNum++;
      classif.inc(60);
        if(flush++ == 5000) { base_histo::flush(); flush = 0;}
      int ErrorSum = 0;
        howManyB= B.howManyFired();
        howManyG= G.howManyFired();
      if(ErrorSum == 0 || ErrorSum == 10) {
          if(event.isBad()) { 
            classif.inc(100);
            for(int kk=0;kk<event.classification.size();kk++) {
              if(event.classification[kk])
                ErrorSum += kk;
            }
          }
          if(howManyG > 0) {            // 1 Ge minimum
            classif.inc(110+howManyB);        // count detectors firing
            classif.inc(120+howManyG);
            classif.inc(150);            // and All events
// BGO BALL
          if(howManyB == 1) {
            Bgo_Mult.inc(B(0).bgo_count);
            Bgo_Esum.inc(B(0).bgo_ener);
            kbval = B(0).bgo_count/16;
            hbval = B(0).bgo_ener/16;
            Mat_Bgo.inc(kbval,hbval);
          }
          printf("Original BGO Values : %d %d\n",(short)B(0).bgo_count,(short)B(0).bgo_ener);
          
          printf("Scaled BGO Values : %d %d\n",kbval,hbval);
          printf("BananaValue : %d\n",BGO_Ban.Banana(kbval,hbval));
          if(BGO_Ban.Banana(kbval,hbval)) {
            BGO_K_Gate.inc(B(0).bgo_count);
            BGO_H_Gate.inc(B(0).bgo_ener);
            for(j=0;j<howManyG;j++)
              Ge_BGO_Gated.inc((long)G(j).ener);
//            Mat_Bgo_Ban.inc(kbval,hbval);
          }

          if(recoil_vc >= 0.0) {    // if v/c has changed
            recoil_vc = recoil_vc*0.01;      // new recoil velocity
            if(recoil_vc > 0.5) {         // a bit too fast
              recoil_vc = 0.;          // disable doppler
              for(i=0;i<40;i++)                      
                doppler[i] = 1./(1.0 + recoil_vc*costheta[i]);  
              recoil_vc = -1.0;          // wait for a new value
            }
          }
// Germani
// Fa tutto qui dentro, Calibrazione, Doppler, spettri e Rings 
          for(i=0;i<G.howManyFired();i++) {  // For each Ge in the event
            detKey = G(i).keyGet();
            geval = G(i).ener;
            if(geval > 0) {          // update pseudo random number
              rnd_val = rnd_val+rnd_step;
              if(rnd_val >= rnd_mod)
                rnd_val = rnd_val-rnd_mod;
              rnd_num = rnd_val*rnd_rnd;
              rtmp = geval + rnd_num + 1.0;  // and apply it
              
              geval = ge_ecal[detKey][0]+
                (ge_ecal[detKey][1]+(ge_ecal[detKey][2]+
                ge_ecal[detKey][3]*rtmp)*rtmp)*rtmp; // apply Calibration
                
              geval *= doppler[detKey]*ge_egain-1.0;   //  ,doppler & gain
              geval = geval < 0 ? 0 : (geval > 4095 ? 4095 : geval);
              G(i).ener = (short) geval;
            }    
            gtval = G(i).time;
            if(gtval > 0) {
              gtval += ge_tshift[detKey];  // Shift time
              tval = gtval < 0 ? 0 : (gtval > 4095 ? 4095 : gtval);
              G(i).time = (short)gtval;
            }                // Increment:
            Ge_Ener[detKey].inc(geval);      // each calibrated Energy
            Ge_E_All.inc(geval);        // the sum of all Ge Energy
            Ge_Rings[ring[detKey]].inc(geval);  // the rings
            Ge_Time[detKey].inc(gtval);      // each shifted Time   
            Ge_T_All.inc(gtval);        // the sum of all Ge Time
          }
          
// This gate could be performed earlier          
          int  geTimeIsInGate;
          if(geTimeIsInGate = Ge_Time_Gate.Gate(gtval)) {  // check Time Gate
            Ge_T_Gate.inc(gtval);        // increment gate window spectrum
//  Ge_Rings[ring[detKey]].inc(geval);  // and time gated rings
          }else{
            G(i).invalidate();
          }
          howManyG = G.howManyFired();  // recalculate remaining Ge's
          if(howManyG > 1) {      // If remains more than 1 Ge     
            for(k=0;k<4;k++) {    // Perform Single and Double gamma-Gate
              for(i=0;i<howManyG;i++) {
                if(Ge_Ener_Gate[k].Gate((long)G(i).ener)) { // if Gate ok
                  Ge_E_Gate[k].inc(G(i).ener);      // increment gate
                  for(j=0;j<howManyG;j++) {
                    if(i != j) {
                      Ge_E_SGated[k].inc(G(j).ener);  // single gated spectra
                      if(Ge_Ener_Gate[k].Gate((long)G(j).ener)) {
                        if(howManyG> 2) {
                          for(m=0;m<howManyG;m++) {
                            if(m != i && m != j) { // double Gated spectra
                              Ge_E_DGated[k].inc(G(m).ener);
                            }
                          }
                        }
                      }
                    }
                  }
                }
              }
            }
          }
        } 
      }
    }catch (RWxmsg & x) { cerr<<x.why()<<endl; }
  }
}

// Subroutines used to change the value of some parameters following a USR Signal

// exit = 1 -> File not found

int  read_file()      // Read the usrAction.txt file descibing what action is requested
{
  FILE    *fd_act;
  char    string[30];
  int    exit = 0;
  
  int    par;
  int    i;
  
  fd_act = fopen("/Gasp_Users/Config/usrAction.txt","r");
  if(fd_act) {
    printf("Reading action file !!\n");
// insert here the read_file code
// add ge_gain and V/C from file
// add ISIS calib and shift
    while(fscanf(fd_act,"%s",string) == 1) {
      if(strstr(string,"Ge_co_cal")) {
        read_calib(0);
        cout<<"Out of Ge_co_cal"<<endl;
      }
      else if(strstr(string,"Ge_df_cal")) {
        read_calib(1);
        cout<<"Out of Ge_df_cal"<<endl;
      }
      else if(strstr(string,"Ge_shift")) {
        read_shift(0);
        cout<<"Out of Ge_shift"<<endl;
      }  
      else if(strstr(string,"Doppler"))  {
        fscanf(fd_act,"%f",&recoil_vc); 
        cout<<"Out of Doppler"<<endl;
      }
      else if(strstr(string,"Gain"))     {
        fscanf(fd_act,"%f",&ge_egain);
        cout<<"Out of Gain"<<endl;
      }
      else if(strstr(string,"BGO_Ban"))  {
        read_banana(0);
        cout<<"Out of Banana"<<endl;
      }
      else if(strstr(string,"Ge_Gate"))  {
        fscanf(fd_act,"%d",&par);
        switch (par) {
          case 0:    // Germanium energy gate
            fscanf(fd_act,"%d",&howManySets);
            if(howManySets > 0 && howManySets < 4) {
              for(i=0;i<howManySets;i++) {
                fscanf(fd_act,"%d %d",&low,&high);
                printf("Ge Energy Gate: Set# %d %d %d \n",i,low,high);
                Ge_Ener_Gate[i].set_Zero();
                Ge_Ener_Gate[i].set_Gate(low,high);
              }
            }
            break;
          case 1:    // Germanium time gate
            fscanf(fd_act,"%d %d",&low,&high);
            printf("Ge Time Gate: %d %d \n",low,high);
            Ge_Time_Gate.set_Zero();
            Ge_Time_Gate.set_Gate(low,high);
            break;
        }
      }
    }
  }else
    exit=1;            // File not found
  do_read_file = 0;
  fclose(fd_act);
  return exit;
}

// exit = 1 -> File not found
// exit = 2 -> Missing Valid header
// exit = 3 -> Index out of range
// exit = 4 -> Too many points

int read_calib(int det)    // Read Energy calibration coefficients from specific files
{
  FILE    *fd_in;
  char    string[101];
  float     a,b;
  int      detPos=-1,exit=0;
  
  switch (det) {
    case 0:    // Germanium
      fd_in = fopen("/Gasp_Users/temp/Calibrations/ge_ener.co60","r");
      break;
    case 1:    // Germanium
      fd_in = fopen("/Gasp_Users/temp/Calibrations/ge_ener.cal","r");
      break;
  }  
  if(fd_in){
    printf("File found\n");
    fgets(string,100,fd_in);
    cout<<string<<endl;
    if(strstr(string,"INIT_DATA")) {
      printf("Correct Header\n");
      fgets(string,100,fd_in);
      while (strlen(string) && !exit && ! strstr(string,"END_DATA")) {
      printf("Correct Header, Entering do loop \n");
        if(sscanf(string,"%2d %f %f",&detPos,&a,&b) == 3) {
          printf("record with 3 values %d %f %f\n",detPos,a,b);
          if(detPos >= 0 && detPos <= 40) {
            printf("Correct Index\n");
            switch (det) {
              case 0:    // Germanium
                ge_ecal[detPos][0] = a;
                ge_ecal[detPos][1] = b;
                ge_ecal[detPos][2] = ge_ecal[detPos][3] = 0;
                break;
              case 1:    // Germanium
                ge_ecal[detPos][0] = a;
                ge_ecal[detPos][1] = b;
                ge_ecal[detPos][2] = ge_ecal[detPos][3] = 0;
                break;
            }
          }else{
            exit = 3;    // Index out of Range
          }
        }
        fgets(string,100,fd_in);
      }
      printf("Exiting do loop \n");
    }else{
      exit = 2;          // Missing Valid Header
    }
    fclose(fd_in);
  }else{
    exit = 1;            // File not found
  }
  return exit;
}

// exit = 1 -> File not found
// exit = 2 -> Missing Valid header
// exit = 3 -> Index out of range
// exit = 4 -> Too many points

int read_shift(int det)    // Read Time calibration coefficients from specific files
{
  FILE    *fd_in;
  char    string[101];
  int     a;
  int      detPos,exit=0;
    
  switch (det) {
    case 0:    // Germanium
      fd_in = fopen("/Gasp_Users/temp/Calibrations/ge_time.cal","r");
      break;
  }  
  if(fd_in) {
//    printf("File found\n");
    fgets(string,100,fd_in);
    if(strstr(string,"INIT_DATA")) {
//      printf("Correct Header\n");
      fgets(string,100,fd_in);
      while (strlen(string) && !exit && ! strstr(string,"END_DATA")) {
//        printf("Correct Header, Entering do loop \n");
        if(sscanf(string,"%d %d ",&detPos,&a) == 2) {
//          printf("record with two values %d %d \n",detPos,a);
          if(detPos >= 0 && detPos <= 40) {
//              printf("Correct Index\n");
            switch (det) {
              case 0:    // Germanium
                ge_tshift[detPos] = a;
                break;
            }
          }else
            exit = 3;  // Index out of Range
        }
        fgets(string,100,fd_in);
      }
    }else{
      exit = 2;          // Missing Valid Header
    }
    fclose(fd_in);
  }else{
    exit = 1;            // File not found
  }
  return exit;
}

// exit = 1 -> File not found
// exit = 2 -> Missing Valid header
// exit = 3 -> Index out of range
// exit = 4 -> Too many points

int read_banana(int type)    // Read Banana poligon 
{
  FILE    *fd_in;
  char    string[101];
  int     a,b,ii;
  long    si_polyX[100],si_polyY[100];
  int      exit=0;
  
  switch (type) {
    case 0:    // BGO
      fd_in = fopen("/Gasp_Users/temp/Calibrations/current/BGO_Ban.ban","r");
      break;
  }  
  if(fd_in){
    printf("File found\n");
    fgets(string,100,fd_in);
    if(strstr(string,"INIT_DATA")) {
      printf("Correct Header\n");
      ii = 0;
      fgets(string,100,fd_in);
      while (strlen(string) && !exit && ! strstr(string,"END_DATA")) {
        printf("Correct Header, Entering do loop \n");
        if(sscanf(string,"%d %d",&a,&b) == 2) {
          printf("record with two values %d %d \n",a,b);
          si_polyX[ii]   = a;
          si_polyY[ii++]  = b;
          if(ii > 99)
            exit = 4;    // Too many points;
        }
        fgets(string,100,fd_in);\
      }
      printf("Exiting do loop \n");
    }else{
      printf("%s\n",string);
      exit = 2;          // Missing Valid Header
    }
    fclose(fd_in);
    if(! exit) {
      switch (type) {      // exit = 5 not a closed Polygon
        case 0:    // BGO
          BGO_Ban.set_Zero();
//          BGO_Ban.set_Banana((long)2,(long)255,(long)2,(long)255);
          cout << "Banana Set" << endl;
          exit = 5*BGO_Ban.set_Banana(ii,si_polyX,si_polyY);
// show the Banana region
          break;
      }
    }
  }else{
    exit = 1;            // File not found
  }
  cout << "Error : " << exit << endl;
  return exit;
}

gaspux4:/Gasp_Users/save

Content: Set of Spectrum Files saved by the XmSpecview Program

gaspux4:/Gasp_Users/temp

Content: Set of Folders with temporary information

gaspux4:/Gasp_Users/temp/OP_Files

Content: Set of Folders with Spectrum Operation descriptors

gaspux4:/Gasp_Users/temp/Calibrations

Content: Set of Files with temporary Calibration, Shift and Gain files.
.60Co extension for the 60 Cobalt gamma calibration files
.cal extension for the generic calibration files

gaspux4:/Gasp_Users/temp/Calibration/current

Content: Set of Files with Banana Poligons and other current related information

Go to Index

XmSpecview.It is an updated version of the one written by Mario DePoli and running under Unix in the HP730.


Menu:

FILE

Add to Display Adds one more spectrum to the display layout

from Histoserver... Reads the spectra from produced by the Neo++ programs.
from Disk...Reads the spectra from the saved folder on the Disk.

Create OP Spectra... Creates an operation spectra. See the Help for more details.
Replace Display Config. Load a display Configuration from the display Folder in the current experiment directory.
Save Level 0 Config. Save the current display Configuration into the display Folder in the current experiment directory.
Save Spectra. Saves the spectra into the Gasp_Users/save Folder according to the format selected by the Menu DEFINE/Disk File format...

Selection.
Screen
All Screens

HARDCOPY (not used)
Single Spectra
Screen
Ps Printer Allows selecting the hardcopy printer.

EDIT

Select All. Selects all spectra of the present screen.
Unselect All. Unselects all spectra of the present screen.
New View Level 1: Selected Opens a new view using only the selected spectra.
Previous View Level * Allows descending one level in the view
Remove Selected. Remove the selected spectra from the view. If the view is the level 0, the spectra are removed from the viewer.
Reposition Selected. Performs rearrangement operations.

Swap Spectra
Move First Spectrum

Regions on Selected. Saves or removes regions of Interest.

SCREENS. Allows selecting one of the 16 available screens.

HISTOG

UPDATE. Updates the spectra according with the following criteria.

Selected
Screen
Level 0 Screen

ERASE. Erases the spectra according with the following criteria.

Selection
Screen
All Screens

DISPLAY

Display Behaviour

Indipendent & Frames. Each spectra has its own frame and scale control bars.
Homogeneous. All spectra have the same scale and common scale control bars.
X-Homogeneous. The spectra have indipendent vertical scales and controls and common horizontal scale and controls.

Display Layout. Allows selecting the spectra layout.

Vertical.
Horizontal.
Row Mayor.
Column Mayor.

Drawing Scales. Allows selecting the proper X-Y scale combination

X Axis Scale

Linear

Y Axis Scale

Linear
Logarithmic
Square Root

X Axis Type

Chan zero Included
Chan zero Excluded

Y Axis Type

0-Max
Min-Max
Simmetric

Drawing Options

Show Overlay. Show the spectra and any associated overlay.
No Clear. Do not clear the existing spectra at the next operation.
Show Work Regions. Show the working regions associated with the spectra
Show Saved Regions. Show the Saved Regions associated with the spectra.

DEFINE

Screen Names... Opens a window allowing to type the screen names.
Peak Width... Allows selecting the aproximate peak width to be used in the peak searching routine.
Good Resolution... Allows to establish what is upper limit of the gamma resolution.
Timer Interval ...Selects the timer interval for repeated operations. If 0 no timer is activated. When de timer is active the DEFINE label appears Red coloured.
Disk File Format ... Selects the disk file format of the saved spectra.

ACTIONS. See more details in the Menu Help.

Show Results. Opens a new window where all calculations are summarised.
Calculate. Performs the selected calculations.

Integral. Uses the I, B and J regions
Gauss Fit. Uses the R and Q regions and the G markers.

Peak Search. Performs a peak searching calculation.
Recal_Cob Performs a global energy linear calibration assuming a 60Co source.
Adjust_Time Performs a shift calibration. Used for Timing alignment.
Eval_Doppler Performs a Doppler estimate.
Adjust_Gain Performs a gain calculation. Used to gain adjust particle spectra.

HELP. Text descrption of the main features of the program.

Go to Index


XmMatview. It is an updated version of the one written by Mario DePoli and running under Unix in the HP730.


Go to Index