using option '-treename HLTesdTree' for EsdCollector, adding default parameter for...
[u/mrichter/AliRoot.git] / PYTHIA8 / pythia8130 / xmldoc / ExternalDecays.xml
1 <chapter name="External Decays">
2
3 <h2>External Decays</h2>
4
5 <code>DecayHandler</code> is a base class for the external handling of 
6 decays. It is intended for normal particle decays, primarily 
7 <ei>B</ei> mesons and <ei>tau</ei>, and cannot be used to redirect
8 decays of heavy resonances like <ei>t</ei> or <ei>Z^0</ei>.  
9 The user-written derived class is called if a pointer to it has
10 been given with the <aloc href="ProgramFlow">
11 <code>pythia.decayPtr()</code></aloc> 
12 method, where it also is specified which particles it will be called for. 
13 This particle information is accessible with the
14 <aloc href="ParticleDataScheme"><code>doExternalDecay()</code></aloc>
15 method. 
16
17 <p/>
18 There is only one pure virtual method in <code>DecayHandler</code>, 
19 to do the decay: 
20 <method name="virtual bool decay(vector&lt;int&gt;&amp; idProd, vector&lt;double&gt;&amp; 
21 mProd, vector&lt;Vec4&gt;&amp; pProd, int iDec, const Event&amp; event)">
22 where 
23 <argument name="idProd"> is a list of particle PDG identity codes,
24 </argument>
25 <argument name="mProd"> is a list of their respective masses (in GeV), and
26 </argument>
27 <argument name="pProd"> is a list of their respective four-momenta.
28 </argument>
29 </method>
30
31 <p/>
32 At input, these vectors each have size one, so that <code>idProd[0]</code>, 
33 <code>mProd[0]</code> and <code>pProd[0]</code> contain information on the 
34 particle that is to be decayed. At output, the vectors should have 
35 increased by the addition of all the decay products. Even if initially 
36 defined in the rest frame of the mother, the products should have been 
37 boosted so that their four-momenta add up to the <code>pProd[0]</code> of 
38 the decaying particle. 
39
40 <p/>
41 Should it be of interest to know the prehistory of the decaying 
42 particle, e.g. to set some helicity information affecting the 
43 decay angular distribution, the full event record is available 
44 read-only, with info in which slot <code>iDec</code> the decaying particle 
45 is stored.
46
47 <p/>
48 The routine should return <code>true</code> if it managed the decay and 
49 <code>false</code> otherwise, in which case <code>Pythia</code> will try 
50 to do the decay itself. This e.g. means you can choose to do some decay 
51 channels yourself, and leave others to <code>Pythia</code>. To avoid
52 doublecounting, the channels you want to handle should be switched off
53 in the <code>Pythia</code> particle database. In the beginning of the  
54 external <code>decay</code> method you should then return 
55 <code>false</code> with a probability given by the sum of the branching
56 ratios for those channels you do not want to handle yourself.   
57
58 <p/>
59 Note that the decay vertex is always set by <code>Pythia</code>, and that 
60 <ei>B-Bbar</ei> oscillations have already been taken into account, 
61 if they were switched on. Thus <code>idProd[0]</code> may be the opposite 
62 of <code>event[iDec].id()</code>, where the latter provides the code at 
63 production.
64
65 <p/>
66 A sample test program is available in <code>main23.cc</code>, providing 
67 a simple example of how to use this facility.
68
69 </chapter>
70
71 <!-- Copyright (C) 2008 Torbjorn Sjostrand -->