#include "TROOT.h"
#include "TFile.h"
#include "TTree.h"

#define FADC_MAX_MODULES 21
#define FADC_MAX_CHAN    16
#define FADC_MAX_SAMPLES 100

class fadc_channel_data : public TObject
{
public:
  // Raw window data
  Int_t nsamples;
  Int_t *raw_data;     //[nsamples]
  Int_t *sample;       //[nsamples]

  fadc_channel_data()  
  {
    nsamples=0; raw_data=0; sample=0;
    raw_data = new Int_t[FADC_MAX_SAMPLES];
    sample   = new Int_t[FADC_MAX_SAMPLES];
  }

  ClassDef(fadc_channel_data,1)

 };

class roc : public TObject
{
public:
  fadc_channel_data fmod[FADC_MAX_MODULES];

  roc ()
  {
  }

  ClassDef(roc,1)
};

void fadc()
{

   TFile f("tree2.root","recreate");
   TTree t2("t2","a Tree with data from a fake Geant3");
   roc *froc = new roc;
   roc &rocref = *froc;
   t2.Branch("roc",&froc);
   t2.Branch("fmod0",&froc->fmod[0]);
   t2.Branch("fmod1",&froc->fmod[1]);

   for(int imod=0; imod<FADC_MAX_MODULES; imod++)
     {
       rocref.fmod[imod].raw_data = new Int_t[FADC_MAX_SAMPLES];
       rocref.fmod[imod].sample = new Int_t[FADC_MAX_SAMPLES];
     }

   for(int event=0; event<100; event++)
     {
       
       for(int isample=0; isample<16; isample++)
	 {
	   rocref.fmod[0].sample[isample]=isample;
	   rocref.fmod[1].sample[isample]=isample;

	   rocref.fmod[0].raw_data[isample]=isample*10;
	   rocref.fmod[1].raw_data[isample]=isample*20;
	 }
       rocref.fmod[0].nsamples=16;
       rocref.fmod[1].nsamples=16;
      
       t2.Fill();
     }
  
   t2.Write();
}