#include "Background.hh"
#include "TRandom.h"
#include "TMath.h"
#include "constants.hh"
#include <cmath>
  

double Background::dP2_dsindec_dlogE( double sindec, double logE)
{
  // some toy energy distribution
  double dP_dlogE   = TMath::Gaus(logE, 2, 1, true); 
  double dP_dsindec = 0.5; // uniform betwen -1 and 1
  return dP_dlogE * dP_dsindec;
}

double Background::dP2_dOmega_dlogE( double sindec, double logE)
{
  return 2*pi * dP2_dsindec_dlogE( sindec, logE );
}

SearchEvent Background::generate()
{
  double logE   = gRandom->Gaus(2,1);
  double sindec = gRandom->Uniform(-1,1);

  return SearchEvent( gRandom->Uniform(-pi,pi),
		      asin( sindec ),
		      logE );
  
}