add option to weight pT hard bin cross section to allow spectra flattening

Author: gconesab

MC/config/common/pythia8/utils/mkpy8cfg.py

@@ -29,6 +29,9 @@
 parser.add_argument('--ptHatMax', type=float,
                     help='The maximum invariant pT')
 
+parser.add_argument('--weightPower', type=float,
+                    help='Weight power to pT hard spectrum')
+
 parser.add_argument('--output', default='pythia8.cfg',
                     help='Where to write the configuration')
 
@@ -103,6 +106,10 @@
     fout.write('PhaseSpace:pTHatMin = %f \n' % (args.ptHatMin))
 if args.ptHatMax is not None :
     fout.write('PhaseSpace:pTHatMax = %f \n' % (args.ptHatMax))
+if args.weightPower is not None :
+    fout.write('PhaseSpace:bias2Selection = on \n')
+    fout.write('PhaseSpace:bias2SelectionPow = %f" \n' % (args.weightPower))
+
 fout.write('\n')
 
 fout.write('### <-- generated by mkpy8cfg.py \n')