#include <iostream>
#include <vector>
#include <sstream>
#include <string>

#include "TGaxis.h"
#include "TString.h"
#include "TF1.h"
#include "TH1.h"
#include "TH2.h"
#include "TH3.h"
#include "TMath.h"
#include "TTree.h"
#include "TChain.h"
#include "TFile.h"
#include "TCanvas.h"
#include "TSystem.h"
#include "TROOT.h"
#include "TGraph.h"
#include "TGraphErrors.h"
#include "TGraphAsymmErrors.h"
#include "TMultiGraph.h"
#include "TCanvas.h"
#include "TPad.h"
#include "TLegend.h"
#include "TLatex.h"
#include "TLine.h"
#include "TAxis.h"
#include "TGraph.h"
#include "TGraphErrors.h"
#include "TLorentzVector.h"
#include "TFitResult.h"
#include "TFitResultPtr.h"
#include "TMatrixDSym.h"
#include "TMatrixD.h"
#include "TArrow.h"

void RiceStyle(){

std::cout << "Welcome to Rice Heavy Ion group!! " << std::endl;

}

//make Canvas

TCanvas* makeCanvas(const char* name, const char *title, bool doLogx = false, bool doLogy = false ){
	
	// Start with a canvas
	TCanvas *canvas = new TCanvas(name,title, 1, 1 ,600 ,600 );
	// General overall stuff
	canvas->SetFillColor      (0);
	canvas->SetBorderMode     (0);
	canvas->SetBorderSize     (10);
	// Set margins to reasonable defaults
	canvas->SetLeftMargin     (0.13);
	canvas->SetRightMargin    (0.10);
	canvas->SetTopMargin      (0.10);
	canvas->SetBottomMargin   (0.13);
	// Setup a frame which makes sense
	canvas->SetFrameFillStyle (0);
	canvas->SetFrameLineStyle (0);
	canvas->SetFrameBorderMode(0);
	canvas->SetFrameBorderSize(10);
	canvas->SetFrameFillStyle (0);
	canvas->SetFrameLineStyle (0);
	canvas->SetFrameBorderMode(0);
	canvas->SetFrameBorderSize(10);

	if( doLogy == true ) gPad->SetLogy(1);
	if( doLogx == true ) gPad->SetLogx(1);
	
	gPad->SetTicks();

	return canvas;
}

TCanvas* makeMultiCanvas(const char* name, 
						 const char* title,
						 int nRows,
						 int nColumns
){

	double ratio = nRows/nColumns;

	TCanvas* canvas = new TCanvas( name, title, 1, 1, 400*nRows, 400*nColumns );
	canvas->SetFillColor      (0);
	canvas->SetBorderMode     (0);
	canvas->SetBorderSize     (10);
	// Set margins to reasonable defaults
	canvas->SetLeftMargin     (0.30);
	canvas->SetRightMargin    (0.10);
	canvas->SetTopMargin      (0.10);
	canvas->SetBottomMargin   (0.30);
	// Setup a frame which makes sense
	canvas->SetFrameFillStyle (0);
	canvas->SetFrameLineStyle (0);
	canvas->SetFrameBorderMode(0);
	canvas->SetFrameBorderSize(10);
	canvas->SetFrameFillStyle (0);
	canvas->SetFrameLineStyle (0);
	canvas->SetFrameBorderMode(0);
	canvas->SetFrameBorderSize(10);
 	
 	canvas->Divide(nRows,nColumns,0.01,0.01);

 	gPad->SetLeftMargin(0.3);
 	gPad->SetBottomMargin(0.3);
	return canvas;

}

void saveCanvas(TCanvas* c, TString dir, TString filename)
{
   TDatime* date = new TDatime();
   //c->Print(Form("../%s/%s_%d.eps",dir.Data(),filename.Data(),date->GetDate()));
   //c->Print(Form("../%s/%s_%d.gif",dir.Data(),filename.Data(),date->GetDate()));
   c->Print(Form("../%s/%s_%d.pdf",dir.Data(),filename.Data(),date->GetDate()));
   //c->Print(Form("../%s/%s_%d.C",dir.Data(),filename.Data(),date->GetDate()));
   c->Print(Form("../%s/%s_%d.png",dir.Data(),filename.Data(),date->GetDate()));
}

void initSubPad(TPad* pad, int i)
{
  //printf("Pad: %p, index: %d\n",pad,i);

  pad->cd(i);
  TPad *tmpPad = (TPad*) pad->GetPad(i);
  tmpPad->SetLeftMargin  (0.20);
  tmpPad->SetTopMargin   (0.06);
  tmpPad->SetRightMargin (0.08);
  tmpPad->SetBottomMargin(0.15);
  return;
}

vector<TPad*> makeMultiPad(const int num){//we only have 4,6,8 options for now

	cout << "num: "<< num << endl;
	vector<TPad*> temp;

	TPad* pad[ num ];

	double setting1[] = {0.0, 0.35, 0.56, 1.0};
	double setting2[] = {0.0, 0.35, 0.40, 0.7, 1.0 };
	double setting3[] = {0.0, 0.35, 0.3, 0.533, 0.766, 1.0};

 	if ( num == 4 ){

		pad[0] = new TPad("pad20", "pad20",setting1[0], setting1[1], setting1[2], setting1[3]);
		pad[1] = new TPad("pad21", "pad21",setting1[2], setting1[1], setting1[3], setting1[3]);
		pad[2] = new TPad("pad28", "pad28",setting1[0], setting1[0], setting1[2], setting1[1]);
		pad[3] = new TPad("pad29", "pad29",setting1[2], setting1[0], setting1[3], setting1[1]);

		for(int i = 0; i < num; i++){

			pad[i]->SetLeftMargin(0.0);
			pad[i]->SetRightMargin(0);
			pad[i]->SetTopMargin(0.0);
			pad[i]->SetBottomMargin(0);
			pad[i]->Draw();
			gPad->SetTicks();

		}

		pad[0]->SetLeftMargin(0.265);
		pad[2]->SetLeftMargin(0.265);

		pad[1]->SetRightMargin(0.05);
		pad[3]->SetRightMargin(0.05);

		pad[0]->SetTopMargin(0.02);
		pad[1]->SetTopMargin(0.02);

		pad[2]->SetBottomMargin(0.3);
		pad[3]->SetBottomMargin(0.3);

	}
	else if( num == 6 ){

	  	pad[0] = new TPad("pad10", "pad10",setting2[0], setting2[1], setting2[2], setting2[4]);
	  	pad[1] = new TPad("pad11", "pad11",setting2[2], setting2[1], setting2[3], setting2[4]);
	  	pad[2] = new TPad("pad12", "pad12",setting2[3], setting2[1], setting2[4], setting2[4]);

	  	pad[3] = new TPad("pad18", "pad18",  setting2[0], setting2[0], setting2[2],  setting2[1]);
	  	pad[4] = new TPad("pad19", "pad19",  setting2[2], setting2[0], setting2[3],  setting2[1]);
	  	pad[5] = new TPad("pad110", "pad110",setting2[3], setting2[0], setting2[4],  setting2[1]);

		for(int i = 0; i < num; i++){

			pad[i]->SetLeftMargin(0.0);
			pad[i]->SetRightMargin(0);
			pad[i]->SetTopMargin(0.0);
			pad[i]->SetBottomMargin(0);
			pad[i]->SetTicks();
			pad[i]->Draw();

		}

		pad[0]->SetLeftMargin(0.265);
		pad[3]->SetLeftMargin(0.265);

		pad[2]->SetRightMargin(0.05);
		pad[5]->SetRightMargin(0.05);

		pad[0]->SetTopMargin(0.02);
		pad[1]->SetTopMargin(0.02);
		pad[2]->SetTopMargin(0.02);

		pad[3]->SetBottomMargin(0.30);
		pad[4]->SetBottomMargin(0.30);
		pad[5]->SetBottomMargin(0.30);

	}
	else if( num == 8 ){

		pad[0] = new TPad("pad10", "pad10",setting3[0], setting3[1], setting3[2], setting3[5]);
		pad[1] = new TPad("pad11", "pad11",setting3[2], setting3[1], setting3[3], setting3[5]);		
  		pad[2] = new TPad("pad12", "pad12",setting3[3], setting3[1], setting3[4], setting3[5]);
		pad[3] = new TPad("pad13", "pad13",setting3[4], setting3[1], setting3[5], setting3[5]);

		pad[4] = new TPad("pad18", "pad18",  setting3[0],  setting3[0], setting3[2], setting3[1]);
		pad[5] = new TPad("pad19", "pad19",  setting3[2],  setting3[0], setting3[3], setting3[1]);
		pad[6] = new TPad("pad110", "pad110",setting3[3],  setting3[0], setting3[4], setting3[1]);
		pad[7] = new TPad("pad111", "pad111",setting3[4],  setting3[0], setting3[5], setting3[1]);

		for( int i = 0; i < num; i++ ){

			pad[i]->SetLeftMargin(0.0);
			pad[i]->SetRightMargin(0);
			pad[i]->SetTopMargin(0.0);
			pad[i]->SetBottomMargin(0);
			pad[i]->SetTicks();
			pad[i]->Draw();		
		}

		pad[0]->SetLeftMargin(0.265);
		pad[4]->SetLeftMargin(0.265);

		pad[3]->SetRightMargin(0.05);
		pad[7]->SetRightMargin(0.05);

		pad[0]->SetTopMargin(0.05);
		pad[1]->SetTopMargin(0.05);
		pad[2]->SetTopMargin(0.05);
		pad[3]->SetTopMargin(0.05);

		pad[4]->SetBottomMargin(0.30);
		pad[5]->SetBottomMargin(0.30);
		pad[6]->SetBottomMargin(0.30);
		pad[7]->SetBottomMargin(0.30);
	}

	for( int i = 0; i < num; i++){

		temp.push_back( pad[i] );
	}

	return temp;
}

TH1D* makeHist(const char*name, const char*title, const char*xtit, const char*ytit, const int nBins, const double lower, const double higher, EColor color = kBlack ){

	TH1D* temp = new TH1D(name, title, nBins, lower, higher);
	
	temp->SetMarkerSize(1.0);
	temp->SetMarkerStyle(20);
	temp->SetMarkerColor(color);
	temp->SetLineColor(color);
	temp->SetStats(kFALSE);

	temp->GetXaxis()->SetTitle( xtit );
	temp->GetXaxis()->SetTitleSize(0.05);
	temp->GetXaxis()->SetTitleFont(42);
	temp->GetXaxis()->SetTitleOffset(1.25);
	temp->GetXaxis()->SetLabelSize(0.05);
	temp->GetXaxis()->SetLabelOffset(0.01);
	temp->GetXaxis()->SetLabelFont(42);
	temp->GetXaxis()->SetLabelColor(kBlack);
	temp->GetXaxis()->CenterTitle();

	temp->GetYaxis()->SetTitle( ytit );
	temp->GetYaxis()->SetTitleSize(0.05);
	temp->GetYaxis()->SetTitleFont(42);
	temp->GetYaxis()->SetTitleOffset(1.4);
	temp->GetYaxis()->SetLabelSize(0.05);
	temp->GetYaxis()->SetLabelOffset(0.01);
	temp->GetYaxis()->SetLabelFont(42);
	temp->GetYaxis()->SetLabelColor(kBlack);
	temp->GetYaxis()->CenterTitle();

	return temp;
}

TH1D* makeHistDifferentBins(const char*name, const char*title, const char*xtit, const char*ytit, const int nBins, double bins[], EColor color = kBlack ){

	TH1D* temp = new TH1D(name, title, nBins, bins);
	
	temp->SetMarkerSize(1.0);
	temp->SetMarkerStyle(20);
	temp->SetMarkerColor(color);
	temp->SetLineColor(color);
	temp->SetStats(kFALSE);

	temp->GetXaxis()->SetTitle( xtit );
	temp->GetXaxis()->SetTitleSize(0.05);
	temp->GetXaxis()->SetTitleFont(42);
	temp->GetXaxis()->SetTitleOffset(1.25);
	temp->GetXaxis()->SetLabelSize(0.05);
	temp->GetXaxis()->SetLabelOffset(0.01);
	temp->GetXaxis()->SetLabelFont(42);
	temp->GetXaxis()->SetLabelColor(kBlack);
	temp->GetXaxis()->CenterTitle();

	temp->GetYaxis()->SetTitle( ytit );
	temp->GetYaxis()->SetTitleSize(0.05);
	temp->GetYaxis()->SetTitleFont(42);
	temp->GetYaxis()->SetTitleOffset(1.4);
	temp->GetYaxis()->SetLabelSize(0.05);
	temp->GetYaxis()->SetLabelOffset(0.01);
	temp->GetYaxis()->SetLabelFont(42);
	temp->GetYaxis()->SetLabelColor(kBlack);
	temp->GetYaxis()->CenterTitle();

	return temp;
}

void fixedFontHist1D(TH1 * h, Float_t xoffset=1.5, Float_t yoffset=2.3)
{
  h->SetLabelFont(43,"X");
  h->SetLabelFont(43,"Y");
  //h->SetLabelOffset(0.01);
  h->SetLabelSize(16);
  h->SetTitleFont(43);
  h->SetTitleSize(20);
  h->SetLabelSize(15,"Y");
  h->SetTitleFont(43,"Y");
  h->SetTitleSize(20,"Y");
  h->SetTitleOffset(xoffset,"X");
  h->SetTitleOffset(yoffset,"Y");
  h->GetXaxis()->CenterTitle();
  h->GetYaxis()->CenterTitle();
  
}

TH2D* make2DHist( const char*name, 
			      const char*title, 
			      const char*xtit, 
			      const char*ytit,
			      const int nxbins,
			      const double xlow, 
			      const double xhigh,
			      const int nybins,
			      const double ylow, 
			      const double yhigh
){

	TH2D* temp2D = new TH2D(name, title, nxbins, xlow, xhigh, nybins, ylow, yhigh);

	temp2D->SetMarkerSize(1.0);
	temp2D->SetMarkerStyle(20);
	temp2D->SetMarkerColor(kBlack);
	temp2D->SetLineColor(kBlack);
	temp2D->SetStats(kFALSE);

	temp2D->GetXaxis()->SetTitle( xtit );
	temp2D->GetXaxis()->SetTitleSize(0.04);
	temp2D->GetXaxis()->SetTitleFont(42);
	temp2D->GetXaxis()->SetTitleOffset(1.4);
	temp2D->GetXaxis()->SetLabelSize(0.04);
	temp2D->GetXaxis()->SetLabelOffset(0.01);
	temp2D->GetXaxis()->SetLabelFont(42);
	temp2D->GetXaxis()->SetLabelColor(kBlack);

	temp2D->GetYaxis()->SetTitle( ytit );
	temp2D->GetYaxis()->SetTitleSize(0.04);
	temp2D->GetYaxis()->SetTitleFont(42);
	temp2D->GetYaxis()->SetTitleOffset(1.7);
	temp2D->GetYaxis()->SetLabelSize(0.04);
	temp2D->GetYaxis()->SetLabelOffset(0.01);
	temp2D->GetYaxis()->SetLabelFont(42);
	temp2D->GetYaxis()->SetLabelColor(kBlack);

	return temp2D;

}

void fixedFontHist(TH2D * h, Float_t xoffset=0.9, Float_t yoffset=2.7)
{
  h->SetLabelFont(43,"X");
  h->SetLabelFont(43,"Y");
  //h->SetLabelOffset(0.01);
  h->SetLabelSize(16);
  h->SetTitleFont(43);
  h->SetTitleSize(20);
  h->SetLabelSize(15,"Y");
  h->SetTitleFont(43,"Y");
  h->SetTitleSize(17,"Y");
  h->SetTitleOffset(xoffset,"X");
  h->SetTitleOffset(yoffset,"Y");
  h->GetXaxis()->CenterTitle();
  h->GetYaxis()->CenterTitle();
}
void make_dNdX( TH1D* hist ){

	for(int i=0;i<hist->GetNbinsX();i++){
		double value = hist->GetBinContent(i+1);
		double error = hist->GetBinError(i+1);
		double binwidth = hist->GetBinWidth(i+1);

		hist->SetBinContent(i+1, value / binwidth );
		hist->SetBinError(i+1, error / binwidth );
	}
}

double calColError(double Ea, double Eb, double Sa, double Sb){

	double temp = Ea/Eb;
	double temp2 = (Sa*Sa)/(Ea*Ea) + (Sb*Sb)/(Eb*Eb);
	double temp3 = (2.*Sa*Sa)/(Ea*Eb);

	return temp*(sqrt(TMath::Abs(temp2-temp3)) );
}

TH1D* make_systematicRatio(TH1D* hist1, TH1D* hist2){

	TH1D* hist_ratio = (TH1D*) hist1->Clone("hist_ratio");

	if( hist1->GetNbinsX() != hist2->GetNbinsX() ){
		std::cout << "Not compatible binning, abort!" << std::endl;
		return 0;
	}

	for(int ibin=0;ibin<hist1->GetNbinsX();ibin++){
		double value_a = hist1->GetBinContent(ibin+1);
		double error_a = hist1->GetBinError(ibin+1);

		double value_b = hist2->GetBinContent(ibin+1);
		double error_b = hist2->GetBinError(ibin+1);

		hist_ratio->SetBinContent(ibin+1, value_a / value_b );
		hist_ratio->SetBinError(ibin+1, calColError(value_a, value_b, error_a, error_b) );
	}

	return hist_ratio;
}

TLegend* makeLegend(){

	TLegend *w2 = new TLegend(0.65,0.15,0.90,0.45);
	w2->SetLineColor(kWhite);
	w2->SetFillColor(0);
	return w2;

}

TGraphAsymmErrors* makeEfficiency(TH1D* hist1, TH1D* hist2, const char*Draw = "cp", EColor color = kBlack  ){

	TGraphAsymmErrors* temp = new TGraphAsymmErrors();
	temp->Divide( hist1, hist2, Draw );
	temp->SetMarkerStyle(20);
	temp->SetMarkerColor(color);
	temp->SetLineColor(color);

	return temp;

}

TLatex* makeLatex(const char* txt,  double x, double y){

	TLatex* r = new TLatex(x, y, txt);
	r->SetTextSize(0.05);
	r->SetNDC();
	return r;

}

void drawBox(TH1D* hist1, double sys, bool doPercentage, double xe= 0.05){

	TBox* temp_box[100];
	int bins = hist1->GetNbinsX();

	for(int deta = 0; deta < bins; deta++){

		double value = hist1->GetBinContent(deta+1);
		double bincenter = hist1->GetBinCenter(deta+1);
		double sys_temp = 0.;

		if( doPercentage ) sys_temp = value*sys;
		else sys_temp = sys;

		temp_box[deta] = new TBox(bincenter-xe,value-sys_temp,bincenter+xe,value+sys_temp);
	    temp_box[deta]->SetFillColorAlpha(kGray+2,0.4);
	    temp_box[deta]->SetFillStyle(1001);
		temp_box[deta]->SetLineWidth(0);
    	temp_box[deta]->Draw("same");
    }

}

void drawBoxRatio(TH1D* hist1, TH1D* hist2, double sys, bool doPercentage){

	TBox* temp_box[100];
	double xe = 0.08;

	int bins = hist1->GetNbinsX();

	for(int deta = 0; deta < bins; deta++){

		if(deta > 6) continue;

		double factor = hist2->GetBinContent(deta+1); 
		double value = hist1->GetBinContent(deta+1);
		double bincenter = hist1->GetBinCenter(deta+1);

		if( doPercentage ) sys = sqrt((value*0.045)*(value*0.045));
		else sys = sys;

		double sys_temp;
		sys_temp = sys;
	
		temp_box[deta] = new TBox(bincenter-xe,value-sys_temp,bincenter+xe,value+sys_temp);
	    temp_box[deta]->SetFillColorAlpha(kGray+2,0.4);
	    temp_box[deta]->SetFillStyle(1001);
		temp_box[deta]->SetLineWidth(0);
    	temp_box[deta]->Draw("same");

    }

}

void drawBoxTGraphRatio(TGraphErrors* gr1, int bins, double sys, bool doPercentage){

	double xe[11];
	TBox* box1[11];

    for(int mult = 0; mult < bins; mult++){

    	if( mult < 6 ){

    		xe[mult] = 15*log(1.1*(mult+1));
    		if(mult == 0) xe[mult] = 10;
    	}
  		if( mult ==  6) xe[mult] = 37;
  		if( mult ==  7) xe[mult] = 50;
  		if( mult ==  8) xe[mult] = 50;
  		if( mult ==  9) xe[mult] = 63;
  		if( mult ==  10) xe[mult] = 73;
    	
    	
    	double x1;
    	double value1;
    	gr1->GetPoint(mult, x1, value1);

    	double ye = sys;
  		if( doPercentage ) ye = sqrt((value1*0.045)*(value1*0.045));
  		else ye = sys;

    	box1[mult] = new TBox(x1-xe[mult],value1-ye,x1+xe[mult],value1+ye);
	 	box1[mult]->SetFillColorAlpha(kGray+2,0.4);
        box1[mult]->SetFillStyle(1001);
    	box1[mult]->SetLineWidth(0);
    	box1[mult]->SetLineColor(kRed);
        box1[mult]->Draw("SAME");

    }

}

void drawBoxTGraph(TGraphErrors* gr1, int bins, double sys, bool doPercentage, bool doConstantWidth){

	double xe[100];
	TBox* box1[100];

    for(int mult = 0; mult < bins; mult++){

    	if(!doConstantWidth){
	    	if( mult < 6 ){

	    		xe[mult] = 15*log(1.1*(mult+1));
	    		if(mult == 0) xe[mult] = 10;
	    	}
	  		if( mult ==  6) xe[mult] = 37;
	  		if( mult ==  7) xe[mult] = 50;
	  		if( mult ==  8) xe[mult] = 50;
	  		if( mult ==  9) xe[mult] = 63;
	  		if( mult ==  10) xe[mult] = 73;
  		}
  		else{ xe[mult] = 0.02;}
    	
    	
    	double x1;
    	double value1;
    	gr1->GetPoint(mult, x1, value1);

    	double ye = sys;
  		if( doPercentage ) ye = value1 * sys;
  		else ye = sys;

    	box1[mult] = new TBox(x1-xe[mult],value1-ye,x1+xe[mult],value1+ye);
	 	box1[mult]->SetFillColorAlpha(kGray+2,0.4);
        box1[mult]->SetFillStyle(1001);
    	box1[mult]->SetLineWidth(0);
    	box1[mult]->SetLineColor(kBlack);
        box1[mult]->Draw("LSAME");

    }
}

void drawBoxTGraph_alt(TGraphErrors* gr1, int bins, double sys, bool doPercentage, bool doConstantWidth){

	double xe[100];
	TBox* box1[100];

    for(int mult = 0; mult < bins; mult++){

    	if(!doConstantWidth){
	    	if( mult < 6 ){

	    		xe[mult] = 15*log(1.1*(mult+1));
	    		if(mult == 0) xe[mult] = 10;
	    	}
	  		if( mult ==  6) xe[mult] = 37;
	  		if( mult ==  7) xe[mult] = 50;
	  		if( mult ==  8) xe[mult] = 50;
	  		if( mult ==  9) xe[mult] = 63;
	  		if( mult ==  10) xe[mult] = 73;
  		}
  		else{ xe[mult] = 0.005;}
    	
    	
    	double x1;
    	double value1;
    	gr1->GetPoint(mult, x1, value1);

    	double ye = sys;
  		if( doPercentage ) ye = value1 * sys;
  		else ye = sys;

    	box1[mult] = new TBox(x1-xe[mult],value1-ye,x1+xe[mult],value1+ye);
	 	box1[mult]->SetFillColorAlpha(kGray+2,0.4);
        box1[mult]->SetFillStyle(1001);
    	box1[mult]->SetLineWidth(0);
    	box1[mult]->SetLineColor(kRed);
        box1[mult]->Draw("SAME");

    }
}


void drawBoxTGraphDiff(TGraphErrors* gr1, TGraphErrors* gr2, int bins, double sys, bool doPercentage){

	double xe[11];
	TBox* box1[11];
	TBox* box2[11];

    for(int mult = 0; mult < bins; mult++){

    	xe[mult] = 10*log(1.1*(mult+1));
    	if(mult == 0) xe[mult] = 7;

    	double x1;
    	double value1;
    	gr1->GetPoint(mult, x1, value1);

    	double x2;
    	double value2;
    	gr2->GetPoint(mult, x2, value2);

    	double value = value2 - value1;

    	double ye = sys;
  		if( doPercentage ) ye = sqrt((value*0.045)*(value*0.045)+sys*sys);
  		else ye = sys;

    	box1[mult] = new TBox(x1-xe[mult],value-ye,x1+xe[mult],value+ye);
	 	box1[mult]->SetFillColorAlpha(kGray+2,0.4);
        box1[mult]->SetFillStyle(1001);
    	box1[mult]->SetLineWidth(0);
    	box1[mult]->SetLineColor(kRed);
        box1[mult]->Draw("SAME");

    }

}