[Rivet] UnstableFinalState::particlesByPt() not working?

[Kiran Joshi]

Hi David,

I've attached the file RIVET_TEST_ANA.cc, which I build with:
rivet-buildplugin RivetRIVET_TEST_ANA.so RIVET_TEST_ANA.cc

When I run the analysis on a sample of events generated with Pythia8.145 (using the example files main32.cc & main32.cmnd which come with Pythia) I get the following output:
Rivet.Analysis.RIVET_TEST_ANA: INFO  n(final state particles) = 326
Rivet.Analysis.RIVET_TEST_ANA: INFO  n(unstable final state particles) = 496
Rivet.Analysis.RIVET_TEST_ANA: INFO  n(unstable final state particlesByPt) = 0

Rivet.Analysis.RIVET_TEST_ANA: INFO  n(final state particles) = 157
Rivet.Analysis.RIVET_TEST_ANA: INFO  n(unstable final state particles) = 232
Rivet.Analysis.RIVET_TEST_ANA: INFO  n(unstable final state particlesByPt) = 0

which I hope demonstrates the problem I described in my previous e-mail.

Thanks,
Kiran
-------------- next part --------------
A non-text attachment was scrubbed...
Name: RIVET_TEST_ANA.cc
Type: application/octet-stream
Size: 2417 bytes
Desc: not available
URL: <http://www.hepforge.org/lists-archive/rivet/attachments/20120126/986791de/attachment.obj>
-------------- next part --------------


On 24 Jan 2012, at 15:33, David Grellscheid wrote:

> Hi Kiran,
> 
> can you please send a minimal analysis file that can be compiled with 'rivet-buildplugin' which shows the problem?
> 
> Thanks,
> 
>  David
> 
> 
> On 24/01/12 14:29, Kiran Joshi wrote:
>> Hello,
>> 
>> When I do something like:
>> const ParticleVector&  unstableFS = applyProjection<UnstableFinalState>(event, "UNSTABLE_FS").particlesByPt();
>> I find that the ParticleVector is empty.
>> 
>> However, when I instead do:
>> const ParticleVector&  unstableFS = applyProjection<UnstableFinalState>(event, "UNSTABLE_FS").particles();
>> I can see the (several hundred) unstable final state particles that I expect to find.
>> 
>> I have been defining my UnstableFinalState as follows:
>> UnstableFinalState unstable_fs(-4.5, 4.5, 0.0*GeV);
>> addProjection(unstable_fs, "UNSTABLE_FS");
>> 
>> Is this a bug, or is particlesByPt() not supposed to be used with the UnstableFinalState?
>> (Or am I just implementing it incorrectly?)
>> 
>> Thanks,
>> Kiran
>> 
>> 
>> _______________________________________________
>> Rivet mailing list
>> Rivet at projects.hepforge.org
>> http://www.hepforge.org/lists/listinfo/rivet