etaDistri.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 "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 etaDistri(){

	/*gStyle->SetPadBorderMode(0);
	  gStyle->SetCanvasColor(10);
	  gStyle->SetCanvasBorderMode(0);
	  gStyle->SetCanvasBorderSize(0);
	  gStyle->SetFrameBorderMode(0);
	  gStyle->SetFrameLineColor(0);
	  gStyle->SetFillColor(0);*/



	//DATA
	// 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/cosmic_data_2016_2017/cosmic_data_16_17.root");
	//treedata->Add("ROOT_Files/MC2016/cosmic_MC_2016.root");
	//treedata->Add("ROOT_Files/MC2016/cosmic_MC_2016.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("theta",theta);
	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.){


			//Pt Bin
			//if (unprop_pt[8][1]>100 && unprop_pt[8][1]<500){

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

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

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

	// this is the root file you have input with your tree inside
	MCtreedata->Add("ROOT_Files/MC2016/MC_2016_P100.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.
	MCtreedata->SetBranchAddress("charge",charge);
	//MCtreedata->SetBranchAddress("theta",theta);
	MCtreedata->SetBranchAddress("unprop_charge",unprop_charge);
	MCtreedata->SetBranchAddress("unprop_theta",unprop_theta);
	MCtreedata->SetBranchAddress("pixel_hits",pixel_hits);
	MCtreedata->SetBranchAddress("strip_layers",strip_layers);
	MCtreedata->SetBranchAddress("unprop_pt",unprop_pt);


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

	Long64_t MCnentries = MCtreedata->GetEntries();


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

	MCtreedata->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.){



	etahistoTunePLOW->Fill(eta[8][1]);
	}
	}// end loop p

	Double_t norm = 1;
	Double_t scale_DATA = norm/(etahistoTunePUP->Integral());
	etahistoTunePUP->Scale(scale_DATA);
	Double_t scale_MC = norm/(etahistoTunePLOW->Integral());
	etahistoTunePLOW->Scale(scale_MC);*/
	TCanvas *c = new TCanvas("c", "canvas", 700, 600);

	// Upper plot will be in pad5a
	TPad *pad5a = new TPad("pad3a", "pad3a", 0, 0.18, 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->SetTitle("TuneP-2016+2017 data-P500");
	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,"Upper leg ","l");
	leg5->AddEntry(etahistoTunePLOW,"Lower leg","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.17);
	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 ");
	etahistoTunePUP->SetLineColor(kBlue+1);
	etahistoTunePUP->SetLineWidth(2);
	etahistoTunePUP->GetXaxis()->SetTitleSize(0.05);


	// 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);
	// Ratio plot (h5) settings
	h5->SetTitle(""); // Remove the ratio title
	// Y axis ratio plot settings
	h5->GetYaxis()->SetTitle("up/low");
	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
	h5->GetYaxis()->SetRangeUser(0.8,1.2);






	// 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,"Upper leg ","l");
	leg0->AddEntry(etahistoGlbLOW,"Lower leg","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("up/low");
	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,"Upper leg ","l");
	leg1->AddEntry(etahistoTkonlyLOW,"Lower leg","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("up/low");
	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,"Upper leg ","l");
	leg2->AddEntry(etahistoTpfmsLOW,"Lower leg","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("up/low");
	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,"Upper leg ","l");
	leg3->AddEntry(etahistoPickyLOW,"Lower leg","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("up/low");
	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,"Upper leg ","l");
	leg4->AddEntry(etahistoDYTLOW,"Lower leg","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("up/low");
	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,"Upper leg ","l");
	leg5->AddEntry(etahistoTunePLOW,"Lower leg","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("up/low");
	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