void shTest(){ Double_t par[] = {0.3164, 3., 1.5}; // guesses Double_t points(24); TF1 * sT = new TF1("sT","0.3164*pow(x,4)*exp(-x*1.5)",0., points); sT->Draw(); Double_t shapeInt = sT->Integral(0., points); Double_t shapeMax = sT->GetMaximum(0.,points); Double_t shapePeak= sT->GetMaximumX(0.,points); Double_t shapeRMS = sqrt(sT->Moment(2,0.,points)-sT->Moment(1,0.,points)*sT->Moment(1,0.,points)); cout<<"SHAPE Peak = "<<shapePeak<<" Max = "<<shapeMax<<" RMS = "<<shapeRMS<<" Integral = "<<shapeInt<<endl; // return par[0]*pow(x1,4.)*(exp(-x1*16./fN)); // fN is set to (float)NSAMPLES }