ReadDaqAlone.C

Go to the documentation of this file.

// Macro to read trc files 
// and put in one single trc file
// Ch. Finck, June 09.


#if !defined(__CINT__) || defined(__MAKECINT__)

#include <Riostream.h>
#include "TString.h"
#include "TH1F.h"
#include "TH2F.h"
#include "TFile.h"
#include "TTree.h"

#include "BufferGANIL.h"

#endif

const Int_t MaxValue = 96;
const Int_t dim      = 16384;

UInt_t value[MaxValue];
TTree *oak = new TTree("TreeValues","Parameters tree");
TH1F* histo[MaxValue];

void MakeHisto(UShort_t* event, UShort_t size )
{
   for (Int_t i = 0; i < size/2; ++i) {
          
          UShort_t id = event[i*2];
          value[id]   = event[i*2+1];
          const Char_t* name = BufferGANIL::GetModuleName(id);

          if (histo[id] == 0x0) {
                 histo[id] = new TH1F(Form("h%s",name), name, dim/4, 0, dim/4);
          }
          histo[id]->Fill(value[id]); 

          if (!oak->GetBranch(name)) 
                 oak -> Branch(name, &value[id], Form("%s/i", name));
   }

   oak->Fill();
}

void ReadDaqAlone(TString file = "run_0002.dat")
{
  
   TString rootFileName = file + ".root";   
   
   TFile* fileRoot = new TFile(rootFileName,"recreate");
   TString fileNameIn;  
   UShort_t* event;
   UShort_t size;

   BufferGANIL* bufferGANIL = new BufferGANIL();
 
   fileNameIn  = file;
   
   FILE* in;
   
   Char_t buf[dim];
      
   if ( !(in = fopen (fileNameIn.Data(),"r") ) )  {
      printf("file %s not found\n", fileNameIn.Data());
      return;
   }
   
   Bool_t ok = false;
   // read header   
   while (fread(buf,sizeof(char),dim,in)) {
          memcpy(bufferGANIL->GetBuffer(),buf, dim);  
          ok = bufferGANIL->InitReading();
          if (!ok) continue;
          do {
                 event = bufferGANIL->NextEvent(size);
                 if (event)
                        MakeHisto(event, size);
          }     while (event != 0x0);
          
   } 
   
   fclose(in );
 
   oak->Write();
   
   for (Int_t i = 0; i < MaxValue; ++i) {
        if (histo[i])
           histo[i]->Write();
   }
          
   fileRoot->Close();
}