{

// This script find location of pedestals for E and dE planes for BB


//gROOT->Reset();
//gStyle->SetOptFit(1);


Int_t nbins=1990;
Int_t bmin=10, bmax=2000;
int nbars;
int a, pnum, tr;
Double_t maxR, maxL;
Double_t Rpedest[35],Lpedest[35];


char *left, *right, *left1, * right1, *c;
left=(char *)malloc(1000);
right=(char *)malloc(1000);
left1=(char *)malloc(1000);
right1=(char *)malloc(1000);
c=(char *)malloc(1000);

TH1D *hr=new TH1D("hr","hr",nbins,bmin,bmax);
TH1D *hl=new TH1D("hl","hl",nbins,bmin,bmax);

cout<<"\n"<<endl;
cout<<"Choose Detector: \n E plane (1) \n dE plane (2) \n" <<endl;

cin>>a;

cout<<"\n"<<endl;
cout<<"Enter trigger type (0 for now cut)\n"<<endl;
cin>>tr;
//===================== ADC pedestals for E plane ==============================

if(tr==0) TCut cut="";
else
	{
	sprintf(c,"fEvtHdr.fEvtType==%i",tr);
	TCut cut=c;
	}
if(a==1)
	{

	nbars=24;

	TCanvas c1;
	TCanvas c2;
	c2.Divide(4,6);
	TCanvas c3;
	c3.Divide(4,6);
	
	for(int i=0;i<nbars; i++)
		{
		c1.cd();
		sprintf(left,"BB.tp.e.LA[%i]>>hl",i);
		sprintf(right,"BB.tp.e.RA[%i]>>hr",i);
 		T->Draw(left,cut);
		T->Draw(right,cut);
		maxL=hl.GetMaximumBin();  	                        // find the maximum channel for Left ADC
		maxR=hr.GetMaximumBin();				// find the maximum channel for Right ADC
		
		TH1D *HL=new TH1D("HL","HL",200,maxL-30,maxL+30);	// define histogram around the maximum location
		TH1D *HR=new TH1D("HR","HR",200,maxR-30,maxR+30);
	
	
		sprintf(left1,"BB.tp.e.LA[%i]>>HL",i);
		sprintf(right1,"BB.tp.e.RA[%i]>>HR",i);
		
		c2.cd(i+1);
		T->Draw(left1,cut);
		Lpedest[i]=HL.GetBinCenter(HL.GetMaximumBin());
		
		
		c3.cd(i+1);
		T->Draw(right1,cut);
		Rpedest[i]=HR.GetBinCenter(HR.GetMaximumBin());

		} //for
	} // if

// ====================== ADC pedestals for dE plane ===================================

if(a==2)
	{
	nbars=24;

	TCanvas c1;
	TCanvas c2;
	c2.Divide(4,6);
	TCanvas c3;
	c3.Divide(4,6);
	
	for(int i=0;i<nbars; i++)
		{
		c1.cd();
		sprintf(left,"BB.tp.de.LA[%i]>>hl",i);
		sprintf(right,"BB.tp.de.RA[%i]>>hr",i);
		T->Draw(left,cut);
		T->Draw(right,cut);
		maxL=hl.GetMaximumBin();  	                        // find the maximum channel for Left ADC
		maxR=hr.GetMaximumBin();				// find the maximum channel for Right ADC
		
		TH1D *HL=new TH1D("HL","HL",200,maxL-30,maxL+30);	// define histogram around the maximum location
		TH1D *HR=new TH1D("HR","HR",200,maxR-30,maxR+30);
	
	
		sprintf(left1,"BB.tp.de.LA[%i]>>HL",i);
		sprintf(right1,"BB.tp.de.RA[%i]>>HR",i);

		c2.cd(i+1);
		T->Draw(left1,cut);
		Lpedest[i]=HL.GetBinCenter(HL.GetMaximumBin());

		c3.cd(i+1);		
		T->Draw(right1,cut);
		Rpedest[i]=GetBinCenter(HR.GetMaximumBin());
		} //for
	}//if



// ========================= Print pedestals on screen =======================
	
for(int j=0;j<nbars;j++)
	{
	cout<<"Bar # --- "<<j+1<<"  LA --- "<<Lpedest[j]<<" | "<<"  RA --- "<<Rpedest[j]<<endl;


	}



//========================= Print pedestals to file ===================================
char *lineL, *lineR, *line, *header, *name, *line2, *header2, *name2;
lineL=(char *)malloc(1000);
lineR=(char *)malloc(1000);
line=(char *)malloc(1000);
line2=(char *)malloc(1000);
header=(char *)malloc(1000);
header2=(char *)malloc(1000);
name=(char *)malloc(1000);
name2=(char *)malloc(1000);
	
if(a==1 || a==2)
	{
	ofstream L,R;

	//====================== E plane ========================
	if(a==1)
		{
		L.open("pedestalsE_left.dat");
		R.open("pedestalsE_right.dat");
		L<<"[ BB.tp.e.left_pedestals ]"<<endl;
		R<<"[ BB.tp.e.right_pedestals ]"<<endl;
		L<<"# first and last index in array"<<endl;
		R<<"# first and last index in array"<<endl;
		L<<"0 23"<<endl;
		R<<"0 23"<<endl;
		L<<"# Table start here"<<endl;
		R<<"# Table start here"<<endl;
		}
	// ===================== dE plane =======================
	if(a==2)
		{
		L.open("pedestalsdE_left.dat");
		R.open("pedestalsdE_right.dat");
		L<<"[ BB.tp.de.left_pedestals ]"<<endl;
		R<<"[ BB.tp.de.right_pedestals ]"<<endl;
		L<<"# first and last index in array"<<endl;
		R<<"# first and last index in array"<<endl;
		L<<"0 23"<<endl;
		R<<"0 23"<<endl;
		L<<"# Table start here"<<endl;
		R<<"# Table start here"<<endl;
		}


		for(int k=0;k<nbars;k+=4)
			{
			sprintf(lineL,"%.3f %.3f %.3f %.3f",Lpedest[k],Lpedest[k+1],Lpedest[k+2],Lpedest[k+3]);
			L<<lineL<<endl;
			sprintf(lineR,"%.3f %.3f %.3f %.3f",Rpedest[k],Rpedest[k+1],Rpedest[k+2],Rpedest[k+3]);
			R<<lineR<<endl;

			}
		L.close();	
		R.close();
	}


}