//-------------------------------------
//  2016.4.1  updated
//  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 "Ramo.h"

void plot_Ey_FDM() {
//------------set main parameters and initialization----------------
m_EfieldModel = 1;
m_VD = 70.;
m_VB = 150.;
//m_isSCTDigiModel = true;
initialize();
//------------------------set up plots -----------------
TCanvas *c1 = new TCanvas("c1","A Simple Graph with error bars",0,0,700,500);
c1->SetGrid();
//-----------------------------------------
double Ex, Ey;
double XX[500], YY[8][500];
TGraph *g1;
TLegend *leg1;
sprintf(m_cid,"Flat Diode Model, V_{D}=%dV",int(m_VD));
leg1 = new TLegend(0.20 ,0.58 ,0.45,0.93, m_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("height [#mum]");
  g1->GetXaxis()->SetTitleOffset(1.0);
  g1->GetXaxis()->SetTitleSize(0.04);
  g1->GetXaxis()->SetLabelSize(0.04);
  g1->GetYaxis()->SetTitle("Vertical Electric Field [V/cm]");
  g1->GetYaxis()->SetTitleOffset(1.5);
  g1->GetYaxis()->SetTitleSize(0.04);
  g1->GetYaxis()->SetLabelSize(0.04);
  g1->GetXaxis()->SetLimits(0.,m_bulk_depth*10000.);
  if( k==0 ) g1->Draw("AL");
  g1->Draw("SAME L");
  sprintf(m_cid,"V_{B}=%4.0fV", m_VB);
  leg1->AddEntry(g1,m_cid,"L"); 
}
leg1->SetFillColor(10);leg1->SetBorderSize(2);leg1->SetTextSize(0.03);
leg1->Draw(); 
sprintf(m_cid,"Ey_VD%d_FDM.png",int(m_VD));
c1->Print(m_cid);
return;
}