void test(){                                                                                                                                                                                                
  TH1F* h = new TH1F("h","h",3,0,3);                                                                                                                                                                        
  for(int i = 0; i < 3; i++)h->Fill(0.5);                                                                                                                                                                   
  for(int i = 0; i < 4; i++)h->Fill(1.5);                                                                                                                                                                   
  for(int i = 0; i < 3; i++)h->Fill(2.5);                                                                                                                                                                   
                                                                                                                                                                                                            
  //This works                                                                                                                                                                                              
  std::cout << "Likelihood Fit:" << std::endl;                                                                                                                                                              
  h->Fit("pol0","L");                                                                                                                                                                                       
                                                                                                                                                                                                            
  //This dosn't                                                                                                                                                                                             
  h->Sumw2();                                                                                                                                                                                               
  h->Scale(10);                                                                                                                                                                                             
  std::cout << "Weighted Likelihood Fit:" << std::endl;                                                                                                                                                     
  h->Fit("pol0","WL");                                                                                                                                                                                      
}