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: EvtSVSCP.cc
12 : //
13 : // Description: Routine to decay scalar -> vectors scalar
14 : // and has CP violation.
15 : //
16 : // Modification history:
17 : //
18 : // RYD Febuary 10, 1997 Module created
19 : //
20 : //------------------------------------------------------------------------
21 : //
22 : #include "EvtGenBase/EvtPatches.hh"
23 : #include <stdlib.h>
24 : #include "EvtGenBase/EvtParticle.hh"
25 : #include "EvtGenBase/EvtGenKine.hh"
26 : #include "EvtGenBase/EvtCPUtil.hh"
27 : #include "EvtGenBase/EvtVector4C.hh"
28 : #include "EvtGenBase/EvtPDL.hh"
29 : #include "EvtGenBase/EvtReport.hh"
30 : #include "EvtGenModels/EvtSVSCP.hh"
31 : #include <string>
32 : #include "EvtGenBase/EvtConst.hh"
33 :
34 0 : EvtSVSCP::~EvtSVSCP() {}
35 :
36 : std::string EvtSVSCP::getName(){
37 :
38 0 : return "SVS_CP";
39 :
40 : }
41 :
42 :
43 : EvtDecayBase* EvtSVSCP::clone(){
44 :
45 0 : return new EvtSVSCP;
46 :
47 0 : }
48 :
49 : void EvtSVSCP::init(){
50 :
51 : // check that there are 7 arguments
52 0 : checkNArg(7);
53 0 : checkNDaug(2);
54 :
55 0 : checkSpinParent(EvtSpinType::SCALAR);
56 :
57 0 : checkSpinDaughter(0,EvtSpinType::VECTOR);
58 0 : checkSpinDaughter(1,EvtSpinType::SCALAR);
59 :
60 :
61 0 : }
62 :
63 :
64 : void EvtSVSCP::initProbMax(){
65 :
66 : //This is probably not quite right, but it should do as a start...
67 : //Anders
68 :
69 0 : setProbMax(2*(getArg(3)*getArg(3)+getArg(5)*getArg(5)));
70 :
71 0 : }
72 :
73 :
74 : void EvtSVSCP::decay( EvtParticle *p ){
75 :
76 : //added by Lange Jan4,2000
77 0 : static EvtId B0=EvtPDL::getId("B0");
78 0 : static EvtId B0B=EvtPDL::getId("anti-B0");
79 :
80 : EvtParticle *v;
81 0 : p->initializePhaseSpace(getNDaug(),getDaugs());
82 0 : v = p->getDaug(0);
83 0 : EvtVector4R momv = v->getP4();
84 0 : EvtVector4R moms = p->getDaug(1)->getP4();
85 0 : double massv = v->mass();
86 0 : double t;
87 0 : EvtId other_b;
88 :
89 0 : EvtCPUtil::getInstance()->OtherB(p,t,other_b,0.5);
90 :
91 0 : EvtComplex amp;
92 :
93 0 : EvtComplex A,Abar;
94 :
95 0 : A=EvtComplex(getArg(3)*cos(getArg(4)),getArg(3)*sin(getArg(4)));
96 0 : Abar=EvtComplex(getArg(5)*cos(getArg(6)),getArg(5)*sin(getArg(6)));
97 :
98 0 : if (other_b==B0B){
99 0 : amp=A*cos(getArg(1)*t/(2*EvtConst::c))+
100 0 : EvtComplex(cos(-2.0*getArg(0)),sin(-2.0*getArg(0)))*
101 0 : getArg(2)*EvtComplex(0.0,1.0)*Abar*sin(getArg(1)*t/(2*EvtConst::c));
102 0 : }
103 0 : if (other_b==B0){
104 0 : amp=A*EvtComplex(cos(2.0*getArg(0)),sin(2.0*getArg(0)))*
105 0 : EvtComplex(0.0,1.0)*sin(getArg(1)*t/(2*EvtConst::c))+
106 0 : getArg(2)*Abar*cos(getArg(1)*t/(2*EvtConst::c));
107 0 : }
108 :
109 0 : EvtVector4R p4_parent;
110 :
111 0 : p4_parent=momv+moms;
112 :
113 0 : double norm=massv/(momv.d3mag()*p4_parent.mass());
114 :
115 0 : vertex(0,amp*norm*p4_parent*(v->epsParent(0)));
116 0 : vertex(1,amp*norm*p4_parent*(v->epsParent(1)));
117 0 : vertex(2,amp*norm*p4_parent*(v->epsParent(2)));
118 :
119 : return ;
120 0 : }
121 :
122 : std::string EvtSVSCP::getParamName(int i) {
123 0 : switch(i) {
124 : case 0:
125 0 : return "weakPhase";
126 : case 1:
127 0 : return "deltaM";
128 : case 2:
129 0 : return "finalStateCP";
130 : case 3:
131 0 : return "Af";
132 : case 4:
133 0 : return "AfPhase";
134 : case 5:
135 0 : return "Abarf";
136 : case 6:
137 0 : return "AbarfPhase";
138 : default:
139 0 : return "";
140 : }
141 0 : }
142 :
143 : std::string EvtSVSCP::getParamDefault(int i) {
144 0 : switch(i) {
145 : case 3:
146 0 : return "1.0";
147 : case 4:
148 0 : return "0.0";
149 : case 5:
150 0 : return "1.0";
151 : case 6:
152 0 : return "0.0";
153 : default:
154 0 : return "";
155 : }
156 0 : }
|
