[Rivet] undocumented behaviour in wfinder/zfinder (bug ?)

[Andy Buckley]

Hi Nicholas,

I'm a bit uncertain about exactly how the remainingFS works for WFinder 
and ZFinder... hopefully someone can clarify and check that there is no 
bug: Frank S, can you help?

We're just about to make a new Rivet release, and I just uploaded the 
release tarballs to the HepForge downloads area... but I'll wait for 
confirmation that there's not a bug before making the release official.

Cheers,
Andy


On 29/07/11 15:29, Nicolas Morange wrote:
> Hello,
>
> I'm writing a Z+jets analysis.
> Thus I use a ZFinder whose remainingFinalState() is chained to a FastJet.
> However, the code in ZFinder
>
> 00044     LeptonClusters leptons(fs, bareleptons, dRmax,
> 00045                            clusterPhotons, excludePhotonsFromRFS,
> 00046                            etaRanges, pTmin);
> 00047     addProjection(leptons, "LeptonClusters");
> 00048     InvMassFinalState imfs(leptons, std::make_pair(pid, -pid), m2_min,
> m2_max);
> 00049     addProjection(imfs, "IMFS");
> 00050
> 00051     VetoedFinalState remainingFS;
> 00052
> remainingFS.addVetoOnThisFinalState(leptons.constituentsFinalState());
> 00053     addProjection(remainingFS, "RFS");
>
> gives me the impression that not only the Z leptons are removed from the
> remainingFS, but also the other leptons in the event.
> Is this correct ?
> And how should I do to remove only the Z leptons (and associated photons) from
> the final state ?
>
> Cheers,
> Nicolas
> _______________________________________________
> Rivet mailing list
> Rivet at projects.hepforge.org
> http://www.hepforge.org/lists/listinfo/rivet
>


-- 
Dr Andy Buckley
SUPA Advanced Research Fellow
Particle Physics Experiment Group, University of Edinburgh

The University of Edinburgh is a charitable body, registered in
Scotland, with registration number SC005336.