{
gROOT->LoadMacro("AtlasUtils.C");
//ATLAS->SetPadRightMargin(0.12);
c1 = new TCanvas("c1","A Simple Graph ",400,30,700,500);
// create the arrays for the points
TDatime db(2010,02,15,00,00,00);
gStyle->SetTimeOffset(db.Convert());
int t0 = db.Convert();
TDatime d1(2010,3,01,00,00,00);
int t1 = d1.Convert() - t0;
TDatime d2(2011,4,15,00,00,00);
int t2 = d2.Convert() - t0;
double run[900],bec[900],layer[900],noise[12][900];
//
int becID=2;
int link=1;
int Layer=8;
//
int k = 0;
double ymax = 58.0; double ymin=53.5;
//
ifstream fin("gain_150V.txt"); 
while (fin>>run[k]>>bec[k]>>layer[k]>>noise[0][k]>>noise[1][k]>>noise[2][k]>>noise[3][k]>>noise[4][k]>>noise[5][k]>>noise[6][k]>>noise[7][k]>>noise[8][k]>>noise[9][k]>>noise[10][k]>>noise[11][k]){
if(bec[k]!=becID)continue;
if(layer[k]!=Layer)continue;
cout<<k<<" "<<run[k]<<" "<<bec[k]<< " "<<layer[k]<<" "<<noise[0][k]<<endl;
k++;
}
cout<<"numder of data points="<<k-1<<endl;
fin.close();
//
//
double time[900],timeL[900],run1[900], integL[900];
int year,month,day,hour,min;
int n = 0;
ifstream fin("run_date_intL.txt"); 
while (fin>>run1[n]>>year>>month>>day>>hour>>min>>integL[n]) {
if( day> 2000) {
   int temp=day;
   day = month;
   month = year;
   year = temp-2000;
}
if(n<20)cout<<n<<" "<< run1[n]<<" "<<year<<" "<<month<<" "<<day<<" "<<hour<<" "<<min<<" "<<integL[n]<<endl;
TDatime da(2000+year,month,day,hour, min,0);
timeL[n] = da.Convert() - t0;
integL[n] = (ymax-ymin)* 0.8/80000.* integL[n] + ymin;
n++;
}
fin.close();
//
int nrunDate = n;
int krun= k;
//
//
for (k=0; k<krun; k++) {
  for (n=0; n<nrunDate; n++) {
     if( run1[n] <= run[k] ) break;
  }
  time[k]=timeL[n];
}
//
TGraph *gr0 = new TGraph(k,time,noise[0+link*6]);
TGraph *grL = new TGraph(nrunDate,timeL,integL);
gr0->GetXaxis()->SetTitle("run number");
gr0->GetYaxis()->SetTitle("< gain >");
gr0->GetXaxis()->SetNdivisions(207);
gr0->GetXaxis()->SetTitle("date");
gr0->GetXaxis()->SetLabelOffset(0.03);
gr0->GetXaxis()->SetTitleOffset(1.3);
gr0->GetYaxis()->SetLabelOffset(0.010);
gr0->GetYaxis()->SetLabelSize(0.04);
gr0->GetYaxis()->SetTitleSize(0.05);
gr0->GetXaxis()->SetLabelSize(0.035);
gr0->SetMaximum(ymax);
gr0->SetMinimum(ymin);
gr0->SetMarkerColor(1);
gr0->SetMarkerSize(0.8);
gr0->GetXaxis()->SetTimeDisplay(1);
gr0->GetXaxis()->SetTimeFormat("#splitline{ %Y}{%b %d}");
//gr0->GetXaxis()->SetLimits(0,t2);
gr0->Draw("ALP");
cid=new char[29];
if(becID==1)sprintf(cid,"Barrel B%d , link %d",Layer,link);
if(becID==0)sprintf(cid,"Endcap A , eta %d , link %d",Layer,link);
if(becID==2)sprintf(cid,"Endcap C , eta %d , link %d",Layer,link);
cid[28]=0;
leg = new TLegend(0.18 ,0.18,0.46,0.45, cid);
cid = new char[9];sprintf(cid,"chip %d",link*6);cid[8]=0;
leg->AddEntry(gr0,cid,"pl");
for(int j=1;j<6;j++){
TGraph *gr1 = new TGraph(k,time,noise[j+link*6]);
int ic=j+1;
if(j>=4) ic = j+2;
gr1->SetMarkerColor(ic);
//gr1->SetLineColor(ic);
gr1->SetMarkerSize(0.8);
cid = new char[9];sprintf(cid,"chip %d",j+link*6);cid[8]=0;
leg->AddEntry(gr1,cid,"pl");
gr1->Draw("SAME LP");
}
//
leg->SetFillColor(10); leg->SetBorderSize(0); leg->Draw();
//
grL->SetLineStyle(2);
grL->SetLineColor(2);
grL->Draw("SAME L");
leg = new TLegend(0.70 ,0.85,0.94,0.93, "SCT calb. nptgain, <gain>");
leg->AddEntry(grL,"int. luminosity (del.)","l");
leg->SetFillColor(10); leg->SetBorderSize(0); leg->Draw();
//
cid = new char[29];
sprintf(cid,"gainPlot_%d_%d_%d.png",becID,Layer,link);
cid[28]=0;
c1->Print(cid);
}