//




{



gROOT->Reset();


  char line[5010];
  //read in cells and LAs
  short LAcell[600];
  FILE* LAfile;
  LAfile = fopen ("jc2_cells.csv","r");

  while ( !feof(LAfile) )
  {
    if ( fgets(&line[0], 5000, LAfile) == NULL)
    {
      cout<<"breaking from LA file"<<endl;
      break;
    }

    if (line[0] == NULL)
    {
      cout<<"breaking"<<endl;
      break;
    }

    int cell  = atoi( strtok( line, ";" ) );
    int la  = atoi( strtok( NULL, ";" ) );

    LAcell[cell] = la;

  }

  for (int i = 0; i<600; i++)
  {
    if (LAcell[i]!=0)
    {
      //cout<<i<<"   "<<LAcell[i]<<endl;
    }
  }
  //read in cells and LAs


  //read in cell positions
  char lat[20];
  char lon[20];

  float CELLlatlon[600][2];
  FILE* LATfile;
  LATfile = fopen ("jc2_lat2.csv","r");

  for (int i=0; i<600; i++) 
  {
    CELLlatlon[i][0]=0;
    CELLlatlon[i][1]=0;
  }

  while ( !feof(LATfile) )
  {
    if ( fgets(&line[0], 5000, LATfile) == NULL)
    {
      cout<<"breaking from LAT file"<<endl;
      break;
    }

    if (line[0] == NULL)
    {
      cout<<"breaking"<<endl;
      break;
    }

    int cell  = atoi( strtok( line, ";" ) );

    strcpy( lat, strtok( NULL, ";" ) );
    strcpy( lon, strtok( NULL, ";" ) );

    int lat_d = atoi( strtok( lat, " " ) );
    int lat_m = atoi( strtok( NULL, " " ) );
    int lat_s = atoi( strtok( NULL, " " ) );

    int lon_d = atoi( strtok( lon, " " ) );
    int lon_m = atoi( strtok( NULL, " " ) );
    int lon_s = atoi( strtok( NULL, " " ) );

    //cout<<cell<<"   "<<lat<<"    "<<lon<<"   "<<lat_d<<"   "<<lat_m<<"   "<<lat_s<<endl;

    CELLlatlon[cell][0] =  (lat_d+ lat_m/60. + lat_s/3600.);
    CELLlatlon[cell][1] = -(lon_d+ lon_m/60. + lon_s/3600.);



  }

  //read in cells positions




  TFile f("histo5.root","READ");
  //TFile f("histo5_all.root","READ");
//TFile f("histo5_all1.root","READ");

 c1 = new TCanvas("c1","sub data",200,10,700,500);

 char c[10]; int n = 1;
 char msc[50]; int n = 1;

 TLatex l; l.SetTextSize(0.03); l.SetTextAlign(22);

gStyle->SetOptStat(111110);

 gStyle->SetStatH(.25);
 gStyle->SetStatW(.25);

//////////////////////////////////////////////////////////////////////
nh_calltime->Draw();
TLine l1; l1->SetLineColor(11);
 TLine l2; l2->SetLineColor(13);
 for (int i = 24; i<720; i+=24) {l1->DrawLine(i, 0, i, 16300);}
 for (int i = 72; i<720; i+=168) {l2->DrawLine(i, 0, i, 16300);}

c1->Print("histo5.ps(","ps");  c1.Clear();


 gStyle->SetOptStat(kFALSE);


// TH2F* huserspercell = new TH2F("userspercell","users per cell", 301, -.5, 300.5, 20, 0., 24.);

 TH1F* husers5 = new TH1F("users5", "users per cell at 5pm", 14, 0., 7000.);

 TH1F* hcell[601];
 for (int i = 0; i<601; i++)
 {
   char name[20];
   char name1[20];
   char title[20];

   sprintf(name, "%s%d%s", "cell",i, NULL);
   sprintf(title, "%s%d%s", "cell ",i,NULL);

   hcell[i] = new TH1F(name , title, 20, 0., 24.);

 }



 for (int i = 0; i<301; i++)
 {

   for (int j = 0; j<20; j++)
   {
     hcell[i]->Fill(.5+(j*72)/60., 5.4 * userspercell->GetBinContent(i,j) );
   }



 }

 int pcnt = 1;
 c1.Divide(4,4);
 for (int i = 0; i<601; i++)
 {
   if (pcnt > 16)
   {
     c1->Print("histo5.ps(","ps");  c1.Clear();
     c1.Divide(4,4);
     pcnt = 1;
   }

   int total=0;
   for (int j = 0; j<20; j++){ total+=hcell[i]->GetBinContent(j);}
   
   if ( total !=0 )
   {
     husers5->Fill(hcell[i]->GetBinContent(15));

     c1.cd(pcnt);
     hcell[i]->Draw();
     pcnt++;    
   }
   
 }
c1->Print("histo5.ps(","ps");  c1.Clear();


gStyle->SetOptStat(111110);
 husers5->Draw();c1->Print("histo5.ps(","ps");  c1.Clear();



c1->Print("histo5.ps)","ps");  c1.Clear();



}