{    

  gROOT->Reset();

  gROOT->ProcessLine(".L fit.C"); 

  
  TH1F* h1= new TH1F("h1","m(dimu)",90,1.9,4.1);
  
  t101->Draw("dimu.mass>>h1");
 
/*
 TH1F* h2= new TH1F("h2","m(mu mu) pt>5",50,2,4);
  
  t101->Draw("dimu.mass>>h2","dimu.pt>5&&dimu.pt<=6");

*/


  TF1* g2p = new TF1("g2p",gauss_gauss_poly, 2.55,3.5,11);
  
  g2p->SetParameters(2.55,5000,-8000,0,0,4000,3.069,0.09,0.15,0,2);
  
  g2p->FixParameter(0,2.55);
  // g2p->FixParameter(3,0);
  g2p->FixParameter(4,0);

  //  g1p->FixParameter(6,3.069);
  //  g1p->FixParameter(7,0.09);
  
  g2p->FixParameter(9,0);

  h1->Fit("g2p","R");
 
  Double_t *para0 = g2p->GetParameters(); 
  
    TF1* g2p = new TF1("g2p",gauss_gauss_poly, 2.55,3.5,11);
  
  g2p->SetParameters(2.55,5000,-8000,0,0,4000,3.069,0.09,0.15,0,2);
  
  g2p->FixParameter(0,2.55);
  // g2p->FixParameter(3,0);
  g2p->FixParameter(4,0);

  //  g1p->FixParameter(6,3.069);
  //  g1p->FixParameter(7,0.09);
  
  g2p->FixParameter(9,0);

  h1->Fit("g2p","R");
 
  Double_t *para0 = g2p->GetParameters(); 
  
 
 /*
  TF1* g1p = new TF1("g1p",gauss_poly, 2.55,3.5,8);
  
  g1p->SetParameters(2.55,5000,-8000,0,0,4000,3.069,0.09);

  h1->Fit("g1p","R");
 
  Double_t *para0 = g1p->GetParameters(); 
 */ 

 /* 

  TF1* g1e = new TF1("g1e",gauss_exp, 2.55,3.5,6);
  
  g1e->SetParameters(2.5,5000,0.1,4000,3.069,0.1);

  h1->Fit("g1e","R");
 
  Double_t *para0 = g1e->GetParameters(); 
 
 
  
  TF1* g1l = new TF1("g1l",gauss_linear, 2.55,3.5,5);
  
  g1l->SetParameters(6000,-0.1,4000,3.069,0.1);

  h1->Fit("g1l","R");
 
  Double_t *para0 = g1l->GetParameters(); 
  
 */  

 }