=== modified file 'Template/LO/SubProcesses/cuts.f' --- Template/LO/SubProcesses/cuts.f 2016-08-30 23:56:34 +0000 +++ Template/LO/SubProcesses/cuts.f 2017-11-09 15:48:50 +0000 @@ -1238,7 +1238,7 @@ c Note the following condition is the first line of setclscales c if(xqcut.gt.0d0.or.ickkw.gt.0.or.scale.eq.0.or.q2fact(1).eq.0)then c Do not duplicate it since some variable are set for syscalc in the fct - if(.not.setclscales(p))then + if(.not.setclscales(p,.false.))then cutsdone=.false. cutspassed=.false. passcuts = .false. === modified file 'Template/LO/SubProcesses/reweight.f' --- Template/LO/SubProcesses/reweight.f 2016-08-29 10:00:16 +0000 +++ Template/LO/SubProcesses/reweight.f 2017-11-09 15:48:24 +0000 @@ -465,13 +465,15 @@ return end - logical function setclscales(p) + logical function setclscales(p, keepq2bck) c************************************************** c Calculate dynamic scales based on clustering c Also perform xqcut and xmtc cuts +c keepq2bck allow to not reset the parameter q2bck c************************************************** implicit none + logical keepq2bck include 'message.inc' include 'genps.inc' include 'maxconfigs.inc' @@ -952,8 +954,10 @@ if(.not. fixed_fac_scale) then q2fact(1)=scalefact**2*q2fact(1) q2fact(2)=scalefact**2*q2fact(2) - q2bck(1)=q2fact(1) - q2bck(2)=q2fact(2) + if (.not.keepq2bck)then + q2bck(1)=q2fact(1) + q2bck(2)=q2fact(2) + endif if (btest(mlevel,3)) $ write(*,*) 'Set central fact scales to ',sqrt(q2bck(1)),sqrt(q2bck(2)) endif @@ -1259,7 +1263,7 @@ endif - if(.not.setclscales(p)) then ! assign the correct id information. + if(.not.setclscales(p,.true.)) then ! assign the correct id information.(preserve q2bck) write(*,*) "Fail to cluster the events from the rewgt function" stop 1 c rewgt = 0d0