void fit_waveform05(void)
{
//=========Macro generated from canvas: c1/c1
//=========  (Wed Dec 17 14:49:01 2014) by ROOT version5.34/19
   TCanvas *c1 = new TCanvas("c1", "c1",1505,246,779,672);

gROOT->Reset();
//TTree *Bfield = (TTree *) gROOT->FindObject("Bfield");
gStyle->SetPalette(1,0);
gStyle->SetOptStat(kTRUE);
gStyle->SetOptFit(kTRUE);
gStyle->SetOptFit(1111);
gStyle->SetPadRightMargin(0.15);
gStyle->SetPadLeftMargin(0.15);
gStyle->SetPadBottomMargin(0.15);
gStyle->SetFillColor(0);
//

   c1->Range(-21.79487,-16.8,106.4103,95.2);
   c1->SetFillColor(0);
   c1->SetBorderMode(0);
   c1->SetBorderSize(2);
   c1->SetLeftMargin(0.17);
   c1->SetRightMargin(0.05);
   c1->SetBottomMargin(0.15);
   c1->SetFrameBorderMode(0);
   c1->SetFrameBorderMode(0);
   
   TH1F *hist05 = new TH1F("hist05","Ch=38,10,15",100,0,100);
   hist05->SetBinContent(1,4);
   hist05->SetBinContent(2,3);
   hist05->SetBinContent(3,4);
   hist05->SetBinContent(4,3);
   hist05->SetBinContent(5,5);
   hist05->SetBinContent(6,5);
   hist05->SetBinContent(7,4);
   hist05->SetBinContent(8,4);
   hist05->SetBinContent(9,2);
   hist05->SetBinContent(10,4);
   hist05->SetBinContent(11,1);
   hist05->SetBinContent(12,3);
   hist05->SetBinContent(13,3);
   hist05->SetBinContent(14,2);
   hist05->SetBinContent(15,2);
   hist05->SetBinContent(16,4);
   hist05->SetBinContent(17,2);
   hist05->SetBinContent(18,2);
   hist05->SetBinContent(19,1);
   hist05->SetBinContent(20,3);
   hist05->SetBinContent(21,3);
   hist05->SetBinContent(22,5);
   hist05->SetBinContent(23,3);
   hist05->SetBinContent(24,2);
   hist05->SetBinContent(25,1);
   hist05->SetBinContent(26,2);
   hist05->SetBinContent(27,2);
   hist05->SetBinContent(28,3);
   hist05->SetBinContent(29,3);
   hist05->SetBinContent(30,1);
   hist05->SetBinContent(31,4);
   hist05->SetBinContent(32,5);
   hist05->SetBinContent(33,2);
   hist05->SetBinContent(34,4);
   hist05->SetBinContent(35,4);
   hist05->SetBinContent(36,20);
   hist05->SetBinContent(37,46);
   hist05->SetBinContent(38,51);
   hist05->SetBinContent(39,52);
   hist05->SetBinContent(40,46);
   hist05->SetBinContent(41,45);
   hist05->SetBinContent(42,42);
   hist05->SetBinContent(43,36);
   hist05->SetBinContent(44,34);
   hist05->SetBinContent(45,32);
   hist05->SetBinContent(46,28);
   hist05->SetBinContent(47,28);
   hist05->SetBinContent(48,25);
   hist05->SetBinContent(49,24);
   hist05->SetBinContent(50,21);
   hist05->SetBinContent(51,23);
   hist05->SetBinContent(52,21);
   hist05->SetBinContent(53,20);
   hist05->SetBinContent(54,19);
   hist05->SetBinContent(55,18);
   hist05->SetBinContent(56,18);
   hist05->SetBinContent(57,17);
   hist05->SetBinContent(58,16);
   hist05->SetBinContent(59,13);
   hist05->SetBinContent(60,13);
   hist05->SetBinContent(61,13);
   hist05->SetBinContent(62,14);
   hist05->SetBinContent(63,12);
   hist05->SetBinContent(64,12);
   hist05->SetBinContent(65,11);
   hist05->SetBinContent(66,12);
   hist05->SetBinContent(67,9);
   hist05->SetBinContent(68,8);
   hist05->SetBinContent(69,8);
   hist05->SetBinContent(70,6);
   hist05->SetBinContent(71,6);
   hist05->SetBinContent(72,7);
   hist05->SetBinContent(73,7);
   hist05->SetBinContent(74,7);
   hist05->SetBinContent(75,6);
   hist05->SetBinContent(76,3);
   hist05->SetBinContent(77,5);
   hist05->SetBinContent(78,5);
   hist05->SetBinContent(79,6);
   hist05->SetBinContent(80,5);
   hist05->SetBinContent(81,4);
   hist05->SetBinContent(82,4);
   hist05->SetBinContent(83,4);
   hist05->SetBinContent(84,4);
   hist05->SetBinContent(85,1);
   hist05->SetBinContent(86,3);
   hist05->SetBinContent(87,2);
   hist05->SetBinContent(88,2);
   hist05->SetBinContent(89,3);
   hist05->SetBinContent(90,2);
   hist05->SetBinContent(91,1);
   hist05->SetBinContent(92,5);
   hist05->SetBinContent(93,4);
   hist05->SetBinContent(95,3);
   hist05->SetBinContent(96,3);
   hist05->SetBinContent(97,2);
   hist05->SetBinContent(98,1);
   hist05->SetBinContent(99,3);
   hist05->SetBinContent(100,4);
   hist05->SetMaximum(84);
   hist05->SetEntries(1100);
   // hist05->SetStats(0);
   hist05->SetLineColor(2);
   hist05->GetXaxis()->SetTitle("sample   (4 ns)");
   hist05->GetXaxis()->SetNdivisions(505);
   hist05->GetXaxis()->SetLabelFont(42);
   hist05->GetXaxis()->SetLabelSize(0.07);
   hist05->GetXaxis()->SetTitleSize(0.07);
   hist05->GetXaxis()->SetTitleFont(42);
   hist05->GetYaxis()->SetTitle("ADC value");
   hist05->GetYaxis()->SetNdivisions(505);
   hist05->GetYaxis()->SetLabelFont(42);
   hist05->GetYaxis()->SetLabelSize(0.07);
   hist05->GetYaxis()->SetTitleSize(0.07);
   hist05->GetYaxis()->SetTitleOffset(1.2);
   hist05->GetYaxis()->SetTitleFont(42);
   hist05->GetZaxis()->SetLabelFont(42);
   hist05->GetZaxis()->SetLabelSize(0.07);
   hist05->GetZaxis()->SetTitleSize(0.07);
   hist05->GetZaxis()->SetTitleFont(42);
   hist05->Draw("");
   
   TPaveText *pt = new TPaveText(0,0.9018323,0.4148387,0.995,"blNDC");
   pt->SetName("title");
   pt->SetBorderSize(0);
   pt->SetFillColor(0);
   pt->SetFillStyle(0);
   pt->SetTextFont(42);
   TText *text = pt->AddText("Ch=38,10,15");
   pt->Draw();
   c1->Modified();
   c1->cd();
   c1->SetSelected(c1);

   // define pulse function

   Double_t xmin=0;
   Double_t xmax=100;

   TF1 *pulse = new TF1("pulse_func",pulse_func,xmin,xmax,5);
   Double_t mu=1/8.;    // unis of 1/x
   Double_t sigma=5;     // units of x 
   Double_t gain=500;   
   Double_t t0=40;
   Double_t pedestal=2;
   pulse->SetParameters(mu,sigma,gain,t0,pedestal);
   pulse->SetParNames("mu","sigma","gain","t0","pedestal");
  
   hist05->Fit(pulse);
   c1->SaveAs("fit_waveform05.pdf");

}
Double_t pulse_func (Double_t *x, Double_t *par) {

   Double_t mu=par[0];
   Double_t sigma=par[1];
   Double_t gain=par[2];
   Double_t t0=par[3];
   Double_t pedestal=par[4];
   Double_t x1=x[0]-t0;

   char string[256];

   Double_t func;
   Double_t amplitude;

   if (x1 < -20) {
           func = pedestal;
           }
   else if {

         Double_t arg = (mu*sigma*sigma-x1)/(sqrt(2)*sigma);
         Double_t arg2 = -mu*x1 +(mu*sigma)*(mu*sigma)/2;
         amplitude = (mu/2)* exp(arg2) *  TMath::Erfc(arg);
          func = gain * amplitude + pedestal;
  }
   
   /*sprintf (string,"x1=%f amplitude=%f mu=%f sigma=%f t0=%f pedestal=%f, func=%f\n",x1,amplitude,mu,sigma,pedestal,func);
   printf ("string=%s",string);*/

   return func;
}