<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=us-ascii">
</head>
<body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">
Hi Rivet Experts,Â
<div class=""><br class="">
</div>
<div class="">I think there might be something happening with the WFinder in the latest version of Rivet. Ive attached two invariant mass plots here and [1] is using the WFinder and [2] is what i estimate by hand. You can see a big difference in the distributions
 even for pp at 2.76 TeV. When we do our heavy ion events we lose a whole lot of Z's which are recovered by estimating by hand. Yes, there are a few differences in the two methods, but between 0.1 and 0.2 in delR and missing ET cut, i dont think it should be
 such a huge difference? </div>
<div class=""><br class="">
</div>
<div class="">Cheers</div>
<div class="">Raghav</div>
<div class=""><br class="">
</div>
<div class="">[1]       WFinder wfinder(fs, Cuts::abseta < 3.0 && Cuts::pT > 20*GeV, PID::MUON, 0.0*GeV, 1000.0*GeV, 25.0*GeV, 0.2);</div>
<div class=""><br class="">
</div>
<div class=""><img apple-inline="yes" id="86C8DB03-6735-43A2-B098-3AFD71E54D83" height="422" width="467" apple-width="yes" apple-height="yes" src="cid:84252C3A-9047-425F-90A1-0527877B491B@cern.ch" class=""></div>
<div class="">[2] in the init() method: </div>
<div class="">      IdentifiedFinalState allleptons;</div>
<div class="">      // allleptons.acceptIdPair(PID::ELECTRON);</div>
<div class="">      allleptons.acceptIdPair(PID::MUON);</div>
<div class="">      Cut cuts = Cuts::abseta < 2.5 && Cuts::pT > 20*GeV;</div>
<div class="">      DressedLeptons leptons(fs, allleptons, 0.1, cuts);</div>
<div class="">      addProjection(leptons, "leptons");</div>
<div class=""><br class="">
</div>
<div class="">      // Leading neutrinos for Etmiss</div>
<div class="">      LeadingParticlesFinalState neutrinos(fs);</div>
<div class="">      // neutrinos.addParticleIdPair(PID::NU_E);</div>
<div class="">      neutrinos.addParticleIdPair(PID::NU_MU);</div>
<div class="">      neutrinos.setLeadingOnly(true);</div>
<div class="">      addProjection(neutrinos, "neutrinos");</div>
<div class=""><br class="">
</div>
<div class="">in the analyze() method:</div>
<div class="">
<div class="">      // implementation of simple Wfinder </div>
<div class="">      const vector<DressedLepton>& leptons = applyProjection<DressedLeptons>(event, "leptons").dressedLeptons();</div>
<div class="">      Particles neutrinos = applyProjection<FinalState>(event, "neutrinos").particlesByPt();</div>
<div class="">      if (leptons.size() != 1 || (neutrinos.size() == 0)) {</div>
<div class=""><span class="Apple-tab-span" style="white-space:pre"></span>vetoEvent;</div>
<div class="">      }</div>
<div class=""><br class="">
</div>
<div class="">      FourMomentum lepton = leptons[0].momentum();</div>
<div class="">      FourMomentum p_miss = neutrinos[0].momentum();</div>
<div class="">      if (p_miss.Et() < 25.0*GeV) {</div>
<div class=""><span class="Apple-tab-span" style="white-space:pre"></span>vetoEvent;</div>
<div class="">      }</div>
<div class=""><br class="">
</div>
<div class="">      double invMass = sqrt(2.0 * lepton.pT()/GeV * p_miss.pT()/GeV * (cosh(lepton.eta()-p_miss.eta()) - cos( lepton.phi()-p_miss.phi()) ) );</div>
<div class="">      _h_W_massCheck->fill(invMass, weight);</div>
</div>
<div class=""><br class="">
</div>
<div class=""><img apple-inline="yes" id="87797FE5-27C4-46FD-9A2F-6CB9E342F5F7" height="422" width="467" apple-width="yes" apple-height="yes" src="cid:375AF03C-519A-4CB4-B220-45CBE3297EAA@cern.ch" class=""></div>
</body>
</html>