Prototype2/hcalLab/shTest.C

0001 void shTest(){
0002   Double_t par[] = {0.3164, 3., 1.5};   //  guesses
0003   Double_t points(24);
0004   TF1 * sT = new TF1("sT","0.3164*pow(x,4)*exp(-x*1.5)",0., points);
0005   sT->Draw();
0006   Double_t shapeInt = sT->Integral(0., points);
0007   Double_t shapeMax = sT->GetMaximum(0.,points);
0008   Double_t shapePeak= sT->GetMaximumX(0.,points);
0009   Double_t shapeRMS = sqrt(sT->Moment(2,0.,points)-sT->Moment(1,0.,points)*sT->Moment(1,0.,points));
0010   cout<<"SHAPE  Peak = "<<shapePeak<<"  Max = "<<shapeMax<<"  RMS = "<<shapeRMS<<"  Integral = "<<shapeInt<<endl;
0011   // return par[0]*pow(x1,4.)*(exp(-x1*16./fN));    //  fN is set to (float)NSAMPLES
0012 }
0013
0014
0015 Double_t signalShape(Double_t *x, Double_t * par){
0016   return par[0]*pow((x-par[1]),par[2])*exp(-(x-par[1])*par[3]);
0017 }
0018