From e0138cde79be1e580e35d5b4bfd6e09f93adfe82 Mon Sep 17 00:00:00 2001
From: Sylvester Joosten <sylvester.joosten@gmail.com>
Date: Wed, 6 Apr 2016 14:46:52 -0400
Subject: [PATCH] re-added PYTHIA CKIN cuts to loadtune

---
 pythia6m/gmc/gmc_loadtune.F | 30 ++++++++++++++++++++++++++++++
 1 file changed, 30 insertions(+)

diff --git a/pythia6m/gmc/gmc_loadtune.F b/pythia6m/gmc/gmc_loadtune.F
index f4bb4b2..ebfc5c9 100644
--- a/pythia6m/gmc/gmc_loadtune.F
+++ b/pythia6m/gmc/gmc_loadtune.F
@@ -255,4 +255,34 @@ C minimum number of flavours that may be assumed in alpha_s expression
 C maximum number of flavours that may be assumed in alpha_s expression
       MSTU(114) = 4
 C
+C******************************************************** kinematical cuts
+C
+C range of m^ = sqrt(s^)
+C (D: 2., -1.)
+      CKIN(1) = 1.
+      CKIN(2) =-1.
+C range of p^_t (DIS = Q^2)
+C (*) (D: 0., -1.) (GMC: 0, -1)
+      CKIN(3) = 0.
+C CKIN    3=1.
+      CKIN(4) = -1.
+C p^_t,min for singular processes in limit p^_t -> 0
+      CKIN(5) = 1.00
+C m_0 hard 2->2 process is classified singular for p^_t -> 0 if the mass
+C of at least one of the two outgoing partons is below m_0
+      CKIN(6) = 1.00
+C range for W^2 in DIS processes (W^2 = Q^2(1-x)/x, neglecting M^2 and ISR!)
+C (*) (D: 4.,-1) before (10. -1)
+      CKIN(39) = 4.
+      CKIN(40) = -1.
+C range for Q^2 
+      CKIN(65) = 1e-09
+      CKIN(66) = 100.
+C CKIN   65=0.1
+C CKIN   66=0.01
+C allowed range for W, i.e. either the photon-hadron or photon-photon invariant mass
+      CKIN(77) = 2.0
+      CKIN(78) = -1.
+C
+STOP
       end
-- 
GitLab