<div>Hi Shayma,</div><br><div>You've not actually said what the event selection involves, but I assume you are wanting to veto the <em>prompt</em> 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.</div><br><div>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!]</div><br><div>Andy</div><br><div><signature id="initial"><div><table cellpadding="0" cellspacing="0"><tbody><tr><td colspan="2"><div style="padding-bottom:15px;"><div><strong>Dr Andy Buckley, Lecturer / Royal Society University Research Fellow</strong></div><div>Particle Physics Experiment Group, University of Glasgow</div></div></td></tr><tr><td style="vertical-align:top;"></td><td><div style="font-size:0.9em;white-space:nowrap;border-left:2px solid gray;margin-left:20px;padding-left:20px;"><div><div></div><div></div></div></div></td></tr></tbody></table></div></signature></div><div class="gmail_quote_attribution">On Feb 14 2018, at 8:42 pm, Deepak Kar <deepak.kar@cern.ch> wrote:</div><blockquote><br><div><div><div>Hi Shayma,</div><div><br></div><div>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?</div><div><br></div><div>Cheers,</div><div>Deepak</div></div><div><br><div><div>On Wed, Feb 14, 2018 at 9:26 PM, Shayma Wahdan <<a href="mailto:shayma.wahdan@cern.ch" title="mailto:shayma.wahdan@cern.ch">shayma.wahdan@cern.ch</a>> wrote:</div><blockquote><div><div><div><div><font style="font-size:10pt"><font style="font-family:Tahoma">Dear Rivet community,</font></font></div><br><div><font style="font-size:10pt"><font style="font-family:Tahoma">I am trying to exclude the charged leptons from the charged particle selection</font></font></div><div><font style="font-size:10pt"><font style="font-family:Tahoma">what I do exactly:</font></font></div><br><div><font style="font-size:10pt"><font style="font-family:Tahoma">in the </font></font><font style="font-size:10pt"><font style="font-family:Tahoma"><span style="color:FF0000">int</span></font></font><font style="font-size:10pt"><font style="font-family:Tahoma"> part I add:</font></font></div><div><font style="font-size:10pt"><font style="font-family:Tahoma"> </font></font></div><div><font style="font-size:10pt"><font style="font-family:Tahoma">     </font></font><font style="font-size:10pt"><font style="font-family:Tahoma"><span style="color:3366FF"> ChargedFinalState cfs(-2.5, 2.5, 1.0*GeV);</span></font></font></div><div><font style="font-size:10pt"><font style="font-family:Tahoma"><span style="color:3366FF">      addProjection(cfs, "CFS");</span></font></font></div><br><div><font style="font-size:10pt"><font style="font-family:Tahoma"><span style="color:3366FF">    // to exclude the charged leptons</span></font></font></div><div><font style="font-size:10pt"><font style="font-family:Tahoma"><span style="color:3366FF">    </span></font></font></div><div><font style="font-size:10pt"><font style="font-family:Tahoma"><span style="color:3366FF">     VetoedFinalState vfs(cfs);</span></font></font></div><div><font style="font-size:10pt"><font style="font-family:Tahoma"><span style="color:3366FF">      vfs.addVetoPairId(PID::ELECTRON);</span></font></font></div><div><font style="font-size:10pt"><font style="font-family:Tahoma"><span style="color:3366FF">      vfs.addVetoPairId(PID::MUON);</span></font></font></div><div><font style="font-size:10pt"><font style="font-family:Tahoma"><span style="color:3366FF">      vfs.addVetoPairId(PID::TAU);</span></font></font></div><div><font style="font-size:10pt"><font style="font-family:Tahoma"><span style="color:3366FF">      addProjection(vfs,"VFS");</span></font></font></div><br><br><div><font style="font-size:10pt"><font style="font-family:Tahoma"><span style="color:3366FF"> </span></font></font><font style="font-size:10pt"><font style="font-family:Tahoma"><span style="color:3366FF"><span style="color:000000">and in the</span></span></font></font><font style="font-size:10pt"><font style="font-family:Tahoma"><span style="color:3366FF"> </span></font></font><font style="font-size:10pt"><font style="font-family:Tahoma"><span style="color:3366FF"><span style="color:FF0000">analyze</span></span></font></font></div><br><div><font style="font-size:10pt"><font style="font-family:Tahoma">    </font></font><font style="font-size:10pt"><font style="font-family:Tahoma"><span style="color:3366FF"> const ChargedFinalState& ChargedParticleCollection = applyProjection<VetoedFinalState>(event, "VFS");</span></font></font></div><br><div><font style="font-size:10pt"><font style="font-family:Tahoma">The compliling goes well, while in the running almost none of the events get selected.</font></font></div><br><div><font style="font-size:10pt"><font style="font-family:Tahoma">Any suggestions please??</font></font></div><br><div><font style="font-size:10pt"><font style="font-family:Tahoma">Note: I got the idea from this Rivet analysis page </font></font><font style="font-size:10pt"><font style="font-family:Tahoma"><a href="http://rivet.hepforge.org/code/2.2.1/a00496_source.html" title="http://rivet.hepforge.org/code/2.2.1/a00496_source.html">http://rivet.hepforge.org/code/2.2.1/a00496_source.html</a></font></font></div><br><div><font style="font-size:10pt"><font style="font-family:Tahoma">Best Regards,</font></font></div><div><font style="font-size:10pt"><font style="font-family:Tahoma">Shayma'</font></font></div><br><br><br><br><br></div></div><div>(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']")) })()</div></div><br><div>_______________________________________________</div><div>Rivet mailing list</div><div><a href="mailto:Rivet@projects.hepforge.org" title="mailto:Rivet@projects.hepforge.org">Rivet@projects.hepforge.org</a></div><div><a href="https://www.hepforge.org/lists/listinfo/rivet" title="https://www.hepforge.org/lists/listinfo/rivet">https://www.hepforge.org/lists/listinfo/rivet</a></div><div><br></div><div><a href="https://www.hepforge.org/lists/listinfo/rivet" title="https://www.hepforge.org/lists/listinfo/rivet">--</a></div><div><div><a href="https://www.hepforge.org/lists/listinfo/rivet" title="https://www.hepforge.org/lists/listinfo/rivet">Deepak Kar</a></div><div><a href="https://www.hepforge.org/lists/listinfo/rivet" title="https://www.hepforge.org/lists/listinfo/rivet">University of Witwatersrand</a></div><div><a href="https://www.hepforge.org/lists/listinfo/rivet" title="https://www.hepforge.org/lists/listinfo/rivet">Room PM15, School of Physics</a></div><div><a href="https://www.hepforge.org/lists/listinfo/rivet" title="https://www.hepforge.org/lists/listinfo/rivet">(0027) 011-7176958 (office) (0027) 0736944181 (mobile)</a></div><br><div><a href="https://www.hepforge.org/lists/listinfo/rivet" title="https://www.hepforge.org/lists/listinfo/rivet">While at CERN:</a></div><div><a href="https://www.hepforge.org/lists/listinfo/rivet" title="https://www.hepforge.org/lists/listinfo/rivet">Building 1, R-016</a></div><div><a href="https://www.hepforge.org/lists/listinfo/rivet" title="https://www.hepforge.org/lists/listinfo/rivet">(0041) 0767321349 (mobile)</a></div><br><div><a href="https://www.hepforge.org/lists/listinfo/rivet" title="https://www.hepforge.org/lists/listinfo/rivet">While at USA:</a></div><div><a href="https://www.hepforge.org/lists/listinfo/rivet" title="https://www.hepforge.org/lists/listinfo/rivet">(001) 330-998-1500 (mobile)</a></div></div></blockquote></div></div><div>_______________________________________________</div><div>Rivet mailing list</div><div>Rivet@projects.hepforge.org</div><div>https://www.hepforge.org/lists/listinfo/rivet</div></div></blockquote>