{
gROOT->LoadMacro("AtlasUtils.C");
//ATLAS->SetPadRightMargin(0.12);
c1 = new TCanvas("c1","A Simple Graph ",300,0,700,700);
// create the arrays for the points
int Bec=2; // 0 for C side(becDB=-2) 2 for A side(becDB=2) 
int Layer=8;
int Short=0;
int chipNO=6;
//------compare two runs--------------
int kmax = 2;
//int runs[kmax]={158393,172621};
//int runs[kmax]={158393,168330};
int runs[kmax]={160510,168330};
int y[kmax]={2010,2010};
int m[kmax]={8,11};
int d[kmax]={1,4};
double chipNumber[12]={0,1,2,3,4,5,6,7,8,9,10,11};
double ymax =1700.;
double ymin=1200.;
if(Short==1){ymax-=500.; ymin-=500.;}
if(Bec!=1 && Layer==6) {ymax+=200.; ymin+=0.;}
if(Layer==8) {ymax-=400.; ymin-=400.;}
int nevt[2]={0,0};
double chipID[2][2000],noiseData[2][2000][12];
//
int r,bec,layer;
double chip_id,noise[12];
//-------------------read the module database-----------------------
int iDB=0,becDB[3300],diskDB[3300],phiDB[3300],etaDB[3300];
double chip_idDB[3300], dummy;
ifstream fin("module_DB.txt"); 
while (fin>>chip_idDB[iDB]>>becDB[iDB]>>diskDB[iDB]>>phiDB[iDB]>>etaDB[iDB]>>dummy>>dummy>>dummy>>dummy>>dummy){
iDB++;
}
fin.close();
cout<<"iDB="<<iDB<<endl;
//-----------------------------------------------------------------
//---------------------------read noise data-----------------------
ifstream fin("chip_150V_2runs.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!=Bec)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;
       for (int j=0;j<12;j++) noiseData[i][nn][j] = noise[j];
       nevt[i]++;
  }
}
fin.close();
cout<<"0, 1="<<nevt[0]<<" "<<nevt[1]<<endl;
//
for (int idisk=0;idisk<9;idisk++){
  if(Layer==7 && idisk==7){ymax-=500.; ymin-=500.;}
  h= new TH2F("h","h title",2000,ymin,ymax,2000,ymin,ymax);
  h1= new TH2F("h1","h1 title",2000,ymin,ymax,2000,ymin,ymax);
  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.; 
// check if this is required disk
    int phi;
    int isModule = 0;
    for (int k=0; k<iDB; k++){     
      if( diskDB[k] != idisk) continue;
      int k1=20220;
      int k2=(chip_idDB[k]-float(k1)*1.e9)/1000.;
      int k3=chip_idDB[k]-float(k1)*1.e9-k2*1000.; 
      if( k2!=i2 || k3!=i3) continue;
      isModule = 1;
      phi=phiDB[k];
    }
    if(isModule==1) {
      for (int j=0; j<nevt[1]; j++) {
      int j1=20220;
      int j2=(chipID[1][j]-float(j1)*1.e9)/1000.;
      int j3=chipID[1][j]-float(j1)*1.e9-j2*1000.;
      if(i2!=j2 || i3!=j3) continue;
      double ratio = noiseData[1][j][6]/noiseData[0][i][6];
      if(fmod(phi,2.0)==0)h->Fill(noiseData[0][i][chipNO],noiseData[1][j][chipNO]);
      if(fmod(phi,2.0)!=0)h1->Fill(noiseData[0][i][chipNO],noiseData[1][j][chipNO]);
    cout<<"idisk,i="<<idisk<<" "<<i<<" "<<k<<" "<<j<<" "<<phi<<" "<<noiseData[0][i][chipNO]<<" "<<ratio<<endl;
      }
    }
  }
//
  h->SetMarkerSize(0.7);
  h1->SetMarkerSize(0.7);
  h->GetXaxis()->SetNdivisions(505); 
  h->GetYaxis()->SetNdivisions(505); 
  h->GetXaxis()->SetLabelSize(0.04);
  h->GetYaxis()->SetLabelSize(0.04);
  h->GetYaxis()->SetLabelOffset(0.02);
  h->GetXaxis()->SetLabelOffset(0.02);
  cid = new char[29];cid[28]=0;
  sprintf(cid,"noise_{0}(%d.%d.%d)",y[0],m[0],d[0]);
  h->GetXaxis()->SetTitle(cid);
  cid = new char[29];cid[28]=0;
  sprintf(cid,"noise(%d.%d.%d)",y[1],m[1],d[1]);
  h->GetYaxis()->SetTitle(cid);
  h->GetXaxis()->SetTitleSize(0.045);
  h->GetYaxis()->SetTitleSize(0.045);
  h->GetYaxis()->SetTitleOffset(1.7);
  h->Draw();
  h1->SetMarkerColor(2);
  h1->Draw("SAME");
//
  leg=new TLegend(0.2,0.7,0.6,0.85);
  leg->SetFillColor(0);
  leg->SetBorderSize(0);
  leg->AddEntry(h,"even #phi numbers","p");
  leg->AddEntry(h1,"odd #phi numbers","p");
  double x1 = ymin/0.93;
  double y2 = ymax*0.93;
  z=new TLine(ymin,ymin,ymax,ymax);
  z->Draw();
  //leg->AddEntry(z,"  y = x","l");
  z=new TLine(x1,ymin,ymax,y2);
  z->SetLineColor(2); z->Draw();
  //leg->AddEntry(z,"  y = 0.93 * x","l");
  leg->Draw();

//
  cid = new char[35];cid[34]=0;
  if(Bec==1) sprintf(cid,"Barrel B%d, disk=%d,  chip=6",Layer, idisk);
  else {
    if(Bec==0){
     if(Layer==6) sprintf(cid,"Endcap C outer, disk=%d chip=6",idisk);
     if(Layer==7) sprintf(cid,"Endcap C middle, disk=%d chip=6",idisk);
     if(Layer==8) sprintf(cid,"Endcap C inner, disk=%d chip=6",idisk);
    } else {
     if(Layer==6) sprintf(cid,"Endcap A outer, disk=%d chip=6",idisk);
     if(Layer==7) sprintf(cid,"Endcap A middle, disk=%d chip=6",idisk);
     if(Layer==8) sprintf(cid,"Endcap A inner, disk=%d chip=6",idisk);
    }
  }
  l=new TLatex(ymin+50., ymax-50., cid);
  l->SetTextSize(0.040);l->Draw();
//
  cid = new char[29];cid[28]=0;
  sprintf(cid,"noise_2D_%d_%d_d%d.png",Bec,Layer,idisk);
  if(Short!=0)sprintf(cid,"noise_2D_%d_%d_d%d_s.png",Bec,Layer,idisk);
  c1->Print(cid);
//
}
}