void fit_etapippim_3pi(){
TF1 *f3=new TF1("f3","gaus(0)+gaus(3)+gaus(6)+[9]*TMath::Voigt(x-[10],[11],[12])+[13]*TMath::Voigt(x-[14],[15],[16])+gaus(17)+gaus(20)",0.8,2.0);
 f3->SetParameters(350,0.958,0.01,100,0.95,0.02,150,2.,1.);
f3->SetParameter(9,21);
f3->SetParameter(10,1.28);
 f3->SetParameter(11,0.01); 
f3->SetParameter(12,0.024);
f3->SetParameter(13,20);
f3->SetParameter(14,1.5);
 f3->SetParameter(15,0.01); 
f3->SetParameter(16,0.04);
f3->SetParameter(17,600);
f3->SetParameter(18,2.);
f3->SetParameter(19,0.5);
f3->SetParameter(20,200);
f3->SetParameter(21,1.3);
f3->SetParameter(22,0.05);
PipPimEtaMass_Eta_to_3pi->Fit("f3","R");




}