remove unwanted comments

higgsfodptN.cpp

@@ -145,18 +145,18 @@ int main(int argc, char* argv[]) {
         if (sectype == "hadronic") {
             results = higgshard.higgsdpt(pt,y1,y2);
         } else if (sectype == "partonic") {
-            results = higgspart.partonichiggsdpt(pt,std::pow(10,nn));
+            results = higgspart.partonichiggsdpt(pt,nn);
         } else {
             std::cout << "Error in entry sectype!" << std::endl;
             exit(EXIT_FAILURE);
         }
 
         // Generate some output logs & write to output file
-        printf("N=%e: dHdpt = %e +- %e. \n", std::pow(10,nn), results[0], results[1]);
+        printf("N=%e: dHdpt = %e +- %e. \n", nn, results[0], results[1]);
         output_file.setf(std::ios_base::scientific);
-        output_file << std::pow(10,nn) << std::setw(space)
-                    << results[0] << std::setw(space) << "\n";
-                    // << results[1] << "\n";
+        output_file << nn << std::setw(space)
+                    << results[0] << std::setw(space)
+                    << results[1] << "\n";
         output_file.flush();
 
         nn += Nbin;

src/higgsptpartonic.cpp

@@ -711,8 +711,6 @@ double HiggsDpTpartonic::deltapartonic(double pt, double nn, double zz) {
   // toghether with setting QQ2=0 is needed to properly handle the delta(Q²)
   // term.
   double jac = 1. / sqrt(std::pow(sh + MH2 - QQ2, 2) - 4. * sh * mt2);
-  // std::cout << gg0(sh, th, uh, MH2) *jac /sh << std::endl;
-
 
   // compute terms proportional to delta(Q^2)
   if (ORD >= 0) {