<html>
<head>
<meta content="text/html; charset=ISO-8859-1"
http-equiv="Content-Type">
</head>
<body bgcolor="#FFFFFF" text="#000000">
Hi again,<br>
<br>
I narrowed it down to the UnstableFinalState Projection,<br>
source code, line 34:<br>
<br>
<blockquote>!isZero(p->momentum().perp())<br>
</blockquote>
If a gun produces tau decays, the taus are produced at rest,
resulting in them being<br>
rejected from the accepted particles.<br>
<br>
Any ideas how to cope with this? Simply dropping the cut on the
transverse momentum<br>
makes my code work as expected but I have no feeling what it would
break in other cases.<br>
<br>
Alternatively, I could put some reconstruction code in the
taufinder.<br>
<br>
Holger<br>
<br>
<br>
<br>
<br>
On 02/08/14 22:16, Holger Schulz wrote:
<blockquote cite="mid:53DD5537.6080900@durham.ac.uk" type="cite">
<meta http-equiv="content-type" content="text/html;
charset=ISO-8859-1">
Hi,<br>
<br>
I am doing some tau stuff and came across a problem I do not know
how to tackle.<br>
<br>
When I look for taus coming from, e.g. an electron-positron
collision, everything is fine.<br>
However, when trying to do the same with a tau-gun, the taus
themselves are <br>
to be found in the incoming particles of a vertex and the
TauFinder is unable to find<br>
them. Can someone with some more experience with HepMC give me
some pointers<br>
how to do this properly? It might also be a generator issue
(Sherpa in DECAYER mode).<br>
<br>
Thanks,<br>
Holger<br>
<tt><br>
</tt>
<blockquote><tt>GenEvent: #1 ID=0 SignalProcessGenVertex
BarcodisZero(e: 0<br>
Momenutm units: GEV Position units: MM<br>
Cross Section: 0 +/- 0<br>
Entries this event: 2 vertices, 5 particles.<br>
Beam Particles are not defined.<br>
RndmState(0)=<br>
Wgts(1)=(0,1) <br>
EventScale -1 [energy] alphaQCD=-1 alphaQED=-1<br>
GenParticle Legend<br>
Barcode PDG ID ( Px, Py, Pz, E
) Stat DecayVtx<br>
________________________________________________________________________________<br>
GenVertex: -1 ID: 6 (X,cT):0<br>
Wgts(1)=1 <br>
I: 1 10001 15
+0.00e+00,+0.00e+00,+0.00e+00,+1.78e+00 2 -1<br>
O: 1 10002 15
+0.00e+00,+0.00e+00,+0.00e+00,+1.78e+00 2 -2<br>
Vertex: -2 ID: 6
(X,cT)=+0.00e+00,+0.00e+00,+0.00e+00,+8.17e-02<br>
Wgts(1)=1 <br>
I: 1 10002 15
+0.00e+00,+0.00e+00,+0.00e+00,+1.78e+00 2 -2<br>
O: 3 10003 -211
-4.88e-01,+1.27e-01,-1.85e-01,+5.55e-01 1<br>
10004 130
+1.86e-02,+3.71e-02,-2.53e-01,+5.60e-01 1<br>
10005 16
+4.69e-01,-1.64e-01,+4.38e-01,+6.62e-01 1</tt><br>
<br>
</blockquote>
<br>
<br>
<fieldset class="mimeAttachmentHeader"></fieldset>
<br>
<pre wrap="">_______________________________________________
Rivet mailing list
<a class="moz-txt-link-abbreviated" href="mailto:Rivet@projects.hepforge.org">Rivet@projects.hepforge.org</a>
<a class="moz-txt-link-freetext" href="https://www.hepforge.org/lists/listinfo/rivet">https://www.hepforge.org/lists/listinfo/rivet</a></pre>
</blockquote>
<br>
</body>
</html>