//------------------------------------------
//reading Ramo_table.txt(FEM solutions) and make PotentialArray and PotentialHVArray
//2019.3.22,5.22 Taka Kondo (KEK)
//------------------------------------------------

void make_EyOnly(const int VFD, const int HV)
{
//  Neff[] like P10E11, M08E11 for Neff=+10e11/cm**3 and -8e11/cm**3
//  HV = bias Voltage applied to SCT backplane

ifstream fin("FEM_dummy.txt");
ofstream fout("Ey_dummy.txt"); 
//fout<<setprecision(2)<<setiosflags(ios::fixed)<<endl;
float X[17][115], Y[17][115], Ey[17][115], Ex[17][115]; 
double x,y,phi,ex,ey,e,dx,dy,d,eng;
int nn=0,nreg, k=0, i=0;
char line[255]; 

fin.getline(line,sizeof(line)); 
while (fin>>x>>y>>nreg>>phi>>ex>>ey>>e>>dx>>dy>>d>>eng) {
nn++;
cout<<k<<" "<<i<<" x,y="<<x<<", "<<y<<" um,  "<<nreg<<" "<<phi<<" "<<ey/100.<<endl;
X[k][i] = x;
Y[k][i] = y;
//----comments----------
//  SATE cordinate starts from strip center
// but ICM coordinate starts from strip edge
int kkk = 16 - k;
Ey[kkk][i] = ey/100.;  // V/m-->V/cm
Ex[kkk][i] = - ex/100.;  // V/m-->V/cm
i++;
if(i<115) continue;
k++;
i = 0;
}

fin.close();


int m = 17;//
int n = 115;//
//--------------------------------------
char cid[20];
sprintf(cid,"%d",VFD);
if(VFD<0) sprintf(cid,"M%d",-VFD);
fout<<"float EyArray_"<<cid<<"_"<<HV<<"[17][115] = {"<<endl;
//
for (int k=0; k<m; k++){
  double x = k*2.5;
//fout<<"//--- PotentialArray["<<k<<"][115] , x = "<<x<<" um --------"<<endl;
  fout<<"//--- EyArray_"<<cid<<"_"<<HV<<"["<<k<<"][115] , x = "<<x<<" um --------"<<endl;
  fout << "{"<<std::endl;
  for(int j=0;j<n;j+=5) { 
   fout<<Form("%9.2f,%9.2f,%9.2f,%9.2f,%9.2f,",
     Ey[k][j],Ey[k][j+1],Ey[k][j+2],Ey[k][j+3],Ey[k][j+4]);
   if(j%2==1) fout<<std::endl;
  }
  fout<<"},"<<std::endl;
}
fout<<" };" << endl;
//fout<<"if(VFD=="<<VFD<<" && HV=="<<HV<<") Ey = EyArray_"<<VFD<<"_"<<HV<<"[ix][iy];"<<endl;
fout<<" "<<endl;
fout<<" "<<endl;

fout.close();
return;

}