void MakePlotsPath(const char* data_path, const char* mc_path){

	printf("data = %s, mc = %s\n",data_path,mc_path);

	gROOT->SetStyle("Plain");
	gStyle->SetPalette(1);
	gStyle->SetOptStat(0);
	gStyle->SetPalette(1);
	gStyle->SetCanvasColor(0);
	gStyle->SetFrameFillColor(0);
	gStyle->SetOptTitle(0);
	
	TString mc_fileName(Form("%s/histograms.root",mc_path));
	TString data_fileName(Form("%s/histograms.root",data_path));
	TFile* mc_file = new TFile(mc_fileName.Data());
	TFile* data_file = new TFile(data_fileName.Data());

	TString mc_fileNameEvStat(Form("%s/event_stat.root",mc_path));
	TString data_fileNameEvStat(Form("%s/event_stat.root",data_path));
	TFile* mc_fileEvSt = new TFile(mc_fileNameEvStat.Data());
	TFile* data_fileEvSt = new TFile(data_fileNameEvStat.Data());
	TDirectoryFile* mc_dir = (TDirectoryFile*)mc_fileEvSt->Get("trigger_histograms_");

	// WARNING!!!! check the trigger classes for your run!!!!
	TDirectoryFile* data_dirCINT1B = (TDirectoryFile*)data_fileEvSt->Get("trigger_histograms_+CINT1B-ABCE-NOPF-ALL");
	TDirectoryFile* data_dirCINT1A = (TDirectoryFile*)data_fileEvSt->Get("trigger_histograms_+CINT1A-ABCE-NOPF-ALL");
	TDirectoryFile* data_dirCINT1C = (TDirectoryFile*)data_fileEvSt->Get("trigger_histograms_+CINT1C-ABCE-NOPF-ALL");
	TDirectoryFile* data_dirCINT1E = (TDirectoryFile*)data_fileEvSt->Get("trigger_histograms_+CINT1-E-NOPF-ALL");

	// mc
	TH2D* mc_etaPhi = (TH2D*)mc_file->Get("fEtaPhi");
	TH1D* mc_eta = (TH1D*)mc_etaPhi->ProjectionX();
	TH1D* mc_phi = (TH1D*)mc_etaPhi->ProjectionY();
	TH1D* mc_phiEta1 = (TH1D*)mc_etaPhi->ProjectionY("mc_phiEta1",mc_eta->FindBin(-1.4),mc_eta->FindBin(1.4)-1);
	TH1D* mc_phiEtaNeg = (TH1D*)mc_etaPhi->ProjectionY("mc_phiEtaNeg",mc_eta->FindBin(-1.4),mc_eta->FindBin(0)-1);
	TH1D* mc_phiEtaPos = (TH1D*)mc_etaPhi->ProjectionY("mc_phiEtaPos",mc_eta->FindBin(0),mc_eta->FindBin(1.4)-1);
	TH1D* mc_deltaPhi = (TH1D*)mc_file->Get("fDeltaPhi");
	TH1D* mc_deltaTheta = (TH1D*)mc_file->Get("fDeltaTheta");
	TH1D* mc_vtxX = (TH1D*)mc_file->Get("fVtxX");
	TH1D* mc_vtxY = (TH1D*)mc_file->Get("fVtxY");
	TH1D* mc_vtxZ = (TH1D*)mc_file->Get("fVtxZ");
	TH3D* mc_vtxContributors = (TH3D*)mc_file->Get("fVtxContributors");
	TH1D* mc_vtxXContributors = (TH1D*)mc_vtxContributors->ProjectionX();
	TH1D* mc_vtxYContributors = (TH1D*)mc_vtxContributors->ProjectionY();
	TH1D* mc_vtxZContributors = (TH1D*)mc_vtxContributors->ProjectionZ();
	TH1F* mc_HistFiredBitsSPD = (TH1F*)mc_dir->Get("fHistFiredBitsSPD");
	mc_etaPhi->SetName("mc_etaPhi");
	mc_eta->SetName("mc_eta");
	mc_phi->SetName("mc_phi");
	mc_deltaPhi->SetName("mc_deltaPhi");
	mc_deltaTheta->SetName("mc_deltaTheta");
	mc_vtxX->SetName("mc_vtxX");
	mc_vtxY->SetName("mc_vtxY");
	mc_vtxZ->SetName("mc_vtxZ");
	mc_vtxXContributors->SetName("mc_vtxXContributors");
	mc_vtxYContributors->SetName("mc_vtxYContributors");
	mc_vtxZContributors->SetName("mc_vtxZContributors");

	// data
	TH2D* data_etaPhi = (TH2D*)data_file->Get("fEtaPhi");
	TH1D* data_eta = (TH1D*)data_etaPhi->ProjectionX();
	TH1D* data_phi = (TH1D*)data_etaPhi->ProjectionY();
	TH1D* data_phiEta1 = (TH1D*)data_etaPhi->ProjectionY("data_phiEta1",data_eta->FindBin(-1.4),data_eta->FindBin(1.4)-1);
	TH1D* data_phiEtaNeg = (TH1D*)data_etaPhi->ProjectionY("data_phiEtaNeg",data_eta->FindBin(-1.4),data_eta->FindBin(0)-1);
	TH1D* data_phiEtaPos = (TH1D*)data_etaPhi->ProjectionY("data_phiEtaPos",data_eta->FindBin(0),data_eta->FindBin(1.4)-1);
	TH1D* data_deltaPhi = (TH1D*)data_file->Get("fDeltaPhi");
	TH1D* data_deltaTheta = (TH1D*)data_file->Get("fDeltaTheta");
	TH1D* data_vtxX = (TH1D*)data_file->Get("fVtxX");
	TH1D* data_vtxY = (TH1D*)data_file->Get("fVtxY");
	TH1D* data_vtxZ = (TH1D*)data_file->Get("fVtxZ");
	TH3D* data_vtxContributors = (TH3D*)data_file->Get("fVtxContributors");
	TH1D* data_vtxXContributors = (TH1D*)data_vtxContributors->ProjectionX();
	TH1D* data_vtxYContributors = (TH1D*)data_vtxContributors->ProjectionY();
	TH1D* data_vtxZContributors = (TH1D*)data_vtxContributors->ProjectionZ();
	TH1F* data_HistFiredBitsSPDCINT1B = (TH1F*)data_dirCINT1B->Get("fHistFiredBitsSPD");
	TH1F* data_HistFiredBitsSPDCINT1A = (TH1F*)data_dirCINT1A->Get("fHistFiredBitsSPD");
	TH1F* data_HistFiredBitsSPDCINT1C = (TH1F*)data_dirCINT1C->Get("fHistFiredBitsSPD");
	TH1F* data_HistFiredBitsSPDCINT1E = (TH1F*)data_dirCINT1E->Get("fHistFiredBitsSPD");
	data_etaPhi->SetName("data_etaPhi");
	data_eta->SetName("data_eta");
	data_phi->SetName("data_phi");
	data_deltaPhi->SetName("data_deltaPhi");
	data_deltaTheta->SetName("data_deltaTheta");
	data_vtxX->SetName("data_vtxX");
	data_vtxY->SetName("data_vtxY");
	data_vtxZ->SetName("data_vtxZ");
	data_vtxXContributors->SetName("data_vtxXContributors");
	data_vtxYContributors->SetName("data_vtxYContributors");
	data_vtxZContributors->SetName("data_vtxZContributors");
	mc_HistFiredBitsSPD->SetName("mc_HistFiredBitsSPD");
	data_HistFiredBitsSPDCINT1B->SetName("data_HistFiredBitsSPDCINT1B");
	data_HistFiredBitsSPDCINT1A->SetName("data_HistFiredBitsSPDCINT1A");
	data_HistFiredBitsSPDCINT1C->SetName("data_HistFiredBitsSPDCINT1C");
	data_HistFiredBitsSPDCINT1E->SetName("data_HistFiredBitsSPDCINT1E");

	// ratios
	TH1D* ratio_phi = new TH1D("ratio_phi","#phi; #phi [rad]; Ratio", 18*5,0,2*TMath::Pi());
	TH1D* ratio_phiEta1 = new TH1D("ratio_phiEta1","#phi; #phi (|#eta|<1.4) [rad]; Ratio", 18*5,0,2*TMath::Pi() );
	TH1D* ratio_phiEtaNeg = new TH1D("ratio_phiEtaNeg","#phi; #phi (#eta<0) [rad]; Ratio", 18*5,0,2*TMath::Pi() );
	TH1D* ratio_phiEtaPos = new TH1D("ratio_phiEtaPos","#phi; #phi (#eta>0) [rad]; Ratio", 18*5,0,2*TMath::Pi() );
	TH1D* ratio_deltaPhi = new TH1D("ratio_deltaPhi","#Delta#phi; #Delta #phi;Ratio", 500, -1, 1);
	TH1D* ratio_deltaTheta = new TH1D("ratio_deltaTheta","#Delta#theta; #Delta #theta;Ratio", 500, -0.2, 0.2);
	TH1D* ratio_vtxX = new TH1D("ratio_vtxX", "x SPD primary vertex; x [cm]; Ratio",100,-1,1);
	TH1D* ratio_vtxY = new TH1D("ratio_vtxY", "y SPD primary vertex; y [cm]; Ratio",100,-1,1);
	TH1D* ratio_vtxZ = new TH1D("ratio_vtxZ", "z SPD primary vertex; z [cm]; Ratio",100,-30,30);
	TH1D* ratio_vtxXContributors = new TH1D("ratio_vtxXContributors", "x SPD primary vertex; x_{SPDvtx (N. Contributors > 0)} [cm]; Ratio",100,-1,1);
	TH1D* ratio_vtxYContributors = new TH1D("ratio_vtxYContributors", "y SPD primary vertex; y_{SPDvtx (N. Contributors > 0)} [cm]; Ratio",100,-1,1);
	TH1D* ratio_vtxZContributors = new TH1D("ratio_vtxZContributors", "z SPD primary vertex; z_{SPDvtx (N. Contributors > 0)} [cm]; Ratio",100,-30,30);
	TH1D* ratio_HistFiredBitsSPD = new TH1D("ratio_HistFiredBitsSPD", "SPD GFO; chip number; Ratio",1200,-0.5,1199.5);

	TCanvas* cEta = new TCanvas("cEta", "#eta", 50,50, 550, 550);
	TCanvas* cPhi = new TCanvas("cPhi", "#phi", 50,50,550, 550);
	TCanvas* cPhiEta1 = new TCanvas("cPhiEta1", "#phi, |#eta| < 1.4", 50,50,550, 550);
	TCanvas* cPhiEtaNeg = new TCanvas("cPhiEtaNeg", "#phi, #eta < 0", 50,50,550, 550);
	TCanvas* cPhiEtaPos = new TCanvas("cPhiEtaPos", "#phi, #eta > 0", 50,50,550, 550);
	TCanvas* cDeltaPhi = new TCanvas("cDeltaPhi", "#Delta#phi", 50,50, 550, 550);
	TCanvas* cDeltaTheta = new TCanvas("cDeltaTheta", "#Delta#theta", 50,50, 550, 550);
	TCanvas* cvtxX = new TCanvas("cvtxX", "vtxX", 50,50,550, 550);
	TCanvas* cvtxY = new TCanvas("cvtxY", "vtxY", 50,50,550, 550);
	TCanvas* cvtxZ = new TCanvas("cvtxZ", "vtxZ", 50,50,550, 550);
	TCanvas* cvtxXContributors = new TCanvas("cvtxXContributors", "vtxXContributors", 50,50,550, 550);
	TCanvas* cvtxYContributors = new TCanvas("cvtxYContributors", "vtxYContributors", 50,50,550, 550);
	TCanvas* cvtxZContributors = new TCanvas("cvtxZContributors", "vtxZContributors", 50,50,550, 550);
	TCanvas* cSPDGFO = new TCanvas("cSPDGFO", "SPD GFO", 50,50,950, 550);
	TCanvas* cratioPhi = new TCanvas("cratioPhi", "#phi", 50,50,550, 550);
	TCanvas* cratioPhiEta1 = new TCanvas("cratioPhiEta1", "#phi, |#eta| < 1.4", 50,50,550, 550);
	TCanvas* cratioPhiEtaNeg = new TCanvas("cratioPhiNeg", "#phi, #eta < 0", 50,50,550, 550);
	TCanvas* cratioPhiEtaPos = new TCanvas("cratioPhiPos", "#phi, #eta > 0", 50,50,550, 550);
	TCanvas* cratioDeltaPhi = new TCanvas("cratioDeltaPhi", "#Delta#phi", 50,50, 550, 550);
	TCanvas* cratioDeltaTheta = new TCanvas("cratioDeltaTheta", "#Delta#theta", 50,50, 550, 550);
	TCanvas* cratiovtxX = new TCanvas("cratiovtxX", "vtxX", 50,50,550, 550);
	TCanvas* cratiovtxY = new TCanvas("cratiovtxY", "vtxY", 50,50,550, 550);
	TCanvas* cratiovtxZ = new TCanvas("cratiovtxZ", "vtxZ", 50,50,550, 550);
	TCanvas* cratiovtxXContributors = new TCanvas("cratiovtxXContributors", "vtxXContributors", 50,50,550, 550);
	TCanvas* cratiovtxYContributors = new TCanvas("cratiovtxYContributors", "vtxYContributors", 50,50,550, 550);
	TCanvas* cratiovtxZContributors = new TCanvas("cratiovtxZContributors", "vtxZContributors", 50,50,550, 550);
	TCanvas* cratioSPDGFO = new TCanvas("cratioSPDGFO", "SPD GFO", 50,50,950, 550);

	cEta->SetLeftMargin(0.15);
	cEta->SetRightMargin(0.05);
	cPhi->SetLeftMargin(0.15);
	cPhi->SetRightMargin(0.05);
	cPhiEta1->SetLeftMargin(0.15);
	cPhiEta1->SetRightMargin(0.05);
	cPhiEtaNeg->SetLeftMargin(0.15);
	cPhiEtaNeg->SetRightMargin(0.05);
	cPhiEtaPos->SetLeftMargin(0.15);
	cPhiEtaPos->SetRightMargin(0.05);
	cDeltaPhi->SetLeftMargin(0.15);
	cDeltaPhi->SetRightMargin(0.05);
	cDeltaTheta->SetLeftMargin(0.15);
	cDeltaTheta->SetRightMargin(0.05);
	cvtxX->SetLeftMargin(0.15);
	cvtxX->SetRightMargin(0.05);
	cvtxY->SetLeftMargin(0.15);
	cvtxY->SetRightMargin(0.05);
	cvtxZ->SetLeftMargin(0.15);
	cvtxZ->SetRightMargin(0.05);
	cvtxXContributors->SetLeftMargin(0.15);
	cvtxXContributors->SetRightMargin(0.05);
	cvtxYContributors->SetLeftMargin(0.15);
	cvtxYContributors->SetRightMargin(0.05);
	cvtxZContributors->SetLeftMargin(0.15);
	cvtxZContributors->SetRightMargin(0.05);
	cSPDGFO->SetLeftMargin(0.15);
	cSPDGFO->SetRightMargin(0.05);

	cratioPhi->SetLeftMargin(0.15);
	cratioPhi->SetRightMargin(0.05);
	cratioPhiEta1->SetLeftMargin(0.15);
	cratioPhiEta1->SetRightMargin(0.05);
	cratioPhiEtaNeg->SetLeftMargin(0.15);
	cratioPhiEtaNeg->SetRightMargin(0.05);
	cratioPhiEtaPos->SetLeftMargin(0.15);
	cratioPhiEtaPos->SetRightMargin(0.05);
	cratioDeltaPhi->SetLeftMargin(0.15);
	cratioDeltaPhi->SetRightMargin(0.05);
	cratioDeltaTheta->SetLeftMargin(0.15);
	cratioDeltaTheta->SetRightMargin(0.05);
	cratiovtxX->SetLeftMargin(0.15);
	cratiovtxX->SetRightMargin(0.05);
	cratiovtxY->SetLeftMargin(0.15);
	cratiovtxY->SetRightMargin(0.05);
	cratiovtxZ->SetLeftMargin(0.15);
	cratiovtxZ->SetRightMargin(0.05);
	cratiovtxXContributors->SetLeftMargin(0.15);
	cratiovtxXContributors->SetRightMargin(0.05);
	cratiovtxYContributors->SetLeftMargin(0.15);
	cratiovtxYContributors->SetRightMargin(0.05);
	cratiovtxZContributors->SetLeftMargin(0.15);
	cratiovtxZContributors->SetRightMargin(0.05);
	cratioSPDGFO->SetLeftMargin(0.15);
	cratioSPDGFO->SetRightMargin(0.05);

	cDeltaPhi->SetLogy();
	cDeltaTheta->SetLogy();
	cvtxX->SetLogy();
	cvtxY->SetLogy();
	cvtxZ->SetLogy();
	cvtxXContributors->SetLogy();
	cvtxYContributors->SetLogy();
	cvtxZContributors->SetLogy();

	mc_eta->SetLineColor(2);
	mc_eta->SetLineWidth(2);
	mc_eta->GetYaxis()->SetTitleOffset(1.6);
	mc_eta->GetYaxis()->SetTitle("Entries");
	mc_eta->Scale(data_eta->Integral()/mc_eta->Integral());
	data_eta->SetLineColor(4);
	data_eta->SetLineWidth(2);
	data_eta->SetLineStyle(2);
	mc_eta->SetMaximum(TMath::Max(mc_eta->GetMaximum(),data_eta->GetMaximum())*1.2);
	cEta->cd();
	mc_eta->Draw("hist");
	data_eta->Draw("hist same");

	mc_phi->SetLineColor(2);
	mc_phi->SetLineWidth(2);
	mc_phi->Scale(data_phi->Integral()/mc_phi->Integral());
	mc_phi->GetYaxis()->SetTitleOffset(1.9);
	mc_phi->GetYaxis()->SetTitle("Entries");
	data_phi->SetLineColor(4);
	data_phi->SetLineWidth(2);
	data_phi->SetLineStyle(2);
	mc_phi->SetMaximum(TMath::Max(mc_phi->GetMaximum(),data_phi->GetMaximum())*1.2);
	cPhi->cd();
	mc_phi->Draw("hist err");
	data_phi->Draw("hist err same");

	TLegend* leg = new TLegend(0.6,0.88,0.9,0.75);
	leg->SetFillColor(0);
	leg->AddEntry(mc_phi,Form("%s",mc_path),"L");
	leg->AddEntry(data_phi,Form("%s",data_path),"L");
	leg->Draw("same");

	mc_phiEta1->SetLineColor(2);
	mc_phiEta1->SetLineWidth(2);
	mc_phiEta1->Scale(data_phiEta1->Integral()/mc_phiEta1->Integral());
	mc_phiEta1->GetYaxis()->SetTitleOffset(1.9);
	mc_phiEta1->GetYaxis()->SetTitle("Entries");
	mc_phiEta1->GetXaxis()->SetTitle("#phi (|#eta|<1.4) [rad]");
	data_phiEta1->SetLineColor(4);
	data_phiEta1->SetLineWidth(2);
	data_phiEta1->SetLineStyle(2);
	mc_phiEta1->SetMaximum(TMath::Max(mc_phiEta1->GetMaximum(),data_phiEta1->GetMaximum())*1.2);
	cPhiEta1->cd();
	mc_phiEta1->Draw("hist err");
	data_phiEta1->Draw("hist err same");

	mc_phiEtaNeg->SetLineColor(2);
	mc_phiEtaNeg->SetLineWidth(2);
	mc_phiEtaNeg->Scale(data_phiEtaNeg->Integral()/mc_phiEtaNeg->Integral());
	mc_phiEtaNeg->GetYaxis()->SetTitleOffset(1.9);
	mc_phiEtaNeg->GetYaxis()->SetTitle("Entries");
	mc_phiEtaNeg->GetXaxis()->SetTitle("#phi (-1.4<#eta<0) [rad]");
	data_phiEtaNeg->SetLineColor(4);
	data_phiEtaNeg->SetLineWidth(2);
	data_phiEtaNeg->SetLineStyle(2);
	mc_phiEtaNeg->SetMaximum(TMath::Max(mc_phiEtaNeg->GetMaximum(),data_phiEtaNeg->GetMaximum())*1.2);
	cPhiEtaNeg->cd();
	mc_phiEtaNeg->Draw("hist err");
	data_phiEtaNeg->Draw("hist err same");

	mc_phiEtaPos->SetLineColor(2);
	mc_phiEtaPos->SetLineWidth(2);
	mc_phiEtaPos->Scale(data_phiEtaPos->Integral()/mc_phiEtaPos->Integral());
	mc_phiEtaPos->GetYaxis()->SetTitleOffset(1.9);
	mc_phiEtaPos->GetYaxis()->SetTitle("Entries");
	mc_phiEtaPos->GetXaxis()->SetTitle("#phi (0<#eta<1.4) [rad]");
	data_phiEtaPos->SetLineColor(4);
	data_phiEtaPos->SetLineWidth(2);
	data_phiEtaPos->SetLineStyle(2);
	mc_phiEtaPos->SetMaximum(TMath::Max(mc_phiEtaPos->GetMaximum(),data_phiEtaPos->GetMaximum())*1.2);
	cPhiEtaPos->cd();
	mc_phiEtaPos->Draw("hist err");
	data_phiEtaPos->Draw("hist err same");

	mc_deltaPhi->SetLineColor(2);
	mc_deltaPhi->SetLineWidth(2);
	mc_deltaPhi->GetYaxis()->SetTitleOffset(1.8);
	mc_deltaPhi->Scale(data_deltaPhi->Integral()/mc_deltaPhi->Integral());
	mc_deltaPhi->GetXaxis()->SetRangeUser(-0.15,0.15);
	data_deltaPhi->SetLineColor(4);
	data_deltaPhi->SetLineWidth(2);
	data_deltaPhi->SetLineStyle(2);
	mc_deltaPhi->SetMaximum(TMath::Max(mc_deltaPhi->GetMaximum(),data_deltaPhi->GetMaximum())*1.2);
	cDeltaPhi->cd();
	mc_deltaPhi->Draw("hist err");
	data_deltaPhi->Draw("hist err same");

	mc_deltaTheta->SetLineColor(2);
	mc_deltaTheta->SetLineWidth(2);
	mc_deltaTheta->GetYaxis()->SetTitleOffset(1.8);
       	mc_deltaTheta->Scale(data_deltaTheta->Integral()/mc_deltaTheta->Integral());
	mc_deltaTheta->GetXaxis()->SetRangeUser(-0.05,0.05);
	data_deltaTheta->SetLineColor(4);
	data_deltaTheta->SetLineWidth(2);
	data_deltaTheta->SetLineStyle(2);
	mc_deltaTheta->SetMaximum(TMath::Max(mc_deltaTheta->GetMaximum(),data_deltaTheta->GetMaximum())*1.2);
	cDeltaTheta->cd();
	mc_deltaTheta->Draw("hist err");
	data_deltaTheta->Draw("hist err same");

	mc_vtxX->SetLineColor(2);
	mc_vtxX->SetLineWidth(2);
	mc_vtxX->GetYaxis()->SetTitleOffset(1.6);
	mc_vtxX->GetYaxis()->SetTitle("Entries");
	mc_vtxX->GetXaxis()->SetTitle("x_{SPDvtx} [cm]");
       	mc_vtxX->Scale(data_vtxX->Integral()/mc_vtxX->Integral());
	data_vtxX->SetLineColor(4);
	data_vtxX->SetLineWidth(2);
	data_vtxX->SetLineStyle(2);
	mc_vtxX->SetMaximum(TMath::Max(mc_vtxX->GetMaximum(),data_vtxX->GetMaximum())*1.2);
	cvtxX->cd();
	mc_vtxX->Draw("hist err");
	data_vtxX->Draw("hist err same");

	mc_vtxY->SetLineColor(2);
	mc_vtxY->SetLineWidth(2);
	mc_vtxY->GetYaxis()->SetTitleOffset(1.6);
	mc_vtxY->GetYaxis()->SetTitle("Entries");
	mc_vtxY->GetXaxis()->SetTitle("y_{SPDvtx} [cm]");
	mc_vtxY->Scale(data_vtxY->Integral()/mc_vtxY->Integral());
	data_vtxY->SetLineColor(4);
	data_vtxY->SetLineWidth(2);
	data_vtxY->SetLineStyle(2);
	mc_vtxY->SetMaximum(TMath::Max(mc_vtxY->GetMaximum(),data_vtxY->GetMaximum())*1.2);
	cvtxY->cd();
	mc_vtxY->Draw("hist err");
	data_vtxY->Draw("hist err same");

	mc_vtxZ->SetLineColor(2);
	mc_vtxZ->SetLineWidth(2);
	mc_vtxZ->GetYaxis()->SetTitleOffset(1.6);
	mc_vtxZ->GetYaxis()->SetTitle("Entries");
	mc_vtxZ->GetXaxis()->SetTitle("z_{SPDvtx} [cm]");
       	mc_vtxZ->Scale(data_vtxZ->Integral()/mc_vtxZ->Integral());
	data_vtxZ->SetLineColor(4);
	data_vtxZ->SetLineWidth(2);
	data_vtxZ->SetLineStyle(2);
	mc_vtxZ->SetMaximum(TMath::Max(mc_vtxZ->GetMaximum(),data_vtxZ->GetMaximum())*1.2);
	cvtxZ->cd();
	mc_vtxZ->Draw("hist err");
	data_vtxZ->Draw("hist err same");

	mc_vtxXContributors->SetLineColor(2);
	mc_vtxXContributors->SetLineWidth(2);
	mc_vtxXContributors->GetYaxis()->SetTitleOffset(1.6);
	mc_vtxXContributors->GetYaxis()->SetTitle("Entries");
	mc_vtxXContributors->GetXaxis()->SetTitle("x_{SPDvtx (N. Contributors > 0)} [cm]");
       	mc_vtxXContributors->Scale(data_vtxXContributors->Integral()/mc_vtxXContributors->Integral());
	data_vtxXContributors->SetLineColor(4);
	data_vtxXContributors->SetLineWidth(2);
	data_vtxXContributors->SetLineStyle(2);
	mc_vtxXContributors->SetMaximum(TMath::Max(mc_vtxXContributors->GetMaximum(),data_vtxXContributors->GetMaximum())*1.2);
	cvtxXContributors->cd();
	mc_vtxXContributors->Draw("hist err");
	data_vtxXContributors->Draw("hist err same");

	mc_vtxYContributors->SetLineColor(2);
	mc_vtxYContributors->SetLineWidth(2);
	mc_vtxYContributors->GetYaxis()->SetTitleOffset(1.6);
	mc_vtxYContributors->GetYaxis()->SetTitle("Entries");
	mc_vtxYContributors->GetXaxis()->SetTitle("y_{SPDvtx (N. Contributors > 0)} [cm]");
	mc_vtxYContributors->Scale(data_vtxYContributors->Integral()/mc_vtxYContributors->Integral());
	data_vtxYContributors->SetLineColor(4);
	data_vtxYContributors->SetLineWidth(2);
	data_vtxYContributors->SetLineStyle(2);
	mc_vtxYContributors->SetMaximum(TMath::Max(mc_vtxYContributors->GetMaximum(),data_vtxYContributors->GetMaximum())*1.2);
	cvtxYContributors->cd();
	mc_vtxYContributors->Draw("hist err");
	data_vtxYContributors->Draw("hist err same");

	mc_vtxZContributors->SetLineColor(2);
	mc_vtxZContributors->SetLineWidth(2);
	mc_vtxZContributors->GetYaxis()->SetTitleOffset(1.6);
	mc_vtxZContributors->GetYaxis()->SetTitle("Entries");
	mc_vtxZContributors->GetXaxis()->SetTitle("z_{SPDvtx (N. Contributors > 0)} [cm]");
       	mc_vtxZContributors->Scale(data_vtxZContributors->Integral()/mc_vtxZContributors->Integral());
	data_vtxZContributors->SetLineColor(4);
	data_vtxZContributors->SetLineWidth(2);
	data_vtxZContributors->SetLineStyle(2);
	mc_vtxZContributors->SetMaximum(TMath::Max(mc_vtxZContributors->GetMaximum(),data_vtxZContributors->GetMaximum())*1.2);
	cvtxZContributors->cd();
	mc_vtxZContributors->Draw("hist err");
	data_vtxZContributors->Draw("hist err same");

	mc_HistFiredBitsSPD->SetLineColor(2);
	mc_HistFiredBitsSPD->SetLineWidth(2);
	mc_HistFiredBitsSPD->GetYaxis()->SetTitleOffset(1.2);
	mc_HistFiredBitsSPD->GetYaxis()->SetTitle("Entries");
       	mc_HistFiredBitsSPD->Scale(data_HistFiredBitsSPDCINT1B->Integral()/mc_HistFiredBitsSPD->Integral());
	data_HistFiredBitsSPDCINT1B->SetLineColor(4);
	data_HistFiredBitsSPDCINT1B->SetLineWidth(2);
	data_HistFiredBitsSPDCINT1B->SetLineStyle(2);
	mc_HistFiredBitsSPD->SetMaximum(TMath::Max(mc_HistFiredBitsSPD->GetMaximum(),data_HistFiredBitsSPDCINT1B->GetMaximum())*1.2);
	cSPDGFO->cd();
	mc_HistFiredBitsSPD->Draw("hist err");
	data_HistFiredBitsSPDCINT1B->Draw("hist err same");

	ratio_phi->SetLineWidth(2);
	ratio_phi->Divide(mc_phi,data_phi);
	ratio_phi->GetYaxis()->SetRangeUser(0.5,1.5);
	ratio_phi->GetYaxis()->SetTitleOffset(1.9);
	ratio_phi->GetYaxis()->SetTitle("Ratio");
	cratioPhi->cd();
	ratio_phi->Draw("hist");

	TPaveText* text = new TPaveText(0.2,1,0.8,0.9,"NDC");
	text->SetFillColor(0);
	text->AddText(Form("%s/%s",mc_path,data_path));
	text->Draw("same");

       	ratio_phiEta1->SetLineWidth(2);
	ratio_phiEta1->Divide(mc_phiEta1,data_phiEta1);
	ratio_phiEta1->GetYaxis()->SetRangeUser(0.5,1.5);
	ratio_phiEta1->GetYaxis()->SetTitleOffset(1.9);
	ratio_phiEta1->GetYaxis()->SetTitle("Ratio");
	cratioPhiEta1->cd();
	ratio_phiEta1->Draw("hist");

       	ratio_phiEtaNeg->SetLineWidth(2);
	ratio_phiEtaNeg->Divide(mc_phiEtaNeg,data_phiEtaNeg);
	ratio_phiEtaNeg->GetYaxis()->SetRangeUser(0.5,1.5);
	ratio_phiEtaNeg->GetYaxis()->SetTitleOffset(1.9);
	ratio_phiEtaNeg->GetYaxis()->SetTitle("Ratio");
	cratioPhiEtaNeg->cd();
	ratio_phiEtaNeg->Draw("hist");

       	ratio_phiEtaPos->SetLineWidth(2);
	ratio_phiEtaPos->Divide(mc_phiEtaPos,data_phiEtaPos);
	ratio_phiEtaPos->GetYaxis()->SetRangeUser(0.5,1.5);
	ratio_phiEtaPos->GetYaxis()->SetTitleOffset(1.9);
	ratio_phiEtaPos->GetYaxis()->SetTitle("Ratio");
	cratioPhiEtaPos->cd();
	ratio_phiEtaPos->Draw("hist");

	ratio_deltaPhi->SetLineWidth(2);
	ratio_deltaPhi->Divide(mc_deltaPhi,data_deltaPhi);
	ratio_deltaPhi->GetYaxis()->SetRangeUser(0.5,1.5);
	ratio_deltaPhi->GetXaxis()->SetRangeUser(-0.15,0.15);
	ratio_deltaPhi->GetYaxis()->SetTitleOffset(1.9);
	ratio_deltaPhi->GetYaxis()->SetTitle("Ratio");
	cratioDeltaPhi->cd();
	ratio_deltaPhi->Draw("hist");

	ratio_deltaTheta->SetLineWidth(2);
	ratio_deltaTheta->Divide(mc_deltaTheta,data_deltaTheta);
	ratio_deltaTheta->GetXaxis()->SetRangeUser(-0.05,0.05);
	ratio_deltaTheta->GetYaxis()->SetRangeUser(0.5,1.5);
	ratio_deltaTheta->GetYaxis()->SetTitleOffset(1.9);
	ratio_deltaTheta->GetYaxis()->SetTitle("Ratio");
	cratioDeltaTheta->cd();
	ratio_deltaTheta->Draw("hist");

	ratio_vtxX->SetLineWidth(2);
	ratio_vtxX->Divide(mc_vtxX,data_vtxX);
	ratio_vtxX->GetYaxis()->SetRangeUser(0.5,1.5);
	ratio_vtxX->GetYaxis()->SetTitleOffset(1.9);
	ratio_vtxX->GetYaxis()->SetTitle("Ratio");
	cratiovtxX->cd();
	ratio_vtxX->Draw("hist");

	ratio_vtxY->SetLineWidth(2);
	ratio_vtxY->Divide(mc_vtxY,data_vtxY);
	ratio_vtxY->GetYaxis()->SetRangeUser(0.5,1.5);
	ratio_vtxY->GetYaxis()->SetTitleOffset(1.9);
	ratio_vtxY->GetYaxis()->SetTitle("Ratio");
	cratiovtxY->cd();
	ratio_vtxY->Draw("hist");

	ratio_vtxZ->SetLineWidth(2);
	ratio_vtxZ->Divide(mc_vtxZ,data_vtxZ);
	ratio_vtxZ->GetYaxis()->SetRangeUser(0.5,1.5);
	ratio_vtxZ->GetYaxis()->SetTitleOffset(1.9);
	ratio_vtxZ->GetYaxis()->SetTitle("Ratio");
	cratiovtxZ->cd();
	ratio_vtxZ->Draw("hist");

	ratio_vtxXContributors->SetLineWidth(2);
	ratio_vtxXContributors->Divide(mc_vtxXContributors,data_vtxXContributors);
	ratio_vtxXContributors->GetYaxis()->SetRangeUser(0.5,1.5);
	ratio_vtxXContributors->GetYaxis()->SetTitleOffset(1.9);
	ratio_vtxXContributors->GetYaxis()->SetTitle("Ratio");
	cratiovtxXContributors->cd();
	ratio_vtxXContributors->Draw("hist");

	ratio_vtxYContributors->SetLineWidth(2);
	ratio_vtxYContributors->Divide(mc_vtxYContributors,data_vtxYContributors);
	ratio_vtxYContributors->GetYaxis()->SetRangeUser(0.5,1.5);
	ratio_vtxYContributors->GetYaxis()->SetTitleOffset(1.9);
	ratio_vtxYContributors->GetYaxis()->SetTitle("Ratio");
	cratiovtxYContributors->cd();
	ratio_vtxYContributors->Draw("hist");

	ratio_vtxZContributors->SetLineWidth(2);
	ratio_vtxZContributors->Divide(mc_vtxZContributors,data_vtxZContributors);
	ratio_vtxZContributors->GetYaxis()->SetRangeUser(0.5,1.5);
	ratio_vtxZContributors->GetYaxis()->SetTitleOffset(1.9);
	ratio_vtxZContributors->GetYaxis()->SetTitle("Ratio");
	cratiovtxZContributors->cd();
	ratio_vtxZContributors->Draw("hist");

	ratio_HistFiredBitsSPD->SetLineWidth(2);
	ratio_HistFiredBitsSPD->Divide(mc_HistFiredBitsSPD,data_HistFiredBitsSPDCINT1B);
	ratio_HistFiredBitsSPD->GetYaxis()->SetRangeUser(0.5,1.5);
	ratio_HistFiredBitsSPD->GetYaxis()->SetTitleOffset(1.9);
	ratio_HistFiredBitsSPD->GetYaxis()->SetTitle("Ratio");
	cratioSPDGFO->cd();
	ratio_HistFiredBitsSPD->Draw("hist");

	cEta->cd();
	cEta->Print(Form("Plots/Plots_%s/eta.gif",mc_path));
	cEta->Print(Form("Plots/Plots_%s/eta.eps",mc_path));
	cPhi->cd();
	cPhi->Print(Form("Plots/Plots_%s/phi.gif",mc_path));
	cPhi->Print(Form("Plots/Plots_%s/phi.eps",mc_path));
	cPhiEta1->cd();
	cPhiEta1->Print(Form("Plots/Plots_%s/phiEta1.gif",mc_path));
	cPhiEta1->Print(Form("Plots/Plots_%s/phiEta1.eps",mc_path));
	cPhiEtaNeg->cd();
	cPhiEtaNeg->Print(Form("Plots/Plots_%s/phiEtaNeg.gif",mc_path));
	cPhiEtaNeg->Print(Form("Plots/Plots_%s/phiEtaNeg.eps",mc_path));
	cPhiEtaPos->cd();
	cPhiEtaPos->Print(Form("Plots/Plots_%s/phiEtaPos.gif",mc_path));
	cPhiEtaPos->Print(Form("Plots/Plots_%s/phiEtaPos.eps",mc_path));
	cDeltaPhi->cd();
	cDeltaPhi->Print(Form("Plots/Plots_%s/deltaPhi.gif",mc_path));
	cDeltaPhi->Print(Form("Plots/Plots_%s/deltaPhi.eps",mc_path));
	cDeltaTheta->cd();
	cDeltaTheta->Print(Form("Plots/Plots_%s/deltaTheta.gif",mc_path));
	cDeltaTheta->Print(Form("Plots/Plots_%s/deltaTheta.eps",mc_path));
	cvtxX->cd();
	cvtxX->Print(Form("Plots/Plots_%s/vtxX.gif",mc_path));
	cvtxX->Print(Form("Plots/Plots_%s/vtxX.eps",mc_path));
	cvtxY->cd();
	cvtxY->Print(Form("Plots/Plots_%s/vtxY.gif",mc_path));
	cvtxY->Print(Form("Plots/Plots_%s/vtxY.eps",mc_path));
	cvtxZ->cd();
	cvtxZ->Print(Form("Plots/Plots_%s/vtxZ.gif",mc_path));
	cvtxZ->Print(Form("Plots/Plots_%s/vtxZ.eps",mc_path));
	cvtxXContributors->cd();
	cvtxXContributors->Print(Form("Plots/Plots_%s/vtxXContributors.gif",mc_path));
	cvtxXContributors->Print(Form("Plots/Plots_%s/vtxXContributors.eps",mc_path));
	cvtxYContributors->cd();
	cvtxYContributors->Print(Form("Plots/Plots_%s/vtxYContributors.gif",mc_path));
	cvtxYContributors->Print(Form("Plots/Plots_%s/vtxYContributors.eps",mc_path));
	cvtxZContributors->cd();
	cvtxZContributors->Print(Form("Plots/Plots_%s/vtxZContributors.gif",mc_path));
	cvtxZContributors->Print(Form("Plots/Plots_%s/vtxZContributors.eps",mc_path));
	cSPDGFO->cd();
	cSPDGFO->Print(Form("Plots/Plots_%s/SPDGFO.gif",mc_path));
	cSPDGFO->Print(Form("Plots/Plots_%s/SPDGFO.eps",mc_path));

	cratioPhi->cd();
	cratioPhi->Print(Form("Plots/Plots_%s/ratiophi.gif",mc_path));
	cratioPhi->Print(Form("Plots/Plots_%s/ratiophi.eps",mc_path));
	cratioPhiEta1->cd();
	cratioPhiEta1->Print(Form("Plots/Plots_%s/ratiophiEta1.gif",mc_path));
	cratioPhiEta1->Print(Form("Plots/Plots_%s/ratiophiEta1.eps",mc_path));
	cratioPhiEtaNeg->cd();
	cratioPhiEtaNeg->Print(Form("Plots/Plots_%s/ratiophiEtaNeg.gif",mc_path));
	cratioPhiEtaNeg->Print(Form("Plots/Plots_%s/ratiophiEtaNeg.eps",mc_path));
	cratioPhiEtaPos->cd();
	cratioPhiEtaPos->Print(Form("Plots/Plots_%s/ratiophiEtaPos.gif",mc_path));
	cratioPhiEtaPos->Print(Form("Plots/Plots_%s/ratiophiEtaPos.eps",mc_path));
	cratioDeltaPhi->cd();
	cratioDeltaPhi->Print(Form("Plots/Plots_%s/ratiodeltaPhi.gif",mc_path));
	cratioDeltaPhi->Print(Form("Plots/Plots_%s/ratiodeltaPhi.eps",mc_path));
	cratioDeltaTheta->cd();
	cratioDeltaTheta->Print(Form("Plots/Plots_%s/ratiodeltaTheta.gif",mc_path));
	cratioDeltaTheta->Print(Form("Plots/Plots_%s/ratiodeltaTheta.eps",mc_path));
	cratiovtxX->cd();
	cratiovtxX->Print(Form("Plots/Plots_%s/ratiovtxX.gif",mc_path));
	cratiovtxX->Print(Form("Plots/Plots_%s/ratiovtxX.eps",mc_path));
	cratiovtxY->cd();
	cratiovtxY->Print(Form("Plots/Plots_%s/ratiovtxY.gif",mc_path));
	cratiovtxY->Print(Form("Plots/Plots_%s/ratiovtxY.eps",mc_path));
	cratiovtxZ->cd();
	cratiovtxZ->Print(Form("Plots/Plots_%s/ratiovtxZ.gif",mc_path));
	cratiovtxZ->Print(Form("Plots/Plots_%s/ratiovtxZ.eps",mc_path));
	cratiovtxXContributors->cd();
	cratiovtxXContributors->Print(Form("Plots/Plots_%s/ratiovtxXContributors.gif",mc_path));
	cratiovtxXContributors->Print(Form("Plots/Plots_%s/ratiovtxXContributors.eps",mc_path));
	cratiovtxYContributors->cd();
	cratiovtxYContributors->Print(Form("Plots/Plots_%s/ratiovtxYContributors.gif",mc_path));
	cratiovtxYContributors->Print(Form("Plots/Plots_%s/ratiovtxYContributors.eps",mc_path));
	cratiovtxZContributors->cd();
	cratiovtxZContributors->Print(Form("Plots/Plots_%s/ratiovtxZContributors.gif",mc_path));
	cratiovtxZContributors->Print(Form("Plots/Plots_%s/ratiovtxZContributors.eps",mc_path));
	cratioSPDGFO->cd();
	cratioSPDGFO->Print(Form("Plots/Plots_%s/ratioSPDGFO.gif",mc_path));
	cratioSPDGFO->Print(Form("Plots/Plots_%s/ratioSPDGFO.eps",mc_path));

	TFile* plotsfile = new TFile(Form("Plots/Plots_%s/plots.root",mc_path),"RECREATE");
	mc_eta->Write();
	mc_phi->Write();
	mc_phiEta1->Write();
	mc_phiEtaNeg->Write();
	mc_phiEtaPos->Write();
	mc_deltaPhi->Write();
	mc_deltaTheta->Write();
	mc_vtxX->Write();
	mc_vtxY->Write();
	mc_vtxZ->Write();
	mc_vtxXContributors->Write();
	mc_vtxYContributors->Write();
	mc_vtxZContributors->Write();
	mc_HistFiredBitsSPD->Write();
	data_eta->Write();
	data_phi->Write();
	data_phiEta1->Write();
	data_phiEtaNeg->Write();
	data_phiEtaPos->Write();
	data_deltaPhi->Write();
	data_deltaTheta->Write();
	data_vtxX->Write();
	data_vtxY->Write();
	data_vtxZ->Write();
	data_vtxXContributors->Write();
	data_vtxYContributors->Write();
	data_vtxZContributors->Write();
	data_HistFiredBitsSPDCINT1B->Write();
	data_HistFiredBitsSPDCINT1A->Write();
	data_HistFiredBitsSPDCINT1C->Write();
	data_HistFiredBitsSPDCINT1E->Write();
	ratio_phi->Write();
	ratio_phiEta1->Write();
	ratio_phiEtaNeg->Write();
	ratio_phiEtaPos->Write();
	ratio_deltaPhi->Write();
	ratio_deltaTheta->Write();
	ratio_vtxX->Write();
	ratio_vtxY->Write();
	ratio_vtxZ->Write();
	ratio_vtxXContributors->Write();
	ratio_vtxYContributors->Write();
	ratio_vtxZContributors->Write();
	ratio_HistFiredBitsSPD->Write();
	plotsfile->Close();
}