|
[Rivet] Question concering PromptFinalState projectionDominic Hirschbuehl dhirsch at mail.cern.chWed Jan 25 16:55:42 GMT 2017
Dear all, I have a question concerning the PromptFinalState projection. I did a study on acceptences of lep+jets ttbar events, where I found an unexpected differences between Powheg+Pythia6 and Powheg+Herwig samples. The differences occured already on the level of the lepton acceptance, which includes a dilepton veto. After some investigation I have the suspicion, that leptons from hadron decays are counted as prompt leptons. Looking to the actual implementation of isPrompt I found this line of code: if (ancestor->status() != 2) continue; // no non-standard statuses or beams to be used in decision making Is there a reason to restrict the status codes to 2? If I print out all ancestors of a "prompt lepton" from a Herwig showered event I get: pdg = 2212 status = 101 0 1 pdg = 2212 status = 102 0 1 pdg = 21 status = 121 1 0 pdg = 21 status = 122 1 0 pdg = 6 status = 123 1 0 pdg = 6 status = 143 1 0 pdg = 6 status = 3 1 0 pdg = 6 status = 155 1 0 pdg = 24 status = 123 0 0 pdg = 24 status = 3 0 0 pdg = 24 status = 155 0 0 pdg = 4 status = 124 1 0 pdg = 4 status = 144 1 0 pdg = 4 status = 2 1 0 pdg = 21 status = 2 1 0 pdg = 4 status = 158 1 0 pdg = -2 status = 158 1 0 pdg = 91 status = 183 0 0 pdg = 411 status = 197 0 1 which would indicate that the identified lepton actually comes fomr a hadron. Do I miss something? Thanks in advance Dominic -- /---------------------------------------------------------------------\ | Dr. Dominic Hirschbuehl | \ Bergische Universitaet Wuppertal - Exp. Elementarteilchenphysik / / hirsch at physik.uni-wuppertal.de / dominic.hirschbuehl at cern.ch \ | office : D.09.22 phone : 0049 - 202 - 439 - 3751 | \---------------------------------------------------------------------/ -------------- next part -------------- A non-text attachment was scrubbed... Name: not available Type: application/pgp-signature Size: 230 bytes Desc: not available URL: <https://www.hepforge.org/lists-archive/rivet/attachments/20170125/c3ec5717/attachment.pgp>
More information about the Rivet mailing list |