#include <iostream>
#include <cassert>

#include "TFile.h"
#include "TH1F.h"

void make_histograms() {

  TH1F* component_000 = new TH1F("component_000","",10,0,10);
  TH1F* component_001 = new TH1F("component_001","",10,0,10);

  component_000->Fill(3);
  component_000->Fill(4);
  component_000->Fill(4);
  component_000->Fill(5);
  component_000->Fill(5);
  component_000->Fill(5);
  component_000->Fill(5);
  component_000->Fill(6);
  component_000->Fill(6);
  component_000->Fill(7);

  component_001->Fill(5);
  component_001->Fill(5);
  component_001->Fill(6);
  component_001->Fill(7);
  component_001->Fill(7);
  component_001->Fill(7);
  component_001->Fill(7);
  component_001->Fill(7);
  component_001->Fill(7);
  component_001->Fill(8);
  for (int i=0; i<100; i++) { // for unequal number of entries
    component_001->Fill(8);
  }

  TH1F* experimental_histogram = new TH1F("experimental_histogram","",10,0,10);
  experimental_histogram->Add(component_000, +1);
  experimental_histogram->Add(component_001, +2);

  // Perform some type of tranformation so that the parameters used by MINUIT will already be normalized to unity ...
  //double experiment_integral = experimental_histogram->Integral();
  //double component_000_integral = component_000->Integral();
  //double component_001_integral = component_001->Integral();
  //double total_simulation_integral = component_000_integral + component_001_integral;
  //std::cout << experiment_integral << std::endl;
  //std::cout << component_000_integral << std::endl;
  //std::cout << component_001_integral << std::endl;
  //std::cout << total_simulation_integral << std::endl;
  
  TFile* fExp = new TFile("experiment_histograms.root","RECREATE");
  experimental_histogram->Write();
  fExp->Close();
  
  TFile* fSim = new TFile("simulation_histograms.root","RECREATE");
  component_000->Write();
  component_001->Write();
  fSim->Close();
  
}