#include #include #include #include "APFEL/APFEL.h" using namespace std; int main() { APFEL::SetPerturbativeOrder(1); APFEL::SetNumberOfGrids(3); APFEL::SetGridParameters(1,130,3,1e-9); APFEL::SetGridParameters(2,60,5,1e-1); APFEL::SetGridParameters(3,20,5,8e-1); // Initializes integrals on the grids APFEL::InitializeAPFEL(); double S = pow(8e3,2); double Q0 = sqrt(2e0); double Qmin = 10; double Qmax = 6e3; int NMX = 30; cout << "Luminosities at sqrt(S) = " << sqrt(S) << " GeV:\n" << endl; cout << "- gg luminosity " << endl; for (int i=1; i < 30; i++) { double Q = Qmin * pow(Qmax/Qmin,(double)(i-1)/(double)(NMX-1)); APFEL::EvolveAPFEL(Q0,Q); cout << "MX = " << Q << " LUMI = " << APFEL::LUMI(0,0,S) << endl; } cout << endl; return 0; }