#include "Background.hh" #include "TRandom.h" #include "TMath.h" #include "constants.hh" #include 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 ); }