// This macro will make a ROOT file containing histograms
// of the charged tracking resolutions suitable for use
// with the src/programs/Simulation/plugins/hdparsim
// plugin for doing parameteric simulations

#include "filelist.h"

//--------------------
// mk_hd_charged_res
//--------------------
void mk_hd_res_charged(void)
{
	gROOT->Reset();
	gStyle->SetPalette(1);

#if 0
	// Create chain of all the input trees
	TChain *trkeff = new TChain("TRACKING/trkeff2","trkeff");
	AddNoTruthFilesToChain(trkeff);
	
	// Open output file to hold histograms
	TFile *f = new TFile("hd_res_charged.root", "RECREATE");
	
	// All of the real work is done in the TSelector derived
	// class in hd_res_charged_selector.cc
	trkeff->Process("hd_res_charged_selector.C");

	// Close resolutions file
	f->Write();
	delete f;
#endif	
	// --- Now do the same thing for the "with truth" case

	// Create chain of all the input trees
	TChain *trkeff2 = new TChain("TRACKING/trkeff2","trkeff");
	AddWithTruthFilesToChain(trkeff2);
	
	// Open output file to hold histograms
	TFile *f = new TFile("hd_res_charged_with_LRtruth.root", "RECREATE");
	
	// All of the real work is done in the TSelector derived
	// class in hd_res_charged_selector.cc
	trkeff2->Process("hd_res_charged_selector.C");

	// Close resolutions file
	f->Write();
	
}