set scale

Hello Marcel,

First check which are the positions (numbers) of the particles you want. It sounds like you only have two final-state particles, in which case the positions are 3 and 4. You can simply add the line (at the end of subroutine set_ren_scale):

      rscale=dsqrt(SumDot(P(0,3),P(0,4),1d0))

and then uncomment the lines (in subroutine set_fac_scale):

      call set_ren_scale(P,q2fact(1))
      q2fact(1)=q2fact(1)**2
      q2fact(2)=q2fact(1) !factorization scale**2 for pdf2

All the best,
Johan

Thanks Johan Alwall, that solved my question.

I think that solved my problem...but one last question...do I have to change sth for that in my run_card or is it set by default that it'll use a dynamic scale?

Hi Marcel,

Yes by default the scale are dynamical.

Cheers,

Olivier

thanks...do u know how to convert barn/bin to dbarn/dpt for a pt distribution and the same for invariant mass?:D..I know..wrong topic now:D

Hi,
shouldn't one also uncomment the calls to set_ren_scale / set_fac_scale in cuts.f?
I cannot see these functions called anywhere else in the code.
Cheers,

Marco