{
#include <fstream>
gROOT->LoadMacro("AtlasUtils.C");
c1 = new TCanvas("c1","A Simple Graph",0,0,280,450);
c1->SetFillStyle(4100);
//c1->SetGrid();
Double_t t[200], y[5][200], dummy;
double ymax=0.40, ymin=-0.20;
//---------------charge-------------------------------------
ifstream fin("case_A_charge.txt"); 
//ifstream fin("case_B_charge.txt"); 
Int_t i=0;
while(fin>>t[i]>>y[0][i]>>y[1][i]>>y[2][i]>>y[3][i]>>y[4][i]>>dummy)i++; 
fin.close();
i++;
//
for(int k=0; k<5; k++){
cout<<"i="<<i<<endl;
h = new TGraph(i,t,y[k]);
h->SetTitle();
h->GetXaxis()->SetTitle("time [ns]");
h->GetXaxis()->SetTitleOffset(1.0);
h->GetYaxis()->SetTitle("induced charge");
h->GetYaxis()->SetTitleOffset(1.5);
h->SetMinimum(ymin);
h->SetMaximum(ymax);
h->SetLineColor(1);
h->GetXaxis()->SetRangeUser(0.,25.);
h->Draw("AL");
//
leg=new TLegend(0.30,0.85,0.93,0.90);
cid=new char[25];sprintf(cid,"Strip %d :(x,y)=( 5, 50)",k-2);cid[24]=0;
//cid=new char[25];sprintf(cid,"Strip %d :(x,y)=(40, 50)",k-2);cid[24]=0;
leg->SetHeader(cid);
leg->SetBorderSize(0);leg->SetFillColor(10); 
leg->SetTextSize(0.06); leg->Draw();
cid=new char[25];sprintf(cid,"Charge_%d_A.png",k-2);cid[24]=0;
//cid=new char[25];sprintf(cid,"Charge_%d_B.png",k-2);cid[24]=0;
c1.Print(cid);
}
//---------------pulse-------------------------------------
double ymax=1.20, ymin=-0.2;
ifstream fin("case_A_pulse.txt"); 
//ifstream fin("case_B_pulse.txt"); 
int i=0;
while(fin>>t[i]>>y[0][i]>>y[1][i]>>y[2][i]>>y[3][i]>>y[4][i]>>dummy)i++; 
fin.close();
i++;
//
for(int k=0; k<5; k++){
cout<<"i="<<i<<endl;
h = new TGraph(i,t,y[k]);
h->SetTitle();
h->GetXaxis()->SetTitle("time [ns]");
h->GetXaxis()->SetTitleOffset(1.0);
h->GetYaxis()->SetTitle("pulse");
h->GetYaxis()->SetTitleOffset(1.5);
h->SetMinimum(ymin);
h->SetMaximum(ymax);
h->SetLineColor(1);
h->GetXaxis()->SetRangeUser(0.,50.);
h->Draw("AL");
//
leg=new TLegend(0.30,0.85,0.93,0.90);
cid=new char[25];sprintf(cid,"Strip %d :(x,y)=( 5, 50)",k-2);cid[24]=0;
//cid=new char[25];sprintf(cid,"Strip %d :(x,y)=(40, 50)",k-2);cid[24]=0;
leg->SetHeader(cid);
leg->SetBorderSize(0);leg->SetFillColor(10); 
leg->SetTextSize(0.06); leg->Draw();
cid=new char[25];sprintf(cid,"Pulse_%d_A.png",k-2);cid[24]=0;
//cid=new char[25];sprintf(cid,"Pulse_%d_B.png",k-2);cid[24]=0;
c1.Print(cid);
}
}