#include <TFile.h>
#include <TString.h>
#include <iostream>
#include <rcpicker.h>

TString xrootd_url("root://nod25.phys.uconn.edu");
TString data_dir("/Gluex/beamline/diamonds/");
TString run_period("chess-5-2013");
TString rootURL(xrootd_url + data_dir + run_period + "/results");

void rockmap(char* file, int scan)
{
   TFile *_file0;
   TString rootfile;
   rootfile.Form("%s/%s_%3.3d_results.root",rootURL.Data(),file,scan);
   _file0 = TFile::Open(rootfile.Data());
   if (! _file0->IsOpen()) {
      std::cerr << "Error - unable to open input root file \"" << rootfile
                << "\"" << std::endl;
   }
   TString rcfile;
   rcfile.Form("%s/%s_%3.3d_rocking_curves.root",rootURL.Data(),file,scan);
   TString logfile;
   logfile.Form("%s/%s/%s",run_period.Data(),file,file);
   rcpicker *rcp = new rcpicker(file, scan, rcfile.Data(), logfile.Data());
   _file0->cd();
   TH2D* hmu;
   _file0->GetObject("hmu",hmu);
   rcp->PickOn(hmu);
}

void cutnoise(double thresh)
{
   std::cout << "cuts go in here" << std::endl;
}