Most of the calculations in lattice QCD has to do with finding averages over set of gauge configurations for given temperature or \(\beta=\frac{2N}{g^2}\).
So, the first process of any calculation is to generate the gauge configurations at given temperature for required dynamical system. On the basis of dynamical system we choose

1. RHMC (for dynamical fermions in background of dynamical system)
2. generateQuenched (for quenched case i.e. infinite/large mass quarks in background)

Let's say we are aiming to calculate the string tension for static quark anti quarks in (2+1) flavour QCD (i.e. 2 light quarks {u,d} and one strange quark(light) in background)

We can proceed with Following steps:

1. The first step is to generate the enough number of gauge configurations to reach equilibrium and take random number file and any one configuration after the point of thermal equilibrium (which will be required to generate any number of gauge configuration at equilibrium).
2. Generate more gauge configurations after equilibrium and measure the quantity on all of these new gauge configurations and average them and do some error analysis.

For this case we will be using RHMC where we will change the parameter file of RHMC.cpp source code or ./rhmc executable (which we build).
The changes in parameter file are crucial for the calculation and the procedure to change these parameters will be discussed in RHMC_Paramfile.