|
[Rivet] Exclude the charged leptons from the charged particle selectionAndy Buckley andy.buckley at cern.chThu Feb 15 13:59:56 GMT 2018
Hi Shayma, You've not actually said what the event selection involves, but I assume you are wanting to veto the prompt charged leptons, i.e. those from the hard process. But this configuration will also veto soft leptons from hadron decays: there are more of those than you might think. You might find that the PromptFinalState projection helps. These days the Identified and Vetoed FinalStates are not needed so much, since you can use a constructor like PromptFinalState(Cuts::abspid != PID::ELECTRON && Cuts::abspid != PID::MUON) [there are probably/hopefully no stable taus in your events!] Andy Dr Andy Buckley, Lecturer / Royal Society University Research Fellow Particle Physics Experiment Group, University of Glasgow On Feb 14 2018, at 8:42 pm, Deepak Kar <deepak.kar at cern.ch> wrote: > > Hi Shayma, > > May be I am being naive here, but would not that depend on the process you are generating? I mean if you for example generate dileptonic ttbar events at leading order, may be you dont have enough pions? > > Cheers, > Deepak > > > On Wed, Feb 14, 2018 at 9:26 PM, Shayma Wahdan <shayma.wahdan at cern.ch (mailto:shayma.wahdan at cern.ch)> wrote: > > Dear Rivet community, > > > > I am trying to exclude the charged leptons from the charged particle selection > > what I do exactly: > > > > in the int part I add: > > > > ChargedFinalState cfs(-2.5, 2.5, 1.0*GeV); > > addProjection(cfs, "CFS"); > > > > // to exclude the charged leptons > > > > VetoedFinalState vfs(cfs); > > vfs.addVetoPairId(PID::ELECTRON); > > vfs.addVetoPairId(PID::MUON); > > vfs.addVetoPairId(PID::TAU); > > addProjection(vfs,"VFS"); > > > > > > and in the analyze > > const ChargedFinalState& ChargedParticleCollection = applyProjection<VetoedFinalState>(event, "VFS"); > > The compliling goes well, while in the running almost none of the events get selected. > > Any suggestions please?? > > Note: I got the idea from this Rivet analysis page http://rivet.hepforge.org/code/2.2.1/a00496_source.html > > Best Regards, > > Shayma' > > > > > > > > > > > > (function(){(function(el) { function removeFromDom(alert) { if (alert.parentNode) if (alert.childNodes.length > 1) { for (var f = document.createDocumentFragment(); alert.childNodes.length > 0; ) f.appendChild(alert.childNodes[0]); alert.parentNode.replaceChild(f, alert); } else alert.firstChild ? alert.parentNode.replaceChild(alert.firstChild, alert) : alert.parentNode.removeChild(alert); } function clean(_el) { if (_el) try { for (var alerts = _el.querySelectorAll(".gr_"), len = alerts.length, i = 0; i < len; i++) removeFromDom(alerts[i]); } catch (e) {} } function redefineInnerHTML(el) { try { Object.defineProperty(el, "innerHTML", { get: function() { try { var r = el.ownerDocument.createRange(); r.selectNodeContents(el); var cnt = r.cloneContents(), d = document.createElement("div"); return d.appendChild(cnt), clean(d), d.innerHTML; } catch (e) { return ""; } }, set: function(value) { try { var r = el.ownerDocument.createRange(); r.selectNodeContents(el), r.deleteContents(); var df = r.createContextualFragment(value); el.appendChild(df); } catch (e) {} } }); } catch (e) {} } if (el) { var nativeClone = el.cloneNode; el.cloneNode = function(val) { var n = nativeClone.call(el, val); if (el.classList.contains("mceContentBody")) n.innerHTML = el.innerHTML, clean(n); else try { redefineInnerHTML(n); } catch (e) {} return n; }, redefineInnerHTML(el); } })(document.querySelector("[data-gramm_id='d0a779df-0022-acfa-f632-b9909d937b5b']")) })() > > _______________________________________________ > > Rivet mailing list > > Rivet at projects.hepforge.org (mailto:Rivet at projects.hepforge.org) > > https://www.hepforge.org/lists/listinfo/rivet > > > > -- (https://www.hepforge.org/lists/listinfo/rivet) > > Deepak Kar (https://www.hepforge.org/lists/listinfo/rivet) > > University of Witwatersrand (https://www.hepforge.org/lists/listinfo/rivet) > > Room PM15, School of Physics (https://www.hepforge.org/lists/listinfo/rivet) > > (0027) 011-7176958 (office) (0027) 0736944181 (mobile) (https://www.hepforge.org/lists/listinfo/rivet) > > > > While at CERN: (https://www.hepforge.org/lists/listinfo/rivet) > > Building 1, R-016 (https://www.hepforge.org/lists/listinfo/rivet) > > (0041) 0767321349 (mobile) (https://www.hepforge.org/lists/listinfo/rivet) > > > > While at USA: (https://www.hepforge.org/lists/listinfo/rivet) > > (001) 330-998-1500 (mobile) (https://www.hepforge.org/lists/listinfo/rivet) > > > > > > _______________________________________________ > Rivet mailing list > Rivet at projects.hepforge.org > https://www.hepforge.org/lists/listinfo/rivet > -------------- next part -------------- An HTML attachment was scrubbed... URL: <https://www.hepforge.org/lists-archive/rivet/attachments/20180215/2019bd75/attachment.html>
More information about the Rivet mailing list |