Line data Source code
1 : //--------------------------------------------------------------------------
2 : //
3 : // Environment:
4 : // This software is part of the EvtGen package developed jointly
5 : // for the BaBar and CLEO collaborations. If you use all or part
6 : // of it, please give an appropriate acknowledgement.
7 : //
8 : // Copyright Information: See EvtGen/COPYRIGHT
9 : // Copyright (C) 1998 Caltech, UCSB
10 : //
11 : // Module: EvtVll.cc
12 : //
13 : // Description: The decay of a vector meson to two leptons,
14 : // or generally, two spin 1/2 particles.
15 : // E.g., J/psi -> e+ e-
16 : //
17 : // Modification history:
18 : //
19 : // RYD January 17, 1997 Module created
20 : //
21 : //------------------------------------------------------------------------
22 : //
23 : #include "EvtGenBase/EvtPatches.hh"
24 : #include <stdlib.h>
25 : #include <iostream>
26 : #include <string>
27 : #include "EvtGenBase/EvtParticle.hh"
28 : #include "EvtGenBase/EvtPDL.hh"
29 : #include "EvtGenBase/EvtGenKine.hh"
30 : #include "EvtGenModels/EvtVll.hh"
31 : #include "EvtGenBase/EvtDiracSpinor.hh"
32 : #include "EvtGenBase/EvtReport.hh"
33 : #include "EvtGenBase/EvtVector4C.hh"
34 :
35 0 : EvtVll::~EvtVll() {}
36 :
37 : std::string EvtVll::getName(){
38 :
39 0 : return "VLL";
40 :
41 : }
42 :
43 :
44 : EvtDecayBase* EvtVll::clone(){
45 :
46 0 : return new EvtVll;
47 :
48 0 : }
49 :
50 : void EvtVll::init(){
51 :
52 : // check that there are 0 arguments
53 0 : checkNArg(0);
54 0 : checkNDaug(2);
55 :
56 0 : checkSpinParent(EvtSpinType::VECTOR);
57 :
58 0 : checkSpinDaughter(0,EvtSpinType::DIRAC);
59 0 : checkSpinDaughter(1,EvtSpinType::DIRAC);
60 :
61 0 : }
62 :
63 : void EvtVll::initProbMax(){
64 :
65 0 : setProbMax(1.0);
66 :
67 0 : }
68 :
69 : void EvtVll::decay(EvtParticle *p){
70 :
71 0 : p->initializePhaseSpace(getNDaug(),getDaugs());
72 :
73 : EvtParticle *l1, *l2;
74 0 : l1 = p->getDaug(0);
75 0 : l2 = p->getDaug(1);
76 :
77 0 : EvtVector4C l11, l12, l21, l22;
78 0 : l11=EvtLeptonVCurrent(l1->spParent(0),l2->spParent(0));
79 0 : l12=EvtLeptonVCurrent(l1->spParent(0),l2->spParent(1));
80 0 : l21=EvtLeptonVCurrent(l1->spParent(1),l2->spParent(0));
81 0 : l22=EvtLeptonVCurrent(l1->spParent(1),l2->spParent(1));
82 :
83 0 : EvtVector4C eps0=p->eps(0);
84 0 : EvtVector4C eps1=p->eps(1);
85 0 : EvtVector4C eps2=p->eps(2);
86 :
87 0 : double M2=p->mass();
88 0 : M2*=M2;
89 0 : double m2=l1->mass();
90 0 : m2*=m2;
91 :
92 0 : double norm=1.0/sqrt(2*M2+4*m2-4*m2*m2/M2);
93 :
94 0 : vertex(0,0,0,norm*(eps0*l11));
95 0 : vertex(0,0,1,norm*(eps0*l12));
96 0 : vertex(0,1,0,norm*(eps0*l21));
97 0 : vertex(0,1,1,norm*(eps0*l22));
98 :
99 0 : vertex(1,0,0,norm*(eps1*l11));
100 0 : vertex(1,0,1,norm*(eps1*l12));
101 0 : vertex(1,1,0,norm*(eps1*l21));
102 0 : vertex(1,1,1,norm*(eps1*l22));
103 :
104 0 : vertex(2,0,0,norm*(eps2*l11));
105 0 : vertex(2,0,1,norm*(eps2*l12));
106 0 : vertex(2,1,0,norm*(eps2*l21));
107 0 : vertex(2,1,1,norm*(eps2*l22));
108 :
109 : return;
110 :
111 0 : }
112 :
113 :
114 :
115 :
116 :
117 :
118 :
119 :
120 :
121 :
122 :
123 :
124 :
|
