/**
 * \file draw_analysis.C
 * \author Firstname Lastname <mail@mail.com>
 * \date 2012
 */

#include "TH1.h"
#include "TFile.h"
#include "TString.h"

// Jpsi file name and file
const TString jpsiFileName = TString(gSystem->Getenv("VMCWORKDIR")) + TString("/macro/tutorial/reco/events_jpsi/analysis.0000.root");
const TFile* jpsiFile = new TFile(jpsiFileName);
// UrQMD file name and file
const TString urqmdFileName = TString(gSystem->Getenv("VMCWORKDIR")) + TString("/macro/tutorial/reco/events_urqmd/analysis.0000.root");
const TFile* urqmdFile = new TFile(urqmdFileName);

void draw_analysis()
{
	TCanvas* canvas = new TCanvas("draw_analysis","draw_analysis", 600, 600);
	canvas->SetGrid();

	// Get the histograms from file
	TH1* jpsiHist = (TH1*) jpsiFile->Get("hStsParamFirstMom");
	TH1* urqmdHist = (TH1*) urqmdFile->Get("hStsParamFirstMom");

	// Set color for histograms
	jpsiHist->SetLineColor(kRed);
	urqmdHist->SetLineColor(kBlue);

	// Draw both histograms on the same canvas
	jpsiHist->Draw();
	urqmdHist->Draw("same");
}