Scale in LHE file and scale definition for NLO

Hi all,

when looking at the scale value given in the LHE files for each event, I was wondering how that value is related to the dynamical scale I added to setscales.f (the process is at NLO in case that is of importance). There I added a custom scale definition that I run via setting dynamical_scale_choice to the corresponding value in the run card. In the LHE header I can see the correct string stored in temp_scale_id that I defined for my scale appear, so I take this as confirmation that my scale definition was actually used.

Now when looking at the particle momenta printed in the LHE file and manually calculating the corresponding scale value I get a different result than what is given as the scale in the LHE file (4th entry of the first line for each event, according to the LHE spec). The scale is calculated using the particle numbers associated to the different particles as given in the Feynman diagrams that can be found in the various SubProcesses/P* folders (alternatively: index.html -> Process Information). I looked at a number of them (Born, Real, Virtual, different initial states) and found the relevant particles used in the scale definition to always have the same ID. So what exactly would the relation of the hard process scales set via dynamical_scale_choice and the scale value printed in the LHE file be?

The background for this is that I would really like to obtain numerical values of the scale that is actually used in the hard process. I tried to put "write(*,*)'scale=',tmp" statements in the part of setscales.f where my scale definition is included - but get no output unfortunately.

Many thanks,
Daniel