{
gROOT->LoadMacro("AtlasUtils.C");
//ATLAS->SetPadRightMargin(0.12);
c1 = new TCanvas("c1","A Simple Graph ",400,0,650,500);
// create the arrays for the points
int Bec=0;
int Layer=7;
int Short=1;
int chipNO=11;
int kmax = 4;
int plotCiS = 1;
//int runs[kmax]={158393,172632};
int runs[kmax]={160510,166544,167165,168330};
int y[kmax]={2010,2010,2010,2010};
int m[kmax]={8,10,10,11};
int d[kmax]={1,9,21,4};
double chipNumber[12]={0,1,2,3,4,5,6,7,8,9,10,11};
double ymax =1800.;
double ymin=1200.;
if(Short==1){ymax-=500.; ymin-=500.;}
if(Bec==0 && Layer==6) {ymax+=100.; ymin+=100.;}
if(Layer==8) {ymax-=400.; ymin-=400.;}
int nevt[kmax]={0,0,0,0};
long double chipID[kmax][2000];
double noiseData[kmax][2000];
h0 = new TH2F("h0","h0 title",2000,ymin,ymax,2000,0.85,1.15);
h1 = new TH2F("h1","h1 title",2000,ymin,ymax,2000,0.85,1.15);
h2 = new TH2F("h2","h2 title",2000,ymin,ymax,2000,0.85,1.15);
h3 = new TH2F("h3","h3 title",2000,ymin,ymax,2000,0.85,1.15);
int r,bec,layer;
long double chip_id;
double noise[12];
//
ifstream fin("chip_150V.txt"); 
while (fin>>r>>bec>>layer>>chip_id>>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<kmax;i++){
    if(r!=runs[i]) continue;
       int nn= nevt[i];
       chipID[i][nn] = chip_id;
       noiseData[i][nn] = noise[chipNO];
       nevt[i]++;
  }
}
fin.close();
//
ifstream fin("CiS.txt"); 
int dummy, iCiS=0, CiSID2[1000], CiSID3[1000];
while (fin>>chip_id>>dummy>>dummy>>dummy>>dummy>>dummy){
   int i1=20220;
   int i2=(chip_id-float(i1)*1.e9)/1000.;
   int i3=chip_id-float(i1)*1.e9-i2*1000.;
   CiSID2[iCiS] = i2;
   CiSID3[iCiS] = i3;
//cout<<iCiS<<"  20220 "<<i2<<" "<<i3<<endl;
iCiS++;
}
fin.close();
//
cout<<nevt[0]<<" "<<nevt[1]<<" "<<nevt[2]<<" "<<nevt[3]<<endl;
int CiSchip = 0;
for(int i=0;i<nevt[0];i++){
   int i1=20220;
   int i2=(chipID[0][i]-float(i1)*1.e9)/1000.;
   int i3=chipID[0][i]-float(i1)*1.e9-i2*1000.;
// remove CiS sensors
   int CiS =0;
   if(Layer>6) {
   for (int j=0; j<iCiS; j++){
      if(i2!=CiSID2[j] || i3!=  CiSID3[j]) continue;
      CiS = 1;
      CiSchip++;
      break;
   }
//   cout<<i<<" "<<i2<<" "<<i3<<" -> "<<CiS<<" "<<CiSchip<<endl;
   }
   if(CiS != plotCiS) continue;
   for (int k=0; k<kmax; k++){
      for (int j=0; j<nevt[k]; j++) {
        int j1=20220;
        int j2=(chipID[k][j]-float(j1)*1.e9)/1000.;
        int j3=chipID[k][j]-float(j1)*1.e9-j2*1000.;
         if(i2!=j2 || i3!=j3) continue;
         if(k==0)h0->Fill(noiseData[0][i],noiseData[k][j]/noiseData[0][i]);
         if(k==1)h1->Fill(noiseData[0][i],noiseData[k][j]/noiseData[0][i]);
         if(k==2)h2->Fill(noiseData[0][i],noiseData[k][j]/noiseData[0][i]);
         if(k==3)h3->Fill(noiseData[0][i],noiseData[k][j]/noiseData[0][i]);
      }
   }
}
cout<<"total="<<nevt[0]<<"  No. of CiS chip identified ="<<CiSchip<<endl;
//
h0->SetMarkerColor(1);
h0->SetMarkerSize(0.6);
h1->SetMarkerSize(0.6);
h2->SetMarkerSize(0.6);
h3->SetMarkerSize(0.6);
h0->GetXaxis()->SetNdivisions(505); 
h0->GetYaxis()->SetNdivisions(510); 
h0->GetXaxis()->SetLabelSize(0.04);
h0->GetYaxis()->SetLabelSize(0.04);
h0->GetYaxis()->SetLabelOffset(0.02);
h0->GetXaxis()->SetLabelOffset(0.02);
cid = new char[29];cid[28]=0;
sprintf(cid,"noise_{0} (%d.%d.%d)",y[0],m[0],d[0]);
h0->GetXaxis()->SetTitle(cid);
cid = new char[29];cid[28]=0;
h0->GetYaxis()->SetTitle("noise / noise_{0}");
h0->GetXaxis()->SetTitleSize(0.045);
h0->GetYaxis()->SetTitleSize(0.045);
h0->GetYaxis()->SetTitleOffset(1.5);
h0->Draw();
h1->SetMarkerColor(2);
h1->Draw("SAME P");
h2->SetMarkerColor(3);
h2->Draw("SAME P");
h3->SetMarkerColor(4);
h3->Draw("SAME P");
//
cid = new char[35];cid[34]=0;
if(plotCiS ==0) {
  if(Bec==1) sprintf(cid,"Barrel B%d, chip=%d, Hamamatsu",Layer, chipNO);
  else {
  if(Layer==6) sprintf(cid,"Endcap outer, chip=%d, Hamamatsu",chipNO);
  if(Layer==7) sprintf(cid,"Endcap middle, chip=%d Hamamatsu",chipNO);
  if(Layer==8) sprintf(cid,"Endcap inner, chip=%d Hamamatsu",chipNO);
  }
} else{
  if(Bec==1) sprintf(cid,"Barrel B%d, chip=%d, CiS",Layer, chipNO);
  else {
  if(Layer==6) sprintf(cid,"Endcap outer, chip=%d, CiS",chipNO);
  if(Layer==7) sprintf(cid,"Endcap middle, chip=%d, CiS",chipNO);
  if(Layer==8) sprintf(cid,"Endcap inner, chip=%d, CiS",chipNO);
  }
}
leg=new TLegend(0.45,0.68,0.92,0.93,cid);
leg->SetFillColor(0); leg->SetBorderSize(1);
cid = new char[29];cid[28]=0;
sprintf(cid,"noise_{%d.%d.%d} / noise_{0}",y[0],m[0],d[0]);
leg->AddEntry(h0,cid,"p");
cid = new char[29];cid[28]=0;
sprintf(cid,"noise_{%d.%d.%d} / noise_{0}",y[1],m[1],d[1]);
leg->AddEntry(h1,cid,"p");
cid = new char[29];cid[28]=0;
sprintf(cid,"noise_{%d.%d.%d} / noise_{0}",y[2],m[2],d[2]);
leg->AddEntry(h2,cid,"p");
cid = new char[29];cid[28]=0;
sprintf(cid,"noise_{%d.%d.%d} / noise_{0}",y[3],m[3],d[3]);
leg->AddEntry(h3,cid,"p");
leg->SetTextSize(0.04);
leg->Draw();
//
line=new TLine(ymin,1., ymax,1.);
line->Draw();
//
cid = new char[29];cid[28]=0;
sprintf(cid,"noise_2Dr_%d_%d_%d.png",Bec,Layer,chipNO);
if(plotCiS!=0)sprintf(cid,"noise_2Dr_%d_%d_%d_cis.png",Bec,Layer,chipNO);
if(Short==1) {
sprintf(cid,"noise_2Dr_%d_%d_%d_s.png",Bec,Layer,chipNO);
if(plotCiS!=0)sprintf(cid,"noise_2Dr_%d_%d_%d_cis_s.png",Bec,Layer,chipNO);
}
c1->Print(cid);
}