#!/usr/bin/env python

import apfel

# activate some options
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();

S = 8e3**2
Q0 = 2.0**0.5
Qmin = 10
Qmax = 6e3
NMX  = 30

print "Luminosities at sqrt(S) = ", S**0.5," GeV:\n"
print "- gg luminosity "
    
for i in range(1,31):
    Q = Qmin * pow(Qmax/Qmin,(float)(i-1)/(float)(NMX-1));
    apfel.EvolveAPFEL(Q0,Q);
    print "MX = ", Q,"LUMI =",apfel.LUMI(0,0,S)

print ""