Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Example hadronizer card for HH decay #873

Merged
merged 5 commits into from May 23, 2016

Conversation

OlivierBondu
Copy link
Contributor

Based on cms-sw/cmssw#12902

NB: need cms-sw/cmssw#14555 to be integrated in CMSSW first before injection in MCM

@bendavid
Copy link
Contributor

cms-sw/cmssw#12902 is already in CMSSW_7_1_22 in fact.

@OlivierBondu
Copy link
Contributor Author

that's why the nota bene was only on the second link 😉

@bendavid
Copy link
Contributor

Yes, sorry, read too fast.

@OlivierBondu
Copy link
Contributor Author

(but since you're around, many thanks again for that feature, very much appreciated)

'ResonanceDecayFilter:eMuTauAsEquivalent = on', #on: treat electrons, muons , and taus as equivalent
'ResonanceDecayFilter:allNuAsEquivalent = on', #on: treat all three neutrino flavours as equivalent
#'ResonanceDecayFilter:mothers =', #list of mothers not specified -> count all particles in hard process+resonance decays (better to avoid specifying mothers when including leptons from the lhe in counting, since intermediate resonances are not gauranteed to appear in general
'ResonanceDecayFilter:daughters = 5,5,11,11,12,12',
Copy link
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

This is not safe in case you decay HH events produced with jet matching and extra ME jets, since this could allow b-quarks from QCD rather than the higgs decay. You should probably explicitly add
'ResonanceDecayFilter:mothers = 25' to protect against this (even if initially the fragment will be used with LO tree level HH events where this won't happen)

Copy link
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Sorry, should actually be
'ResonanceDecayFilter:mothers = 25,23,24'

since the leptons and neutrinos will actually come from W's or Z's in the event history

Copy link
Contributor Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

mmm, but in your comment you mention not specifying W's and Z's in case of leptonic decays ?

Copy link
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

This is a special case. Since the W"s and Z's would come from pythia decays, they ARE guaranteed to be included the event history.

My comment is a warning in case the hard process includes lepton production at the lhe level.

Copy link
Contributor Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

ok, I understand

actually I realized I should also have the exclusive = on if I don't want HH -> Z(bb)Z(ll)Z(ll)Z(nunu) decays possibly laying around...

Copy link
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Yes that is also true.

Copy link
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

(but you can and should turn off the Z->qqbar decay modes for this config anyways)

@OlivierBondu
Copy link
Contributor Author

@bendavid: card updated

@OlivierBondu
Copy link
Contributor Author

@bendavid : card updated

@OlivierBondu
Copy link
Contributor Author

actually re-reading your comment I'm not sure if you were talking about 23:onIfMatch or ResonanceDecayFilter:daughters, if it's the later then I'm not sure how to ensure 6 final objects for the hard process and having all lepton flavours...

@bendavid
Copy link
Contributor

I was indeed talking about 23:onIfMatch, so this looks ok now (for really maximum efficiency could do similar to switch on only the needed higgs and W decay modes, but ok probably this is pretty irrelevant, since it will only change the number of retries needed for the decays and should have a pretty small impact on cpu time)

@OlivierBondu
Copy link
Contributor Author

@bendavid card updated:

  • reshuffling the ordering of the forced decays for easier reading
  • W now forced to decay leptonically
  • open only leptonic tau decays (not quite sure of the interplay with the ResonanceDecayFilter here, feedback is welcome !)
  • onIfMatch -> onIfAny for W, Z, tau

@bendavid
Copy link
Contributor

This looks ok to me.

Note that the ResonanceDecayFilter checks the decay state after the Higgs/W/Z/Top/etc decays but BEFORE the tau decays. So there is no interplay here, and the taus are counted by the ResonanceDecayFilter if so configured, but the decay products of the tau never are.

@OlivierBondu
Copy link
Contributor Author

Thanks !
@covarell @perrozzi @govoni this is ready to be merged

@covarell covarell merged commit 921bfd1 into cms-sw:master May 23, 2016
@OlivierBondu OlivierBondu deleted the HH_decay_example branch May 23, 2016 12:17
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

Successfully merging this pull request may close these issues.

None yet

3 participants