{

  double gallong;
  double gallat;
  double x40d, y40d, x70d, y70d, x100d, y100d;
  
 
 
 vector<double> x40;
 vector<double> y40;
 vector<double> x70;
 vector<double> y70;
 vector<double> x100;
 vector<double> y100;

 vector<double> gallongv;
 vector<double> gallatv;
 double temp1, temp2, temp3,temp4,temp5;
  
 ifstream in;
 in.open("/home/macolino/Work_stuff/work0/anisotropy/autocorrelation/eventi/evs.nov08.sort.dat");

 //while((in.good())
 while(1)
{
   if (! in.good()) break;
  
  in>>temp1>>temp2>>temp3>>temp4>>temp5;
  in>>gallong>>gallat;
  
  //cout<<temp1<<"  "<<temp2<<"  "<<temp3<<"  "<<temp4<<"  "<<temp5<<"  "<<gallong<<"  "<<gallat<<endl;
  //cout<<endl;
  
  
  gallongv.push_back(gallong);
  gallatv.push_back(gallat);
  
  if(temp3>40.0) {
  x40.push_back(gallong);
  y40.push_back(gallat);
  }
  if(temp3>70.0) 
    {
      x70.push_back(gallong);
      y70.push_back(gallat);
    }
  if(temp3>100.0){ 
        x100.push_back(gallong);
	y100.push_back(gallat);
  }
  


 
}

in.close();

int n=gallongv.size();
cout<<n<<endl;

int n4=x40.size();
int n7=x70.size();
int n10=x100.size();

cout<< n4 <<" "<<n7<<" "<<n10<<endl;

double x[n];
double y[n];

double x4x[n4];
double y4y[n4];

double x7x[n7];
double y7y[n7];

double x10x[n10];
double y10y[n10];


for(int i=0;i<gallongv.size();i++){
  
  x[i]=-gallongv[i];
  y[i]=gallatv[i];
}

for(int i=0;i<n4;i++){
  cout<<i<<endl;
  x4x[i]=-x40[i];
  y4y[i]=y40[i];

}

for(int i=0;i<n7;i++){
  cout<<i<<endl;
  x7x[i]=-x70[i];
  y7y[i]=y70[i];
}

for(int i=0;i<n10;i++){
  x10x[i]=-x100[i];
  y10y[i]=y100[i];

}
  //  cout<<"gallatv["<<i<<"]="<<gallatv[i]<<endl; 

TGraph *plotevents=new TGraph(n,x,y);
plotevents->SetMarkerColor(2);
plotevents->SetMarkerStyle(6);
plotevents->GetYaxis()->SetTitle("Galactic Latitude [#circ]");
plotevents->GetXaxis()->SetTitle("Galactic Longitude [#circ]");
plotevents->GetYaxis()->SetRangeUser(-90.,90.);
plotevents->Draw("AP");

// Remove the current axis

plotevents->GetXaxis()->SetLabelOffset(999);
plotevents->GetXaxis()->SetTickLength(0);

// Redraw the new axis
gPad->Update();
/*TGaxis *newaxis = new TGaxis(gPad->GetUxmax(),
			     gPad->GetUymin(),
			     gPad->GetUxmin(),     
			     gPad->GetUymin(),
			     plotevents->GetXaxis()->GetXmin(),
			     plotevents->GetXaxis()->GetXmax(),
			     510,"-SDH");*/

TGaxis *newaxis = new TGaxis(gPad->GetUxmax(),
			     gPad->GetUymin(),
			     gPad->GetUxmin(),     
			     gPad->GetUymin(),
			     -180,
			     180,
			     510,"-SDH");
 
newaxis->SetLabelOffset(-0.03);   

newaxis->Draw();


TGraph *plotevents2=new TGraph(n4, x4x, y4y);
plotevents2->SetMarkerColor(4);
plotevents2->SetMarkerStyle(8);
plotevents2->Draw("P");

TGraph *plotevents3=new TGraph(n7, x7x, y7y);
plotevents3->SetMarkerColor(3);
plotevents3->SetMarkerStyle(8);
plotevents3->Draw("P");

TGraph *plotevents3=new TGraph(n10, x10x, y10y);
plotevents3->SetMarkerColor(6);
plotevents3->SetMarkerStyle(8);
plotevents3->Draw("P");

}