Skip to content

Commit c2a97f3

Browse files
committed
Give the QED interaction spec the size of the QED event pool
This fixes a problem in the QED specification handed to the collision context tool. - The MC number string was 10000000:NEventsQED, so the round robin wrapped the QED event IDs at 10000000 instead of at the number of events the workflow simulates, and the MC labels named QED events which do not exist. - Both numbers are NEventsQED now, which also repairs the labels on the O2 versions already deployed. https://its.cern.ch/jira/browse/O2-7132
1 parent dfaedea commit c2a97f3

1 file changed

Lines changed: 5 additions & 1 deletion

File tree

MC/bin/o2dpg_sim_workflow.py

Lines changed: 5 additions & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -659,7 +659,11 @@ def getDPL_global_options(bigshm=False, ccdbbackend=True, runcommand=True):
659659
includeQED = False
660660
else:
661661
qedrate = INTRATE * QEDXSecExpected[COLTYPE] / XSecSys[COLTYPE] # hadronic interaction rate * cross_section_ratio
662-
qedspec = 'qed' + ',' + str(qedrate) + ',10000000:' + str(NEventsQED)
662+
# the QED events are reused in a round robin, so both numbers are the size of the QED pool.
663+
# Older O2 versions wrap the QED event IDs at the first number instead of the second one and
664+
# would otherwise put event IDs into the MC labels which are not in the QED kinematics.
665+
# See https://its.cern.ch/jira/browse/O2-7132
666+
qedspec = 'qed' + ',' + str(qedrate) + ',' + str(NEventsQED) + ':' + str(NEventsQED)
663667

664668
PreCollContextTask['cmd'] = task_finalizer([
665669
'${O2_ROOT}/bin/o2-steer-colcontexttool',

0 commit comments

Comments
 (0)