#include <iostream>
#include "read_gain.C"
#include "read_voltages.C"
#include "read_hv_mapping.C"


void adjust()	{


	int ly,a,b;
	cout<<"Enter layer #"<<endl;
	cin>>ly;

	if(ly==1)	{a=10; const b=30;}
	if(ly==2)	{a=6;  const b=24;} 
	if(ly==3)	{a=6;  const b=22;}
	if(ly==4)	{a=4;  const b=12;}
	if(ly==5)	{a=4;  const b=12;}
	if(ly==6)	{a=4;  const b=12;}


	int d;

	double meas_ADC_L[30]={0};
	double meas_ADC_R[30]={0};
	double V_l[30]={0};
	double V_r[30]={0};


	double trg=1500;

	int ind;


	double ADC_L[30]={0};
	double ADC_R[30]={0};
	double DL[30]={0};
	double DL1[30]={0};
	double DR[30]={0};
	double d_vL[30]={0};
	double d_vR[30]={0};
	double N_VL[30]={0};
	double N_VR[30]={0};
	double dADC_R[30]={0};
	double dADC_L[30]={0};




	//=====================================================

	double aL[112]={0};
	double aR[112]={0};
	double bL[112]={0},bR[112]={0};
	double pL[112],pR[112];

	//=====================================================
	int bars_p[6]={0,30,54,76,88,100};



	//=====================================================

	int CrateL[112]={-1}, SlotL[112]={-1}, ChanL[112]={-1}; 
	int CrateR[112]={-1}, SlotR[112]={-1}, ChanR[112]={-1};
	int Layer[112]={-1},Bar[112]={-1};
	double HV1[16][12]={0}, HV2[16][12]={0};	// variables for read_voltages.C

	read_hv_mapping(CrateL,SlotL,ChanL,CrateR,SlotR,ChanR,Layer,Bar);
	read_gain(aL,bL,pL,aR,bR,pR);
	read_voltages(HV1,HV2);

//================================= Read The measured ADCs and Voltages from file. This file created by ADC.C ============
ifstream fL,fR;
char val[200];
int indx=0,j=0;
fL.open("ADC_left.dat");
for(int i=0;i<=2*b;i++)
	{
	if(indx==0)
		{
		fL>>val;
		meas_ADC_L[j]=atof(val);
		j++;
		}
	if(indx==1)
		{
		fL>>val;
		V_l[j]=atof(val);
		j++;
		}
	if(!strcmp(val,"HV"))
		{
		if(i!=b)	cout<<"NON SYMMETRIC FILE, DATA IS MISSING in ADC_left"<<endl;
		indx=1;
		j=0;
		}
	
	}
fL.close();

fR.open("ADC_right.dat");
j=0;
indx=0;
for(int i=0;i<=2*b;i++)
	{
	if(indx==0)
		{
		fR>>val;
		meas_ADC_R[j]=atof(val);
		j++;
		}
	if(indx==1)
		{
		fR>>val;
		V_r[j]=atof(val);
		j++;
		}
	if(!strcmp(val,"HV"))
		{
		if(i!=b)	cout<<"NON SYMMETRIC FILE, DATA IS MISSING in ADC_right"<<endl;
		indx=1;
		j=0;
		}
	
	}
fR.close();

//====================================================================================================================




for(i=0;i<b;i++)
	{

	ind=bars_p[ly-1]+i;
	ADC_L[i]=(aL[ind])*(V_l[i])**(pL[ind])+bL[ind];
	dADC_L[i]=-meas_ADC_L[i]+ADC_L[i];
	if((ADC_L[i]-bL[ind]+dADC_L[i])/aL[ind]>0)
		{
		N_VL[i]=((ADC_L[i]-bL[ind]+dADC_L[i])/aL[ind])**(1./pL[ind]);
		}
	else N_VL[i]=V_l[i];

cout<<" =  = == = =  = = = = = = = = = = = = = ==  = ="<<endl;
cout<<" aL[i]        = "<<aL[ind]<<endl;
cout<<" pL[i]        = "<<pL[ind]<<endl;
cout<<" bL[i]        = "<<bL[ind]<<endl;
cout<<" V 	     = "<<V_l[i]<<endl;
cout<<" ADC_L[i]     = "<<ADC_L[i]<<endl;
cout<<" ADC_meas[i]  = "<<meas_ADC_L[i]<<endl;
cout<<"++++++++++++++++++++++++++++++++++++++++++++++"<<endl;

	ADC_R[i]=(aR[ind])*(V_r[i])**(pR[ind])+bR[ind];

	V_r[i]=((ADC_R[i]-bR[ind])/aR[ind])**(1./pR[ind]);
	dADC_R[i]=-meas_ADC_R[i]+ADC_R[i];
	if((ADC_R[i]-bR[ind]+dADC_R[i])/aR[ind]>0)
		{
		N_VR[i]=((ADC_R[i]-bR[ind]+dADC_R[i])/aR[ind])**(1./pR[ind]);
		}
	else N_VL[i]=V_r[i];
	}


int crL[30]={0},slL[30]={0},chL[30]={0};
int crR[30]={0},slR[30]={0},chR[30]={0};
//ofstream file;
//file.open("volt.dat", ios::ate);
//file
for(i=0;i<112;i++)
	{
	if(Layer[i]==ly)
		{
		for(j=0;j<b;j++)
			{
			crL[j]=CrateL[i];
			slL[j]=SlotL[i];
			chL[j]=ChanL[i];
			
			crR[j]=CrateR[i];
			slR[j]=SlotR[i];
			chR[j]=ChanR[i];
			}
		}
	}

/*for(i=0;i<b;i++)
	{
	if(crL[i]==1)
		{
		if(
		sprintf(line,"bbps1:2001 S%d",slL[i]);

*/	



for(int j=0;j<12;j++)
	{
	cout<<"Left Bar = "<<j<<"  Old Voltage = "<<V_l[j]<<"  New Voltage = "<<N_VL[j]<<endl;
 	cout<<"Right Bar = "<<j<<" Old Voltage = "<<V_r[j]<<"  New Voltage = "<<N_VR[j]<<endl;
	cout<<endl;
	}
}