//-------------------------------------
//  2016.2.18 Ey for FDM (flat diode model)
//  2016.2.17 recovere from old ROOT program by Taka Kondo
//            using ROOT 6.02/12
//  2009.11.27 Electric field calculation,  Taka Kondo (KEK) 
//------------------------------------
#include <iostream>
#include <math.h>
#include "sub_functions_v7.C"
#include "induced_charge.h"
#include "GetExEy_150V.cxx"
#include "GetPotential.cxx"

void plot_Ey_FDM() {
//-----seit initial parameters----------------------
m_model = 1;
m_VD = 70. ; // [Volt]
int VD = m_VD;
m_VB = 30. ;  // [Volt]
initialize();
// unit: nsec, cm, sec, K, 
TCanvas *c1 = new TCanvas("c1","A Simple Graph with error bars",0,0,700,500);
c1->SetFillStyle(4100);
c1->SetGrid();
c1->GetFrame()->SetBorderSize(12);
char cid[50];
double Ex, Ey;
double d = m_bulk_depth;  // [
double XX[500], YY[8][500];
//-----------------------------------------
TGraph *g1;
TLegend *leg1, *leg2;
sprintf(cid,"Flat Diode Model, Vd=%d V",VD);
leg1 = new TLegend(0.20 ,0.55 ,0.45,0.92, cid);
for (int k=0; k<8; k++){
  m_VB = 150. - k * 20.;
  m_depletion_depth = sqrt(m_VB/m_VD) * m_bulk_depth ;
  int npoint = 0;
  for (double y=0.; y<m_bulk_depth; y+= 0.0002) {
     XX[npoint] = y * 10000. ;
     EField(40./10000., y , Ex, Ey);    // [V/cm]
     YY[k][npoint] = Ey;
     if(npoint < 10) cout<<m_VB<<" V,  "<<y*10000.<<" "<<Ey<<endl;
     npoint++;
  }
  g1= new TGraph(npoint,XX,YY[k]);
  g1->SetLineColor(k+1);
  g1->SetLineWidth(2);
  g1->SetMaximum(10000.);
  g1->SetMinimum(0.);
  g1->GetXaxis()->SetTitle("depth [#mum]");
  g1->GetXaxis()->SetTitleOffset(1.0);
  g1->GetXaxis()->SetTitleSize(0.05);
  g1->GetXaxis()->SetLabelSize(0.05);
  g1->GetYaxis()->SetTitle("Vertical Electric Field [V/cm]");
  g1->GetYaxis()->SetTitleOffset(1.5);
  g1->GetYaxis()->SetTitleSize(0.05);
  g1->GetYaxis()->SetLabelSize(0.05);
  g1->GetXaxis()->SetLimits(0.,m_bulk_depth*10000.);
  if( k==0 ) g1->Draw("AL");
  g1->Draw("SAME L");
  sprintf(cid,"V_{B}=%4.0fV", m_VB);
  leg1->AddEntry(g1,cid,"L"); 
}
leg1->SetFillColor(10);leg1->SetBorderSize(1);leg1->SetTextSize(0.03);
leg1->Draw(); 
sprintf(cid,"Ey_Vd%d_FDM.png",VD);
c1->Print(cid);
return;
}