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) 2000 Caltech, UCSB
10 : //
11 : // Module: EvtGen/EvtRaritaSchwing.hh
12 : //
13 : // Description:Class to handle spin 3/2
14 : //
15 : // Modification history:
16 : //
17 : // RYD April 23, 2000 Module created
18 : //
19 : //------------------------------------------------------------------------
20 :
21 : #ifndef EVTRARITASCHWINGER_HH
22 : #define EVTRARITASCHWINGER_HH
23 :
24 : #include "EvtGenBase/EvtComplex.hh"
25 : #include "EvtGenBase/EvtVector4C.hh"
26 : #include "EvtGenBase/EvtVector4R.hh"
27 : #include "EvtGenBase/EvtDiracSpinor.hh"
28 :
29 : class EvtRaritaSchwinger ;
30 : EvtRaritaSchwinger rotateEuler(const EvtRaritaSchwinger& rs,
31 : double alpha,double beta,double gamma);
32 : EvtRaritaSchwinger boostTo(const EvtRaritaSchwinger& rs,
33 : const EvtVector4R p4);
34 : EvtRaritaSchwinger boostTo(const EvtRaritaSchwinger& rs,
35 : const EvtVector3R boost);
36 : EvtRaritaSchwinger dirProd(EvtVector4R v,EvtDiracSpinor u);
37 : EvtRaritaSchwinger dirProd(EvtVector4C v,EvtDiracSpinor u);
38 : EvtRaritaSchwinger operator+(const EvtRaritaSchwinger& u1,
39 : const EvtRaritaSchwinger& u2);
40 : EvtRaritaSchwinger operator-(const EvtRaritaSchwinger& u1,
41 : const EvtRaritaSchwinger& u2);
42 : EvtComplex operator*(const EvtRaritaSchwinger& u1,
43 : const EvtRaritaSchwinger& u2);
44 :
45 : class EvtRaritaSchwinger{
46 :
47 : friend EvtRaritaSchwinger rotateEuler(const EvtRaritaSchwinger& rs,
48 : double alpha,double beta,double gamma);
49 : friend EvtRaritaSchwinger boostTo(const EvtRaritaSchwinger& rs,
50 : const EvtVector4R p4);
51 : friend EvtRaritaSchwinger boostTo(const EvtRaritaSchwinger& rs,
52 : const EvtVector3R boost);
53 :
54 : friend EvtRaritaSchwinger dirProd(EvtVector4R v,EvtDiracSpinor u);
55 : friend EvtRaritaSchwinger dirProd(EvtVector4C v,EvtDiracSpinor u);
56 :
57 : friend EvtRaritaSchwinger operator+(const EvtRaritaSchwinger& u1,
58 : const EvtRaritaSchwinger& u2);
59 : friend EvtRaritaSchwinger operator-(const EvtRaritaSchwinger& u1,
60 : const EvtRaritaSchwinger& u2);
61 :
62 : friend EvtComplex operator*(const EvtRaritaSchwinger& u1,
63 : const EvtRaritaSchwinger& u2);
64 :
65 : public:
66 :
67 : inline EvtRaritaSchwinger();
68 : virtual ~EvtRaritaSchwinger();
69 : inline EvtRaritaSchwinger(const EvtRaritaSchwinger& rs);
70 : inline EvtRaritaSchwinger& operator=(const EvtRaritaSchwinger& rs);
71 :
72 : void set(int i,int j,const EvtComplex& sp);
73 :
74 : void applyRotateEuler(double alpha,double beta,double gamma);
75 : void applyBoostTo(const EvtVector4R p4);
76 : void applyBoostTo(const EvtVector3R boost);
77 :
78 : EvtRaritaSchwinger& operator+=(const EvtRaritaSchwinger& u2);
79 : EvtRaritaSchwinger& operator-=(const EvtRaritaSchwinger& u2);
80 :
81 : EvtComplex get(int i,int j) const;
82 : friend std::ostream& operator<<(std::ostream& s, const EvtRaritaSchwinger& rs);
83 :
84 : EvtVector4C getVector(int i) const;
85 : EvtDiracSpinor getSpinor(int i) const;
86 :
87 : void setVector(int i,const EvtVector4C& v);
88 : void setSpinor(int i,const EvtDiracSpinor& sp);
89 :
90 :
91 :
92 : private:
93 :
94 : //First index in spinor index, second is Lorentz index.
95 : EvtComplex _rs[4][4];
96 :
97 : };
98 :
99 0 : EvtRaritaSchwinger::EvtRaritaSchwinger(){
100 :
101 : int i,j;
102 0 : for(i=0;i<4;i++){
103 0 : for(j=0;j<4;j++){
104 0 : _rs[i][j]=0.0;
105 : }
106 : }
107 :
108 0 : }
109 :
110 0 : EvtRaritaSchwinger::EvtRaritaSchwinger(const EvtRaritaSchwinger& rs){
111 :
112 : int i,j;
113 0 : for(i=0;i<4;i++){
114 0 : for(j=0;j<4;j++){
115 0 : _rs[i][j]=rs._rs[i][j];
116 : }
117 : }
118 :
119 0 : }
120 :
121 : EvtRaritaSchwinger& EvtRaritaSchwinger::operator=(const EvtRaritaSchwinger& rs){
122 :
123 : int i,j;
124 0 : for(i=0;i<4;i++){
125 0 : for(j=0;j<4;j++){
126 0 : _rs[i][j]=rs._rs[i][j];
127 : }
128 : }
129 :
130 0 : return *this;
131 :
132 : }
133 :
134 : #endif
135 :
136 :
|
