{
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=0;
int Layer=6;
int Short=0;
int plotCiS = 0;
//------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==0 && 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];
//-------------------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){
if(Bec==1 && becDB[iDB]!=0) continue;
if(Bec!=1 && becDB[iDB]==0) continue;
iDB++;
}
fin.close();
cout<<"iDB="<<iDB<<endl;
//----------------------read noise data-------------------------------
int r,bec,layer;
double chip_id,noise[12];
//
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==1 && bec!=Bec)continue;
  if(Bec==0 && bec==1)continue;
  if(layer!=Layer)continue;
  if(Layer==7 && Short==0 && noise[1]<1100.) continue;
  if(Layer==7 && Short==1 && noise[1]>=1100.) 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();
//---------------------read CiS modules---------------------------
ifstream fin("CiS.txt"); 
int idummy, iCiS=0, CiSID2[1000], CiSID3[1000];
while (fin>>chip_id>>idummy>>idummy>>idummy>>idummy>>idummy){
   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();
//-----------------------------------------------------------------
for (chipNO=0;chipNO<12;chipNO++){
int total=0;
h= new TH2F("h","h title",2000,ymin,ymax,2000,ymin,ymax);
h2= new TH2F("h2","h2 title",2000,ymin,ymax,2000,ymin,ymax);
cout<<nevt[0]<<" "<<nevt[1]<<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.;
//------------------find phi number--------------------------------
    int disk, bec, phi, eta;
    int isModule = 0;
    for (int k=0; k<iDB; k++){     
      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];
      disk = diskDB[k];
      bec = becDB[k];
      eta = etaDB[k];
//    cout<<k2<<" & "<<k3<<" "<<i<<" "<<k<<" "<<phi<<" "<<disk<<endl;
    }
    if(isModule==0) continue;
//------------------remove or take 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++;
      }
   }
   if(CiS != plotCiS) continue;
//----------------------------------------------------------------
   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;
        
//cout<<i2<<"="<<j2<<" "<<i3<<"="<<j3<<":"<<i<<" "<<k<<" "<<j<<" phi="<<phi<<" "<<disk<<" "<<noiseData[1][j][chipNO]/noiseData[0][i][chipNO]<<endl;
        total++;
        if((bec!=-2 || (disk!=2 && disk!=8)) && !(bec==2 && eta==1)) {
          if(fmod(phi,2.0)==0) 
             h->Fill(noiseData[0][i][chipNO],noiseData[1][j][chipNO]);
          else 
             h2->Fill(noiseData[0][i][chipNO],noiseData[1][j][chipNO]);
        } else {
           if(fmod(phi,2.0)==0) 
              h2->Fill(noiseData[0][i][chipNO],noiseData[1][j][chipNO]);
            else 
              h->Fill(noiseData[0][i][chipNO],noiseData[1][j][chipNO]);
        }
   }
}
//----------------------------------------------------------------
cout<<"total="<<total<<endl;
h->SetMarkerSize(0.6);
h2->SetMarkerSize(0.6);
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();
h2->SetMarkerColor(2);
h2->Draw("SAME");
//---------------------------------------------------------------
cid = new char[35];cid[34]=0;
if(plotCiS ==0) {
  if(Bec==1) sprintf(cid,"Barrel B%d, chip=%d",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==7 && Short==1) sprintf(cid,"Endcap middle(short), chip=%d, CiS",chipNO);
  if(Layer==8) sprintf(cid,"Endcap inner, chip=%d, CiS",chipNO);
  }
}
l=new TLatex(ymin+50., ymax-50., cid);
l->SetTextSize(0.040);l->Draw();
//--------------------------print legend etc---------------------
leg=new TLegend(0.17,0.75,0.65,0.86);
leg->SetFillColor(0);leg->SetBorderSize(0);
if(Bec!=1){
leg->AddEntry(h,"even #phi (w/exceptions)","p");
leg->AddEntry(h2,"odd #phi (w/exceptions)","p");
} else {
leg->AddEntry(h,"even #phi","p");
leg->AddEntry(h2,"odd #phi","p");
}
leg->Draw();
//-----------------------------------------------------
leg2=new TLegend(0.70,0.20,0.93,0.30);
leg2->SetFillColor(0);leg2->SetBorderSize(0);
z=new TLine(ymin,ymin,ymax,ymax);
z->Draw();
leg2->AddEntry(z," y=x","l");
z=new TLine(ymin/0.93,ymin,ymax,ymax*0.93);
z->SetLineColor(2);z->Draw();
leg2->AddEntry(z," y=0.93*x","l");
leg2->Draw();
//----------------------------------------------------
cid = new char[29];cid[28]=0;
sprintf(cid,"noise_2D_%d_%d_%d.png",Bec,Layer,chipNO);
if(Short!=0)sprintf(cid,"noise_2D_%d_%d_%d_s.png",Bec,Layer,chipNO);
if(plotCiS==1){
sprintf(cid,"noise_2D_%d_%d_%d_cis.png",Bec,Layer,chipNO);
if(Short!=0)sprintf(cid,"noise_2D_%d_%d_%d_cis_s.png",Bec,Layer,chipNO);
}
c1->Print(cid);
}
//
}