{
#include <iostream>
#include <math.h>

// Electric field calculation,  Taka Kondo (KEK) 2009.11.27
// unit: nsec, cm, sec, K, 
gROOT->LoadMacro("AtlasUtils.C");
ATLAS->SetPadRightMargin(0.20);
c1 = new TCanvas("c1","A Simple Graph with error bars",0,0,700,500);
c1->GetFrame()->SetBorderSize(12);
//
TH2F *hE2 = new TH2F("hE2","Electric field2",80,0.,80.,114,0.,285.); 
 
double x,y,Ex,Ey;
ifstream fin("E_150V.DAT"); //150V=Vb
for (int k=0; k<40; k++){
//for (int k=0; k<2; k++){
for (int i=0; i<117; i++){
x = 1.*(k+0.5);
y = 2.5*(i + 0.5);
EField( x/10000., y/10000., Ex, Ey );
if(i==0) cout<<"x,y,Ex,Ey="<<x<<" "<<y<<" "<<Ex<<" "<<Ey<<endl;
hE2->Fill(x,y,Ey);
hE2->Fill(80.-x,y,Ey);
}
}

hE2->SetXTitle("X [#mum]");
hE2->GetXaxis()->SetTitleOffset(1.0);
hE2->SetYTitle("Y [#mum]");
hE2->GetYaxis()->SetTitleOffset(1.2);
gStyle->SetPalette(1);
hE2->Draw("COLZ");
TLatex *t=new TLatex(103.,100.,"E_{y}[V/cm]"); 
t->SetTextAngle(90.);
t->SetTextSize(0.05);
t->Draw();
c1->Print("E_2D.pdf");
}