#include "TH1F.h"
#include "TMath.h"
#include <vector> 

void CopiaHisto2(){
    
    TFile *f = new TFile("output.root","RECREATE");

    vector< TH1F* > sample; 

    Int_t nsamples=0;
    for(Int_t i=1;i<=31; i++){

	if(i==4) continue;
	if(i==19) continue;
	if(i==21) continue;
	//if(i==27) continue;
	TString istring=TString::Itoa(i,10);
	TString filename="MassHistogram_mclasseventv2_D0_"+istring+".root"; 
	TFile *infile=new TFile(filename);
	TH1F *hnew = (TH1F*)hist->Clone("hnew");
	sample.push_back(hnew);
	nsamples++;
    }

//

    TH1F *hns=new TH1F("hns","Numero di samples",30,0,31);

    sample[0]->SetTitle("");
    sample[0]->SetXTitle("Massa K-#pi [MeV/c^{2}]");
    sample[0]->Draw();
    c1->Print("plot"+TString::Itoa(0,10)+".gif");
    Int_t N=sample[0]->GetEntries();

    cout <<"Numero eventi sample "<< "0" <<" " << sample[0]->GetEntries() << endl;

    hns->Fill(sample[0]->GetEntries());
        
    for(Int_t j=1;j<nsamples;j++){

	/*
	if(j==4) continue;
	if(j==19) continue;
	if(j==21) continue;
	if(j==27) continue;
	*/

	cout << "Numero eventi sample " << j << " " << sample[j]->GetEntries() << endl;
	hns->Fill(sample[j]->GetEntries());
	
	sample[0]->Add(sample[j]);
	sample[0]->Draw("e");

	Float_t max=sample[0]->GetMaximum();
	N=sample[0]->GetEntries();

	cout << max << endl;

	TPaveText *pt = new TPaveText(1880,0.98*(max+sqrt(max)),1925,0.75*(max+sqrt(max)));
//	pt->SetFillColor(18);
	pt->SetTextFont(62); //20 times new roman
	pt->SetTextColor(2);
	pt->AddText("Numero di Candidati ");
	pt->AddText(TString::Itoa(N,10));
	pt->Draw();

	c1->Update();
        c1->Print("plot"+TString::Itoa(j,10)+".gif");
	
	c1->WaitPrimitive();	
	
//	gets(s);
//	c1->Clear();
    }

    hns->Draw();
    c1->Print("EventsPerSample.gif");

}