Very slow when using custom differential flux function

[xpzhang]

Hello, When I use GenerateFromCustomJ() with a user-defined differential flux function, same with the example code in the paper:

double PrimaryGeneratorAction::J_Gaisser(double p, double theta)
{
  double A = 0.14*pow(p, -2.7);
  double B = 1. / (1. + 1.1*p*cos(theta)/115.);
  double C = 0.054 / (1. + 1.1*p*cos(theta)/850.);
  return A*(B+C);
}

The generation gets very slow (several tenth slower compared to the pre-defined J. What can I do to improve the performance? Thanks.

[dr4kan]

Version 2 will improve this. It will be released in few weeks.

[dr4kan]

I checked and the code is not slow. I suspect that you did not set a reasonable range for the momuntum. Please note that, in the Gaisser parametrization if p is small J diverges...