{
gROOT->LoadMacro("AtlasUtils.C");
//ATLAS->SetPadRightMargin(0.12);
c1 = new TCanvas("c1","A Simple Graph ",300,0,650,650);
// create the arrays for the points
int Bec=0;
int Layer=8;
int Short=0;
int plotCiS = 0;
int kmax=4;
int runs[kmax]={155271,155410,175603,175710};//for 150V, 50V, 150V, 50V
int y[kmax]={2010,2010,2011,2011};
int m[kmax]={5,5,2,2};
int d[kmax]={18,20,18,19};
double chipNumber[12]={0,1,2,3,4,5,6,7,8,9,10,11};
double xmax =1800.,xmin=1200.,ymin=0.98,ymax=1.20;
if(Short==1){xmax-=500.; xmin-=500.;}
int nevt[kmax]={0,0,0,0};
double chipID[kmax][2000],noiseData[kmax][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];
//
for (int irun=0; irun<kmax; irun++) {
  if(irun==0 || irun==2) ifstream fin("chip_150V.txt"); 
  if(irun==1 || irun==3) ifstream fin("chip_50V.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;
    if(r!=runs[irun]) continue;
       int nn= nevt[irun];
       chipID[irun][nn] = chip_id;
       for (int j=0; j<12; j++) noiseData[irun][nn][j] = noise[j];
       nevt[irun]++;
  }
}
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 (int i=0;i<kmax;i++) {cout<<i<<" "<<nevt[i]<<endl;}
//
for(int ichip=0; ichip<12; ichip++){
h = new TH2F("h","h title",2000,ymin,ymax,2000,ymin,ymax);
h2 = new TH2F("h2","h2 title",2000,ymin,ymax,2000,ymin,ymax);//
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[2]; j++) {
    int j1=20220;
    int j2=(chipID[2][j]-float(j1)*1.e9)/1000.;
    int j3= chipID[2][j]-float(j1)*1.e9-j2*1000.;
    if( i2!=j2 || i3!=j3) continue;
    for (int k2=0; k2<nevt[1];k2++) {
     if(chipID[0][i] != chipID[1][k2] ) continue;
     for( int k3=0;k3<nevt[3];k3++) {
       if(chipID[0][i] != chipID[3][k3] ) continue;
       if((bec!=-2 || (disk!=2 && disk!=8)) && !(bec==2 && eta==1)) {
            if((Bec==1 && ((eta<0 && fmod(-eta,2.0)==0)||(eta>0 && fmod(eta,2.0)!=0)))
               || (Bec!=1 && fmod(phi,2.0)==0))
              h->Fill(noiseData[1][k2][ichip]/noiseData[0][i][ichip],             
                 noiseData[3][k3][ichip]/noiseData[2][j][ichip]);
            else h2->Fill(noiseData[1][k2][ichip]/noiseData[0][i][ichip],
                 noiseData[3][k3][ichip]/noiseData[2][j][ichip]);
        } else {
            if((Bec==1 && ((eta<0 && fmod(-eta,2.0)==0)||(eta>0 && fmod(eta,2.0)!=0)))
               || (Bec!=1 && fmod(phi,2.0)==0))
              h2->Fill(noiseData[1][k2][ichip]/noiseData[0][i][ichip],
                 noiseData[3][k3][ichip]/noiseData[2][j][ichip]);
            else h->Fill(noiseData[1][k2][ichip]/noiseData[0][i][ichip],
                 noiseData[3][k3][ichip]/noiseData[2][j][ichip]);
              
        }    
      }
    }
   }
}
//
h->SetMarkerSize(0.6);
h2->SetMarkerSize(0.6);
h2->SetMarkerColor(2);
h->GetXaxis()->SetNdivisions(510); 
h->GetYaxis()->SetNdivisions(510); 
h->GetXaxis()->SetLabelSize(0.04);
h->GetYaxis()->SetLabelSize(0.04);
h->GetYaxis()->SetLabelOffset(0.02);
h->GetXaxis()->SetLabelOffset(0.02);
h->GetXaxis()->SetTitle("noise@50V/noise@150V, May 2010");
h->GetYaxis()->SetTitle("noise@50V/noise@150V, Feb 2011");
h->GetXaxis()->SetTitleSize(0.045);
h->GetYaxis()->SetTitleSize(0.045);
h->GetYaxis()->SetTitleOffset(1.7);
h->Draw();
h2->Draw("SAME");
//-----------------------------------------------------
z=new TLine(ymin,ymin,ymax,ymax);
z->Draw();
z=new TLine(ymin,1.,ymax,1.);
z->Draw();
z=new TLine(1.,ymin,1.,ymax);
z->Draw();
//--------------------print legend etc---------------------
cid = new char[35];cid[34]=0;
  if(Bec==1) sprintf(cid,"Barrel B%d, chip=%d",Layer, ichip);
  else {
  if(Layer==6) sprintf(cid,"Endcap outer, chip=%d",ichip);
  if(Layer==7) sprintf(cid,"Endcap middle, chip=%d",ichip);
  if(Layer==8) sprintf(cid,"Endcap inner, chip=%d",ichip);
  }
//
leg=new TLegend(0.25,0.78,0.65,0.93,cid);
leg->SetFillColor(10);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,"eta=-6,-4,-2,1,3,5","p");
leg->AddEntry(h2,"eta=-5,-3,-1,2,4,6","p");
}
leg->Draw();
//-----------------------------------------------------
//
cid = new char[39];
sprintf(cid,"noise_50R150_%d_%d_%d.png",Bec,Layer,ichip);
cid[38]=0;
c1->Print(cid);
}
}