Small adjustment to make collinear evolution with APFEL++ correctly work also...

Small adjustment to make collinear evolution with APFEL++ correctly work also with an ActiveFlavorsThresholdsConstant object

src/partons/modules/evolution/collinear_distribution/CollinearDistributionEvolutionApfel.cpp

@@ -89,8 +89,13 @@ void CollinearDistributionEvolutionApfel::prepareSubModules(const std::map<std::
 
     // Get thresholds. Set to zero whatever is below one.
     std::vector<double> thresholds;
-    for (ActiveFlavorsThresholds aft : m_pActiveFlavorsModule->getNfIntervals())
-      thresholds.push_back(aft.getLowerBound() < 1 ? 0 : sqrt(aft.getLowerBound()));
+    std::vector<ActiveFlavorsThresholds> afts = m_pActiveFlavorsModule->getNfIntervals();
+    if (afts.size() == 1)
+      for (int i = 0; i < afts[0].getNf(); i++)
+	thresholds.push_back(0);
+    else
+      for (ActiveFlavorsThresholds aft : afts)
+	thresholds.push_back(aft.getLowerBound() < 1 ? 0 : sqrt(aft.getLowerBound()));
 
     // Initialize QCD evolution objects
     m_dglapobj = InitializeDglapObjectsQCD(*m_g, thresholds);

src/partons/modules/evolution/collinear_distribution/CollinearDistributionEvolutionModule.cpp

@@ -16,8 +16,8 @@ const std::string CollinearDistributionEvolutionModule::COLLINEAR_DISTRIBUTION_E
 CollinearDistributionEvolutionModule::CollinearDistributionEvolutionModule(const std::string &className) :
         ModuleObject(className, ChannelType::UNDEFINED),
 	m_x(0), m_MuF2(0), m_MuR2(0),
-	m_pertOrder(PerturbativeQCDOrderType::UNDEFINED),
-	m_type(CollinearDistributionType::UNDEFINED),
+	m_pertOrder(PerturbativeQCDOrderType::NNLO),
+	m_type(CollinearDistributionType::UnpolPDF),
 	m_pRunningAlphaStrongModule(0),
 	m_pActiveFlavorsModule(0) {
 }