{
// T. Kondo, 2010.4.19
c1=new TCanvas("c1","canvas",0,0,700,500);
c1->SetGrid();
int n=8;
int m=1;
double X[n], Y[n], Potential[n], E[n],j[n], PDens[n], Ey[n], Ex[n], Jx[n], Jy[n], Epsi[n], Rho[n]; 
ifstream fin("dopeDependence_70V.txt"); //150V=Vb
for (int i=0; i<n; i++){
fin >> X[i]>>Y[i]>>Potential[i]>>E[i]>>j[i]>>PDens[i]>>Ex[i]>>Ey[i]>>Jx[i]>>Jy[i]>>Epsi[i]>>Rho[i]; 
  E[i] = E[i]/100.;
  Ex[i] = Ex[i]/100.;
  Ey[i] = Ey[i]/100.;
  Rho[i] = Rho[i];
  cout <<i<<" x,y="<<X[i] << " "<<Y[i]<<" "<<E[i]<<" Ex="<<Ex[i]<< " Ey=" <<Ey[i]<<" Rho="<< Rho[i]<<endl; 
}
h1 = new TGraph(n,Rho,Ey);
h1->GetXaxis()->SetTitle("Doping concentration [ Coulomb/m^{3} ]");
h1->GetXaxis()->SetTitleOffset(1.0);
h1->GetYaxis()->SetTitle("Electric field at y = 0 [V/cm]");
h1->Draw("ALP");
b = new TPave(0.158,35,0.184,410); 
b->SetFillColor(10); 
b->SetBorderSize(1);
b->Draw(); 
z1=new TLatex(0.16,340., "SCT Barrel geometry");
z1->Draw();
z1=new TLatex(0.16, 270, "bulk depth = 285 #mum");
z1->Draw();
z1=new TLatex(0.16, 200, "strip pitch = 80 #mum");
z1->Draw();
z1=new TLatex(0.16, 130, "strip width = 20 #mum");
z1->Draw();
z1=new TLatex(0.16,  60., "V_{B} = 70 V, #epsilon_{Si} = 11.6");
z1->Draw();
}