// Sara Vanini: main program muon radiography for EM reconstruction
// June, 10th 2008

#include <stdio.h>
#include <iostream.h>
#include <math.h>
#include <vector.h>
#include <stdlib.h>
#include <fstream.h>
#include <string.h>
using std::string;  

// root classes
#include "TH1F.h"
#include "TFile.h"
#include "TRandom.h"

// algorithm parameters
int iterNum = 10;		// number of iterations
bool dump_root = true;

// global variables
TFile * fout;

// useful histograms
TH1F * hlambda[10]; 

int main()
{
  fout = new TFile("EM.root","recreate");

  // init histograms
  for(int i=0; i<iterNum; i++)
  {
        string name = "lambda_";
        name += i +'0';	
cout << name<<endl;
  	hlambda[i] = new TH1F(name.c_str(),"lambda values",1000,0.,1.);
	hlambda[i]->SetDirectory(fout); 
  }

  // iterations  
  TRandom * rand = new TRandom; 
  for(int it=0; it<iterNum; it++) 
  {
	for(int m=0; m<100; m++)
	{
		float lambda_j_new = rand->Rndm();
		hlambda[it]->Fill(lambda_j_new);
	}
  }// end iterations
 
  if(dump_root)
  {
  	// *** Draw histograms just for sanity check and save them to file
  	for(int i=0; i<iterNum; i++)
  	{
  		hlambda[i]->Write();
  		fout->Close();
  	}
  } 

  return 0;
}