{
gROOT->Reset();


TChain *ch5 =new TChain("h101");
ch5->Add("/home/hanouf/ROOT/geminihf5_cl_1000_fe_100k_Input.root");
TChain *ch4 =new TChain("tree");
ch4->Add("/home/hanouf/ROOT/INCL4clusters_GEMINI_p56Fe_output.root");


TH1F *H1[28];
char name[28];
char title[28];
char func[28];
char select[28];
TH1F *H2[28];
char n2[28];
char t2[28];
char f2[28];
char s2[28];



        
   sprintf(func,"Avv","Zvv==%d");
         printf("func=%s\n",func);

     for (Int_t i=0;i<=28;i++){

         
        sprintf(title,"A{Z=%d}",i);
	sprintf(name,"mass%d",i);
        sprintf(select,"Zvv==%d",i);
	printf("select=%s\n",select);
	sprintf(name,"mass%d",i);
        H1[i] = new TH1F(name,title,60,0,60);
H1[i]->SetLineColor(60);

        H1[i]->GetXaxis()->SetTitle("Avv");

        h101->Project(name,func,select);





   sprintf(f2,"A","Z==%d");
         printf("f2=%s\n",f2);

        sprintf(t2,"A{Z=%d}",i);
	sprintf(n2,"MASS%d",i);
        sprintf(s2,"Z==%d",i);
	printf("s2=%s\n",s2);
	sprintf(n2,"MASS%d",i);
        H2[i] = new TH1F(n2,t2,60,0,60);
        H2[i]->GetXaxis()->SetTitle("A");
H2[i]->SetLineColor(90);
        tree->Project(n2,f2,s2);

}



    TFile *f=new TFile("comparison6.root","recreate");   
 for (Int_t i=0;i<=28;i++){
H1[i]->Write();
H1[i]->Draw();
H2[i]->Write();
H2[i]->Draw("same");
f->Write();
}
gPad->BuildLegend();











   
   f->Close();
 return;

   }