{
gROOT->LoadMacro("AtlasUtils.C");
//ATLAS->SetPadRightMargin(0.12);
c1 = new TCanvas("c1","A Simple Graph ",300,0,700,500);
// create the arrays for the points
int bec,layer;
int Bec=0;
int Layer=7;
int Short=1;
double chipID;
int k = 8;
int runs[k]={158393,165169,165899,166544,167165,168330,170137,172621},r;
int y[k]={2010,2010,2010,2010,2010,2010,2010,2011};
int m[k]={7,9,9,10,10,11,12,1};
int d[k]={1,16,29,9,21,4,2,21};
double chipNumber[12]={0,1,2,3,4,5,6,7,8,9,10,11};
double ymax =1.1;
double ymin=0.85;
double noise[12],nevt[k][12],nsum[k][12],nsum2[k][12];
for(int i=0;i<k;i++){
  for(int j=0;j<12;j++) {nevt[i][j]=0.; nsum[i][j]=0.; nsum2[i][j]=0.;}
}
//
ifstream fin("chip_150V.txt"); 
while (fin>>r>>bec>>layer>>chipID>>noise[0]>>noise[1]>>noise[2]>>noise[3]>>noise[4]>>noise[5]>>noise[6]>>noise[7]>>noise[8]>>noise[9]>>noise[10]>>noise[11]){
  if(Bec==1 && bec!=Bec)continue;
  if(Bec==0 && bec==1)continue;
  if(layer!=Layer)continue;
  if(Layer==7 && Short==0 && noise[1]<1000.) continue;
  if(Layer==7 && Short==1 && noise[1]>=1000.) continue;
  for (int i=0;i<k;i++){
    if(r!=runs[i]) continue;
    for(int j=0;j<12;j++) {
       nevt[i][j]++;
       nsum[i][j]+=noise[j];
       nsum2[i][j]+=noise[j]*noise[j];
    }
  }
}
fin.close();
//--- do averaging-----
for(int i=0;i<k;i++){
   for(int j=0;j<12;j++){
      nsum[i][j]/=nevt[i][j];
      nsum2[i][j]=sqrt(nsum2[i][j]/nevt[i][j]-nsum[i][j]*nsum[i][j]);
   //cout<<i<<" "<<j<<"="<<nevt[i][j]<<" "<<nsum[i][j]<<" "<<nsum2[i][j]<<endl;
   }
}
//
double y1=0.72;
if(Layer==8) y1=0.22;
leg = new TLegend(0.19 , y1 ,0.55, y1+0.22);
leg2 = new TLegend(0.39, y1 ,0.75, y1+0.22);
//
for (int i=k-1;i>=0;i--) {
   for(int j=0;j<12;j++) {nsum[i][j] = nsum[i][j]/nsum[0][j];}
}
//
for (int i=0;i<k;i++) {
gr0 = new TGraph(12,chipNumber,nsum[i]);
gr0->GetYaxis()->SetTitle("<noise> / <noise>_{0}");
gr0->GetXaxis()->SetTitle("chip number");
gr0->GetXaxis()->SetTitleOffset(1.0);
gr0->GetYaxis()->SetLabelOffset(0.015);
gr0->GetYaxis()->SetLabelSize(0.045);
gr0->GetXaxis()->SetLabelSize(0.045);
gr0->SetMaximum(ymax);
gr0->SetMinimum(ymin);
gr0->SetMarkerColor(1);
gr0->SetMarkerSize(1);
gr0->SetMarkerColor(i+1);
gr0->GetXaxis()->SetLimits(-0.5,11.5);
if(i==0) gr0->Draw("ALP");
if(i!=0) gr0->Draw("SAME LP");
cid = new char[25];sprintf(cid,"%d-%d-%d",y[i],m[i],d[i]);cid[24]=0;
if(i<4)leg->AddEntry(gr0,cid,"pl");
if(i>=4)leg2->AddEntry(gr0,cid,"pl");
}
//
leg->SetFillColor(0);leg->SetBorderSize(0);
leg->Draw();
leg2->SetFillColor(0);leg2->SetBorderSize(0);
leg2->Draw();
cid = new char[39];sprintf(cid,"Barrel B%d",Layer);cid[38]=0;
if(Bec!=1 && Layer==6)sprintf(cid,"Endcap outer",Layer);
if(Layer==7)sprintf(cid,"Endcap middle",Layer);
if(Layer==8)sprintf(cid,"Endcap inner",Layer);
z=new TLatex(7.5,1.075,cid);
z->SetTextSize(0.05);z->Draw();
if(Short==1) {
z=new TLatex(9.3,1.055,"(short)");
z->SetTextSize(0.05);z->Draw();
}
//
cid = new char[29];
sprintf(cid,"noise_ratio_%d_%d.png",Bec,Layer);
if(Layer==7 && Short==1) sprintf(cid,"noise_ratio_%d_%d_s.png",Bec,Layer);
cid[28]=0;
c1->Print(cid);
}