#include <CobremsGeneration.hh>

#include <iostream>

#include <TROOT.h>
#include <TH1D.h>
#include <TH2D.h>
#include <TF1.h>
#include <TFile.h>

int main() 
{
   double Ebeam = 12.;
   double Eedge = 9.;

   TFile fout("cobrems.root", "recreate");
   CobremsGeneration cobrems(Ebeam, Eedge);

   int nbins = 200;
   TH1D incoh("incoh", "incoherent bremsstrahlung spectrum",
              nbins, 1.2, 12.2);
   for (int i=1; i < nbins; ++i) {
      double x = incoh->GetXaxis()->GetBinCenter(i) / Ebeam;
      incoh->SetBinContent(i, cobrems.Rate_dNidx(x));
   }
   incoh.Write();
   return 0;
}