{
#include <fstream.h>
gROOT->LoadMacro("AtlasUtils.C");
ATLAS->SetPadTickX(1);
ATLAS->SetPadTickY(0);
ATLAS->SetPadRightMargin(0.12);
ATLAS->SetPadLeftMargin(0.12);
c1 = new TCanvas("c1","A Simple Graph",500,0,700,500);
//c1->SetFillStyle(4100);
TDatime db(2010,08,01,00,00,00);
int t1 = db.Convert();
TDatime db(2011,06,01,00,00,00);
int t2 = db.Convert();
gStyle->SetOptStat(0);
c1->SetGrid();
// create the arrays for the points
Double_t ymax= 1.0;
double Tscale=ymax*0.1/10., Toffset = 0.3;
Int_t n = 300;
int i=0, iT=0;
double intL[n],T3[n],T4[n],T5[n],T6[n],leak[4][n],leakop[4][n],leak3[5][n],time[n],
Ttime[n],weeklyL, T3temp,T4temp,T5temp, T6temp; 
int week,year1,month1,day1,year2,month2,day2; 
//
ifstream fin("leak_prediction_Barrel_17.5C.txt");
while(fin>>year1>>month1>>day1>> year2>>month2>>day2>>intL[i]>>T3temp>>T4temp>>
T5temp>>T6temp>> leak[0][i]>>leak[1][i]>>leak[2][i]>>leak[3][i]>>leak3[0][i]>>leak3[1][i]>>leak3[2][i]>>leak3[3][i]>>leak3[4][i]>>leakop[0][i]>>leakop[1][i]>>leakop[2][i]>>leakop[3][i] )
{
intL[i]= intL[i]/1000.; 
TDatime d1(year1,month1,day1,12, 0,0);
int time1=d1.Convert();
TDatime da(year2,month2,day2,12, 0,0);
time[i]=da.Convert();
Ttime[iT]=time1;
T3[iT]=T3temp*Tscale+Toffset; 
T4[iT]=T4temp*Tscale+Toffset;
T5[iT]=T5temp*Tscale+Toffset; 
T6[iT]=T6temp*Tscale+Toffset;
iT++;
Ttime[iT]=time[i];
T3[iT]=T3temp*Tscale+Toffset; 
T4[iT]=T4temp*Tscale+Toffset;
T5[iT]=T5temp*Tscale+Toffset; 
T6[iT]=T6temp*Tscale+Toffset;
iT++;
//cout<<i<<" "<<time[i]<<" "<<intL[i]<<" "<<leak[0][i]<<endl;
i++;
}
fin.close();
//
int j=0,idata[4]={0,0,0,0};
int bec, layer;
double tempdata,temptime,dtime[4][100],data[4][100],error;
ifstream fin("leak_data.txt");
while(fin>>bec>>layer>>year2>>month2>>day2>>tempdata>>error)
{
//cout<<bec<<" "<<layer<<" "<<year2<<" "<<month2<<" "<<day2<<" "<<tempdata<<endl;
if(bec!=1)continue;
int L=layer-3;
if(L<0 || L>3) continue;
int m = idata[L];
TDatime da(year2,month2,day2,12, 0,0);
dtime[L][m]=da.Convert();
data[L][m]=tempdata;
cout<<L<<" "<<m<<" "<<dtime[L][m]<<" "<<data[L][m]<<endl;
idata[L]++;
}
fin.close();
//
leg1 = new TLegend(0.15 ,0.60 ,0.48,0.94,"B3 data and prediction");
for (int k=0; k<6; k++){
for (int j=0; j<i; j++) {
 if(k==0)leak[k][j] /= 1000.;
 if(k!=0)leak3[k-1][j] /= 1000.;
}
if(k==0) g1 = new TGraph(i,time,leak[k]);
if(k!=0) g1 = new TGraph(i,time,leak3[k-1]);
g1->SetMaximum(ymax);
g1->GetYaxis()->SetTitle("I _{leak} [#muA] (-10#circC)");
//g1->GetXaxis()->SetTitle("DATE");
g1->SetLineColor(k+1);
g1->GetXaxis()->SetTitleOffset(1.4);
g1->GetXaxis()->SetLabelOffset(0.03);
g1->GetXaxis()->SetLabelSize(0.04);
g1->GetYaxis()->SetTitleSize(0.05);
g1->GetYaxis()->SetTitleOffset(1.0);
g1->GetXaxis()->SetRangeUser(13,45);
g1->GetXaxis()->SetTimeDisplay(1);
g1->GetXaxis()->SetTimeFormat("#splitline{%Y}{%b %d}");
g1->GetXaxis()->SetLimits(t1,t2);
g1->GetXaxis()->SetNdivisions(206);
if(k==0)g1->Draw("AL");
if(k!=0)g1->Draw("SAME L");
if(k==0)leg1->AddEntry(g1,"total: data & prediction","LP");
if(k==1)leg1->AddEntry(g1,"i_{1} (0.42 833 days)","L");
if(k==2)leg1->AddEntry(g1,"i_{2} (0.10 28.5 days)","L");
if(k==3)leg1->AddEntry(g1,"i_{3} (0.23  2.6 days)","L");
if(k==4)leg1->AddEntry(g1,"i_{4} (0.21 0.09 days)","L");
if(k==5)leg1->AddEntry(g1,"i_{5} (0.04 .006 days)","L");
if(k!=0) continue;
for (int j=0; j<idata[k]; j++) data[k][j] /= 1000.;
g2 = new TGraph(idata[k],dtime[k],data[k]);
g2->SetMarkerColor(k+1);
g2->SetMarkerSize(1.0);
g2->SetLineColor(k+1);
g2->Draw("SAME P");
}
gT = new TGraph(iT,Ttime,T3);
gT->SetLineStyle(2); gT->SetLineColor(1);gT->Draw("SAME L");
leg1->AddEntry(gT,"sensor temperature (right scale)","L");
//
leg1->SetFillColor(10); leg1->SetBorderSize(0); leg1->SetTextSize(0.04);
leg1->Draw();
//draw an axis on the right side
TGaxis *axis = new TGaxis(t2,0.,t2, ymax,-Toffset/Tscale,(ymax-Toffset)/Tscale,505,"+L");
axis->SetTitle("sensor temperature (#circC)");
axis->SetLabelOffset(0.01);
axis->SetTitleOffset(1.2);
axis->SetLabelFont(42);
axis->SetTitleFont(42);
axis->Draw();
//
TDatime db(2010,6,29, 0,0,0); t1=db.Convert();
z=new TLatex(t1,62,"T_{shutdown}=17.5#circC");
//z->Draw();
c1->Print("leak_current_B3.png");
}