#include "chain_fine.h"
#include "chain_course.h"

void mk_hists(const string &who)
{

	TChain *chain = NULL;
	
	if(who=="fine"){
		chain = chain_fine("reconPhoton");
	}else if(who=="course"){
		chain = chain_course("reconPhoton");
	}

	if(!chain){
		cout<<endl<<"INVALID DATA SET! Use \"fine\" or \"course\"!"<<endl<<endl;
		return;
	}

	// Open File
	char fname[256];
	sprintf(fname, "reconPhoton_%s.root", who.c_str());
	TFile *f = new TFile(fname, "RECREATE");
	
	// Call processor
	chain->Process("bcal_calib_selector.C");
	
	// Flush file and close it cleanly
	f->Write();
	f->Close();
	delete f;
}