Index: /branches/pap_branch_080617/ppStack/src/ppStackMatch.c
===================================================================
--- /branches/pap_branch_080617/ppStack/src/ppStackMatch.c	(revision 18290)
+++ /branches/pap_branch_080617/ppStack/src/ppStackMatch.c	(revision 18291)
@@ -37,4 +37,5 @@
     int renormWidth = psMetadataLookupS32(&mdok, config->arguments,
                                           "RENORM.WIDTH"); // Width for renormalisation box
+    float penalty = psMetadataLookupF32(NULL, config->arguments, "PENALTY"); // Penalty for wideness
 
     if (!pmReadoutMaskNonfinite(readout, maskBad)) {
@@ -61,5 +62,4 @@
         int ringsOrder = psMetadataLookupS32(NULL, recipe, "RINGS.ORDER"); // RINGS polynomial order
         int binning = psMetadataLookupS32(NULL, recipe, "SPAM.BINNING"); // Binning for SPAM kernel
-        float penalty = psMetadataLookupF32(NULL, config->arguments, "PENALTY"); // Penalty for wideness
         psMaskType maskBlank = pmConfigMask(psMetadataLookupStr(NULL, recipe, "MASK.BLANK"),
                                             config); // Mask for blank reg.
@@ -140,5 +140,6 @@
                          PS_DATA_REGION | PS_META_DUPLICATE_OK, "Fake subtraction region", region);
         psFree(region);
-        pmSubtractionKernels *kernels = pmSubtractionKernelsPOIS(FAKE_SIZE, 0, PM_SUBTRACTION_MODE_1);
+        pmSubtractionKernels *kernels = pmSubtractionKernelsPOIS(FAKE_SIZE, 0, penalty,
+                                                                 PM_SUBTRACTION_MODE_1);
         // Set solution to delta function
         kernels->solution1 = psVectorAlloc(kernels->num + 2, PS_TYPE_F64);