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: EvtVVSPwave.cc
12 : //
13 : // Description: Routine to decay vector-> vector pi pi where the
14 : // decay is S-wave dominated.
15 : //
16 : // Modification history:
17 : //
18 : // Jim Hunt June 4, 2008 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/EvtPDL.hh"
27 : #include "EvtGenBase/EvtVector4C.hh"
28 : #include "EvtGenBase/EvtTensor4C.hh"
29 : #include "EvtGenBase/EvtReport.hh"
30 : #include "EvtGenModels/EvtVVPIPI_WEIGHTED.hh"
31 : #include <string>
32 : using std::endl;
33 :
34 0 : EvtVVPIPI_WEIGHTED::~EvtVVPIPI_WEIGHTED() {}
35 :
36 : std::string EvtVVPIPI_WEIGHTED::getName(){
37 :
38 0 : return "VVPIPI_WEIGHTED";
39 :
40 : }
41 :
42 :
43 : EvtDecayBase* EvtVVPIPI_WEIGHTED::clone(){
44 :
45 0 : return new EvtVVPIPI_WEIGHTED;
46 :
47 0 : }
48 :
49 : void EvtVVPIPI_WEIGHTED::init(){
50 :
51 0 : static EvtId PIP=EvtPDL::getId("pi+");
52 0 : static EvtId PIM=EvtPDL::getId("pi-");
53 0 : static EvtId PI0=EvtPDL::getId("pi0");
54 :
55 : // check that there are 0 arguments
56 0 : checkNArg(0);
57 0 : checkNDaug(3);
58 :
59 0 : checkSpinParent(EvtSpinType::VECTOR);
60 0 : checkSpinDaughter(0,EvtSpinType::VECTOR);
61 :
62 :
63 :
64 0 : if ((!(getDaug(1)==PIP&&getDaug(2)==PIM))&&
65 0 : (!(getDaug(1)==PI0&&getDaug(2)==PI0))) {
66 0 : report(Severity::Error,"EvtGen") << "EvtVVPIPI_WEIGHTED generator expected "
67 0 : << " pi+ and pi- (or pi0 and pi0) "
68 0 : << "as 2nd and 3rd daughter. "<<endl;
69 0 : report(Severity::Error,"EvtGen") << "Will terminate execution!"<<endl;
70 0 : ::abort();
71 : }
72 :
73 0 : }
74 :
75 : void EvtVVPIPI_WEIGHTED::initProbMax() {
76 : //Hard coded... should not be hard to calculate...
77 0 : setProbMax(0.08* 1.13 );
78 0 : }
79 :
80 : double reweight_event(double pipi_mass)
81 : {
82 0 : pipi_mass *= 1000.0;
83 0 : return sqrt(
84 0 : -3.6911336508223251 +
85 0 : 0.019119831948029617 * pipi_mass +
86 0 : -1.8962883732377376e-05 * pipi_mass * pipi_mass
87 : );
88 : }
89 :
90 : void EvtVVPIPI_WEIGHTED::decay( EvtParticle *psi_prime){
91 :
92 0 : psi_prime->initializePhaseSpace(getNDaug(),getDaugs());
93 :
94 : EvtParticle *jpsi,*pi1,*pi2;
95 :
96 0 : jpsi=psi_prime->getDaug(0);
97 0 : pi1=psi_prime->getDaug(1);
98 0 : pi2=psi_prime->getDaug(2);
99 :
100 : // Put phase space results into the daughters.
101 :
102 0 : EvtVector4C ep0,ep1,ep2;
103 :
104 0 : ep0=psi_prime->eps(0);
105 0 : ep1=psi_prime->eps(1);
106 0 : ep2=psi_prime->eps(2);
107 :
108 0 : EvtVector4C e0,e1,e2;
109 :
110 0 : e0 = jpsi->epsParent(0);
111 0 : e1 = jpsi->epsParent(1);
112 0 : e2 = jpsi->epsParent(2);
113 :
114 0 : double mass2 = (pi1->getP4()+pi2->getP4()).mass2();
115 :
116 0 : double fac = mass2-4*pi1->mass()*pi2->mass();
117 :
118 0 : fac *= reweight_event(sqrt(mass2));
119 :
120 0 : vertex(0,0,fac*(ep0*e0.conj()));
121 0 : vertex(0,1,fac*(ep0*e1.conj()));
122 0 : vertex(0,2,fac*(ep0*e2.conj()));
123 :
124 0 : vertex(1,0,fac*(ep1*e0.conj()));
125 0 : vertex(1,1,fac*(ep1*e1.conj()));
126 0 : vertex(1,2,fac*(ep1*e2.conj()));
127 :
128 0 : vertex(2,0,fac*(ep2*e0.conj()));
129 0 : vertex(2,1,fac*(ep2*e1.conj()));
130 0 : vertex(2,2,fac*(ep2*e2.conj()));
131 :
132 : return ;
133 :
134 0 : }
135 :
|
