{
// version light de analysis2.C
// pour etude de l'energie reconstruite.

 gROOT->Reset();

 Char_t DATAFILE[50];
 
 Int_t Nfile=14;
 TFile* file[Nfile];
 Double_t Corr[Nfile]={0.,50.,10.,20.,30.,40.,60.,70.,80.,90.,100.,150.,200.,250.};
 
 //Ouverture des fichiers
 for (Int_t i=0;i<Nfile;i++) {
	memset(DATAFILE,0,sizeof(DATAFILE));
	sprintf(DATAFILE,"~/scratch0/RecData/ntuple_Gamma_0000%2d_1_20_NNnoSshapeCorr.root",60+i);
	file[i]=new TFile(DATAFILE);
	cout<<"File: "<<DATAFILE<<" opened "<<endl;
 }

TCanvas c1("c1","c1",5,5,600,400);
c1->Divide(1,2);

TH1F* h1;
Double_t gMean[Nfile];
Double_t gRMS[Nfile];

//loop on files
for(Int_t k=0;k<Nfile;k++) {
 cout<<"k: "<<k<<endl;
 file[k]->cd();

 h1 = new TH1F("h1","h1",100,-20,10); //create new histo
 cout<<"coucou1"<<endl;
 TTree *t1 =(TTree*) gDirectory->Get("CBNT/3333;1");
 if (t1==NULL) cout << " Ntuples 3333 not found." << endl;

 t1->Draw("ecluster[0]-PtGen[0]*cosh(EtaGen[0])>>h1","","goff");
 
 gMean[k]=h1->GetMean();
 gRMS[k]=h1->GetRMS();
 
 cout<<gMean[k]<<" "<<gRMS[k]<<endl;

 h1->Delete();
 t1->Delete();
} //end k loop on file

// create and draw Graphs
TGraph* gr1= new TGraph(Nfile,Corr,gMean);
TGraph* gr2= new TGraph(Nfile,Corr,gRMS);


c1->cd(1);
gr1->SetTitle("Mean");
gr1->SetMarkerStyle(21);
gr1->SetMarkerColor(kBlue);
gr1->Draw("AP");
c1->cd(2);
gr2->SetTitle("RMS");
gr2->SetMarkerStyle(21);
gr2->SetMarkerColor(kBlue);
gr2->Draw("AP");

gStyle->SetOptStat(0);

}