#include "../../Ramo/Ramo_v3.h"

void plot_Ey(const int VFD, const int VB) {
//
m_isSCTDigiModel = true;
m_EfieldModel = 1; 
char sign[2], Minus[2]="M";
int NEFF;
//double neff;
//sscanf(Neff,"%c%2dE11",sign,&NEFF);
//neff = NEFF*1.e11;
//if(sign[0]==Minus[0]) neff *= -1.;
//double VD = neff/(2.*11.6*8.854e-14)*m_e*m_bulk_depth*m_bulk_depth;
//cout<<Neff<<"  -->"<<sign[0]<<" , neff="<<neff<<" VD="<<VD<<"V"<<endl;  

m_VD=VFD;
m_VB=VB;
//
initialize();
//
const int n=115;
const int m=17;
double X[m][n], Y[m][n], Potential[m][n], E[m][n],j[m][n], PDens[m][n], Ey[m][n], Ex[m][n], Jx[m][n], Jy[m][n], Epsi[m][n], Rho[m][n], dummy; 
char line[200];
sprintf(m_cid, "SCT_%d_%d.txt",VFD,VB);
sprintf(m_cid, "FEM_dummy.txt");
ifstream fin(m_cid);
fin.getline(line, sizeof(line));
for (int k=0; k<m; k++){
  for (int i=0; i<n; i++){
    fin >> X[k][i]>>Y[k][i]>>dummy>>dummy>>Ex[k][i]>>Ey[k][i]>>E[k][i]>>dummy>>dummy>>dummy>>dummy; 
    E[k][i] = E[k][i]/100.;
    Ex[k][i] = Ex[k][i]/100.;
    Ey[k][i] = Ey[k][i]/100.;
//    cout <<"k,i="<<k<<" "<<i<<" x,y="<<X[k][i] << " "<<Y[k][i]<<" Ey="<<Ey[k][i]<<endl; 
  }
} 
//-----------------------------------------
double d = 0.0285;           // [cm]
double xmin=0, xmax=285.;
double ymin=-5000., ymax = 9000.;
TCanvas *c1;
c1 = new TCanvas("c1","A Simple Graph",500,30,1400,1000);
c1->SetGrid();
TGraph *g1;
TLegend *leg1, *leg2, *leg;
leg = new TLegend(0.58 ,0.88,0.81,0.93);
leg->SetHeader("distance from strip center");leg->SetBorderSize(0);
leg->SetFillColor(10);leg->SetTextSize(0.035);leg->Draw();
leg1 = new TLegend(0.58 ,0.67 ,0.71,0.88);
leg1->SetFillColor(10);leg1->SetBorderSize(0);leg1->SetTextFont(132);
leg2 = new TLegend(0.71 ,0.71 ,0.84,0.88);
leg2->SetFillColor(10);leg2->SetBorderSize(0);leg2->SetTextFont(132);
char cid[100];
for (int k=0; k<m; k++){
  int kk=k/2.0;
  if(fmod(k,2.0)>0) continue;
  g1= new TGraph(n,&Y[k][0],&Ey[k][0]);
  g1->GetXaxis()->SetTitle("depth [#mum]");
  g1->GetYaxis()->SetTitle("Vertical Electric Field [V/cm]");
  g1->GetYaxis()->SetTitleOffset(1.5);
  g1->GetXaxis()->SetTitleOffset(1.0);
  g1->SetMarkerStyle(20);
  g1->SetMarkerSize(0.8);
  g1->SetMarkerColor(kk+1);
  g1->SetMinimum(ymin);
  g1->SetMaximum(ymax);
  g1->GetXaxis()->SetLimits(xmin, xmax);
  if(k==0) g1->Draw("AP");
  else g1->Draw("SAME P");
  sprintf(cid," %d #mum",kk*5); 
  if (kk<5)leg1->AddEntry(g1,cid,"P"); 
  else leg2->AddEntry(g1,cid,"P"); 
}
leg->Draw(); 
leg1->Draw(); 
leg2->Draw();
//-----------------flat diode model------------------------
double yy[n],Exx, Eyy[n];
for (int iy=0; iy<n; iy++) {
   yy[iy] = iy*2.5;
   EField( 0.004, iy*2.5*0.0001, Exx, Eyy[iy]) ;
}
g1= new TGraph(n,yy,Eyy);
g1->SetMarkerSize(0.8);
g1->SetMarkerColor(1);
g1->SetMarkerStyle(20);
g1->Draw("SAME P");
//
//---------------------------------
TPave *b;
TLatex *t;
double x1=(xmax-xmin)*22./285 + xmin;
double x2=(xmax-xmin)*150./285 + xmin;
double y1=(ymax-ymin)*0.73 + ymin;
double y2=(ymax-ymin)*0.96 + ymin;
b = new TPave(x1,y1, x2, y2); b->SetFillColor(10);
b->SetBorderSize(2); b->SetLineColor(1); b->Draw();
x1=(xmax-xmin)*28./285 + xmin; y1=(ymax-ymin)*0.90 + ymin;
sprintf(cid,"V_{FD} = %dV", VFD);
//t=new TLatex(x1, y1,cid); t->SetTextSize(0.04); t->Draw();
y1=(ymax-ymin)*0.90 + ymin;
sprintf(cid,"V_{FD}= %5.1fV, HV=%4.1fV", m_VD, m_VB);
t=new TLatex(x1, y1,cid); t->SetTextSize(0.04); t->Draw();
y1=(ymax-ymin)*0.83 + ymin;
t=new TLatex(x1, y1,"colored dots : FEM solutions");
t->SetTextColor(4); t->SetTextSize(0.04); t->Draw();
y1=(ymax-ymin)*0.76 + ymin;
t=new TLatex(x1, y1, "black dots : Flat diode model");
t->SetTextSize(0.04); t->Draw();
//
//---------------------------plot printout------------
sprintf(cid,"Ey_%d_%d.png",int(m_VD), int(m_VB));
c1->Print(cid); 
}