/********************* PROBAB.C *********************/
#include  <stdlib.h>
#include  <math.h>
#include  "randefs.h"
#include  "probab.h"

#define  NSIGMA  3

index SumNPicks(RAND *dat)
{ int    i, n = dat->nobjs;
  index  trial, event = 0;

  for (i = 0; i < dat->npick; i++)
  {  while (1)
     {  trial = random(n);
        if (random(101) <= (*dat->prob)(trial, n) )
        {   event += trial;
            break;
        }
     }
  }
  return(event);
}

int Dummy(index dummy1, int dummy2)
{
  return(100);
}

int LoadZero(index trial, int dummy)
{
  if (trial == 0)
     return(100);
  return(50);
}

int Normal(index trial, int n)
{ double  x;

  x = NSIGMA*( (2.0/(n-1))*trial-1.0);
  return( 100 * exp(-(x*x)/2.0) );
}
