BOSS 7.0.9
BESIII Offline Software System
Loading...
Searching...
No Matches
BesEvtGen-00-04-08/src/EvtGen/EvtGenModels/EvtChi1BB2.cc
Go to the documentation of this file.
1//--------------------------------------------------------------------------
2//
3// Environment:
4// This software is part of models developed at BES collaboration
5// based on the EvtGen framework. If you use all or part
6// of it, please give an appropriate acknowledgement.
7//
8// Copyright Information: See EvtGen/BesCopyright
9// Copyright (A) 2006 Ping Rong-Gang @IHEP
10//
11// Module: EvtChi1BB2.cc
12//
13// Description: Routine to decay Chi1 to B10 B10bar.
14//
15// Modification history:
16//
17// Pang C.Y. and Ping R.-G. April, 2007 Module created
18//
19//------------------------------------------------------------------------
20//
21#include "EvtGenBase/EvtPatches.hh"
22#include <stdlib.h>
23#include "EvtGenBase/EvtParticle.hh"
24#include "EvtGenBase/EvtGenKine.hh"
25#include "EvtGenBase/EvtPDL.hh"
26#include "EvtGenBase/EvtVector4C.hh"
27#include "EvtGenBase/EvtVector4R.hh"
28#include "EvtGenBase/EvtTensor4C.hh"
29#include "EvtGenBase/EvtReport.hh"
30#include "EvtGenBase/EvtdFunction.hh"
31#include "EvtGenBase/EvtHelSys.hh"
32#include "EvtChi1BB2.hh"
33#include <string>
34using std::endl;
35using std::cout;
36EvtChi1BB2::~EvtChi1BB2() {}
37
38void EvtChi1BB2::getName(std::string& model_name){
39
40 model_name="Chi1BB2";
41
42}
43
44
45EvtDecayBase* EvtChi1BB2::clone(){
46
47 return new EvtChi1BB2;
48
49}
50
51void EvtChi1BB2::init(){
52// check that there are 1 arguments
53 checkNArg(4);
54 checkNDaug(2);
55 checkSpinParent(EvtSpinType::VECTOR);
56 checkSpinDaughter(0,EvtSpinType::RARITASCHWINGER);
57 checkSpinDaughter(1,EvtSpinType::RARITASCHWINGER);
58
59 }
60
61/*
62void EvtChi1BB2::initProbMax() {
63
64 //Hard coded... should not be hard to calculate...
65 setProbMax(50.0);
66
67}
68*/
69
70void EvtChi1BB2::decay( EvtParticle *p){
71
72 p->initializePhaseSpace(getNDaug(),getDaugs());
73
74 EvtParticle *v,*s1;
75 EvtVector4R pv,ps,ppr;
76
77 v =p->getDaug(0);
78 s1=p->getDaug(1);
79 pv=v->getP4();
80 ps=s1->getP4();
81 ppr=p->getP4();
82
83// Put phase space results into the daughters.
84 EvtHelSys angles(ppr,pv); //using helicity sys. angles
85 double theta =angles.getHelAng(1);
86 double phi =angles.getHelAng(2);
87 double gamma=0;
88
89/* double alpha=getArg(0);
90 double F01=sqrt((1+alpha)/2);
91 double F00=sqrt((1-alpha)/4);
92*/
93 double F00=getArg(0);
94 double F01=getArg(1);
95 double F03=getArg(2);
96 double F11=getArg(3);
97 double F33=-F00; double F32=-F01; double F30=-F03;
98 double F10=F01; double F23=-F10; double F22=-F11;
99//cout<<"********F00,F01,F03,F10,F11,F22,F23,F30,F32,F33="<<F00<<F01<<F03<<F10<<F11<<F22<<F23<<F30<<F32<<F33<<endl;
100/*cout<<"********Djmn(1, 1, 0,phi,theta,gamma)="<<Djmn(1, 1, 0,phi,theta,gamma);
101cout<<"********Djmn(1, 1,-1,phi,theta,gamma)="<<Djmn(1, 1,-1,phi,theta,gamma);
102cout<<"********Djmn(1, 1, 1,phi,theta,gamma)="<<Djmn(1, 1, 1,phi,theta,gamma);
103*/
104 int j1,j2,j3;
105 for(j1=0;j1<=2;){
106 for(j2=0;j2<=3;){
107 for(j3=0;j3<=3;){
108 vertex(j1,j2,j3,0.0);j3++;}
109j2++;}
110j1++;}
111
112 vertex(0,0,0,Djmn(1, 1, 0,phi,theta,gamma)*F00);//Chi1 helicity =1 corresponding index=0
113 vertex(0,0,1,Djmn(1, 1,-1,phi,theta,gamma)*F01);
114 vertex(0,0,3,Djmn(1, 1, 1,phi,theta,gamma)*F03);
115 vertex(0,1,0,Djmn(1, 1, 1,phi,theta,gamma)*F10);
116 vertex(0,1,1,Djmn(1, 1, 0,phi,theta,gamma)*F11);
117 vertex(0,2,2,Djmn(1, 1, 0,phi,theta,gamma)*F22);
118 vertex(0,2,3,Djmn(1, 1,-1,phi,theta,gamma)*F23);
119 vertex(0,3,0,Djmn(1, 1,-1,phi,theta,gamma)*F30);
120 vertex(0,3,2,Djmn(1, 1, 1,phi,theta,gamma)*F32);
121 vertex(0,3,3,Djmn(1, 1, 0,phi,theta,gamma)*F33);
122
123 vertex(1,0,0,Djmn(1,-1, 0,phi,theta,gamma)*F00);//Chi1 helicity =-1 corresponding index=1
124 vertex(1,0,1,Djmn(1,-1,-1,phi,theta,gamma)*F01);
125 vertex(1,0,3,Djmn(1,-1, 1,phi,theta,gamma)*F03);
126 vertex(1,1,0,Djmn(1,-1, 1,phi,theta,gamma)*F10);
127 vertex(1,1,1,Djmn(1,-1, 0,phi,theta,gamma)*F11);
128 vertex(1,2,2,Djmn(1,-1, 0,phi,theta,gamma)*F22);
129 vertex(1,2,3,Djmn(1,-1,-1,phi,theta,gamma)*F23);
130 vertex(1,3,0,Djmn(1,-1,-1,phi,theta,gamma)*F30);
131 vertex(1,3,2,Djmn(1,-1, 1,phi,theta,gamma)*F32);
132 vertex(1,3,3,Djmn(1,-1, 0,phi,theta,gamma)*F33);
133
134 vertex(2,0,0,Djmn(1, 0, 0,phi,theta,gamma)*F00);//Chi1 helicity =0 corresponding index=2
135 vertex(2,0,1,Djmn(1, 0,-1,phi,theta,gamma)*F01);
136 vertex(2,0,3,Djmn(1, 0, 1,phi,theta,gamma)*F03);
137 vertex(2,1,0,Djmn(1, 0, 1,phi,theta,gamma)*F10);
138 vertex(2,1,1,Djmn(1, 0, 0,phi,theta,gamma)*F11);
139 vertex(2,2,2,Djmn(1, 0, 0,phi,theta,gamma)*F22);
140 vertex(2,2,3,Djmn(1, 0,-1,phi,theta,gamma)*F23);
141 vertex(2,3,0,Djmn(1, 0,-1,phi,theta,gamma)*F30);
142 vertex(2,3,2,Djmn(1, 0, 1,phi,theta,gamma)*F32);
143 vertex(2,3,3,Djmn(1, 0, 0,phi,theta,gamma)*F33);
144
145 return ;
146
147}
148
v
**********Class see also m_nmax DOUBLE PRECISION m_amel DOUBLE PRECISION m_x2 DOUBLE PRECISION m_alfinv DOUBLE PRECISION m_Xenph INTEGER m_KeyWtm INTEGER m_idyfs DOUBLE PRECISION m_zini DOUBLE PRECISION m_q2 DOUBLE PRECISION m_Wt_KF DOUBLE PRECISION m_WtCut INTEGER m_KFfin *COMMON c_KarLud $ !Input CMS energy[GeV] $ !CMS energy after beam spread beam strahlung[GeV] $ !Beam energy spread[GeV] $ !z boost due to beam spread $ !electron beam mass *ff pair spectrum $ !minimum v
Definition: KarLud.h:35
Djmn
EvtComplex Djmn(int j, int m, int n, double phi, double theta, double gamma)
Definition: bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtHelSys.cc:151
EvtDecayBase::checkSpinDaughter
void checkSpinDaughter(int d1, EvtSpinType::spintype sp)
Definition: bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtDecayBase.cc:533
EvtDecayBase::checkNArg
void checkNArg(int a1, int a2=-1, int a3=-1, int a4=-1)
Definition: bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtDecayBase.cc:482
EvtParticle::initializePhaseSpace
double initializePhaseSpace(int numdaughter, EvtId *daughters, double poleSize=-1., int whichTwo1=0, int whichTwo2=1)
Definition: bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtParticle.cc:1071