/*

This script is reading the HV mapping of BB scintillators planes dE and E.

*/



#include <cstring>
#include <string>


void read_hv_mapping(int *NaL_layer, int *NaL_bar,  int *NaL_crate, int *NaL_slot, int *NaL_ch, int *NaR_layer, int *NaR_bar, int *NaR_crate, int *NaR_slot, int *NaR_ch)
{

char line1[200],line2[200],s[200],s0[200],s1[200],s2[200],s3[200],*file,*dd,*side,*strL,*strR,*detector;
file=(char*)malloc(1000);
dd=(char*)malloc(1000);
side=(char*)malloc(1000);
strL=(char*)malloc(1000);
strR=(char*)malloc(1000);
detector=(char*)malloc(1000);

int ind1,ind2,ind3,ind4,ind5,indx;
int k1=0,k2=0;

sprintf(file,"./HV-mapping.dat");
FILE *dbFile3 = fopen(file,"r");
int t=0;
while(fgets(line1, sizeof(line1), dbFile3) !=NULL) 
	{
	sscanf(line1,"%s %s",&s1, &s2);
	sprintf(side,"%s",s2);
	sprintf(strL,"left");
	sprintf(strR,"right");

	if(!strcmp(side,strL)) t=1;
	if(!strcmp(side,strR)) t=2;

	if(s1[0]=='#' && t==1)	indx=1;	// left E plane scintillators
	if(s1[0]=='#' && t==2)	indx=2;	// right E plane scintillators

	if(indx==1 && s1[0]!='#')
		{
		sscanf(line1," %d %d %d %d %d", &ind1, &ind2, &ind3, &ind4, &ind5);
//		cout<<" ind 1 = "<<ind1<<" ind 2 = "<<ind2<<" ind 3 = "<<ind3<<" ind 4 = "<<ind4<<endl;

		*(NaL_layer+k1)=ind1;
		*(NaL_bar+k1)=ind2;
		*(NaL_crate+k1)=ind3;
		*(NaL_slot+k1)=ind4;
		*(NaL_ch+k1)=ind5;
//cout<<ind1<<"    "<<ind5<<endl;
		k1++;
		}
	if(indx==2 && s1[0]!='#')
		{
		sscanf(line1,"%d %d %d %d %d",&ind1,&ind2,&ind3, &ind4,&ind5);

		*(NaR_layer+k2)=ind1;
		*(NaR_bar+k2)=ind2;
		*(NaR_crate+k2)=ind3;
		*(NaR_slot+k2)=ind4;
		*(NaR_ch+k2)=ind5;
		k2++;
		}

	}
}