void fit_eta2pi0(){
TF1 *f4=new TF1("f4","gaus(0)+gaus(3)+gaus(6)+[9]*TMath::Voigt(x-[10],[11],[12])",0.8,1.6);
 f4->SetParameters(350,0.958,0.01,100,0.95,0.02,150,2.,1.);
f4->SetParameter(9,10);
f4->SetParameter(10,1.285);
 f4->SetParameter(11,0.01); 
f4->SetParameter(12,0.024);
TwoPi0EtaMass_kf->Fit("f4","R");




}