phiDistri.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 phi[n_tracks][n_track_pos];
short charge[n_tracks][n_track_pos];

float unprop_phi[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 *phihistoGlbUP =new TH1F("phihistoGlb","",15,-3,0);
TH1F *phihistoGlbLOW =new TH1F("phihistoGlbLOW","",15,-3,0);

TH1F *phihistoTkonlyUP =new TH1F("phihistoTkonly","",15,-3,0);
  TH1F *phihistoTkonlyLOW =new TH1F("phihistoTkonlyLOW","",15,-3,0);

  TH1F *phihistoTpfmsUP =new TH1F("phihistoTpfms","",15,-3,0);
  TH1F *phihistoTpfmsLOW =new TH1F("phihistoTpfmsLOW","",15,-3,0);

  TH1F *phihistoPickyUP =new TH1F("phihistoPicky","",15,-3,0);
  TH1F *phihistoPickyLOW =new TH1F("phihistoPickyLOW","",15,-3,0);

  TH1F *phihistoDYTUP =new TH1F("phihistoDYT","",15,-3,0);
  TH1F *phihistoDYTLOW =new TH1F("phihistoDYTLOW","",15,-3,0);

TH1F *phihistoTunePUP =new TH1F("phihistoTuneP","",15,-3,0);
TH1F *phihistoTunePLOW =new TH1F("phihistoTunePLOW","",15,-3,0);


//Begining of teh function
void phiDistri(){



	// 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/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, dilephion mass in my case for DY events.
	treedata->SetBranchAddress("charge",charge);
	treedata->SetBranchAddress("phi",phi);
	treedata->SetBranchAddress("unprop_charge",unprop_charge);
	treedata->SetBranchAddress("unprop_phi",unprop_phi);
	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
		//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]>500 /*&& unprop_pt[8][1]<500*/){

		//
		//Algorithm choice
		//
		//
		//		//global
					phihistoGlbUP->Fill(unprop_phi[0][0]);
					phihistoGlbLOW->Fill(unprop_phi[0][1]);
				//Tkonly
					phihistoTkonlyUP->Fill(unprop_phi[2][0]);
					phihistoTkonlyLOW->Fill(unprop_phi[2][1]);
				//Tpfms
					phihistoTpfmsUP->Fill(unprop_phi[3][0]);
					phihistoTpfmsLOW->Fill(unprop_phi[3][1]);
				//Picky
					phihistoPickyUP->Fill(unprop_phi[4][0]);
					phihistoPickyLOW->Fill(unprop_phi[4][1]);
				//DYT
					phihistoDYTUP->Fill(unprop_phi[5][0]);
					phihistoDYTLOW->Fill(unprop_phi[5][1]);
				//TuneP
					phihistoTunePUP->Fill(unprop_phi[8][0]);
					phihistoTunePLOW->Fill(unprop_phi[8][1]);
	}// end loop p

//}	
}

	//Define the Canvas
	TCanvas *c = new TCanvas("c", "canvas", 700, 600);

			// Upper plot will be in pad5a
			TPad *pad5a = new TPad("pad5a", "pad5a", 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
			phihistoTunePUP->SetStats(111);          // No statistics on upper plot
			phihistoTunePUP->Draw("E1");               // Draw h1
			phihistoTunePLOW->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(phihistoTunePUP,"Upper leg ","l");
			leg5->AddEntry(phihistoTunePLOW,"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*)phihistoTunePUP->Clone("h5");
	h5->SetLineColor(kBlack);
	h5->Sumw2();
			h5->SetStats(0);      // No statistics on lower plot
			h5->Divide(phihistoTunePUP,phihistoTunePLOW ,1,1,"B");
			h5->SetMarkerStyle(21);
			h5->Draw("E1");       // Draw the ratio plot
			// h1 settings
			phihistoTunePUP->SetXTitle("#phi-TuneP ");
			phihistoTunePUP->SetLineColor(kBlue+1);
			phihistoTunePUP->SetLineWidth(2);
			phihistoTunePUP->GetXaxis()->SetTitleSize(0.05);

			// Y axis h1 plot settings
			phihistoTunePUP->GetYaxis()->SetTitleSize(20);
			phihistoTunePUP->GetYaxis()->SetTitleFont(43);
			phihistoTunePUP->GetYaxis()->SetTitleOffset(1.);
			// h5 settings
			phihistoTunePLOW->SetLineColor(kRed);
			phihistoTunePLOW->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
			//phihistoGlbUP->SetStats(0);          // No statistics on upper plot
			phihistoGlbUP->Draw("E1");               // Draw h1
			phihistoGlbLOW->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(phihistoGlbUP,"Upper leg ","l");
			leg0->AddEntry(phihistoGlbLOW,"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*)phihistoGlbUP->Clone("h0");
			h0->SetLineColor(kBlack);
			h0->Sumw2();
			h0->SetStats(0);      // No statistics on lower plot
			h0->Divide(phihistoGlbUP,phihistoGlbLOW ,1,1,"B");
			h0->SetMarkerStyle(21);
			h0->Draw("E1");       // Draw the ratio plot
			// h1 settings
			phihistoGlbUP->SetXTitle("#phi ");
			phihistoGlbUP->SetLineColor(kBlue+1);
			phihistoGlbUP->SetLineWidth(2);
			// Y axis h1 plot settings
			phihistoGlbUP->GetYaxis()->SetTitleSize(20);
			phihistoGlbUP->GetYaxis()->SetTitleFont(43);
			phihistoGlbUP->GetYaxis()->SetTitleOffset(1.55);
			// h2 settings
			phihistoGlbLOW->SetLineColor(kRed);
			phihistoGlbLOW->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
			//phihistoTkonlyUP->SetStats(0);          // No statistics on upper plot
			phihistoTkonlyUP->Draw("E1");               // Draw h1
			phihistoTkonlyLOW->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(phihistoTkonlyUP,"Upper leg ","l");
			leg1->AddEntry(phihistoTkonlyLOW,"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*)phihistoTkonlyUP->Clone("h1");
			h1->SetLineColor(kBlack);
			h1->Sumw2();
			h1->SetStats(0);      // No statistics on lower plot
			h1->Divide(phihistoTkonlyUP,phihistoTkonlyLOW ,1,1,"B");
			h1->SetMarkerStyle(21);
			h1->Draw("E1");       // Draw the ratio plot
			// h1 settings
			phihistoTkonlyUP->SetXTitle("#phi ");
			phihistoTkonlyUP->SetLineColor(kBlue+1);
			phihistoTkonlyUP->SetLineWidth(2);
			// Y axis h1 plot settings
			phihistoTkonlyUP->GetYaxis()->SetTitleSize(20);
			phihistoTkonlyUP->GetYaxis()->SetTitleFont(43);
			phihistoTkonlyUP->GetYaxis()->SetTitleOffset(1.55);
			// h2 settings
			phihistoTkonlyLOW->SetLineColor(kRed);
			phihistoTkonlyLOW->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
			//phihistoTpfmsUP->SetStats(0);          // No statistics on upper plot
			phihistoTpfmsUP->Draw("E1");               // Draw h1
			phihistoTpfmsLOW->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(phihistoTpfmsUP,"Upper leg ","l");
			leg2->AddEntry(phihistoTpfmsLOW,"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*)phihistoTpfmsUP->Clone("h2");
			h2->SetLineColor(kBlack);
			h2->Sumw2();
			h2->SetStats(0);      // No statistics on lower plot
			h2->Divide(phihistoTpfmsUP,phihistoTpfmsLOW ,1,1,"B");
			h2->SetMarkerStyle(21);
			h2->Draw("E1");       // Draw the ratio plot
			// h1 settings
			phihistoTpfmsUP->SetXTitle("#phi ");
			phihistoTpfmsUP->SetLineColor(kBlue+1);
			phihistoTpfmsUP->SetLineWidth(2);
			// Y axis h1 plot settings
			phihistoTpfmsUP->GetYaxis()->SetTitleSize(20);
			phihistoTpfmsUP->GetYaxis()->SetTitleFont(43);
			phihistoTpfmsUP->GetYaxis()->SetTitleOffset(1.55);
			// h2 settings
			phihistoTpfmsLOW->SetLineColor(kRed);
			phihistoTpfmsLOW->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
			//phihistoPickyUP->SetStats(0);          // No statistics on upper plot
			phihistoPickyUP->Draw("E1");               // Draw h1
			phihistoPickyLOW->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(phihistoPickyUP,"Upper leg ","l");
			leg3->AddEntry(phihistoPickyLOW,"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*)phihistoPickyUP->Clone("h3");
			h3->SetLineColor(kBlack);
			h3->Sumw2();
			h3->SetStats(0);      // No statistics on lower plot
			h3->Divide(phihistoPickyUP,phihistoPickyLOW ,1,1,"B");
			h3->SetMarkerStyle(21);
			h3->Draw("E1");       // Draw the ratio plot
			// h1 settings
			phihistoPickyUP->SetXTitle("#phi ");
			phihistoPickyUP->SetLineColor(kBlue+1);
			phihistoPickyUP->SetLineWidth(2);
			// Y axis h1 plot settings
			phihistoPickyUP->GetYaxis()->SetTitleSize(20);
			phihistoPickyUP->GetYaxis()->SetTitleFont(43);
			phihistoPickyUP->GetYaxis()->SetTitleOffset(1.55);
			// h3 settings
			phihistoPickyLOW->SetLineColor(kRed);
			phihistoPickyLOW->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
			//phihistoDYTUP->SetStats(0);          // No statistics on upper plot
			phihistoDYTUP->Draw("E1");               // Draw h1
			phihistoDYTLOW->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(phihistoDYTUP,"Upper leg ","l");
			leg4->AddEntry(phihistoDYTLOW,"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*)phihistoDYTUP->Clone("h4");
			h4->SetLineColor(kBlack);
			h4->Sumw2();
			h4->SetStats(0);      // No statistics on lower plot
			h4->Divide(phihistoDYTUP,phihistoDYTLOW ,1,1,"B");
			h4->SetMarkerStyle(21);
			h4->Draw("E1");       // Draw the ratio plot
			// h1 settings
			phihistoDYTUP->SetXTitle("#phi ");
			phihistoDYTUP->SetLineColor(kBlue+1);
			phihistoDYTUP->SetLineWidth(2);
			// Y axis h1 plot settings
			phihistoDYTUP->GetYaxis()->SetTitleSize(20);
			phihistoDYTUP->GetYaxis()->SetTitleFont(43);
			phihistoDYTUP->GetYaxis()->SetTitleOffset(1.55);
			// h4 settings
			phihistoDYTLOW->SetLineColor(kRed);
			phihistoDYTLOW->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
			//phihistoTunePUP->SetStats(0);          // No statistics on upper plot
			phihistoTunePUP->Draw("E1");               // Draw h1
			phihistoTunePLOW->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(phihistoTunePUP,"Upper leg ","l");
			leg5->AddEntry(phihistoTunePLOW,"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*)phihistoTunePUP->Clone("h5");
			h5->SetLineColor(kBlack);
			h5->Sumw2();
			h5->SetStats(0);      // No statistics on lower plot
			h5->Divide(phihistoTunePUP,phihistoTunePLOW ,1,1,"B");
			h5->SetMarkerStyle(21);
			h5->Draw("E1");       // Draw the ratio plot
			// h1 settings
			phihistoTunePUP->SetXTitle("#phi ");
			phihistoTunePUP->SetLineColor(kBlue+1);
			phihistoTunePUP->SetLineWidth(2);
			// Y axis h1 plot settings
			phihistoTunePUP->GetYaxis()->SetTitleSize(20);
			phihistoTunePUP->GetYaxis()->SetTitleFont(43);
			phihistoTunePUP->GetYaxis()->SetTitleOffset(1.55);
			// h5 settings
			phihistoTunePLOW->SetLineColor(kRed);
			phihistoTunePLOW->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("testphi.root ","RECREATE");

	  phihistoGlb->Write();
	  phihistoStalone->Write();

	  fout->Close();*/

}// end of function