#include <iostream>#include <stdlib.h>#include <sstream>#include <math.h>#include <map>#include <stdexcept>#include <gsl/gsl_errno.h>#include <gsl/gsl_math.h>#include <gsl/gsl_roots.h>#include <gsl/gsl_sf.h>#include <TF1.h>#include <Math/WrappedTF1.h>#include <Math/BrentRootFinder.h>#include "QCD.h"