// Hit and Cluster Profile	saved as :	profile_plots.pdf
// Offset for RUN0133		saved as :	Offset_Run0133.pdf


//	Formula for resolution
//	DX= track.m * zOff + trackx.q - sCMSNS2LC2.geoposX
//
//	track.m & trackx.q  will we available in RECO*.root file
//	zOff
//	Call tree like : T_AN59		T_RECO59	T_CRC59
	
#include "Loaod_all_root_files.h"
#include "Riostream.h"
#include "TLegend.h"

int alignment()
{
////////////////////////////////////////////////////////////////////////////////
//		MUON TRACK PARAMETERS
////////////////////////////////////////////////////////////////////////////////
/*
const Int_t nbranch = 9;
char *bnames[nbranch] = {"trackx.m","trackx.q","trackx.chi2",
			"tracky.m","tracky.q","tracky.chi2",
			"sCMSNS2LC1.geoposX", "sCMSNS2LC2.geoposX","sCMSNS2LC3.geoposX" 
			};

char *fnames[nbranch] = {"trackx_m","trackx_q","trackx_chi2",
			"tracky_m","tracky_q","tracky_chi2",
			"Hit Profile for LC1", "Hit Profile for LC2", "Hit Profile for LC3"
			};

const Float_t bins[nbranch] = {20,	20,	20,
				20,	20,	20,
				50,	50,	50
			      };

const Float_t range[2*nbranch] ={-0.003,0.003,	-10,100,	-1,10,
				-0.003,0.003,	-10,100,	-1,10
				0,120,		0,120,		0,120
				};

TH1F * hist[nbranch];
//TProfile * hprofile[nbranch];
TCanvas *Canvas[nbranch];
//TF1 *function[nbranch];

Char_t message[90];
vector<int> vecm,vecn; 
Double_t fit_Low, fit_High;
Int_t k=0;

//ofstream o_file;
//o_file.open("fit_detail.txt");

//o_file<<"NAME_OF_PROFILE_PLOT\t\t\tChi-sqr/NDF\tINTERCEPT(P0)\tINTERCEPT_ERROR\t\tSLOPE(P1)\tSLOPE_ERROR"<<endl;
for(Int_t i=0;i<nbranch;i++){
	k=i*2;
	Canvas[i] = new TCanvas(Form("Canvas%d",i));
	Canvas[i]->Divide(3,2);
	Canvas[i]->cd(1);

	hist[i] = new TH1F("H1_"+TString(fnames[i]),TString(bnames[i])+"_R139",bins[i],range[k],range[k+1]);
	T_RECO139->Draw(TString(bnames[i])+">>H1_"+TString(fnames[i]));
	hist[i]->Draw();
	//delete gDirectory->FindObject("H_"+TString(fnames[i]));

	Canvas[i]->cd(2);

	hist[i] = new TH1F("H2_"+TString(fnames[i]),TString(bnames[i])+"_R110",bins[i],range[k],range[k+1]);
	T_RECO110->Draw(TString(bnames[i])+">>H2_"+TString(fnames[i]));
	hist[i]->Draw();

	Canvas[i]->cd(3);

	hist[i] = new TH1F("H3_"+TString(fnames[i]),TString(bnames[i])+"_R115",bins[i],range[k],range[k+1]);
	T_RECO115->Draw(TString(bnames[i])+">>H3_"+TString(fnames[i]));
	hist[i]->Draw();

	Canvas[i]->cd(4);

	hist[i] = new TH1F("H4_"+TString(fnames[i]),TString(bnames[i])+"_R121",bins[i],range[k],range[k+1]);
	T_RECO121->Draw(TString(bnames[i])+">>H4_"+TString(fnames[i]));
	hist[i]->Draw();

	Canvas[i]->cd(5);

	hist[i] = new TH1F("H5_"+TString(fnames[i]),TString(bnames[i])+"_R127",bins[i],range[k],range[k+1]);
	T_RECO127->Draw(TString(bnames[i])+">>H5_"+TString(fnames[i]));
	hist[i]->Draw();

	Canvas[i]->cd(6);

	hist[i] = new TH1F("H6_"+TString(fnames[i]),TString(bnames[i])+"_R133",bins[i],range[k],range[k+1]);
	T_RECO133->Draw(TString(bnames[i])+">>H6_"+TString(fnames[i]));
	hist[i]->Draw();

	Canvas[i]->SaveAs("parameter1_"+TString(fnames[i])+".pdf");
	Canvas[i]->Clear();

	Canvas[i]->Divide(3,2);
	Canvas[i]->cd(1);

	hist[i] = new TH1F("H7_"+TString(fnames[i]),TString(bnames[i])+"_R144",bins[i],range[k],range[k+1]);
	T_RECO144->Draw(TString(bnames[i])+">>H7_"+TString(fnames[i]));
	hist[i]->Draw();

	Canvas[i]->cd(2);

	hist[i] = new TH1F("H8_"+TString(fnames[i]),TString(bnames[i])+"_R150",bins[i],range[k],range[k+1]);
	T_RECO150->Draw(TString(bnames[i])+">>H8_"+TString(fnames[i]));
	hist[i]->Draw();

	Canvas[i]->cd(3);

	hist[i] = new TH1F("H9_"+TString(fnames[i]),TString(bnames[i])+"_R156",bins[i],range[k],range[k+1]);
	T_RECO156->Draw(TString(bnames[i])+">>H9_"+TString(fnames[i]));
	hist[i]->Draw();

	Canvas[i]->cd(4);

	hist[i] = new TH1F("H10_"+TString(fnames[i]),TString(bnames[i])+"_R162",bins[i],range[k],range[k+1]);
	T_RECO162->Draw(TString(bnames[i])+">>H10_"+TString(fnames[i]));
	hist[i]->Draw();

	Canvas[i]->cd(5);

	hist[i] = new TH1F("H11_"+TString(fnames[i]),TString(bnames[i])+"_R168",bins[i],range[k],range[k+1]);
	T_RECO168->Draw(TString(bnames[i])+">>H11_"+TString(fnames[i]));
	hist[i]->Draw();

	Canvas[i]->cd(6);

	hist[i] = new TH1F("H12_"+TString(fnames[i]),TString(bnames[i])+"_R133",bins[i],range[k],range[k+1]);
	T_RECO133->Draw(TString(bnames[i])+">>H12_"+TString(fnames[i]));
	hist[i]->Draw();

	Canvas[i]->SaveAs("parameter2_"+TString(fnames[i])+".pdf");
	Canvas[i]->Clear();


}
*/

const Int_t nbranch1 = 3;
char *bnames1[nbranch1] = {
			"sCMSNS2LC1.geoposX", "sCMSNS2LC2.geoposX","sCMSNS2LC3.geoposX" 
			};

char *fnames1[nbranch1] = {
			"Hit_Profile_for_LC1", "Hit_Profile_for_LC2", "Hit_Profile_for_LC3"
			};

const Float_t bins1[nbranch1] = {
				30,	30,	30
			      };

const Float_t range1[2*nbranch1] ={
				-20,70,		-20,70,		-20,70
				};

TH1F * hist1[nbranch1];
TCanvas *Canvas1[nbranch1];

k=0;

for(Int_t i=0;i<nbranch1;i++){
	k=i*2;
	Canvas1[i] = new TCanvas(Form("Canvas1%d",i));
	Canvas1[i]->Divide(3,2);
	Canvas1[i]->cd(1);

	hist1[i] = new TH1F("H1_"+TString(fnames1[i]),TString(bnames1[i])+"_R139",bins1[i],range1[k],range1[k+1]);
	T_CRC139->Draw(TString(bnames1[i])+">>H1_"+TString(fnames1[i]));
	hist1[i]->Draw();
	//delete gDirectory->FindObject("H_"+TString(fnames1[i]));

	Canvas1[i]->cd(2);

	hist1[i] = new TH1F("H2_"+TString(fnames1[i]),TString(bnames1[i])+"_R110",bins1[i],range1[k],range1[k+1]);
	T_CRC110->Draw(TString(bnames1[i])+">>H2_"+TString(fnames1[i]));
	hist1[i]->Draw();

	Canvas1[i]->cd(3);

	hist1[i] = new TH1F("H3_"+TString(fnames1[i]),TString(bnames1[i])+"_R115",bins1[i],range1[k],range1[k+1]);
	T_CRC115->Draw(TString(bnames1[i])+">>H3_"+TString(fnames1[i]));
	hist1[i]->Draw();

	Canvas1[i]->cd(4);

	hist1[i] = new TH1F("H4_"+TString(fnames1[i]),TString(bnames1[i])+"_R121",bins1[i],range1[k],range1[k+1]);
	T_CRC121->Draw(TString(bnames1[i])+">>H4_"+TString(fnames1[i]));
	hist1[i]->Draw();

	Canvas1[i]->cd(5);

	hist1[i] = new TH1F("H5_"+TString(fnames1[i]),TString(bnames1[i])+"_R127",bins1[i],range1[k],range1[k+1]);
	T_CRC127->Draw(TString(bnames1[i])+">>H5_"+TString(fnames1[i]));
	hist1[i]->Draw();

	Canvas1[i]->cd(6);

	hist1[i] = new TH1F("H6_"+TString(fnames1[i]),TString(bnames1[i])+"_R133",bins1[i],range1[k],range1[k+1]);
	T_CRC133->Draw(TString(bnames1[i])+">>H6_"+TString(fnames1[i]));
	hist1[i]->Draw();

	Canvas1[i]->SaveAs("parameter1_"+TString(fnames1[i])+".pdf");
	Canvas1[i]->Clear();

	Canvas1[i]->Divide(3,2);
	Canvas1[i]->cd(1);

	hist1[i] = new TH1F("H7_"+TString(fnames1[i]),TString(bnames1[i])+"_R144",bins1[i],range1[k],range1[k+1]);
	T_CRC144->Draw(TString(bnames1[i])+">>H7_"+TString(fnames1[i]));
	hist1[i]->Draw();

	Canvas1[i]->cd(2);

	hist1[i] = new TH1F("H8_"+TString(fnames1[i]),TString(bnames1[i])+"_R150",bins1[i],range1[k],range1[k+1]);
	T_CRC150->Draw(TString(bnames1[i])+">>H8_"+TString(fnames1[i]));
	hist1[i]->Draw();

	Canvas1[i]->cd(3);

	hist1[i] = new TH1F("H9_"+TString(fnames1[i]),TString(bnames1[i])+"_R156",bins1[i],range1[k],range1[k+1]);
	T_CRC156->Draw(TString(bnames1[i])+">>H9_"+TString(fnames1[i]));
	hist1[i]->Draw();

	Canvas1[i]->cd(4);

	hist1[i] = new TH1F("H10_"+TString(fnames1[i]),TString(bnames1[i])+"_R162",bins1[i],range1[k],range1[k+1]);
	T_CRC162->Draw(TString(bnames1[i])+">>H10_"+TString(fnames1[i]));
	hist1[i]->Draw();

	Canvas1[i]->cd(5);

	hist1[i] = new TH1F("H11_"+TString(fnames1[i]),TString(bnames1[i])+"_R168",bins1[i],range1[k],range1[k+1]);
	T_CRC168->Draw(TString(bnames1[i])+">>H11_"+TString(fnames1[i]));
	hist1[i]->Draw();

	Canvas1[i]->cd(6);

	hist1[i] = new TH1F("H12_"+TString(fnames1[i]),TString(bnames1[i])+"_R133",bins1[i],range1[k],range1[k+1]);
	T_CRC133->Draw(TString(bnames1[i])+">>H12_"+TString(fnames1[i]));
	hist1[i]->Draw();

	Canvas1[i]->SaveAs("parameter2_"+TString(fnames1[i])+".pdf");
}
}