h61->Fit("gaus")
   h61->GetFunction("gaus")->GetChisquare()
   h61->GetFunction("gaus")->GetNDF()