#include <stdio.h>
#include <math.h>
//#include <fstream.h>
#include <string.h>
#include <stdlib.h>

int main(int argc, char *argv[])

{

  long double ra1[750],dec1[750],hx[750],errhx[750],ra2[750],dec2[750],b2[750],r2[750],i2[750];
  long double dx1[750],dx2[750],dy1[750],dy2[750],s1[750],s2[750],rate1[750],errate1[750],errpos1[750],rate2[750],errate2[750],errpos2[750];
  long double cts1[750],errcts1[750],ext1[750],ml1[750],bkg1[750],exp1[750],ontime1[750];
  long double cts2[750],errcts2[750],ext2[750],ml2[750],bkg2[750],exp2[750],ontime2[750];
  long double dist,distref,err_sist,numero_sigma;
  int i,j,k,nx,no;
  char filex[100], fileo[100];                                                  
  FILE *p;

  //puts("\n file src 1 ");
  //scanf("%s",&filex);
  //puts("\n file src 2");
  //scanf("%s",&fileo);
  //puts("\n n src file 1");
  //scanf("%d",&nx);
  //puts("\n n src file 2");
  //scanf("%d",&no);

   if(argc != 7)
   {
      printf("Incorrect number of args\n");
      exit(-1);
   }
   if(sscanf(argv[1],"%s", &filex)!=1)
   {
      printf("input file 1 scanf failed.\n");
      exit(-2);
   }
   if (sscanf(argv[2],"%s", &fileo) != 1)
   {
      printf("input file 2 scanf failed.\n");
      exit(-3);
   }
   if(sscanf(argv[3],"%d", &nx)!=1)
   {
      printf("n src file1 scan failed.\n");
      exit(-2);
   }
   if (sscanf(argv[4],"%d", &no) != 1)
   {
      printf("n src file1 scan failed..\n");
      exit(-3);
   }
   if(sscanf(argv[5],"%Lf", &err_sist) !=1)
   {
      printf("n src file1 scan failed.\n");
      exit(-2);
   }
   if (sscanf(argv[6],"%Lf", &numero_sigma) != 1)
   {
      printf("n src file1 scan failed..\n");
      exit(-3);
   }
   //   if (sscanf(argv[5],"%Lf", &distref) != 1)
   //{
   // printf("dist scan failed..\n");
   // exit(-4);
   //}


	nx=nx-2; // eliminazione prime righe vuote nei file in lettura

	//printf("\nreading first file....");	
  p=fopen(filex,"r"); 

	for(i=0;i<nx;i++)
	{
	ra1[i]=0.0;
	dec1[i]=0.0;
	rate1[i]=0.0;
	errate1[i]=0.0;
	errpos1[i]=0.0;
	fscanf(p,"%Lf %Lf  %Lf %Lf %Lf  %Lf %Lf  %Lf %Lf %Lf %Lf %Lf",&rate1[i],&errate1[i],&ra1[i],&dec1[i],&errpos1[i],&cts1[i],&errcts1[i],&ext1[i],&ml1[i],&bkg1[i],&exp1[i],&ontime1[i]);
	}
	
  fclose(p);

	//printf("done");
  
	//printf("\nreading second file....");
  p=fopen(fileo,"r"); 

	for(i=0;i<no;i++)
	{
	rate2[i]=0.0;
	errate2[i]=0.0;
	errpos2[i]=0.0;
	ra2[i]=0.0;
	dec2[i]=0.0;
	fscanf(p,"%Lf %Lf  %Lf %Lf %Lf  %Lf %Lf %Lf  %Lf %Lf %Lf  %Lf",&rate2[i],&errate2[i],&ra2[i],&dec2[i],&errpos2[i],&cts2[i],&errcts2[i],&ext2[i],&ml2[i],&bkg2[i],&exp2[i],&ontime2[i]);
	}
	
  fclose(p);

	//printf("done");


  k=0;

  p=fopen("out.txt","w"); 

	
  	for(i=0;i<nx;i++)
	  {
		for(j=0;j<no;j++)
		  {
			dist=0;
			distref=0.0;
			// 2= err_sist
			// 6= numero_sigma
        // vecchia distanza    
		//	dist=sqrt(((ra1[i]-ra2[j])*(ra1[i]-ra2[j]))+((dec1[i]-dec2[j])*(dec1[i]-dec2[j])));
        // nuova distanza 1 app per piccoli angoli    
//			dist=((sqrt((((ra1[i]-ra2[j])*(ra1[i]-ra2[j]))*cos(dec1[i]/(57.29577951308232))*cos(dec1[i]/(57.29577951308232)))+((dec1[i]-dec2[j])*(dec1[i]-dec2[j])))));
        // nuova distanza 2    
		    dist=(atan2(sqrt(1-(cos((dec2[j])/57.29577951308232)*cos((dec1[i])/57.29577951308232)*cos(((ra2[j])/57.29577951308232)-((ra1[i])/57.29577951308232))+sin((dec2[j])/57.29577951308232)*sin((dec1[i])/57.29577951308232))*(cos((dec2[j])/57.29577951308232)*cos((dec1[i])/57.29577951308232)*cos(((ra2[j])/57.29577951308232)-((ra1[i])/57.29577951308232))+sin((dec2[j])/57.29577951308232)*sin((dec1[i])/57.29577951308232))),cos((dec2[j])/57.29577951308232)*cos((dec1[i])/57.29577951308232)*cos((ra2[j])/57.29577951308232-(ra1[i])/57.29577951308232)+sin((dec2[j])/57.29577951308232)*sin((dec1[i])/57.29577951308232))*57.29577951308232);
            distref=(err_sist+(numero_sigma*sqrt((errpos1[i]*errpos1[i])+(errpos2[j]*errpos2[j]))))/3600;
            //distref=(err_sist+(numero_sigma*errpos1[i]))/3600;
			//printf("\nerr_sist: %Lf, numero_sigma: %Lf, dist: %Lf, distref: %Lf",err_sist, numero_sigma, dist, distref);
			if(dist<distref)
			  {
			    k=k+1;
			    fprintf(p,"%Lf %Lf %Lf %Lf %Lf %Lf %Lf %Lf %Lf %Lf %Lf  %Lf %Lf  %Lf  %Lf %Lf  %Lf  %Lf %Lf  %Lf %Lf %Lf %Lf %Lf  %Lf  \n",ra1[i],dec1[i],errpos1[i],ra2[j],dec2[j],errpos2[j],dist,rate1[i],errate1[i],rate2[j],errate2[j],cts1[i],errcts1[i],cts2[j],errcts2[j],ext1[i],ext2[j],ml1[i],ml2[j],bkg1[i],bkg2[j],exp1[i],exp2[j],ontime1[i],ontime2[j]);
			  }
		  }
		if(k<1)
		  {
		    fprintf(p,"sorgente assente su src list complessiva %Lf %Lf %Lf %Lf  %Lf %Lf %Lf %Lf %Lf %Lf %Lf  %Lf \n",ra1[i],dec1[i],errpos1[i],rate1[i],errate1[i],cts1[i],errcts1[i],ext1[i],ml1[i],bkg1[i],exp1[i],ontime1[i]);
		  }
		k=0;		  
	  }
  fclose(p);

//printf("\nnx: %d, no: %d\n", nx, no);
	   
printf("\n"); 
}