etaRatio.C

#include "TFile.h"
#include <iostream>
#include <vector>
#include "TChain.h"
#include "TTree.h"
#include "TSystem.h"
#include "TCanvas.h"
#include "TH1.h"
#include "TH2.h"
#include "TF1.h"
#include "TGraphErrors.h"
#include "TGraphAsymmErrors.h"
#include "TLorentzVector.h"
#include "TClonesArray.h"
#include "TROOT.h"
#include "TMath.h"
#include "TStyle.h"
#include "TColor.h"
#include <math.h>
#include "TLegend.h"
#include "TPad.h"
#include "TLine.h"
#include "TH1F.h"
//#include "TRatioPlot.h"
#include "../../plugins/CosmicSplittingResolutionNtuple.h"

//you first decalre all the varaibles you are goingt o pick from the tree and use in this code. Be careful that the type of the varaible (Int, Float etc) have to be the same that when it was filled in the tree maker


//my interesting variables 

float eta[n_tracks][n_track_pos];
short charge[n_tracks][n_track_pos];

float unprop_theta[n_tracks][n_track_pos];
short unprop_charge[n_tracks][n_track_pos];

ushort pixel_hits[n_tracks][n_track_pos];

ushort strip_layers[n_tracks][n_track_pos];

float unprop_pt[n_tracks][n_track_pos];


////////////////////////////////////////////////////////////////////////////////


//Declaration of hitsogramm i plan to plot with this code

TH1F *etahistoGlbUP =new TH1F("etahistoGlb","",20,-3,3);
TH1F *etahistoGlbLOW =new TH1F("etahistoGlbLOW","",20,-3,3);

TH1F *etahistoTkonlyUP =new TH1F("etahistoTkonly","",20,-3,3);
TH1F *etahistoTkonlyLOW =new TH1F("etahistoTkonlyLOW","",20,-3,3);

TH1F *etahistoTpfmsUP =new TH1F("etahistoTpfms","",20,-3,3);
TH1F *etahistoTpfmsLOW =new TH1F("etahistoTpfmsLOW","",20,-3,3);

TH1F *etahistoPickyUP =new TH1F("etahistoPicky","",20,-3,3);
TH1F *etahistoPickyLOW =new TH1F("etahistoPickyLOW","",20,-3,3);

TH1F *etahistoDYTUP =new TH1F("etahistoDYT","",20,-3,3);
TH1F *etahistoDYTLOW =new TH1F("etahistoDYTLOW","",20,-3,3);

TH1F *etahistoTunePUP =new TH1F("etahistoTuneP","",20,-3,3);
TH1F *etahistoTunePLOW =new TH1F("etahistoTunePLOW","",20,-3,3);


//Begining of teh function
void etaRatio(){



	// the next line is to link the tree you have

	TChain *treedata = new TChain("UTpickedTracks/t");

	// this is the root file you have input with your tree inside
	//treedata->Add("ROOT_Files/MC2016/cosmic_MC_2016.root");
	treedata->Add("ROOT_Files/cosmic_data_2016_2017/cosmic_data_16_17.root");

	// the following lines are to extract the branch of the trees, when there is a & it means that there is only one entry per event, like for run, lumi, event, dileetaon mass in my case for DY events.
	treedata->SetBranchAddress("charge",charge);
	//treedata->SetBranchAddress("eta",eta);
	treedata->SetBranchAddress("unprop_charge",unprop_charge);
	treedata->SetBranchAddress("unprop_theta",unprop_theta);
	treedata->SetBranchAddress("pixel_hits",pixel_hits);
	treedata->SetBranchAddress("strip_layers",strip_layers);
	treedata->SetBranchAddress("unprop_pt",unprop_pt);



	// This will get the numbers of entries of the tree.

	Long64_t nentries = treedata->GetEntries();

	//Decalaration of variable not in the tree that i use for the results

	Float_t sameCharge = 0.;
	Float_t goodMuon = 0.;

	for ( int p=0; p<nentries ;p++){
		// this loop over all the events in your tree

		treedata->GetEntry(p);
		//apply some quality cuts to select the entries that we think are interesting
		//
		//
		//Algorithm choice
		//
		for (int i=0;i<n_tracks;i++){
			for (int j=0;j<n_track_pos;j++){
				eta[i][j] = - log(tan(unprop_theta[i][j]/2));
			}
		}

		//selection cuts

		if (pixel_hits[8][0]>=1 && strip_layers[8][0]>=5 && pixel_hits[8][1]>=1 && strip_layers[8][1]>=5 && unprop_pt[8][1]>30.){

		//global
		if (charge[0][1] > 0) {
			etahistoGlbUP->Fill(eta[0][1]);
		}
		if (charge[0][1] < 0) {
			etahistoGlbLOW->Fill(eta[0][1]);
		}
		//Tkonly
		if (charge[2][1] > 0) {
			etahistoTkonlyUP->Fill(eta[2][1]);
		}
		if (charge[2][1] < 0) {
			etahistoTkonlyLOW->Fill(eta[2][1]);
		}
		//Tpfms
		if (charge[3][1] > 0) {
			etahistoTpfmsUP->Fill(eta[3][1]);
		}
		if (charge[3][1] < 0) {
			etahistoTpfmsLOW->Fill(eta[3][1]);
		}
		//Picky
		if (charge[4][1] > 0) {
			etahistoPickyUP->Fill(eta[4][1]);
		}
		if (charge[4][1] < 0) {
			etahistoPickyLOW->Fill(eta[4][1]);
		}
		//DYT
		if (charge[5][1] > 0) {
			etahistoDYTUP->Fill(eta[5][1]);
		}
		if (charge[5][1] < 0) {
			etahistoDYTLOW->Fill(eta[5][1]);
		}
		//TuneP
		if (charge[8][1] > 0) {
			etahistoTunePUP->Fill(eta[8][1]);
		}
		if (charge[8][1] < 0) {
			etahistoTunePLOW->Fill(eta[8][1]);
		}
}
	}// end loop p

	


        TCanvas *c = new TCanvas("c", "canvas", 700, 600);


	// Upper plot will be in pad5a
	TPad *pad5a = new TPad("pad3a", "pad3a", 0, 0.19, 1, 1.0);
	//pad5a->SetBottomMargin(0); // Upper and lower plot are joined
	pad5a->SetGridx();         // Vertical grid
	pad5a->Draw();             // Draw the upper pad: pad1
	pad5a->cd();               // pad1 becomes the current pad
	etahistoTunePUP->SetStats(111);          // No statistics on upper plot
	etahistoTunePUP->Draw("E1");               // Draw h1
	etahistoTunePLOW->Draw("E1 same");         // Draw h5 on top of h1
	// Do not draw the Y axis label on the upper plot and redraw a small
	// axis instead, in order to avoid the first label (0) to be clipped.

	TLegend * leg5 = new TLegend(0.6, 0.60, 0.98, 0.77);
	leg5->SetHeader("TuneP-2016+2017 data");
	leg5->SetFillColor(10);
	leg5->AddEntry(etahistoTunePUP,"#mu^{+} ","l");
	leg5->AddEntry(etahistoTunePLOW,"#mu^{-}","l");
	leg5 -> Draw();


	// lower plot will be in pad
	c->cd();          // Go back to the main canvas before defining pad5b
	TPad *pad5b = new TPad("pad5b", "pad5b", 0, 0, 1, 0.16);
	pad5b->SetTopMargin(0);
	pad5b->SetBottomMargin(0.2);
	pad5b->SetGridx(); // vertical grid
	pad5b->Draw();
	pad5b->cd();       // pad2 becomes the current pad
	// Define the ratio plot
	TH1F *h5 = (TH1F*)etahistoTunePUP->Clone("h5");
	h5->SetLineColor(kBlack);
	h5->Sumw2();
	h5->SetStats(0);      // No statistics on lower plot
	h5->Divide(etahistoTunePUP,etahistoTunePLOW ,1,1,"B");
	h5->SetMarkerStyle(21);
	h5->Draw("E1");       // Draw the ratio plot
	// h1 settings
	etahistoTunePUP->SetXTitle("#eta-TuneP lower leg tag");
	etahistoTunePUP->SetLineColor(kBlue+1);
	etahistoTunePUP->SetLineWidth(2);
	// Y axis h1 plot settings
	etahistoTunePUP->GetYaxis()->SetTitleSize(20);
	etahistoTunePUP->GetYaxis()->SetTitleFont(43);
	etahistoTunePUP->GetYaxis()->SetTitleOffset(1.);
	// h5 settings
	etahistoTunePLOW->SetLineColor(kRed);
	etahistoTunePLOW->SetLineWidth(2);
	etahistoTunePUP->GetXaxis()->SetTitleSize(0.05);


	// Ratio plot (h5) settings
	h5->SetTitle(""); // Remove the ratio title
	// Y axis ratio plot settings
	h5->GetYaxis()->SetTitle("mu+/mu-");
	h5->GetYaxis()->SetNdivisions(505);
	h5->GetYaxis()->SetTitleSize(20);
	h5->GetYaxis()->SetTitleFont(43);
	h5->GetYaxis()->SetTitleOffset(1.);
	h5->GetYaxis()->SetLabelFont(43); // Absolute font size in pixel (precision 3)
	h5->GetYaxis()->SetLabelSize(15);
	// X axis ratio plot settings
	h5->GetXaxis()->SetTitleSize(20);
	h5->GetXaxis()->SetTitleFont(43);
	h5->GetXaxis()->SetTitleOffset(4.);
	h5->GetXaxis()->SetLabelFont(43); // Absolute font size in pixel (precision 3)
	h5->GetXaxis()->SetLabelSize(15);
	h5->SetStats(0);      // No statistics on lower plot

	// Define the Canvas
	/*TCanvas *c = new TCanvas("c", "canvas", 1200, 800);
	c->Divide(3,2);
	c->cd(1);
	// Upper plot will be in pad0a
	TPad *pad0a = new TPad("pad0a", "pad0a", 0, 0.3, 1, 1.0);
	//pad0a->SetBottomMargin(0); // Upper and lower plot are joined
	pad0a->SetGridx();         // Vertical grid
	pad0a->Draw();             // Draw the upper pad: pad1
	pad0a->cd();               // pad1 becomes the current pad
	//etahistoGlbUP->SetStats(0);          // No statistics on upper plot
	etahistoGlbUP->Draw("E1");               // Draw h1
	etahistoGlbLOW->Draw("E1 same");         // Draw h2 on top of h1
	// Do not draw the Y axis label on the upper plot and redraw a small
	// axis instead, in order to avoid the first label (0) to be clipped.

	TLegend * leg0 = new TLegend(0.6, 0.60, 0.98, 0.77);
	leg0->SetHeader("Global fit");
	leg0->SetFillColor(10);
	leg0->AddEntry(etahistoGlbUP,"#mu^{+} ","l");
	leg0->AddEntry(etahistoGlbLOW,"#mu^{-}","l");
	leg0 -> Draw();


	// lower plot will be in pad
	c->cd();          // Go back to the main canvas before defining pad0b
	TPad *pad0b = new TPad("pad0b", "pad0b", 0, 0.56, 0.33, 0.66);
	pad0b->SetTopMargin(0);
	pad0b->SetBottomMargin(0.2);
	pad0b->SetGridx(); // vertical grid
	pad0b->Draw();
	pad0b->cd();       // pad2 becomes the current pad
	// Define the ratio plot
	TH1F *h0 = (TH1F*)etahistoGlbUP->Clone("h0");
	h0->SetLineColor(kBlack);
	h0->Sumw2();
	h0->SetStats(0);      // No statistics on lower plot
	h0->Divide(etahistoGlbUP,etahistoGlbLOW ,1,1,"B");
	h0->SetMarkerStyle(21);
	h0->Draw("E1");       // Draw the ratio plot
	// h1 settings
	etahistoGlbUP->SetXTitle("#eta ");
	etahistoGlbUP->SetLineColor(kBlue+1);
	etahistoGlbUP->SetLineWidth(2);
	// Y axis h1 plot settings
	etahistoGlbUP->GetYaxis()->SetTitleSize(20);
	etahistoGlbUP->GetYaxis()->SetTitleFont(43);
	etahistoGlbUP->GetYaxis()->SetTitleOffset(1.55);
	// h2 settings
	etahistoGlbLOW->SetLineColor(kRed);
	etahistoGlbLOW->SetLineWidth(2);
	// Ratio plot (h0) settings
	h0->SetTitle(""); // Remove the ratio title
	// Y axis ratio plot settings
	h0->GetYaxis()->SetTitle("mu+/mu-");
	h0->GetYaxis()->SetNdivisions(505);
	h0->GetYaxis()->SetTitleSize(20);
	h0->GetYaxis()->SetTitleFont(43);
	h0->GetYaxis()->SetTitleOffset(1.55);
	h0->GetYaxis()->SetLabelFont(43); // Absolute font size in pixel (precision 3)
	h0->GetYaxis()->SetLabelSize(15);
	// X axis ratio plot settings
	h0->GetXaxis()->SetTitleSize(20);
	h0->GetXaxis()->SetTitleFont(43);
	h0->GetXaxis()->SetTitleOffset(4.);
	h0->GetXaxis()->SetLabelFont(43); // Absolute font size in pixel (precision 3)
	h0->GetXaxis()->SetLabelSize(15);
	h0->SetStats(0);      // No statistics on lower plot
	//h0->GetYaxis()->SetRangeUser(0.8,1.2);

	c->cd(2);

	// Upper plot will be in pad1a
	TPad *pad1a = new TPad("pad1a", "pad1a", 0, 0.3, 1, 1.0);
	//pad1a->SetBottomMargin(0); // Upper and lower plot are joined
	pad1a->SetGridx();         // Vertical grid
	pad1a->Draw();             // Draw the upper pad: pad1
	pad1a->cd();               // pad1 becomes the current pad
	//etahistoTkonlyUP->SetStats(0);          // No statistics on upper plot
	etahistoTkonlyUP->Draw("E1");               // Draw h1
	etahistoTkonlyLOW->Draw("E1 same");         // Draw h2 on top of h1
	// Do not draw the Y axis label on the upper plot and redraw a small
	// axis instead, in order to avoid the first label (0) to be clipped.

	TLegend * leg1 = new TLegend(0.6, 0.60, 0.98, 0.77);
	leg1->SetHeader("Tkonly");
	leg1->SetFillColor(10);
	leg1->AddEntry(etahistoTkonlyUP,"#mu^{+} ","l");
	leg1->AddEntry(etahistoTkonlyLOW,"#mu^{-}","l");
	leg1 -> Draw();


	// lower plot will be in pad
	c->cd();          // Go back to the main canvas before defining pad1b
	TPad *pad1b = new TPad("pad1b", "pad1b", 0.33, 0.56, 0.66,0.66);
	pad1b->SetTopMargin(0);
	pad1b->SetBottomMargin(0.2);
	pad1b->SetGridx(); // vertical grid
	pad1b->Draw();
	pad1b->cd();       // pad2 becomes the current pad
	// Define the ratio plot
	TH1F *h1 = (TH1F*)etahistoTkonlyUP->Clone("h1");
	h1->SetLineColor(kBlack);
	h1->Sumw2();
	h1->SetStats(0);      // No statistics on lower plot
	h1->Divide(etahistoTkonlyUP,etahistoTkonlyLOW ,1,1,"B");
	h1->SetMarkerStyle(21);
	h1->Draw("E1");       // Draw the ratio plot
	// h1 settings
	etahistoTkonlyUP->SetXTitle("#eta ");
	etahistoTkonlyUP->SetLineColor(kBlue+1);
	etahistoTkonlyUP->SetLineWidth(2);
	// Y axis h1 plot settings
	etahistoTkonlyUP->GetYaxis()->SetTitleSize(20);
	etahistoTkonlyUP->GetYaxis()->SetTitleFont(43);
	etahistoTkonlyUP->GetYaxis()->SetTitleOffset(1.55);
	// h2 settings
	etahistoTkonlyLOW->SetLineColor(kRed);
	etahistoTkonlyLOW->SetLineWidth(2);
	// Ratio plot (h1) settings
	h1->SetTitle(""); // Remove the ratio title
	// Y axis ratio plot settings
	h1->GetYaxis()->SetTitle("mu+/mu-");
	h1->GetYaxis()->SetNdivisions(505);
	h1->GetYaxis()->SetTitleSize(20);
	h1->GetYaxis()->SetTitleFont(43);
	h1->GetYaxis()->SetTitleOffset(1.55);
	h1->GetYaxis()->SetLabelFont(43); // Absolute font size in pixel (precision 3)
	h1->GetYaxis()->SetLabelSize(15);
	// X axis ratio plot settings
	h1->GetXaxis()->SetTitleSize(20);
	h1->GetXaxis()->SetTitleFont(43);
	h1->GetXaxis()->SetTitleOffset(4.);
	h1->GetXaxis()->SetLabelFont(43); // Absolute font size in pixel (precision 3)
	h1->GetXaxis()->SetLabelSize(15);
	h1->SetStats(0);      // No statistics on lower plot
	//h1->GetYaxis()->SetRangeUser(0.8,1.2);

	c->cd(3);
	// Upper plot will be in pad2ba
	TPad *pad2a = new TPad("pad2a", "pad2a", 0, 0.3, 1, 1.0);
	//pad2a->SetBottomMargin(0); // Upper and lower plot are joined
	pad2a->SetGridx();         // Vertical grid
	pad2a->Draw();             // Draw the upper pad: pad1
	pad2a->cd();               // pad1 becomes the current pad
	//etahistoTpfmsUP->SetStats(0);          // No statistics on upper plot
	etahistoTpfmsUP->Draw("E1");               // Draw h1
	etahistoTpfmsLOW->Draw("E1 same");         // Draw h2 on top of h1
	// Do not draw the Y axis label on the upper plot and redraw a small
	// axis instead, in order to avoid the first label (0) to be clipped.

	TLegend * leg2 = new TLegend(0.6, 0.60, 0.98, 0.77);
	leg2->SetHeader("TPFMS");
	leg2->SetFillColor(10);
	leg2->AddEntry(etahistoTpfmsUP,"#mu^{+} ","l");
	leg2->AddEntry(etahistoTpfmsLOW,"#mu^{-}","l");
	leg2 -> Draw();


	// lower plot will be in pad
	c->cd();          // Go back to the main canvas before defining pad2b
	TPad *pad2b = new TPad("pad2b", "pad2b", 0.66, 0.56, 1, 0.66);
	pad2b->SetTopMargin(0);
	pad2b->SetBottomMargin(0.2);
	pad2b->SetGridx(); // vertical grid
	pad2b->Draw();
	pad2b->cd();       // pad2 becomes the current pad
	// Define the ratio plot
	TH1F *h2 = (TH1F*)etahistoTpfmsUP->Clone("h2");
	h2->SetLineColor(kBlack);
	h2->Sumw2();
	h2->SetStats(0);      // No statistics on lower plot
	h2->Divide(etahistoTpfmsUP,etahistoTpfmsLOW ,1,1,"B");
	h2->SetMarkerStyle(21);
	h2->Draw("E1");       // Draw the ratio plot
	// h1 settings
	etahistoTpfmsUP->SetXTitle("#eta ");
	etahistoTpfmsUP->SetLineColor(kBlue+1);
	etahistoTpfmsUP->SetLineWidth(2);
	// Y axis h1 plot settings
	etahistoTpfmsUP->GetYaxis()->SetTitleSize(20);
	etahistoTpfmsUP->GetYaxis()->SetTitleFont(43);
	etahistoTpfmsUP->GetYaxis()->SetTitleOffset(1.55);
	// h2 settings
	etahistoTpfmsLOW->SetLineColor(kRed);
	etahistoTpfmsLOW->SetLineWidth(2);
	// Ratio plot (h2) settings
	h2->SetTitle(""); // Remove the ratio title
	// Y axis ratio plot settings
	h2->GetYaxis()->SetTitle("mu+/mu-");
	h2->GetYaxis()->SetNdivisions(505);
	h2->GetYaxis()->SetTitleSize(20);
	h2->GetYaxis()->SetTitleFont(43);
	h2->GetYaxis()->SetTitleOffset(1.55);
	h2->GetYaxis()->SetLabelFont(43); // Absolute font size in pixel (precision 3)
	h2->GetYaxis()->SetLabelSize(15);
	// X axis ratio plot settings
	h2->GetXaxis()->SetTitleSize(20);
	h2->GetXaxis()->SetTitleFont(43);
	h2->GetXaxis()->SetTitleOffset(4.);
	h2->GetXaxis()->SetLabelFont(43); // Absolute font size in pixel (precision 3)
	h2->GetXaxis()->SetLabelSize(15);
	h2->SetStats(0);      // No statistics on lower plot
	//h2->GetYaxis()->SetRangeUser(0.8,1.2);

	c->cd(4);
	// Upper plot will be in pad3a
	TPad *pad3a = new TPad("pad3a", "pad3a",0, 0.3, 1, 1.0);
	//pad3a->SetBottomMargin(0); // Upper and lower plot are joined
	pad3a->SetGridx();         // Vertical grid
	pad3a->Draw();             // Draw the upper pad: pad1
	pad3a->cd();               // pad1 becomes the current pad
	//etahistoPickyUP->SetStats(0);          // No statistics on upper plot
	etahistoPickyUP->Draw("E1");               // Draw h1
	etahistoPickyLOW->Draw("E1 same");         // Draw h3 on top of h1
	// Do not draw the Y axis label on the upper plot and redraw a small
	// axis instead, in order to avoid the first label (0) to be clipped.

	TLegend * leg3 = new TLegend(0.6, 0.60, 0.98, 0.77);
	leg3->SetHeader("Picky fit");
	leg3->SetFillColor(10);
	leg3->AddEntry(etahistoPickyUP,"#mu^{+} ","l");
	leg3->AddEntry(etahistoPickyLOW,"#mu^{-}","l");
	leg3 -> Draw();


	// lower plot will be in pad
	c->cd();          // Go back to the main canvas before defining pad3b
	TPad *pad3b = new TPad("pad3b", "pad3b", 0, 0.06, 0.33, 0.16);
	pad3b->SetTopMargin(0);
	pad3b->SetBottomMargin(0.2);
	pad3b->SetGridx(); // vertical grid
	pad3b->Draw();
	pad3b->cd();       // pad2 becomes the current pad
	// Define the ratio plot
	TH1F *h3 = (TH1F*)etahistoPickyUP->Clone("h3");
	h3->SetLineColor(kBlack);
	h3->Sumw2();
	h3->SetStats(0);      // No statistics on lower plot
	h3->Divide(etahistoPickyUP,etahistoPickyLOW ,1,1,"B");
	h3->SetMarkerStyle(21);
	h3->Draw("E1");       // Draw the ratio plot
	// h1 settings
	etahistoPickyUP->SetXTitle("#eta ");
	etahistoPickyUP->SetLineColor(kBlue+1);
	etahistoPickyUP->SetLineWidth(2);
	// Y axis h1 plot settings
	etahistoPickyUP->GetYaxis()->SetTitleSize(20);
	etahistoPickyUP->GetYaxis()->SetTitleFont(43);
	etahistoPickyUP->GetYaxis()->SetTitleOffset(1.55);
	// h3 settings
	etahistoPickyLOW->SetLineColor(kRed);
	etahistoPickyLOW->SetLineWidth(2);
	// Ratio plot (h3) settings
	h3->SetTitle(""); // Remove the ratio title
	// Y axis ratio plot settings
	h3->GetYaxis()->SetTitle("mu+/mu-");
	h3->GetYaxis()->SetNdivisions(505);
	h3->GetYaxis()->SetTitleSize(20);
	h3->GetYaxis()->SetTitleFont(43);
	h3->GetYaxis()->SetTitleOffset(1.55);
	h3->GetYaxis()->SetLabelFont(43); // Absolute font size in pixel (precision 3)
	h3->GetYaxis()->SetLabelSize(15);
	// X axis ratio plot settings
	h3->GetXaxis()->SetTitleSize(20);
	h3->GetXaxis()->SetTitleFont(43);
	h3->GetXaxis()->SetTitleOffset(4.);
	h3->GetXaxis()->SetLabelFont(43); // Absolute font size in pixel (precision 3)
	h3->GetXaxis()->SetLabelSize(15);
	h3->SetStats(0);      // No statistics on lower plot
	//h3->GetYaxis()->SetRangeUser(0.8,1.2);

	c->cd(5);
	// Upper plot will be in pad4a
	TPad *pad4a = new TPad("pad4a", "pad4a", 0, 0.3, 1, 1.0);
	//pad4a->SetBottomMargin(0); // Upper and lower plot are joined
	pad4a->SetGridx();         // Vertical grid
	pad4a->Draw();             // Draw the upper pad: pad1
	pad4a->cd();               // pad1 becomes the current pad
	//etahistoDYTUP->SetStats(0);          // No statistics on upper plot
	etahistoDYTUP->Draw("E1");               // Draw h1
	etahistoDYTLOW->Draw("E1 same");         // Draw h4 on top of h1
	// Do not draw the Y axis label on the upper plot and redraw a small
	// axis instead, in order to avoid the first label (0) to be clipped.

	TLegend * leg4 = new TLegend(0.6, 0.60, 0.98, 0.77);
	leg4->SetHeader("DYT");
	leg4->SetFillColor(10);
	leg4->AddEntry(etahistoDYTUP,"#mu^{+} ","l");
	leg4->AddEntry(etahistoDYTLOW,"#mu^{-}","l");
	leg4 -> Draw();


	// lower plot will be in pad
	c->cd();          // Go ba to the main canvas before defining pad4b
	TPad *pad4b = new TPad("pad4b", "pad4b", 0.33, 0.06, 0.66,0.16);
	pad4b->SetTopMargin(0);
	pad4b->SetBottomMargin(0.2);
	pad4b->SetGridx(); // vertical grid
	pad4b->Draw();
	pad4b->cd();       // pad2 becomes the current pad
	// Define the ratio plot
	TH1F *h4 = (TH1F*)etahistoDYTUP->Clone("h4");
	h4->SetLineColor(kBlack);
	h4->Sumw2();
	h4->SetStats(0);      // No statistics on lower plot
	h4->Divide(etahistoDYTUP,etahistoDYTLOW ,1,1,"B");
	h4->SetMarkerStyle(21);
	h4->Draw("E1");       // Draw the ratio plot
	// h1 settings
	etahistoDYTUP->SetXTitle("#eta ");
	etahistoDYTUP->SetLineColor(kBlue+1);
	etahistoDYTUP->SetLineWidth(2);
	// Y axis h1 plot settings
	etahistoDYTUP->GetYaxis()->SetTitleSize(20);
	etahistoDYTUP->GetYaxis()->SetTitleFont(43);
	etahistoDYTUP->GetYaxis()->SetTitleOffset(1.55);
	// h4 settings
	etahistoDYTLOW->SetLineColor(kRed);
	etahistoDYTLOW->SetLineWidth(2);
	// Ratio plot (h4) settings
	h4->SetTitle(""); // Remove the ratio title
	// Y axis ratio plot settings
	h4->GetYaxis()->SetTitle("mu+/mu-");
	h4->GetYaxis()->SetNdivisions(505);
	h4->GetYaxis()->SetTitleSize(20);
	h4->GetYaxis()->SetTitleFont(43);
	h4->GetYaxis()->SetTitleOffset(1.55);
	h4->GetYaxis()->SetLabelFont(43); // Absolute font size in pixel (precision 3)
	h4->GetYaxis()->SetLabelSize(15);
	// X axis ratio plot settings
	h4->GetXaxis()->SetTitleSize(20);
	h4->GetXaxis()->SetTitleFont(43);
	h4->GetXaxis()->SetTitleOffset(4.);
	h4->GetXaxis()->SetLabelFont(43); // Absolute font size in pixel (precision 3)
	h4->GetXaxis()->SetLabelSize(15);
	h4->SetStats(0);      // No statistics on lower plot
	//h4->GetYaxis()->SetRangeUser(0.8,1.2);

	c->cd(6);
	// Upper plot will be in pad5a
	TPad *pad5a = new TPad("pad3a", "pad3a", 0, 0.3, 1, 1.0);
	//pad5a->SetBottomMargin(0); // Upper and lower plot are joined
	pad5a->SetGridx();         // Vertical grid
	pad5a->Draw();             // Draw the upper pad: pad1
	pad5a->cd();               // pad1 becomes the current pad
	//etahistoTunePUP->SetStats(0);          // No statistics on upper plot
	etahistoTunePUP->Draw("E1");               // Draw h1
	etahistoTunePLOW->Draw("E1 same");         // Draw h5 on top of h1
	// Do not draw the Y axis label on the upper plot and redraw a small
	// axis instead, in order to avoid the first label (0) to be clipped.

	TLegend * leg5 = new TLegend(0.6, 0.60, 0.98, 0.77);
	leg5->SetHeader("TuneP");
	leg5->SetFillColor(10);
	leg5->AddEntry(etahistoTunePUP,"#mu^{+} ","l");
	leg5->AddEntry(etahistoTunePLOW,"#mu^{-}","l");
	leg5 -> Draw();


	// lower plot will be in pad
	c->cd();          // Go back to the main canvas before defining pad5b
	TPad *pad5b = new TPad("pad5b", "pad5b", 0.66, 0.06, 1, 0.16);
	pad5b->SetTopMargin(0);
	pad5b->SetBottomMargin(0.2);
	pad5b->SetGridx(); // vertical grid
	pad5b->Draw();
	pad5b->cd();       // pad2 becomes the current pad
	// Define the ratio plot
	TH1F *h5 = (TH1F*)etahistoTunePUP->Clone("h5");
	h5->SetLineColor(kBlack);
	h5->Sumw2();
	h5->SetStats(0);      // No statistics on lower plot
	h5->Divide(etahistoTunePUP,etahistoTunePLOW ,1,1,"B");
	h5->SetMarkerStyle(21);
	h5->Draw("E1");       // Draw the ratio plot
	// h1 settings
	etahistoTunePUP->SetXTitle("#eta ");
	etahistoTunePUP->SetLineColor(kBlue+1);
	etahistoTunePUP->SetLineWidth(2);
	// Y axis h1 plot settings
	etahistoTunePUP->GetYaxis()->SetTitleSize(20);
	etahistoTunePUP->GetYaxis()->SetTitleFont(43);
	etahistoTunePUP->GetYaxis()->SetTitleOffset(1.55);
	// h5 settings
	etahistoTunePLOW->SetLineColor(kRed);
	etahistoTunePLOW->SetLineWidth(2);
	// Ratio plot (h5) settings
	h5->SetTitle(""); // Remove the ratio title
	// Y axis ratio plot settings
	h5->GetYaxis()->SetTitle("mu+/mu-");
	h5->GetYaxis()->SetNdivisions(505);
	h5->GetYaxis()->SetTitleSize(20);
	h5->GetYaxis()->SetTitleFont(43);
	h5->GetYaxis()->SetTitleOffset(1.55);
	h5->GetYaxis()->SetLabelFont(43); // Absolute font size in pixel (precision 3)
	h5->GetYaxis()->SetLabelSize(15);
	// X axis ratio plot settings
	h5->GetXaxis()->SetTitleSize(20);
	h5->GetXaxis()->SetTitleFont(43);
	h5->GetXaxis()->SetTitleOffset(4.);
	h5->GetXaxis()->SetLabelFont(43); // Absolute font size in pixel (precision 3)
	h5->GetXaxis()->SetLabelSize(15);
	h5->SetStats(0);      // No statistics on lower plot
	//h5->GetYaxis()->SetRangeUser(0.8,1.2);*/



	/*TFile *fout = new TFile("testeta.root ","RECREATE");

	  etahistoGlb->Write();
	  etahistoStalone->Write();

	  fout->Close();*/

}// end of function