#include <TH1D.h>
#include <TGraph.h>
#include <string>

TGraph *deq_ambient(const char *input_file);

TH1D *h1dose(const TH1D *hcounts, double height_cm, double duration_s)
{
   double width_cm = hcounts->GetBinWidth(1);
   double bin_area = width_cm * height_cm;
   TGraph *deq_graph = deq_ambient();

   int nbins = hcounts->GetNbinsX();
   double xlow = hcounts->GetXaxis()->GetXlow();
   double xhigh = hcounts->GetXaxis()->GetXhigh();
   std::string name(hcounts->GetName());
   std::string title(hcounts->GetTitle());
   name += "_dose";
   title += " ambient dose equivalent (mrem/hr)";
   TH1D *hdose = new TH1D(name,title,nbins,xlow,xhigh);

   for (int bin=1; bin <= nbins; ++bin) {
      double count = hcounts->GetBinContent(bin);
      double deqf = deq_graph->Eval();
      hdose->SetBinContent(bin,count
   }
}