herwig is hosted by Hepforge, IPPP Durham
close Warning: Can't synchronize with repository "(default)" (Repository path '/hepforge/hg/herwig/public/herwig' does not exist.). Look in the Trac log for more information.

Opened 7 years ago

Closed 7 years ago

#7 closed defect (fixed)

Eta asymmetry in W production, H7

Reported by: Peter Richardson Owned by: Peter Richardson
Priority: critical Milestone: Herwig 7.0.3
Version: 7.0.2 Keywords:
Cc: Miroslav.Myska@…

Description

Dear Herwig experts,
I would like to ask you for a help with understanding of an output of herwig, when generating
very simple single W production with muon decay. (in file attached)

I plot muon pseudorapidity (after pT cut 5 GeV, |eta| < 5.0 ) for two sets of muons:
1) Matrix Element (ME) outgoing muons:
(using: event->primarySubProcess()->outgoing(); )
2) Final state muons - muons at the end of the loop over children of the previous muon - having no further children.

I observe asymmetry, where muons slightly tend to go into the negative side, mean eta = -0.072
It is visible by eye in the plot - red line. Statistics is 10x1M events, all samples have similar asym.
Strange thing is that muons coming directly from ME show much smaller asymmetry, mean =-0.0025 (greem line)

Please, do you see an obvious reason (or a bug/missing thing in the in file) for that behavior?

I tried to turn off the QED radiation by:
set /Herwig/Shower?/Evolver:Interactions QCDOnly
but the result did not really changed (blue and black lines)

I run H7.0.2 now, but I remember seeing this also in older H++ versions.

Thank you!

All the best,
Miroslav

Attachments (2)

LHC_13_Wp.in (849 bytes) - added by Peter Richardson 7 years ago.
input file
Muon_pseudorapidity_asmmetry.gif (8.0 KB) - added by Peter Richardson 7 years ago.
plot

Download all attachments as: .zip

Change History (3)

Changed 7 years ago by Peter Richardson

Attachment: LHC_13_Wp.in added

input file

Changed 7 years ago by Peter Richardson

plot

comment:1 Changed 7 years ago by Peter Richardson

Cc: Miroslav.Myska@… added
Resolution: fixed
Status: newclosed

This is a problem with the handling of events with the quark in the negative z direction in the hard matrix element correction, fixed in rfe6900a23531

Note: See TracTickets for help on using tickets.