Toy code to calculate production cross sections in pp hard scattering collisions to leading order at tree level using Monte Carlo integration over the ELHQ parton distribution functions.

make
./ppmc_test
./ppmc

Let's do our own calculation, in a new directory called src/my_calc. Create a file in src/my_calc/ with a main function. Include any other files you need as normal, and place them in src/my_calc/. Using mc::Integrate2D and the parton distribution functions in elhq:: a variety of tree-level hard scattering diagrams can be computed. See src/ppmc.cpp for an example; calculating the Drell-Yan pp -> μ+μ- cross section for multiple low-energy cutoffs.

Add your file(s) (lets call it src/my_calc/my_calculation.cpp, which #includes src/my_calc/support.hpp, which is implemented in src/my_calc/support.cpp) as a target to the Makefile like so:

MYCALC_DIR         = ./my_calc
MYCALC_FILES       = $(wildcard $(MYCALC_DIR)/*.cpp)
MYCALC_OBJ_FILES   = $(patsubst $(SRC_DIR)/%.cpp,$(OBJ_DIR)/my_calc/%.o,$(MYCALC_FILES))

my_calculation: $(OBJ_FILES_NOPP) $(MYCALC_FILES)
  $(CC) $(LDFLAGS) -o $@ $^ 

And build and run:

make my_calculation
./my_calculation