#include <math.h>
#include <iostream.h>
#include <fstream.h>

void panic (char* str); 

main (int argc, char* argv[])
{
   ifstream info (argv[1]); 
   ifstream in   (argv[2]); 
   ofstream out  (argv[3]); 


   double v1, v2;  
   int nvar = atoi(argv[4]); 
   int idt; 
   int nnoise, nsignal; 

   

   for (double cut = 0.01 ; cut < 4; cut += 0.005) 
   { 
      nnoise = nsignal = 0; 
      do { 
         in >> v1 >> v2 >> idt;  
         if ( 
((4500*exp(-v1*v1)+800*exp(-v1/1.5))*550*v2/v1*exp(-v2*v2*64/(v1*4.4*v1*4.4)))
/((500*exp(-v1*v1)+700*exp(-v1/3))*150*v2/v1*exp(-v2*v2*64/(v1*5*v1*5)))
             < cut         ) 
         { 
             if (idt == 0) nnoise++; 
             else          nsignal++;   
         } 
      }  while (!in.eof()); 
         out << nsignal/44000.0 << " " << (double)nsignal/(nsignal+nnoise)
          << endl; 
         in.close (); 
         in.open (argv[1]); 
   } 
   in.close(); 
   out.close (); 
}


void panic (char* str)
{
   cout << str << endl; 
   exit (0); 
}