---

DispHistogram.cpp

Go to the documentation of this file.

unsigned short GetWord(FILE* file){

  unsigned char  lo_char, hi_char;
  unsigned short word;

  if(file==NULL){
    printf("bad file\n");
    return 0;
  }

  // Read two bytes
  lo_char = fgetc(file);
  hi_char = fgetc(file);

  if(lo_char == EOF || hi_char == EOF) {
    printf("end of file\n");
    return 0;
  }

  // Swap bytes to form 16 bit word
  word = (hi_char<<8)+lo_char;

  return word;
}

unsigned long GetLWord(FILE* file){

  unsigned char  lo_char, hi_char;
  unsigned char  lo_char2, hi_char2;
  unsigned long word;

  if(file==NULL){
    printf("bad file\n");
    return 0;
  }

  // Read four bytes
  lo_char = fgetc(file);
  hi_char = fgetc(file);
  lo_char2 = fgetc(file);
  hi_char2 = fgetc(file);

  if(lo_char == EOF || hi_char == EOF) {
    printf("end of file\n");
    return 0;
  }

  // Swap bytes to form 16 bit word
  word = (hi_char2<<24)+(lo_char2<<16) + (hi_char<<8)+lo_char;

  return word;
}

void DispHistogram(char* infname) { // , float start, float step, int bursts) {

  // Convert ROD histograms to SCTDAQ format

  char htitle[128];
  char hname[32];

  int i,j;
  unsigned short header[8];
  unsigned short trailer[8];
  unsigned long data[1024*2];
  unsigned short chan,chip;

  TH1F* h_scan_evcnt;
  TH2F* h_scan0;
  TH2F* h_scan1;

  FILE* infile;
  TFile* outfile;

  infile = fopen(infname,"rb");

  if(infile==NULL) {
    printf("Could not open file %s\n",infname);
    return;
  }  

  int dontSkip = 0;

  short version = GetWord(infile);          // version
  short length = GetWord(infile);
  dontSkip += 1;

  int run = GetLWord(infile);
  int scan = GetLWord(infile);
  dontSkip += 2;

  char moduleName[16];
  for(int i=0; i<16; i++) {
    moduleName[i] = fgetc(infile);
  }
  dontSkip +=4;

  short scanType = GetWord(infile);           // scanType
  short nBins = GetWord(infile);
  dontSkip ++;

  int size = GetLWord(infile);
  dontSkip ++;
  short type = GetWord(infile);
  short width = GetWord(infile);
  dontSkip ++;

  int bursts = nBins;

  dontSkip +=4;            // Dont skip pointers

  for(int i=0; i<length/4-dontSkip; i++) {
    unsigned long dummy = GetLWord(infile);
  }

  unsigned long pointer = GetLWord(infile);     // Pointer to thresholds
  unsigned long pEvents = GetLWord(infile);
  unsigned long pErrors = GetLWord(infile);
  unsigned long pData = GetLWord(infile);       // Pointer to data

  for(int i=0; i<(pointer-length)/4; i++) {   // Presumably 0
    unsigned long dummy = GetLWord(infile);
  }

  printf("nbins = %d\n", nBins);
  float *threshold = new float[bursts];

  for(int i=0; i<bursts; i++) {
    unsigned long word = GetLWord(infile);
    threshold[i] = *((float*)(&word));
  }

  for(int i=0; i<bursts; i++) {
    unsigned long events = GetLWord(infile);
  }

  for(int i=0; i<bursts; i++) {
    unsigned long errors = GetLWord(infile);
  }

  double hist_start;
  double hist_stop;
  float fist_start;
  float fist_stop;

  //  hist_start = start-(step/2);
  //  hist_stop =  start+(112*step)+(step/2);
  float step = threshold[1] - threshold[0];
  if(threshold[0] < threshold[bursts-1]) {
    hist_start = threshold[0]; // start;
    hist_stop =  threshold[bursts-1] + step; // start+(bursts*step);
  } else {
    hist_start =  threshold[bursts-1] + step; // start+(bursts*step);
    hist_stop = threshold[0]; // start;
  }
  printf("start %3.2f stop %3.2f\n",hist_start, hist_stop);

  fist_start = hist_start;
  fist_stop =  hist_stop;
  printf("start %3.2f stop %3.2f\n",fist_start, fist_stop);

  sprintf(htitle,"ROD Events Requested");
  sprintf(hname,"h_scan_evcnt");
  h_scan_evcnt = new TH1F(hname, htitle, bursts, hist_start, hist_stop);
  h_scan_evcnt->SetFillColor(46);

  sprintf(htitle, "%s Stream 0 THRESHOLD Scan", moduleName);
  sprintf(hname,"h_scan0_0");
  h_scan0 = new TH2F(hname,htitle, 768,0.5,768.5, bursts, hist_start, hist_stop);
  h_scan0->SetXTitle("Channel Number");
  h_scan0->SetYTitle("THRESHOLD (mV)");

  sprintf(htitle, "%s Stream 1 THRESHOLD Scan", moduleName);
  sprintf(hname,"h_scan1_0");
  h_scan1 = new TH2F(hname,htitle, 768,0.5,768.5, bursts, hist_start, hist_stop);
  h_scan1->SetXTitle("Channel Number");
  h_scan1->SetYTitle("THRESHOLD (mV)");

  for(j=0;j<bursts;j++) {
    h_scan_evcnt->Fill(threshold[j],50);  // hardcoded to 50 events
  }

  int slice = 1;
  if(slice) {
    for(j=0;j<bursts;j++){
      int words, skip;
      if(type == 16 || type == 6) {
        words = 768*2; skip = 0;
      } else {
        words = 1024*2; skip = 1024-768;
      }
      for(i=0;i<words;i++){
        // data
        data[i]=GetLWord(infile);

        if(i>0    && i<768)
          h_scan0->Fill(i,threshold[j],data[i]);
        if(i>(768 + skip) && i<(words-skip))
          h_scan1->Fill(i-(768 + skip),threshold[j],data[i]);
      }
//      printf("%s\n", (char *)&data[i-8]);
    }
  } else {
    for(i=0;i<512;i++){

      // For each channel we get 100 hex data words:
      //    
      // Event format
      // Header - 8 16 bit words
      // 00 - Chan Chip
      // 04 - 0000 0000
      // 08 - BBBB BBBB
      // 0C - BBBB BBBB
      //
      // Now the real data
      //
      // Trailer - 8 16 bit words
      // F0 - FFFF FFFF 
      // F4 - FFFF FFFF 
      // F8 - FFFF FFFF 
      // FC - FFFF FFFF

      // header
      for(j=0;j<8;j++){
        header[j]=GetWord(infile);
      }

      // data
      for(j=0;j<bursts;j++){
        data[j]=GetWord(infile);
      }

      // trailer
      for(j=0;j<8;j++){
        trailer[j]=GetWord(infile);
      }

      // Fill SCTDAQ scan Histogram
      chan = (header[1]*128)+header[0];

      for(j=0;j<bursts;j++){
        h_scan0->Fill(chan,threshold[j],data[j]);
      }
    
    } 

  }

  fclose(infile);

  outfile = new TFile("rod.root","RECREATE");

  h_scan0->Write("h_scan0");
  h_scan1->Write("h_scan1");
  h_scan_evcnt->Write("h_scan_evcnt");

  outfile->Close();

  TCanvas *c1 = new TCanvas;
  h_scan_evcnt->Print();
  h_scan0->Print();
  h_scan0->Draw("colz");

  TCanvas *c2 = new TCanvas;
  h_scan1->Print();
  h_scan1->Draw("colz");
}

---